Efficient DEM‐based overland flow routing using integrated recursive algorithms          下载免费PDF全文

Pin‐Chun Huang  Kwan Tun Lee 《水文研究》2017,31(5):1007-1017

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.  相似文献   

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images          下载免费PDF全文

Pei Leng  Xiaoning Song  Si‐Bo Duan  Zhao‐Liang Li 《水文研究》2017,31(6):1398-1407

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.  相似文献   

In‐situ analysis of trace elements and Sr–Pb isotopes of K‐feldspars from Tongshankou Cu–Mo deposit,SE Hubei Province,China: Insights into early potassic alteration of the porphyry mineralization system          下载免费PDF全文

Yao‐Ming Xu  Shao‐Yong Jiang 《地学学报》2017,29(6):343-355

In the Tongshankou porphyry deposit (SE Hubei Province, South China), three types of K‐feldspars are recognized: (I) the phenocryst type in the porphyry that crystalized during the magmatic stage, (II) the megacryst type and (III) the vein type in the altered porphyry and orebody that was produced by hydrothermal fluids. A detailed in‐situ analysis of trace elements and Sr–Pb isotopes was carried out on K‐feldspars in an attempt to unravel their formation processes and to trace the element sources during potassic alteration. The Type III K‐feldspars show lower Sr contents and Sr‐isotope ratios but higher Pb contents and Pb‐isotope ratios than the Type I and II K‐feldspars, possibly reflecting a contribution from the country carbonate rocks with less radiogenic Sr but more radiogenic Pb sources, and indicate that the ore‐forming fluids and materials may have been partially derived from external sources such as the host sedimentary rocks during the early potassic alteration stage.  相似文献   

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications          下载免费PDF全文

I. Baziotis  C.‐H. Tsai  W. G. Ernst  B.‐M. Jahn  Y. Iizuka 《Journal of Metamorphic Geology》2017,35(1):35-54

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.  相似文献   

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin          下载免费PDF全文

Jun‐Zhi Wang  Xiao‐Wei Jiang  Zhi‐Yuan Zhang  Li Wan  Xu‐Sheng Wang  Hailong Li 《水文研究》2017,31(22):4006-4018

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.  相似文献   

Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA     

Bryan G. Moravec  C. Kent Keller  Jeffrey L. Smith  Richelle M. Allen‐King  Angela J. Goodwin  Jerry P. Fairley  Peter B. Larson 《水文研究》2010,24(4):446-460

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Weichselian interstadials at Riipiharju,northern Sweden – interpretation of vegetation and climate from fossil and modern pollen records     

MARTINA HÄTTESTRAND  ANN‐MARIE ROBERTSSON 《Boreas: An International Journal of Quaternary Research》2010,39(2):296-311

Hättestrand, M. & Robertsson, A.‐M. 2010: Weichselian interstadials at Riipiharju, northern Sweden – interpretation of vegetation and climate from fossil and modern pollen records. Boreas, 10.1111/j.1502‐3885.2009.00129.x. ISSN 0300‐9483. The most complete records of Weichselian ice‐free conditions in northern Sweden have been retrieved from kettleholes in the Riipiharju esker. In an earlier study, the Riipiharju I core was described as containing two Weichselian interstadials and Riipiharju was chosen as type site for the second Weichselian interstadial in northern Sweden. Here, we present a palynological investigation of two new sediment cores (Riipiharju II and III) retrieved from Riipiharju. Together, the new cores comprise a late cold part of the first Weichselian interstadial recorded in northeastern Sweden (Tärendö I, earlier correlated with Peräpohjola in Finland) as well as a long sequence of the second Weichselian interstadial (Tärendö II, earlier named Tärendö). The results indicate that the climate during deposition of the Tärendö II sequence was more variable than earlier suggested. According to the present interpretation it was relatively warm in the early part of Tärendö II; thereafter a long cold phase persisted, and finally the climate was warmer again in the late part of Tärendö II. The warm phases are characterized by Betula‐dominant pollen assemblages, while the cold phase is characterized by high percentages of Artemisia and Gramineae pollen. Since there is still no firm chronology established of the interstadials in northeastern Sweden, two possible correlations are discussed; either Tärendö I and II are correlated with Brörup (MIS 5c) and Odderade (MIS 5a), or, perhaps more likely, they are correlated with Odderade and early Middle Weichselian (MIS 3) time.  相似文献   

The data‐driven approach as an operational real‐time flood forecasting model     

Phuoc Khac‐Tien Nguyen  Lloyd Hock‐Chye Chua 《水文研究》2012,26(19):2878-2893

Accurate water level forecasts are essential for flood warning. This study adopts a data‐driven approach based on the adaptive network–based fuzzy inference system (ANFIS) to forecast the daily water levels of the Lower Mekong River at Pakse, Lao People's Democratic Republic. ANFIS is a hybrid system combining fuzzy inference system and artificial neural networks. Five ANFIS models were developed to provide water level forecasts from 1 to 5 days ahead, respectively. The results show that although ANFIS forecasts of water levels up to three lead days satisfied the benchmark, four‐ and five‐lead‐day forecasts were only slightly better in performance compared with the currently adopted operational model. This limitation is imposed by the auto‐ and cross‐correlations of the water level time series. Output updating procedures based on the autoregressive (AR) and recursive AR (RAR) models were used to enhance ANFIS model outputs. The RAR model performed better than the AR model. In addition, a partial recursive procedure that reduced the number of recursive steps when applying the AR or the RAR model for multi‐step‐ahead error prediction was superior to the fully recursive procedure. The RAR‐based partial recursive updating procedure significantly improved three‐, four‐ and five‐lead‐day forecasts. Our study further shows that for long lead times, ANFIS model errors are dominated by lag time errors. Although the ANFIS model with the RAR‐based partial recursive updating procedure provided the best results, this method was able to reduce the lag time errors significantly for the falling limbs only. Improvements for the rising limbs were modest. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Dissolved nitrous oxide (N2O) dynamics in agricultural field drains and headwater streams in an intensive arable catchment          下载免费PDF全文

Zanist Q. Hama‐Aziz  Kevin M. Hiscock  Richard J. Cooper 《水文研究》2017,31(6):1371-1381

Indirect nitrous oxide (N₂O) emissions produced by nitrogen (N) leaching into surface water and groundwater bodies are poorly understood in comparison to direct N₂O emissions from soils. In this study, dissolved N₂O concentrations were measured weekly in both lowland headwater streams and subsurface agricultural field drain discharges over a 2‐year period (2013–2015) in an intensive arable catchment, Norfolk, UK. All field drain and stream water samples were found to have dissolved N₂O concentrations higher than the water–air equilibrium concentration, illustrating that all sites were acting as a net source of N₂O emissions to the atmosphere. Soil texture was found to significantly influence field drain N₂O dynamics, with mean concentrations from drains in clay loam soils (5.3 μg N L^?1) being greater than drains in sandy loam soils (4.0 μg N L^?1). Soil texture also impacted upon the relationships between field drain N₂O concentrations and other water quality parameters (pH, flow rate, and nitrate (NO₃) and nitrite (NO₂) concentrations), highlighting possible differences in N₂O production mechanisms in different soil types. Catchment antecedent moisture conditions influenced the storm event mobilisation of N₂O in both field drains and streams, with the greatest concentration increases recorded during precipitation events preceded by prolonged wet conditions. N₂O concentrations also varied seasonally, with the lowest mean concentrations typically occurring during the summer months (JJA). Nitrogen fertiliser application rates and different soil inversion regimes were found to have no effect on dissolved N₂O concentrations, whereas higher N₂O concentrations recorded in field drains under a winter cover crop compared to fallow fields revealed cover crops are an ineffective greenhouse gas emission mitigation strategy. Overall, this study highlights the complex interactions governing the dynamics of dissolved N₂O concentrations in field drains and headwater streams in a lowland intensive agricultural catchment.  相似文献   

Tropical precipitation anomalies and d‐excess evolution during El Niño 2014‐16          下载免费PDF全文

Ricardo Sánchez‐Murillo  Ana M. Durán‐Quesada  Christian Birkel  Germain Esquivel‐Hernández  Jan Boll 《水文研究》2017,31(4):956-967

The last 2014‐16 El Niño event was among the three strongest episodes on record. El Niño considerably changes annual and seasonal precipitation across the tropics. Here, we present a unique stable isotope data set of daily precipitation collected in Costa Rica prior to, during, and after El Niño 2014‐16, in combination with Lagrangian moisture source and precipitation anomaly diagnostics. δ²H composition ranged from ‐129.4 to +18.1 (‰) while δ¹⁸O ranged from ‐17.3 to +1.0 (‰). No significant difference was observed among δ¹⁸O (P=0.186) and δ²H (P=0.664) mean annual compositions. However, mean annual d‐excess showed a significant decreasing trend (from +13.3 to +8.7 ‰) (P<0.001) with values ranging from +26.6 to ‐13.9 ‰ prior to and during the El Niño evolution. The latter decrease in d‐excess can be partly explained by an enhanced moisture flux convergence across the southeastern Caribbean Sea coupled with moisture transport from northern South America by means of an increased Caribbean Low Level Jet regime. During 2014‐15, precipitation deficit across the Pacific domain averaged 46% resulting in a very severe drought; while a 94% precipitation surplus was observed in the Caribbean domain. Understanding these regional moisture transport mechanisms during a strong El Niño event may contribute to a) better understanding of precipitation anomalies in the tropics and b) re‐evaluate past stable isotope interpretations of ENSO events in paleoclimatic archives within the Central America region.  相似文献