Impact of lateral flow and spatial scaling on the simulation of semi‐arid urban land surfaces in an integrated hydrologic and land surface model

Understanding and representing hydrologic fluxes in the urban environment is challenging because of fine scale land cover heterogeneity and lack of coherent scaling relationships. Here, the impact of urban land cover heterogeneity, scale, and configuration on the hydrologic and surface energy budget (SEB) is assessed using an integrated, coupled land surface/hydrologic model at high spatial resolutions. Archetypes of urban land cover are simulated at varying resolutions using both the National Land Cover Database (NLCD; 30 m) and an ultra high‐resolution land cover dataset (0.6 m). The analysis shows that the impact of highly organized, yet heterogeneous, land cover typical of the urban domain can cause large variations in hydrologic and energy fluxes within areas of similar land cover. The lateral flow processes that occur within each simulation create variations in overland flow of up to ±200% and ±4% in evapotranspiration. The impact on the SEB is smaller and largely restricted to the wet season for our semi‐arid forcing scenarios. Finally, we find that this seasonal bias, predominantly caused by lateral flow, is displaced by a systematic diurnal bias at coarser resolutions caused by deficiencies in the method used for scaling of land surface and hydrologic parameters. As a result of this research, we have produced land surface parameters for the widely used NLCD urban land cover types. This work illustrates the impact of processes that remain unrepresented in traditional high‐resolutions land surface models and how they may affect results and uncertainty in modeling of local water resources and climate. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of two modeling approaches for groundwater–surface water interactions

An assessment of interactions between groundwater and surface water was carried out by applying two different modeling approaches to a small‐scale study area in the municipality of Havelock, Quebec. The first approach involved a commonly used sequential procedure that consists in determining the daily recharge rate using a quasi 2D infiltration model (HELP), applied in the next step as an imposed flux to a 3D finite‐element groundwater flow model. The flow model was calibrated under steady‐state and transient conditions against measured water levels. The second approach was based on a recently developed physically based, 3D fully coupled groundwater–surface water flow model (CATHY) applied to the entire flow domain in an integrated manner. Implementation, calibration, and results of the simulations for both approaches are presented and discussed. For equal annual precipitation (1038 mm/y) and evapotranspiration (556 mm/y), the second approach computed a recharge rate of 233 mm/y (8.9% higher than the first approach) and a net upward flow from the fractured aquifer (the first approach predicted a net downward flow to the rock). The simulated annual discharge was similar for the two approaches (9.6% difference). Both approaches were found to be useful in understanding the interactions between groundwater and surface water, although limitations are apparent in the sequential procedure's inability to account for surface–subsurface feedbacks, for instance near stream reaches where groundwater discharge is prevalent. The decoupled, two‐model approach provides disaggregated surface, vadose, and aquifer flows, and a simple aperçu at the different components of total discharge. The fully coupled model accounts for continuous water exchanges between the land surface, subsurface, and stream channel in a more complex manner, and produces a better match against observed data. Copyright © 2012 John Wiley & Sons, Ltd.     

Stochastic–conceptual analysis of near‐surface hydrological response

The stochastic–conceptual rainfall–runoff simulator (SCRRS) developed by R. A. Freeze in 1980 was used in this study to demonstrate quantitatively the interplay of the factors that control the occurrence of overland flow by the Horton and Dunne mechanisms. The simulation domain and input data for the SCRRS simulations reported here were abstracted from the R‐5 catchment (Chickasha, OK) data sets. The results illustrate that the identification of a dominant hydrological response process may not be as simple as a singular Horton or Dunne characterization. The SCRRS simulations show that the Horton and Dunne processes can (i) occur simultaneously at different locations during a given rainfall event, (ii) change from one process to the other with time depending on the characteristics of the rainfall event, and (iii) be strongly dependent on the initial conditions. Copyright © 2001 John Wiley & Sons, Ltd.     

Recent studies on surface water–groundwater relationships at hydro‐projects in India using environmental isotopes

Environmental isotopes (²H or D, ¹⁸O, ³H), along with geology, hydrochemistry and in situ physicochemical parameters (EC, T, DO, pH) were employed to study surface water (reservoir, lake)–groundwater (spring) relationships at (1) Nagewadi, a minor irrigation project in the State of Maharashtra, Western India; (2) Kanhirapuzha reservoir in the State of Kerala, Southern India and (3) Ghatghar Pumped Storage Hydroelectric Project in the State of Maharashtra, Western India for the purpose of understanding the seepage/leakage and its associated problems. The studies concluded that the springs found downstream of the Nagewadi project originate from the reservoir and not from the abutments or shallow aquifers. The Kanhirapuzha reservoir receives a substantial base‐flow component compared to riverine inputs. The reason for the water‐logging problem at a nearby downstream village during the non‐summer periods is due to the change in the upstream groundwater flow direction under reservoir filling conditions and is not due to reservoir leakage. Most of the springs in the approach tunnel to the underground power house of the Ghatghar Project originate from the lower reservoir and not from the upper reservoir or the overburden rock matrix. The above case studies illustrate the diversity of environmental isotope applications in surface water (reservoir, lake)–groundwater (spring) relationships related to sustainability of hydro‐projects. Copyright © 2008 John Wiley & Sons, Ltd.     

Field assessment of surface water–groundwater connectivity in a semi‐arid river basin (Murray–Darling,Australia)

In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater–surface‐water interactions in a braided river: a tracer‐based assessment

Natural tracers (alkalinity and silica) were used to infer groundwater–surface‐water exchanges in the main braided reach of the River Feshie, Cairngorms, Scotland. Stream‐water samples were collected upstream and downstream of the braided section at fortnightly intervals throughout the 2001–2002 hydrological year and subsequently at finer resolution over two rainfall events. The braided reach was found to exert a significant downstream buffering effect on the alkalinity of these waters, particularly at moderate flows (4–8 m³ s^?1/?Q_30–70). Extensive hydrochemical surveys were undertaken to characterize the different source waters feeding the braids. Shallow groundwater flow systems at the edge of the braided floodplain, recharged by effluent streams and hillslope drainage, appeared to be of particular significance. Deeper groundwater was identified closer to the main channel, upwelling through the hyporheic zone. Both sources contributed to the significant groundwater–surface‐water interactions that promote the buffering effect observed through the braided reach. Their impact was less significant at higher flows (>15 m³ s^?1/>Q₁₀) when acidic storm runoff from the peat‐covered catchment headwaters dominated, as well as under baseflow conditions (<4 m³ s^?1/<Q₇₀), when upstream alkalinity was already buffered owing to headwater groundwater sources assuming dominance. The significant temporally and spatially dynamic influence of these groundwater–surface‐water interactions was therefore seen to have important implications for both catchment functioning and instream ecology. Copyright © 2004 John Wiley & Sons, Ltd.     

Heterogeneity in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream

Spatial and temporal variability in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream was investigated. Four locations in a 150‐m reach of the stream were studied using hydrometric and hydrochemical tracing techniques. A high degree of hydrological connectivity between the riparian hillslope and the stream channel was indicated at two locations, where hydrochemical changes and hydraulic gradients indicated that the hyporheic zone was dominated by upwelling ground water. The chemistry of ground water reflected relatively long residence times and reducing conditions with high levels of alkalinity and conductivity, low dissolved oxygen (DO) and nitrate. At the other locations, connectivity was less evident and, at most times, the hyporheic zone was dominated by downwelling stream water characterized by high DO, low alkalinity and conductivity. Substantial variability in hyporheic chemistry was evident at fine (<10 m) spatial scales and changed rapidly over the course of hydrological events. The nature of the hydrochemical response varied among locations depending on the strength of local ground water influence. It is suggested that greater emphasis on spatial and temporal heterogeneity in ground water–surface water interactions in the hyporheic zone is necessary for a consideration of hydrochemical effects on many aspects of stream ecology. For example, the survival of salmonid eggs in hyporheic gravels varied considerably among the locations studied and was shown to be associated with variation in interstitial chemistry. River restoration schemes and watershed management strategies based only on the surface expression of catchment characteristics risk excluding consideration of potentially critical subsurface processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a lowland watershed: source contribution to stream flow

The lower coastal plain of the Southeast USA is undergoing rapid urbanisation as a result of population growth. Land use change has been shown to affect watershed hydrology by altering stream flow and, ultimately, impairing water quality and ecologic health. However, because few long‐term studies have focused on groundwater–surface water interactions in lowland watersheds, it is difficult to establish what the effect of development might be in the coastal plain region. The objective of this study was to use an innovative improvement to end‐member mixing analysis (EMMA) to identify time sequences of hydrologic processes affecting storm flow. Hydrologic and major ion chemical data from groundwater, soil water, precipitation and stream sites were collected over a 2‐year period at a watershed located in USDA Forest Service's Santee Experimental Forest near Charleston, South Carolina, USA. Stream flow was ephemeral and highly dependent on evapotranspiration rates and rainfall amount and intensity. Hydrograph separation for a series of storm events using EMMA allowed us to identify precipitation, riparian groundwater and streambed groundwater as main sources to stream flow, although source contribution varied as a function of antecedent soil moisture condition. Precipitation, as runoff, dominated stream flow during all storm events while riparian and streambed groundwater contributions varied and were mainly dependent on antecedent soil moisture condition. Sensitivity analyses examined the influence of 10% and 50% increases in analyte concentration on EMMA calculations and found that contribution estimates were very sensitive to changes in chemistry. This study has implications on the type of methodology used in traditional forms of EMMA research, particularly in the recognition and use of median end‐member water chemistry in hydrograph separation techniques. Potential effects of urban development on important hydrologic processes (groundwater recharge, interflow, runoff, etc.) that influence stream flow in these lowland watersheds were qualitatively examined. Copyright © 2011 John Wiley & Sons, Ltd.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

A fundamental dual‐zone continuum theory for dynamic soil–structure interaction

A continuum theory for an improved characterization of dynamic soil–structure interaction in the framework of three‐dimensional elastodynamics is presented. Effective in demonstrating the importance of integrating free‐field and near‐field effects under general soil and foundation conditions, a compact two‐zone delineation of the soil medium is proposed as a quintessential mechanics perspective for this class of problems. Sufficient to deliver a practical resolution of some perennial analytical and experimental conflicts, a fundamental formulation commensurate to a gradated unification of the homogenization approach and any sole free‐field inhomogeneous representation is developed and implemented computationally. Specialized to the problem of a rigid circular footing on sand, a nominal set of dynamic contact stress distributions and related impedance functions by the dual‐zone theory is included for theoretical and engineering evaluation. Through its comparison with benchmark analytical solutions and relevant physical measurements, the usage of the underlying conceptual platform as an advanced yet practical foundation for general dynamic soil–structure interaction is illustrated. Copyright © 2010 John Wiley & Sons, Ltd.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluation of sensitivity of land surface hydrology representations with and without land–atmosphere feedback

This paper compares various ways of quantifying the importance of land–atmosphere feedback. A widely used land surface hydrology model is used in coupled (to a planetary boundary layer model) and uncoupled modes to compare the adequacy of different feedback indices. It is found that existing feedback indices are primarily based on ‘one factor at a time’ sensitivity analysis and cannot adequately capture the interaction between land and atmosphere. A new index is used which combines factorial design concepts and traditional sensitivity analysis. This index is shown to capture and quantify the strength of interaction between land surface parameters and atmosphere. To assess the effects of forcing characteristics on the stand alone model sensitivity, several ways to specify near-surface atmospheric conditions are evaluated. It is found that commonly used forcing conditions (e.g. model generated or observed time-series of near-surface atmospheric variables) may not be adequate to mimic the coupled model environment for evaluating the land surface representations. The partially coupled model sensitivity is shown to capture a major feedback loop related to water holding capacity, surface fluxes and near-surface atmospheric processes. These results suggest that sensitivity from the stand alone model should be interpreted with caution and future evaluations should strive to incorporate land–atmosphere feedback, at least within a partially coupled model. © 1997 John Wiley & Sons, Ltd.     

Modeling surface water–groundwater interaction in New Zealand: Model development and application

Most rivers worldwide have a strong interaction with groundwater when they leave the mountains and flow over alluvial plains before flowing into the seas or disappearing in the deserts, and in New Zealand, typically, rivers lose water to the groundwater in the upper plains and generally gain water from the groundwater in the lower plains. Aiming at simulating surface water–groundwater interaction nationally in New Zealand, we developed a conceptual groundwater module for the national hydrologic model TopNet to simulate surface water–groundwater interaction, groundwater flow, and intercatchment groundwater flow. The developed model was applied to the Pareora catchment in South Island of New Zealand, where there are concurrent spot gauged flows. Results show that the model simulations not only fit quite well to flow measurement but also to concurrent spot gauged flows, and compared to the original TopNet, it has a significant improvement in the low flows. Sensitivity analysis shows river flow is sensitive to the river losing/gaining rate instead of groundwater characteristic, while groundwater storage is sensitive to both river losing/gaining rate and groundwater characteristic. This indicates our conceptual approach is promising for nationwide modeling without the large amount of geology and aquifer data typically required by physically‐based modeling approaches.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Impact of soil moisture heterogeneity length scale and gradients on daytime coupled land‐cloudy boundary layer interactions

Using a coupled large‐eddy simulation–land surface model framework, the impact of two‐dimensional soil moisture heterogeneity on the cloudy boundary layer under varied free‐atmosphere stabilities is investigated. Specifically, the impacts of soil moisture heterogeneity length scale and heterogeneity in terms of soil moisture gradients on micrometeorological states, surface fluxes, boundary layer characteristics, and cloud development are examined. The results show that mesoscale circulations due to surface heterogeneity in soil moisture play an important role in transferring water vapour within the boundary layer and in regulating cloud distribution at the entrainment zone, which, in turn, provides feedbacks on boundary layer/surface energy budgets. The initial domain‐averaged soil moisture is identical for all homogenous and heterogeneous cases; however, the soil moisture heterogeneity in gradient and length scale between dry and wet regions has a significant impact on the estimates of near‐surface micrometeorological properties and surface fluxes, which further affect the boundary layer states and characteristics. Both liquid water potential temperature and liquid water mixing ratio increase with an increasing soil moisture gradient, whereas the amount of specific humidity decreases. Heterogeneity length scale and free atmosphere stability also amplify these impacts on the boundary layer structure and cloud formation. In a low atmospheric stability condition that potentially allows for a deeper boundary layer and a higher entrainment rate, cloud base height and cloud thickness significantly increase as the soil moisture gradient and length scale increase. Analysis to differentiate the influences of surface heterogeneity type (i.e. length scale vs gradient) shows that in general soil moisture gradient provides a larger impact than heterogeneity length scale, although the heterogeneity length scale is large enough to initiate circulation features responsible for differences in the coupled system between homogeneous and heterogeneous soil moisture cases. Copyright © 2012 John Wiley & Sons, Ltd.     

Intermittent chaos in nonlinear wave–wave interactions in space plasmas

There is increasing observational evidence of nonlinear wave–wave interactions in space and astrophysical plasmas. We first review a number of theoretical models of nonlinear wave–wave interactions which our group has developed in the past years. We next describe a nonlinear three-mode truncated model of Alfvén waves, involving resonant interactions of one linearly unstable mode and two linearly damped modes. We construct a bifurcation diagram for this three-wave model and investigate the phenomenon of intermittent chaos. The theoretical results presented in this paper can improve our understanding of intermittent time series frequently observed in space and astrophysical plasmas.     

Implementation of a second‐order paraxial boundary condition for a water‐saturated layered half‐space in plane strain

A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd.