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.     

Land-surface segmentation as a method to create strata for spatial sampling and its potential for digital soil mapping

Sampling efforts are constrained by limited availability of resources. Therefore, methods to reduce the number of samples, while still achieving reasonable accuracy are needed. Land-surface segmentation (LSS) has proven a powerful technique to partition digital elevation models (DEMs) and their derivatives into relatively homogeneous areas, which can be further employed as support in soil sampling. Though topography is one of the main soil forming factors, a robust assessment of the potential of this technique to digital soil mapping (DSM) is still missing. In this study, we aimed at evaluating the potential of LSS in stratifying a landscape into relatively homogeneous areas, which can be used as strata for guiding the selection of sampling points in DSM. The experiments were carried out in two study areas where soil samples were available. Land-surface derivatives were derived from DEMs and segmented with a tool based on the multiresolution segmentation algorithm, into objects considered as homogeneous soil-landscape divisions. Thus, one sample was randomly selected within each segment from the existing sample data, based on which predictions of soil classes/sub-orders and properties, i.e. soil texture and A-horizon thickness, were made. Results were compared with predictions based on simple random sampling (SRS) and conditioned Latin hypercube (cLHS). The segmentation-based sampling (SBS) scheme performed better than SRS and cLHS schemes in predicting the A-horizon thickness, soil texture fractions and soil classes, showing a high potential of LSS in stratifying a landscape for the purposes of DSM. The novelty of this study is in the way strata are constructed, rather than in the sampling design itself. Further research is needed to demonstrate the value of a SBS design for practical use. The analyses presented here further highlight the importance of considering locally adaptive techniques in optimization of sampling schemes and predictions of soil properties.     

Oxygen and Nutrient Exchanges at the Sediment-Water Interface: a Global Synthesis and Critique of Estuarine and Coastal Data

Estuarine and coastal marine sediment-water fluxes are considered to be important ecological features, but a global-scale assessment has yet to be developed. Goals of this work were to assemble a global-scale database of net sediment-water flux measurements, examine measurement techniques, characterize the geographic distribution and magnitude of sediment fluxes, explore the data for controls on sediment flux magnitude, and assess the importance of sediment fluxes in ecosystem-level metabolism and primary production. We examined 480 peer-reviewed sources and found sediment flux data for 167 estuarine and coastal systems. Most measurements were made in North America, Europe, and Australia. Fluxes varied widely among systems, some by several orders of magnitude. Inter-annual variability within sites was less than an order of magnitude but time series flux data to evaluate this were rare. However, limited time series data exhibited large and rapid responses to decreased external nutrient loading rates, climate change effects (possible temperature effects), and variability in trophic conditions. Comparative analyses indicated organic matter supply to sediments set the upper limits of flux magnitude, with other factors playing secondary roles. Two metrics were developed to assess ecosystem-level importance of sediment-water fluxes. Sediments represented 30% or more of depth-integrated rates of aerobic system respiration at depths of <10 m. An annual phytoplankton production data set was used to estimate N and P demand; sediments supplied an average of 15–32% of N and 17–100% of P demand and, in some cases, was as large or larger than external nutrient inputs. The percent of demand supplied by sediments was highest in temperate latitudes and lower in high and tropical latitudes.     

Dimethylsulfide production in Sargasso Sea eddies

Lagrangian time series of dimethylsulfide (DMS) concentrations from a cyclonic and an anticyclonic eddy in the Sargasso Sea were used in conjunction with measured DMS loss rates and a model of vertical mixing to estimate gross DMS production in the upper 60 m during summer 2004. Loss terms included biological consumption, photolysis, and ventilation to the atmosphere. The time- and depth (0–60 m)-averaged gross DMS production was estimated to be 0.73±0.09 nM d⁻¹ in the cyclonic eddy and 0.90±0.15 nM d⁻¹ in the anticyclonic eddy, with respective DMS replacement times of 5±1 and 6±1 d. The higher estimated rate of gross production and lower measured loss rate constants in the anticyclonic eddy were equally responsible for this eddy's 50% higher DMS inventory (0–60 m). When normalized to chlorophyll and total dimethylsulfoniopropionate (DMSP), estimated gross production in the anticyclonic eddy was about twice that in the cyclonic eddy, consistent with the greater fraction of phytoplankton that were DMSP producers in the anticyclonic eddy. Higher rates of gross production were estimated below the mixed layer, contributing to the subsurface DMS maximum found in both eddies. In both eddies, gas exchange, microbial consumption, and photolysis were roughly equal DMS loss terms in the surface mixed layer (0.2–0.4 nM d⁻¹). Vertical mixing was a substantial source of DMS to the surface mixed layer in both eddies (0.2–0.3 nM d⁻¹) owing to the relatively high DMS concentrations below the mixed layer. Estimated net biological DMS production rates (gross production minus microbial consumption) in the mixed layer were substantially lower (by almost a factor of 3) than those estimated in a previous study of the Sargasso Sea, which may explain the relatively low mixed-layer DMS concentrations found here during July 2004 (3 nM) compared to previous summers (4–6 nM).     

A Modified Form of Mild-Slope Equation with Weakly Nonlinear Effect

Nonlinear effect is of importance to waves propagating from deep water to shallow water.Thenon-linearity of waves is widely discussed due to its high precision in application.But there are still someproblems in dealing with the nonlinear waves in practice.In this paper,a modified form of mild-slope equa-tion with weakly nonlinear effect is derived by use of the nonlinear dispersion relation and the steady mild-slope equation containing energy dissipation.The modified form of mild-slope equation is convenient to solvenonlinear effect of waves.The model is tested against the laboratory measurement for the case of a submergedelliptical shoal on a slope beach given by Berkhoff et al,The present numerical results are also comparedwith those obtained through linear wave theory.Better agreement is obtained as the modified mild-slope e-quation is employed.And the modified mild-slope equation can reasonably simulate the weakly nonlinear ef-fect of wave propagation from deep water to coast.     

Chemical variations and genetic relationships between the Hess and Foy offset dikes at the Sudbury impact structure

Offset dikes are found concentrically around—and extending radially outward from—the Sudbury Igneous Complex (SIC), which represents an ~3 km thick differentiated impact melt sheet. The dikes are typically composed of an inclusion‐rich, so‐called quartz diorite (IQD) in the center of the dike, and an inclusion‐poor quartz diorite (QD) along the margins of the dike. New exposures of the intersection between the concentric Hess and radial Foy offset dikes provide an excellent opportunity to understand the relationship between the radial and concentric offset dikes and their internal phases. The goal was to constrain the timing of the dike emplacements relative to the impact and formation of the SIC. Results herein suggest that (1) the timing between the emplacement of the QD and IQD melts was geologically short, (2) the Hess and Foy dikes coexisted as melts at the same time and the intersection between them represents a mixture of the two, (3) the Foy dike has a slightly more evolved chemical composition than the Hess dike, and (4) the IQD melt from the Foy dike underwent some degree of chemical fractionation after its initial emplacement.     

DARCY'S COEFFICIENT OF PERMEABILITY AS SYMMETRIC TENSOR OF SECOND RANK

Synopsis

The dynamic equation of motion that governs the laminar flow of water through soils is the empirical equation of Darcy. According to Darcy's equation the velocity of the flowing water is proportional to the hydraulic gradient under which the water is flowing, with the constant of proportionality being the coefficient of permeability. The interesting question arising is whether or not the coefficient of permeability is a scalar quantity (having only a magnitude) or a vector (having both magnitude and direction). It is proved, in the present paper, that the permeability coefficient is neither a scalar nor a vector but a symmetric tensor of second rank. The fact that the permeability tensor is symmetric gives rise to great simplifications and permits a simple graphical construction of the tensor ellipsoid. Having the tensor ellipsoid, the determination of the direction at which the water will flow under a known imposed hydraulic gradient can be found graphically. In case of isotropic soils (the permeability coefficient has the same value along any direction) the ellipsoid reduces to a sphere and the tensor becomes a scalar. In the general case of anisotropic soils the permeability tensor is an entity with nine elements, six of which are independent representing pure extension or contraction along the three principal coordinate axes, thus transforming the permeability sphere into an ellipsoid and vice versa. It should be noted that in anisotropic soils the only directions along which the flow takes place in the direction of the hydraulic gradient are those of the principal axes of the tensor ellipsoid.

Permeability tests were conducted on anisotropic sandstone samples taken at different directions with respect to rectangular coordinates. The permeability coefficient values plotted on a two-dimensional polar coordinate graph paper give rise to an ellipse substantiating therefore the tensor concept of the permeability coefficient. The graphical construction of the tensor ellipse and the use of it in order to obtain the direction of flow by knowing the direction of the hydraulic gradient is also shown.