Spatially distributed denitrification in a karst springshed

Karst spring measurements assess biogeochemical processes occurring within groundwater contributing areas to springs (springsheds) but can only provide aggregated information. To better understand spatially distributed processes that comprise these aggregated measures, we investigated aquifer denitrification evidence in groundwater wells (n = 16) distributed throughout a springshed in the Upper Floridan aquifer in northern Florida. Aquifer geochemistry, nitrate isotopes, and dissolved gases were compared against similar measurements at the spring outlet to evaluate spatial heterogeneity of denitrification evidence in relation to land surface–aquifer connectivity. Sample locations spanned spatial variation in recharge processes (i.e., diffuse vs. focused recharge) and proximity to sources of denitrification reactants (e.g., wetlands). Although no distinct spatial pattern in denitrification was uncovered, excess dissolved N₂ gas measurements were only above detection in the unconfined springshed, with some evidence of a wetland proximity effect. Measured oxidation–reduction potential and dissolved oxygen poorly predicted denitrification, indicating that measured denitrification may be occurring upgradient from sampled wells. Despite dramatic spatial chemical heterogeneity across wells, mean values for recharge nitrate concentrations (0.02 to 5.56 mg N L^?1) and excess N₂ from aquifer denitrification (below detection to 1.37 mg N L^?1) corresponded reasonably with mean spring outlet measurements for initial nitrate (0.78 to 1.36 mg N L^?1) and excess N₂ (0.15 to 1.04 mg N L^?1). Congruence between groundwater and spring measurements indicates that combining sampling at the spring outlet and across the springshed is useful for understanding spatial aquifer denitrification. However, this approach would be improved with a high‐density sampling network with transects of wells along distinct groundwater flow paths.     

Agricultural land use and hydrology affect variability of shallow groundwater nitrate concentration in South Florida

South Florida's Miami‐Dade agricultural area is located between two protected natural areas, the Biscayne and Everglades National Parks, subject to the costliest environmental restoration project in history. Agriculture, an important economic activity in the region, competes for land and water resources with the restoration efforts and Miami's urban sprawl. The objective of this study, understanding water quality interactions between agricultural land use and the shallow regional aquifer, is critical to the reduction of agriculture's potentially negative impacts. A study was conducted in a 4‐ha square field containing 0·9 ha of corn surrounded by fallow land. The crop rows were oriented NW–SE along the dominant groundwater flow in the area. A network of 18 monitoring wells was distributed across the field. Shallow groundwater nitrate–nitrogen concentration [N‐NO₃^?] was analyzed on samples collected from the wells biweekly for 3 years. Detailed hydrological (water table elevation [WTE] at each well, groundwater flow direction [GwFD], rainfall) and crop (irrigation, fertilization, calendar) data were also recorded in situ. Flow direction is locally affected by seasonal regional drainage through canal management exercised by the local water authority. The data set was analyzed by dynamic factor analysis (DFA), a specialized time series statistical technique only recently applied in hydrology. In a first step, the observed nitrate variation was successfully described by five common trends representing the unexplained variability. By including the measured hydrological series as explanatory variables the trends were reduced to only three. The analysis yields a quantification of the effects of hydrological factors over local groundwater nitrate concentration. Furthermore, a spatial structure across the field, matching land use, was found in the five remaining common trends whereby the groundwater [N‐NO₃^?] in wells within the corn rows could be generally separated from those in fallow land NW and SE of the crop strip. Fertilization, masked by soil/water/plant‐delayed processes, had no discernible effect on groundwater nitrate levels. Copyright © 2006 John Wiley & Sons, Ltd.     

Spatial variation in soil resistance to flowing water erosion along a regional transect in the Loess Plateau

The factors influencing soil erosion may vary with scale. It remains unclear whether the spatial variation in soil erosion resistance is controlled by regional variables (e.g. precipitation, temperature, and vegetation zone) or by local specific variables (e.g. soil properties, root traits, land use, and farming operations) when the study area enlarges from a hillslope or catchment to the regional scale. This study was performed to quantify the spatial variations in soil erosion resistance to flowing water under three typical land uses along a regional transect on the Loess Plateau and to identify whether regional or local specific variables are responsible for these changes. The results indicated that the measured soil detachment capacities (D_c) of cropland exhibited an irregular trend along the regional transect. The D_c of grassland increased with mean annual precipitation, except for two sites (Yijun and Erdos). The measured D_c of woodland displayed an inverted ‘U’ shape. The changes in rill erodibility (K_r) of three land uses were similar to D_c, whereas no distinguishable trend was found for critical shear stress (τ_c). No significant correlation was detected between D_c, K_r and τ_c, and the regional variables. The spatial variation in soil erosion resistance could be explained reasonably by changes in soil properties, root traits, land use, and farming operations, rather than regional variables. The adjustment coefficient of K_r for grassland and woodland could be well simulated by soil cohesion and root mass density (R² = 0.70, P < 0.01), and the adjustment coefficient of critical shear stress could be estimated with aggregate stability (R² = 0.57, P < 0.01). The results are helpful for quantifying the spatial variation in soil detachment processes by overland flow and to develop process‐based erosion model at a regional scale. Copyright © 2015 John Wiley & Sons, Ltd.     

The impact of pipe flow in riparian peat deposits on nitrate transport and removal

The effect of preferential flow in soil pipes on nitrate retention in riparian zones is poorly understood. The characteristics of soil pipes and their influence on patterns of groundwater transport and nitrate dynamics were studied along four transects in a 1‐ to >3‐m deep layer of peat and marl overlying an oxic sand aquifer in a riparian zone in southern Ontario, Canada. The peat‐marl deposit, which consisted of several horizontal layers with large differences in bulk density, contained soil pipes that were generally 0.1 to 0.2 m in diameter and often extended vertically for 1 to >2 m. Springs that produced overland flow across the riparian area occurred at some sites where pipes extended to the peat surface. Concentrations of NO₃^?–N (20–30 mg L^?1) and dissolved oxygen (DO) (4–6 mg L^?1) observed in peat pipe systems and surface springs were similar to values in the underlying sand aquifer, indicating that preferential flow transported groundwater with limited nitrate depletion. Low NO₃^?–N concentrations of <5 mg L^?1 and enriched δ¹⁵N values indicated that denitrification was restricted to small areas of the peat where pipes were absent. Groundwater DO concentrations declined rapidly to <2 mg L^?1 in the peat matrix adjacent to pipes, whereas high NO₃^?–N concentrations of >15 mg L^?1 extended over a larger zone. Low dissolved organic carbon values at these locations suggest that supplies of organic carbon were not sufficient to support high rates of denitrification, despite low DO conditions. These data indicate that it is important to develop a greater understanding of pipes in peat deposits, which function as sites where the transport of large fluxes of water with low biogeochemical reaction rates can limit the nitrate removal capacity of riparian zones. Copyright © 2011 John Wiley & Sons, Ltd.     

Predicting ground water nitrate concentration from land use

Ground water nitrate concentrations on Nantucket Island, Massachusetts, were analyzed to assess the effects of land use on ground water quality. Exploratory data analysis was applied to historic ground water nitrate concentrations to determine spatial and temporal trends. Maximum likelihood Tobit and logistic regression analyses of explanatory variables that characterize land use within a 1000-foot radius of each well were used to develop predictive equations for nitrate concentration at 69 wells. The results demonstrate that historic nitrate concentrations downgradient from agricultural land are significantly higher than nitrate concentrations elsewhere. Tobit regression results demonstrate that the number of septic tanks and the percentages of forest, undeveloped, and high-density residential land within a 1000-foot radius of a well are reliable predictors of nitrate concentration in ground water. Similarly, logistic regression revealed that the percentages of forest, undeveloped, and low-density residential land are good indicators of ground water nitrate concentration > 2 mg/L. The methodology and results outlined here provide a useful tool for land managers in communities with shallow water tables overlain with highly permeable materials to evaluate potential effects of development on ground water quality.     

Hydrogeochemistry and environmental isotopes of ground water in Jeju volcanic island,Korea: implications for nitrate contamination

Ground water from springs and public supply wells was investigated for hydrochemistry and environmental isotopes of ³H, ¹⁸O and D in Jeju volcanic island, Korea. The wells are completed in a basaltic aquifer and the upper part of hydrovolcanic sedimentary formation. Nitrate contamination is conspicuous in the coastal area where most of the samples have nitrate concentrations well above 1 mg NO₃ N/l. Agricultural land use seems to have a strong influence on the distribution of nitrate in ground water. Comparison of stable isotopic compositions of precipitation and ground water show that ground water mostly originates from rainy season precipitation without significant secondary modification and that local recharge is dominant. ³H concentration of ground water ranged from nearly zero to 5 TU and is poorly correlated with vertical location of well screens. The occurrence of the ³H‐free, old ground water is due to the presence of low permeability layers near the boundary of the basaltic aquifer and the hydrovolcanic sedimentary formation, which significantly limits ground water flow from the upper basaltic aquifer. The old ground water exhibited background‐level nitrate concentrations despite high nitrate loadings, whereas young ground water had considerably higher nitrate concentrations. This correlation of ³H and nitrate concentration may be ascribed to the history of fertilizer use that has increased dramatically since the early 1960s in the island. This suggests that ³H can be used as a qualitative indicator for aquifer vulnerability to nitrate contamination. Copyright © 2005 John Wiley & Sons, Ltd.     

Spatio-temporal evolution and time-stable characteristics of soil moisture within remote sensing footprints with varying soil, slope, and vegetation

Air-borne passive microwave remote sensors measure soil moisture at the footprint scale, a scale of several hundred square meters or kilometers that encompasses different characteristic combinations of soil, topography, vegetation, and climate. Studies of within-footprint variability of soil moisture are needed to determine the factors governing hydrologic processes and their relative importance, as well as to test the efficacy of remote sensors. Gridded ground-based impedance probe water content data and aircraft-mounted Electronically Scanned Thinned Array Radiometer (ESTAR) pixel-average soil moisture data were used to investigate the spatio-temporal evolution and time-stable characteristics of soil moisture in three selected (LW03, LW13, LW21) footprints from the Southern Great Plains 1997 (SGP97) Hydrology Experiment. Better time-stable features were observed within a footprint containing sandy loam soil than within two pixels containing silty loam soil. Additionally, flat topography with split wheat/grass land cover produced the largest spatio-temporal variability and the least time stability in soil moisture patterns. A comparison of ground-based and remote sensing data showed that ESTAR footprint-average soil moisture was well calibrated for the LW03 pixel with sandy loam soil, rolling topography, and pasture land cover, but improved calibration is warranted for the LW13 (silty loam soil, rolling topography, pasture land) and LW21 (silty loam soil, flat topography, split vegetation of wheat and grass land with tillage practice) pixels. Footprint-scale variability and associated nonlinear soil moisture dynamics may prove to be critical in the regional-scale hydroclimatic models.     

Modelling soil moisture in a high‐latitude landscape using LiDAR and soil data

Soil moisture has a pronounced effect on earth surface processes. Global soil moisture is strongly driven by climate, whereas at finer scales, the role of non‐climatic drivers becomes more important. We provide insights into the significance of soil and land surface properties in landscape‐scale soil moisture variation by utilizing high‐resolution light detection and ranging (LiDAR) data and extensive field investigations. The data consist of 1200 study plots located in a high‐latitude landscape of mountain tundra in north‐western Finland. We measured the plots three times during growing season 2016 with a hand‐held time‐domain reflectometry sensor. To model soil moisture and its temporal variation, we used four statistical modelling methods: generalized linear models, generalized additive models, boosted regression trees, and random forests. The model fit of the soil moisture models were R² = 0.60 and root mean square error (RMSE) 8.04 VWC% on average, while the temporal variation models showed a lower fit of R² = 0.25 and RMSE 13.11 CV%. The predictive performances for the former were R² = 0.47 and RMSE 9.34 VWC%, and for the latter R² = 0.01 and RMSE 15.29 CV%. Results were similar across the modelling methods, demonstrating a consistent pattern. Soil moisture and its temporal variation showed strong heterogeneity over short distances; therefore, soil moisture modelling benefits from high‐resolution predictors, such as LiDAR based variables. In the soil moisture models, the strongest predictor was SAGA (System for Automated Geoscientific Analyses) wetness index (SWI), based on a 1 m² digital terrain model derived from LiDAR data, which outperformed soil predictors. Thus, our study supports the use of LiDAR based SWI in explaining fine‐scale soil moisture variation. In the temporal variation models, the strongest predictor was the field‐quantified organic layer depth variable. Our results show that spatial soil moisture predictions can be based on soil and land surface properties, yet the temporal models require further investigation. Copyright © 2017 John Wiley & Sons, Ltd.     

Hillslope scale surface runoff,sediment and nutrient losses associated with tramline wheelings

Research on arable sandy loam and silty clay loam soils on 4° slopes in England has shown that tramlines (i.e. the unseeded wheeling areas used to facilitate spraying operations in cereal crops) can represent the most important pathway for phosphorus and sediment loss from moderately sloping fields. Detailed monitoring over the October–March period in winters 2005–2006 and 2006–2007 included event‐based sampling of surface runoff, suspended and particulate sediment, and dissolved and particulate phosphorus from hillslope segments (each ～300–800 m²) established in a randomized block design with four replicates of each treatment at each of two sites on lighter and heavier soils. Experimental treatments assessed losses from the cropped area without tramlines, and from the uncropped tramline area, and were compared to losses from tramlines which had been disrupted once in the autumn with a shallow tine. On the lighter soil, the effects of removal or shallow incorporation of straw residues was also determined. Research on both sandy and silty clay loam soils across two winters showed that tramline wheelings represented the dominant pathway for surface runoff and transport of sediment, phosphorus and nitrogen from cereal crops on moderate slopes. Results indicated 5·5–15·8% of rainfall lost as runoff, and losses of 0·8–2·9 kg TP ha^?1 and 0·3–4·8 t ha^?1 sediment in tramline treatments, compared to only 0·2–1·7% rainfall lost as runoff, and losses of 0·0–0·2 kg TP ha^?1 and 0·003–0·3 t ha^?1 sediment from treatments without tramlines or those where tramlines had been disrupted. The novel shallow disruption of tramline wheelings using a tine once following the autumn spray operation consistently and dramatically reduced (p < 0·001) surface runoff and loads of sediment, total nitrogen and total phosphorus to levels similar to those measured in cropped areas between tramlines. Results suggest that options for managing tramline wheelings warrant further refinement and evaluation with a view to incorporating them into spatially‐targeted farm‐level management planning using national or catchment‐based agri‐environment policy instruments aimed at reducing diffuse pollution from land to surface water systems. Copyright © 2010 John Wiley & Sons, Ltd.     

The distribution of solute processes on an acid hillslope and the delivery of solutes to a stream: I. Exchangeable bases

This paper aims to identify the spatial distribution of exchangeable base cations in soils on an acid hillslope and to investigate possible cation release processes from slope soils to the stream. The basic assumption underlying this research is that the amount of exchangeable cations in soils reflects the nutrient stores and cation leaching processes across the slope where vegetation and parent materials are similar. The distribution of exchangeable Ca²⁺, Mg²⁺, K⁺ and Na⁺ has been investigated on a three-dimensional hillslope on the Quantock Hills, Somerset, UK. A two-way ANOVA shows that soil depth is predominant in explaining the total variance of exchangeable bases, despite the steep slope gradient and clear podzolic catena development. Major nutrient base cations, such as Ca²⁺, Mg²⁺ and K⁺, display homogeneous topsoil storage right across the slope. This spatial pattern may indicate that the spatial distribution of major nutrient cations is tightly controlled by the soil–vegetation system in nutrient-poor heathland environments. Na⁺ is an exception to this vegetation-controlled spatial distribution, because of its small involvement in the soil–vegetation and soil exchangeable systems. In subsurface soils, cations liberated from the soil–vegetation system are subject to redistribution over the slope according to the hydrological flowpaths operating on the slope, with some eventually released into the stream. The saturated wedge developed at the base of the slope plays a key role in the storage and release processes of base cations from slope soils to the stream. Ca²⁺, Mg²⁺ and Na⁺ carried by throughflow are stored in the saturated wedge and gradually released into the stream at times of high flow. K⁺, however, shows an apparently different spatial behaviour, being deficient in the saturated wedge. Copyright © 1999 John Wiley & Sons, Ltd.     

基于非线性剪切梁方法的堤防液化判别分析研究

很多河流的堤防工程修建在人口密集地区,失事后往往会造成严重后果,尤其是细粒土填筑而成的堤防,在地震作用下存在液化破坏的可能。采用合理方法对堤防进行液化判别,确定其在地震作用下的动力稳定性十分必要。针对某实际堤防工程,进行五种拟用筑堤土料动力特性试验研究,并基于试验成果,采用能够考虑地基-堤身相互作用的非线性剪切梁方法,对堤防的液化可能性进行判别。研究表明:提出的非线性剪切梁方法能够考虑地基和堤身间的相互作用,且简单实用、可操作性强;土体的相对密度、空隙比和颗粒粒径大小及其他粒分布特性对其抗地震液化能力影响较大;在设定运行工况下,Ⅷ度地震时除拟选粉质壤土外,其他四种土筑堤均有发生液化的可能,选用粉质壤土筑堤较合适。     

Introducing the recharge–discharge relationship to evaluate the recharge coefficient of karstic aquifers: thematic and indexing approaches

ABSTRACT

Karst aquifers and springs are important with respect to their potential for supplying drinking water in regions suffering from water scarcity in Iran. Accordingly, it is essential to determine the recharge potential of the catchment and the regions with higher obtainability potential. This study provides a road map for the Sheshpeer catchment in southern Iran. A recharge potential (RP) map was produced from which a recharge index (RI) was computed for several selected springs in the catchment. Furthermore, the unit discharge (q) – defined as the average annual discharge for a given catchment area and unit rainfall depth for each spring – was calculated. The plot of q versus RI for the springs showed a linear positive relationship between the two variables (R ² = 0.9). Applying the trend equation of this plot to the whole Sheshpeer karstic catchment reveals that its long-term recharge coefficient is 0.74.     

Electromagnetic induction prediction of soil salinity and groundwater properties in a Tunisian Saharan oasis

Abstract

Electromagnetic induction measurements (EM) were taken in a saline gypsiferous soil of the Saharan-climate Fatnassa oasis (Tunisia) to predict the electrical conductivity of saturated soil extract (ECe) and shallow groundwater properties (depth, Dgw, and electrical conductivity, ECgw) using various models. The soil profile was sampled at 0.2 m depth intervals to 1.2 m for physical and chemical analysis. The best input to predict the log-transformed soil salinity (lnECe) in surface (0–0.2 m) soil was the EMh/EMv ratio. For the 0–0.6 m soil depth interval, the performance of multiple linear regression (MLR) models to predict lnECe was weaker using data collected over various seasons and years (R _a ² = 0.66 and MSE = 0.083 dS m^-1) as compared to those collected during the same period (R _a ² = 0.97, MSE = 0.007 dS m^-1). For similar seasonal conditions, for the Dgw–EMv relationship, R ² was 0.88 and the MSE was 0.02 m for Dgw prediction. For a validation subset, the R ² was 0.85 and the MSE was 0.03 m. Soil salinity was predicted more accurately when groundwater properties were used instead of soil moisture with EM variables as input in the MLR.

Editor D. Koutsoyiannis; Associate editor K. Heal

Citation Bouksila, F., Persson, M., Bahri, A., and Berndtsson, R., 2012. Electromagnetic induction predictions of soil salinity and groundwater properties in a Tunisian Saharan oasis. Hydrological Sciences Journal, 57 (7), 1473–1486.     

Spatial distribution and content of soil organic matter in an agricultural field in eastern Canada,as estimated from geostatistical tools

Soil erosion induces soil redistribution within the landscape and thus contributes to the spatial variability of soil quality. This study complements a previous experimentation initiated by the authors focusing on soil redistribution as a result of soil erosion, as indicated by caesium‐137 (¹³⁷Cs) measurements, in a small agricultural field in Canada. The spatial variability of soil organic matter (SOM) was characterized using geostatistics, which consider the randomized and structured nature of spatial variables and the spatial distribution of the samples. The spatial correlation of SOM (in percentages) patterns in the topsoil was established taking into account the spatial structure present in the data. A significant autocorrelation and reliable variograms were found with a R² ≥ 0·9, thus demonstrating a strong spatial dependence. Ordinary Kriging (OK) interpolation provided the best cross validation (r² = 0·35). OK and inverse distance weighting power two (IDW2) interpolation approaches produced similar estimates of the total SOM content of the topsoil (0–20 cm) of the experimental field, i.e. 211 and 213 tonnes, respectively. However, the two approaches produced differences in the spatial distribution patterns and the relative magnitude of some SOM content classes. The spatialization of SOM and soil redistribution variability – as evidenced by ¹³⁷Cs measurements – is a first step towards the assessment of the impact of soil erosion on SOM losses to recommend conservation measures. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulations of water movement and solute transport through different soil texture configurations under negative‐pressure irrigation

This study examined the effects of different soil texture configurations on water movement and solute transport to provide a reliable scientific basis for the application of negative‐pressure irrigation (NPI) technology. HYDRUS‐2D was used to analyse water movement and solute transport under NPI. The main results are as follows: (a) HYDRUS‐2D can be used to simulate water movement and solute transport under NPI, as there was good agreement between the simulated and measured values for water contents, NaCl concentrations, cumulative water infiltration, and wetting distances in the horizontal and vertical directions; the Nash–Sutcliffe efficiency coefficients were in the range of 0.94–0.97. (b) Layered soils have obvious effects on water movement under NPI. With the emitter position in the loam layer, when a coarse texture of loamy sand was present below the loam layer (namely, L‐LS), irrigation water accumulated in the topsoil, and this led to an increase in evaporation compared with the homogeneous loam profile. However, fine texture silty loam or silty clay loam layers beneath the loam layer (namely, L‐SiL or L‐SiCL, respectively) was more conducive to water infiltration into the lower layer, and this increased the amount of water infiltration and simultaneously reduced the surface evaporation effectively. (c) Layered soils have obvious effects on solute transport under NPI, and salt accumulation will readily occur in the clay‐rich soil layer at the interface. The maximum soil salt accumulation of L‐LS occurred above the soil interface between the two soil layers with a value of 21.80 g/kg; however, for L‐SiCL and L‐SiL, the maximum salt accumulation occurred below the soil interface between the two soil layers, with values of 23.80 g/kg and 20.08 g/kg, respectively. (d) Interlayered soils showed remarkable changes in the water infiltration characteristics and salt‐leaching intensities under NPI, and the properties for the soil profile with a silty loam interlayer were better than those for the soil profile with a silty clay loam interlayer. The soil profile with a loamy sand interlayer had the lowest amount of water infiltration, which resulted in reductions of the salt‐leaching intensities. Thus, NPI is clearly not suitable for loamy sand soil. Overall, the results demonstrated that soil texture configurations affected water movement and solute transport under NPI. Therefore, careful consideration should be given to the use of NPI to achieve target soil water and solution conditions and reduce water loss.     

ESTIMATING INFILTRATION PARAMETERS FROM BASIC SOIL PROPERTIES

Infiltration data were collected on two rectangular grids with 25 sampling points each. Both experimental grids were located in tropical rain forest (Guyana), the first in an Arenosol area and the second in a Ferralsol field. Four different infiltration models were evaluated based on their performance in describing the infiltration data. The model parameters were estimated using non-linear optimization techniques. The infiltration behaviour in the Ferralsol was equally well described by the equations of Philip, Green–Ampt, Kostiakov and Horton. For the Arenosol, the equations of Philip, Green–Ampt and Horton were significantly better than the Kostiakov model. Basic soil properties such as textural composition (percentage sand, silt and clay), organic carbon content, dry bulk density, porosity, initial soil water content and root content were also determined for each sampling point of the two grids. The fitted infiltration parameters were then estimated based on other soil properties using multiple regression. Prior to the regression analysis, all predictor variables were transformed to normality. The regression analysis was performed using two information levels. The first information level contained only three texture fractions for the Ferralsol (sand, silt and clay) and four fractions for the Arenosol (coarse, medium and fine sand, and silt and clay). At the first information level the regression models explained up to 60% of the variability of some of the infiltration parameters for the Ferralsol field plot. At the second information level the complete textural analysis was used (nine fractions for the Ferralsol and six for the Arenosol). At the second information level a principal components analysis (PCA) was performed prior to the regression analysis to overcome the problem of multicollinearity among the predictor variables. Regression analysis was then carried out using the orthogonally transformed soil properties as the independent variables. Results for the Ferralsol data show that the parameters of the Green–Ampt and Kostiakov model were estimated relatively accurately (maximum R² = 0.76). For the Arenosol, use of the second information level together with PCA produced regression models with an R² value ranging from 0.38 to 0.68. For the Ferralsol, most of the variance was explained by the root content and organic matter content. In the Arenosol plot, the fractions medium and fine sand explained most of the observed variance.     

Integration of remotely sensed C factor into SWAT for modelling sediment yield

The C factor, representing the impact of plant and ground cover on soil loss, is one of the important factors of the Modified Universal Soil Loss Equation (MUSLE) in the Soil and Water Assessment Tool (SWAT) to model sediment yield. The daily update of C factors in SWAT was originally determined by land use types and plant growth cycles. This does not reflect the spatial variation of C values that exists within a large land use area. We present a new approach to integrate remotely sensed C factors into SWAT for highlighting the effect of detailed vegetative cover data on soil erosion and sediment yield. First, the C factor was estimated using the abundance of ground components extracted from remote sensing images. Then, the gridding data of the C factor were aggregated to hydrological response units (HRUs), instead of to land use units of SWAT. In the end, the C factor values in HRUs were integrated into SWAT to predict sediment yield by modifying the ysed subroutine. This substitution work not only increases the spatial variation of the C factor in SWAT, but also makes it possible to utilize other sources of C databases rather than those from the United States. The demonstration in the Dage basin shows that the modified SWAT produces reasonable results in water flow simulation and sediment yield prediction using remotely sensed C values. The Nash–Sutcliffe efficiency coefficient (E_NS) and R² for surface runoff range from 0·69 to 0·77 and 0·73 to 0·87, respectively. The coefficients E_NS and R² for sediment yield were generally above 0·70 and 0·60, respectively. The soil erosion risk map based on sediment yield prediction at the HRU level illustrates instructive details on spatial distribution of soil loss. Copyright © 2011 John Wiley & Sons, Ltd.     

Aqueous fluids derived from a subducting slab: Observed high 3He emanation and conductive anomaly in a non-volcanic region,Kii Peninsula southwest Japan

It has long been recognized that the Kii Peninsula in the southwest Japan arc is peculiar in a non-volcanic region, indicated by the presence of high temperature hot springs, high terrestrial heat flow and high ³He content in hot spring gases. Geophysical and geochemical studies were carried out to understand the geotectonic environment in the southern part of the Kii Peninsula. Most of the measured ³He / ⁴He ratios are similar or higher than air, indicating wide spread incorporation of mantle-derived helium into meteoric water. A region with rather high ³He / ⁴He ratios (> 4 R_A) on the west side of the Omine Mountains coincides with the occurrence of high temperature hot springs. A deep crustal resistivity structure across the Omine Mountains was imaged by wide-band magnetotelluric soundings. A 2-D inversion with N–S strike using both TM and TE modes reveals two conductors, one in the upper (3–7 km depth) and the other in lower crust (25–35 km depth) to the west of the Omine Mountains. The distribution of microearthquakes and low-frequency tremors, and the existence of seismic reflectors indicate that the large conductor in the lower crust is related to aqueous fluids derived from the Philippine Sea plate. The upper-crustal conductive zone may also reflect the aqueous fluids trapped in the upper crust, which are presumably derived from the subducting slab. Considering the occurrence of seismic events in the subducting slab beneath the southern Kii Peninsula, the aqueous fluids generated by dehydration of the slab mantle could plausibly include MORB-type helium derived from the residual lithospheric mantle. Therefore, the high temperature hot springs and high ³He emanations in hot spring gases and other geotectonic events in the southern Kii Peninsula may be due to heat flux and mantle-derived helium discharged from aqueous fluid in the upper crust.     

An analytical solution for rapidly predicting post-fire peak streamflow for small watersheds in southern California

Following wildfires, the probability of flooding and debris flows increase, posing risks to human lives, downstream communities, infrastructure, and ecosystems. In southern California (USA), the Rowe, Countryman, and Storey (RCS) 1949 methodology is an empirical method that is used to rapidly estimate post-fire peak streamflow. We re-evaluated the accuracy of RCS for 33 watersheds under current conditions. Pre-fire peak streamflow prediction performance was low, where the average R² was 0.29 and average RMSE was 1.10 cms/km² for the 2- and 10-year recurrence interval events, respectively. Post-fire, RCS performance was also low, with an average R² of 0.26 and RMSE of 15.77 cms/km² for the 2- and 10-year events. We demonstrated that RCS overgeneralizes watershed processes and does not adequately represent the spatial and temporal variability in systems affected by wildfire and extreme weather events and often underpredicted peak streamflow without sediment bulking factors. A novel application of machine learning was used to identify critical watershed characteristics including local physiography, land cover, geology, slope, aspect, rainfall intensity, and soil burn severity, resulting in two random forest models with 45 and five parameters (RF-45 and RF-5, respectively) to predict post-fire peak streamflow. RF-45 and RF-5 performed better than the RCS method; however, they demonstrated the importance and reliance on data availability. The important parameters identified by the machine learning techniques were used to create a three-dimensional polynomial function to calculate post-fire peak streamflow in small catchments in southern California during the first year after fire (R² = 0.82; RMSE = 6.59 cms/km²) which can be used as an interim tool by post-fire risk assessment teams. We conclude that a significant increase in data collection of high temporal and spatial resolution rainfall intensity, streamflow, and sediment loading in channels will help to guide future model development to quantify post-fire flood risk.     

Seasonal and event‐based drivers of runoff and phosphorus export through agricultural tile drains under sandy loam soil in a cool temperate region

Frequent algal blooms in surface water bodies caused by nutrient loading from agricultural lands are an ongoing problem in many regions globally. Tile drains beneath poorly and imperfectly drained agricultural soils have been identified as key pathways for phosphorus (P) transport. Two tile drains in an agricultural field with sandy loam soil in southern Ontario, Canada were monitored over a 28‐month period to quantify discharge and the concentrations and loads of dissolved reactive P (DRP) and total P (TP) in their effluent. This paper characterizes seasonal differences in runoff generation and P export in tile drain effluent and relates hydrologic and biogeochemical responses to precipitation inputs and antecedent soil moisture conditions. The generation of runoff in tile drains was only observed above a clear threshold soil moisture content (~0.49 m³·m^?3 in the top 10 cm of the soil; above field capacity and close to saturation), indicating that tile discharge responses to precipitation inputs were governed by the available soil‐water storage capacity of the soil. Soil moisture content approached this threshold throughout the non‐growing season (October – April), leading to runoff responses to most events. Concentrations of P in effluent were variable throughout the study but were not correlated with discharge (p > 0.05). However, there were significant relationships between discharge volume (mm) and DRP and TP loads (kg ha^?1) for events occurring over the study period (R² ≥ 0.49, p ≤ 0.001). This research has shown that the hydrologic and biogeochemical responses of tile drains in a sandy loam soil can be predicted to within an order of magnitude from simple hydrometric data such as precipitation and soil moisture once baseline conditions at a site have been determined. Copyright © 2016 John Wiley & Sons, Ltd.