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.     

Hyporheic exchange flow around constructed in‐channel structures and implications for restoration design

In‐channel rock vane structures are widely used in stream restoration as a way to reduce stream channel erosion and create pool or riffle features. When these structures change hydraulic gradients they may affect ecological stream functions, such as hyporheic exchange flow (HEF) patterns. A study of constructed in‐channel structure controls on HEF was conducted in the third‐order Batavia Kill, New York using stream and hyporheic temperature amplitude analysis and computational fluid dynamics (CFD) hydraulic simulations. Temperature monitors were installed in the water column and channel bed at six locations around each of seven in‐channel restoration structures (three cross‐vanes and four J‐hooks) at baseflow in 2007. Elevation surveys of the structures were then used to simulate HEF using CFD. The results indicate a pattern of pronounced upwelling in the run section just below the structure, upwelling transitioning to downwelling within the pool, and pronounced downwelling in the glide out of the pool. This pattern is consistent with natural riffle pool sequences. The direction of HEF inferred from the temperature amplitude analysis agreed with the direction of flow simulated with CFD at 80% of the locations, and the few disagreements were expected due to model limitations. CFD simulation demonstrated that increasing stream flows result in changes in HEF spatial patterns and magnitude at each structure. This work illustrates how CFD simulations can guide design of in‐channel restoration structures for HEF function. Copyright © 2009 John Wiley & Sons, Ltd.     

Geostatistical analysis and conditional simulation for estimating the spatial variability of hydraulic conductivity in the Choushui River alluvial fan,Taiwan

This work evaluated the spatial variability and distribution of heterogeneous hydraulic conductivity (K) in the Choushui River alluvial fan in Taiwan, using ordinary kriging (OK) and mean and individual sequential Gaussian simulations (SGS). A baseline flow model constructed by upscaling parameters was inversely calibrated to determine the pumping and recharge rates. Simulated heads using different K realizations were then compared with historically measured heads. A global/local simulated error between simulated and measured heads was analysed to assess the different spatial variabilities of various estimated K distributions. The results of a MODFLOW simulation indicate that the OK realization had the smallest sum of absolute mean simulation errors (SAMSE) and the SGS realizations preserved the spatial variability of the measured K fields. Moreover, the SAMSE increases as the spatial variability of the K field increases. The OK realization yields small local simulation errors in the measured K field of moderate magnitude, whereas the SGS realizations have small local simulation errors in the measured K fields, with high and low values. The OK realization of K can be applied to perform a deterministic inverse calibration. The mean SGS method is suggested for constructing a K field when the application focuses on extreme values of estimated parameters and small calibration errors, such as in a simulation of contaminant transport in heterogeneous aquifers. The individual SGS realization is useful in stochastically assessing the spatial uncertainty of highly heterogeneous aquifers. Copyright © 2004 John Wiley & Sons, Ltd.     

Testing solar forcing of pervasive Holocene climate cycles

The temporal and spatial extent of Holocene climate change is an area of considerable uncertainty, with solar forcing recently proposed to be the origin of cycles identified in the North Atlantic region. To address these issues we have developed an annually resolved record of changes in Irish bog tree populations over the last 7468 years which, together with radiocarbon‐dated bog and lake‐edge populations, extend the dataset back to ～9000 yr ago. The Irish trees underpin the internationally accepted radiocarbon calibration curve, used to derive a proxy of solar activity, and allow us to test solar forcing of Holocene climate change. Tree populations and age structures provide unambiguous evidence of major shifts in Holocene surface moisture, with a dominant cyclicity of 800 yr, similar to marine cycles in the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. The cycles, however, are not coherent with changes in solar activity (both being on the same absolute timescale), indicating that Holocene North Atlantic climate variability at the millennial and centennial scale is not driven by a linear response to changes in solar activity. Copyright © 2005 John Wiley & Sons, Ltd.     

Detecting subtle hydrochemical anomalies with multivariate statistics: an example from ‘homogeneous’ groundwaters in the Great Artesian Basin,Australia

The major ion composition of Great Artesian Basin groundwater in the lower Namoi River valley is relatively homogeneous in chemical composition. Traditional graphical techniques have been combined with multivariate statistical methods to determine whether subtle differences in the chemical composition of these waters can be delineated. Hierarchical cluster analysis and principal components analysis were successful in delineating minor variations within the groundwaters of the study area that were not visually identified in the graphical techniques applied. Hydrochemical interpretation allowed geochemical processes to be identified in each statistically defined water type and illustrated how these groundwaters differ from one another. Three main geochemical processes were identified in the groundwaters: ion exchange, precipitation, and mixing between waters from different sources. Both statistical methods delineated an anomalous sample suspected of being influenced by magmatic CO₂ input. The use of statistical methods to complement traditional graphical techniques for waters appearing homogeneous is emphasized for all investigations of this type. Copyright © 2006 John Wiley & Sons, Ltd.     

Relative sea‐level changes,glacial isostatic modelling and ice‐sheet reconstructions from the British Isles since the Last Glacial Maximum

While contributing <1 m equivalent eustatic sea‐level rise the British Isles ice sheet produced glacio‐isostatic rebound in northern Britain of similar magnitude to eustatic sea‐level change, or global meltwater influx, over the last 18 000 years. The resulting spatially variable relative sea‐level changes combine with observations from far‐field locations to produce a rigorous test for quantitative models of glacial isostatic adjustment, local ice‐sheet history and global meltwater influx. After a review of the attributes of relative sea‐level observations significant for constraining large‐scale models of the isostatic adjustment process we summarise long records of relative sea‐level change from the British Isles and far‐field locations. We give an overview of different global theoretical models of the isostatic adjustment process before presenting intercomparisons of observed and predicted relative sea levels at sites in the British Isles and far‐field for a range of Earth and ice model parameters in order to demonstrate model sensitivity and the resolving power available from using evidence from the British Isles. For the first time we show a good degree of fit between relative sea‐level observations and predictions that are based upon global Earth and ice model parameters, independently derived from analysis of far‐field data, with a terrain‐corrected model of the British Isles ice sheet that includes extensive glaciation of the North Sea and western continental shelf, that does not assume isostatic equilibrium at the Last Glacial Maximum and keeps to trimline constraints of ice surface elevation. We do not attempt to identify a unique solution for the model lithosphere thickness parameter or the local‐scale detail of the ice model in order to provide a fit for all sites, but argue that the next stage should be to incorporate an ice‐sheet model that is based on quantitative, glaciological model simulations. We hope that this paper will stimulate this debate and help to integrate research in glacial geomorphology, glaciology, sea‐level change, Earth rheology and quantitative modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Variations in dissolved CO2 and CH4 in a first‐order stream and catchment: an investigation of soil–stream linkages

Spatial and seasonal variations in CO₂ and CH₄ concentrations in streamwater and adjacent soils were studied at three sites on Brocky Burn, a headwater stream draining a peatland catchment in upland Britain. Concentrations of both gases in the soil atmosphere were significantly higher in peat and riparian soils than in mineral soils. Peat and riparian soil CO₂ concentrations varied seasonally, showing a positive correlation with air and soil temperature. Streamwater CO₂ concentrations at the upper sampling site, which mostly drained deep peats, varied from 2·8 to 9·8 mg l^?1 (2·5 to 11·9 times atmospheric saturation) and decreased markedly downstream. Temperature‐related seasonal variations in peat and riparian soil CO₂ were reflected in the stream at the upper site, where 77% of biweekly variation was explained by an autoregressive model based on: (i) a negative log‐linear relationship with stream flow; (ii) a positive linear relationship with soil CO₂ concentrations in the shallow riparian wells; and (iii) a negative linear relationship with soil CO₂ concentrations in the shallow peat wells, with a significant 2‐week lag term. These relationships changed markedly downstream, with an apparent decrease in the soil–stream linkage and a switch to a positive relationship between stream flow and stream CO₂. Streamwater CH₄ concentrations also declined sharply downstream, but were much lower (<0·01 to 0·12 mg l^?1) than those of CO₂ and showed no seasonal variation, nor any relationship with soil atmospheric CH₄ concentrations. However, stream CH₄ was significantly correlated with stream flow at the upper site, which explained 57% of biweekly variations in dissolved concentrations. We conclude that stream CO₂ can be a useful integrative measure of whole catchment respiration, but only at sites where the soil–stream linkage is strong. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of soil interaction on the performance of liquid column dampers for seismic applications

The effects of soil–structure interaction (SSI) while designing the liquid column damper (LCD) for seismic vibration control of structures have been presented in this study. The formulation for the input–output relation of a flexible‐base structure with attached LCD has been presented. The superstructure has been modelled by a single‐degree‐of‐freedom (SDOF) system. The non‐linearity in the orifice damping of the LCD has been replaced by equivalent linear viscous damping by using equivalent linearization technique. The force–deformation relationships and damping characteristics of the foundation have been described by complex valued impedance functions. Through a numerical stochastic study in the frequency domain, the various aspects of SSI on the functioning of the LCD have been illustrated. A simpler approach for studying the LCD performance considering SSI, using an equivalent SDOF model for the soil–structure system available in literature by Wolf (Dynamic Soil–Structure Interaction. International Series in Civil Engineering and Engineering Mechanics. Prentice‐Hall: Englewood Cliffs, NJ, 1985) has also been presented. Copyright © 2005 John Wiley & Sons, Ltd.