Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of biochar and bio-fertilizer on yield and qualitative properties of soybean and some chemical properties of soil

This study attempted to examine the effects of biochar amendment together with bio-fertilizer on soybean yield and its qualitative properties, as well as a few chemical properties of soil through a factorial randomized complete block design at three replications in east of Golestan Province (Iran) during 2014. The two factors under study included the following: (1) biochar amendment (in four levels of 0, 2.5, 8, and 16 tons per hectare), (2) bio-fertilizer containing phosphorus and sulfur growth-promoting rhizobacteria (in two levels of inoculation and non-inoculation) applied through foliar feeding. The results of analysis of variance indicated that interactions of biochar amendment and bio-fertilizer on harvest index and grain yield were significant (p ≤?0.01). According to the results of this study, the highest harvest index and oil content were 56.9, and 17.7%, respectively, in the treatment of 8 tons per hectare biochar and inoculation with bio-fertilizer. The lowest harvest index and the lowest oil content were in the control treatment. The interaction of biochar and bio-fertilizer on bulk density and cation exchange capacity was significant (p ≤?0.01). The results of this study demonstrated that biochar affected the amount of residual nitrogen in the soil after harvest, cation exchange capacity (CEC), acidity (pH), and electrical conductivity (EC). The highest grain yield (3440 kg/ha) was in the 8-ton biochar treatment with inoculated bio-fertilizer. Our study concludes that the biochar and bio-fertilizers can improve grain yield of soybean till 51% relative to the control.     

Developing a climate model for Iran using GIS

Summary In order to develop a climate model for Iran, monthly mean climatic variables from 117 synoptic stations were obtained from the Iranian Meteorological Organisation. These variables were reduced to six orthogonal factors using factor analysis. The stations were then divided into six groups using cluster analysis. Within each climatic group, the lowest and highest thresholds for each factor were identified. The factor scores of the stations within each factor were interpolated across the country applying Inverse Squared Distance Weight in the ArcGIS environment. Based on the factor scores, six conditional functions were defined to allocate each pixel to a region. In order to simplify the models, one index variable was substituted for each factor. Then, through Discriminant Analysis, the constants and coefficients of the models were determined. The final models were evaluated against some examples, one of which, Yazd, was demonstrated fully. Authors’ address: Bohloul Alijani, Manijeh Ghohroudi, Nahid Arabi, Department of Geography, Teacher Training University, Mofetteh Avenue, Tehran, Iran.     

The Selloum Basin: new element of the Middle Liassic paleogeography in the southern Middle Atlas (Morocco)

A comparative study of the Middle Liassic series in the south western part of the Middle Atlas system leads to the definition of three new formations called Fellat 1, 2 and 3. They mainly involve gravity-flow deposits related to the breakup of the regional carbonate platform in a subsiding area, the Selloum Basin. The creation of this basin provides evidence for a rifting stage during the Carixian–Domerian interval, inducing the separation of the Selloum unit with respect to the Middle Atlasic trough and its possible connection with the High Atlasic sedimentary domain.     

Impact of rehabilitation of Assiut barrage, Nile River, on groundwater rise in urban areas

To make optimum use of the most vital natural resource of Egypt, the River Nile water, a number of regulating structures (in the form of dams and barrages) for control and diversion of the river flow have been constructed in this river since the start of the 20th century. One of these barrages is the Assiut barrage which will require considerable repairs in the near future. The design of the rehabilitation of the barrage includes a headpond with water levels maintained at a level approximately 0.60 m higher than the highest water level in the headpond of the present barrage. This development will cause an increase of the seepage flow from the river towards the adjacent agricultural lands, Assiut Town and villages. The increased head pond level might cause a rise of the groundwater levels and impedance of drainage outflows. The drainage conditions may therefore be adversely affected in the so-called impacted areas which comprise floodplains on both sides of the Nile for about 70 km upstream of the future barrage. A rise in the groundwater table, particularly when high river levels impede drainage, may result in waterlogging and secondary salinization of the soil profile in agricultural areas and increase of groundwater into cellars beneath buildings in the urban areas. In addition, a rise in the groundwater table could have negative impact on existing sanitation facilities, in particular in the areas which are served with septic tanks. The impacts of increasing the headpond level were assessed using a three-dimensional groundwater model. The mechanisms of interactions between the Nile River and the underlying Quaternary aquifer system as they affect the recharge/discharge processes are comprehensively outlined. The model has been calibrated for steady state and transient conditions against historical data from observation wells. The mitigation measures for the groundwater rise in the urban areas have been tested using the calibrated mode.     

Global sensitivity analysis of hydrologic processes in major snow‐dominated mountainous river basins in Colorado

The performance of watershed models in simulating stream discharge depends on the adequate representation of important watershed processes. In snow‐dominated systems, snow, surface and subsurface hydrologic processes comprise a complex network of nonlinear interactions that influence the magnitude and timing of discharge. This study aims to identify critical processes and interactions that control discharge hydrographs in five major mountainous snow‐dominated river basins in Colorado, USA. A comprehensive watershed model (Soil and Water Assessment Tool) and a variance‐based global sensitivity analysis technique (Fourier Amplitude Sensitivity Test) were used in conjunction to identify critical models parameters and processes that they represent. Average monthly streamflow and streamflow root mean square error over a period of 20 years were used as two separate objective functions in this analysis. Examination of the sensitivity of monthly streamflow revealed the influence of parameters on flow volume, whereas the sensitivity of streamflow root mean square error also exposed the influence of parameters on the timing of the hydrographs. A stability analysis was performed to investigate the computational requirements for a robust sensitivity analysis. Results show that streamflow volume is mostly influenced by shallow subsurface processes, whereas interactions between groundwater and snow processes were the key in the timing of streamflows. A large majority of important parameters were common among all study watersheds, which underlies the prospect for regionalization of process‐based hydrologic modelling in headwater river basins in Colorado. Copyright © 2013 John Wiley & Sons, Ltd.     

Practical Use of Computationally Frugal Model Analysis Methods

Three challenges compromise the utility of mathematical models of groundwater and other environmental systems: (1) a dizzying array of model analysis methods and metrics make it difficult to compare evaluations of model adequacy, sensitivity, and uncertainty; (2) the high computational demands of many popular model analysis methods (requiring 1000's, 10,000 s, or more model runs) make them difficult to apply to complex models; and (3) many models are plagued by unrealistic nonlinearities arising from the numerical model formulation and implementation. This study proposes a strategy to address these challenges through a careful combination of model analysis and implementation methods. In this strategy, computationally frugal model analysis methods (often requiring a few dozen parallelizable model runs) play a major role, and computationally demanding methods are used for problems where (relatively) inexpensive diagnostics suggest the frugal methods are unreliable. We also argue in favor of detecting and, where possible, eliminating unrealistic model nonlinearities—this increases the realism of the model itself and facilitates the application of frugal methods. Literature examples are used to demonstrate the use of frugal methods and associated diagnostics. We suggest that the strategy proposed in this paper would allow the environmental sciences community to achieve greater transparency and falsifiability of environmental models, and obtain greater scientific insight from ongoing and future modeling efforts.     

Representation of agricultural conservation practices with SWAT

Results of modelling studies for the evaluation of water quality impacts of agricultural conservation practices depend heavily on the numerical procedure used to represent the practices. Herein, a method for the representation of several agricultural conservation practices with the Soil and Water Assessment Tool (SWAT) is developed and evaluated. The representation procedure entails identifying hydrologic and water quality processes that are affected by practice implementation, selecting SWAT parameters that represent the affected processes, performing a sensitivity analysis to ascertain the sensitivity of model outputs to selected parameters, adjusting the selected parameters based on the function of conservation practices, and verifying the reasonableness of the SWAT results. This representation procedure is demonstrated for a case study of a small agricultural watershed in Indiana in the Midwestern USA. The methods developed in the present work can be applied with other watershed models that employ similar underlying equations to represent hydrologic and water quality processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of climate change on diffuse pollutant fluxes at the watershed scale

This study aims to assess watershed‐scale impacts of changing climate on sediment, phosphorus, nitrogen and pesticide (atrazine) fluxes over the 21st century at the watershed scale. In particular, changes in dissolved and particulate forms of water quality constituents in response to climate change are investigated. The hydrologic model Soil and Water Assessment Tool was calibrated and evaluated in a primarily agricultural watershed in the Midwestern United States to simulate hydrologic and water quality processes on a daily basis over the 2015–2099 time horizon. The model was then driven with 112 distinct statistically downscaled climate projections representing Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC SRES) low, moderate and high greenhouse gas emission scenarios. Projected hydrologic and water quality responses were categorized according to the three IPCC SRES emission scenarios for summarizing and synthesizing results over early‐century (2015–2034), mid‐century (2045–2064) and late‐century (2080–2099) assessment. Results revealed clear warming trends in the study area, whereas small increases in precipitation were predicted. Streamflow, sediment and total nutrient loads did not differ noticeably between assessment periods. However, the proportion of dissolved to total nutrients increased significantly from early‐century to late‐century periods. With the exception of total atrazine in the mid‐century period, predicted pollutant loads for a given assessment period did not differ between emission pathways for a given assessment period. Changes in pollutant fluxes showed pronounced monthly variability. The projected increase in readily available forms of nutrients has important implications for the ecological health of water systems and management of drinking water supplies. Copyright © 2013 John Wiley & Sons, Ltd.     

Combination of Precise Leveling and InSAR Data to Constrain Source Parameters of the Mw = 6.5, 26 December 2003 Bam Earthquake

We used new precise leveling data acquired 40 days after the Bam earthquake in combination with radar interferometry observations from both ascending and descending orbits to investigate static deformation associated with the 2003 Bam earthquake. We invert this geodetic data set to gain insight into the fault geometry and slip distribution of the rupture. The best-fitting dislocation model is a steeply east-dipping right-lateral strike-slip fault that has a size of 11 by 8 km and strikes N2°W. We find that such smooth geometry fits available geodetic data better than previously proposed models for this earthquake. Our distributed slip model indicates a maximum strike slip of 3 m occurring about 3 to 5 km deep. The slip magnitude and depth of faulting taper to the north, where the fault approaches the Bam city. Inclusion of crustal layering increases the amount of maximum slip inferred at depth by about 4%.