Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

GEOCHEMICAL RECOGNITION OF SPILLED SEDIMENTS USED IN NUMERICAL MODEL VALIDATION

A fixed link (tunnel and bridge, in total 16 km) was constructed between Sweden and Denmark during 1995-2000. As part of the work, approximately 16 million tonnes of seabed materials (limestone and clay till) were dredged, and about 0.6 million tonnes of these were spilled in the water. Modelling of the spreading and sedimentation of the spilled sediments took place as part of the environmental monitoring of the construction activities. In order to verify the results of the numerical modelling of sediment spreading and sedimentation, a new method with the purpose of distinguishing between the spilled sediments and the naturally occurring sediments was developed. Because the spilled sediments tend to accumulate at the seabed in areas with natural sediments of the same size, it is difficult to separate these based purely on the physical properties. The new method is based on the geo-chemical differences between the natural sediment in the area and the spill. The basic propertiesused are the higher content of calcium carbonate material in the spill as compared to the natural sediments and the higher Ca/Sr ratio in the spill compared to shell fragments dominating the natural calcium carbonate deposition in the area. The reason for these differences is that carbonate derived from recent shell debris can be discriminated from Danien limestone, which is the material in which the majority of the dredging took place, on the basis of the Ca/Sr ratio being 488 in Danien Limestone and 237 in shell debris. The geochemical recognition of the origin of the sediments proved useful in separating the spilled from the naturally occurring sediments. Without this separation, validation of the modelling of accumulation of spilled sediments would not have been possible. The method has general validity and can be used in many situations where the origin of a given sediment is sought.     

Gravity reduction with three-dimensional atmospheric pressure data for precise ground gravity measurements

The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.     

Effects of pollen leaching and microbial degradation on organic carbon and nutrient availability in lake water

Nutrient fluxes across terrestrial-aquatic boundaries and their subsequent integration into lake nutrient cycles are currently a major topic of aquatic research. Although pollen represents a good substrate for microorganisms, it has been neglected as a terrestrial source of organic matter in lakes. In laboratory experiments, we incubated pollen grains of Pinus sylvestris in water of lakes with different trophy and pH to estimate effects of pollen input and its subsequent microbial degradation on nutrient dynamics. In this ex situ experiment, we measured concentrations of organic carbon, phosphorus and nitrogen in the surrounding water as well as microbial dynamics (bacteria and fungal sporangia) at well-controlled conditions. Besides leaching, chemical and microbial decomposition of pollen was strongest within the first week of incubation. This led to a marked increase of soluble reactive phosphorus and total dissolved nitrogen (up to 0.04 and 1.5 mg L⁻¹, respectively, after 5 days of incubation) in the ambient water. In parallel, pollen grains were rapidly colonized by heterotrophic bacteria and aquatic fungi. Leaching and microbial degradation of pollen accounted for ≥80, ≥40, ≥50% for organic C, N and P, respectively, and did not significantly differ among water samples from the studied lakes. Thus, pollen introduces high amounts of bio-available terrestrial organic matter and nutrients into surface waters within a short time. A rough calculation on P input into oligotrophic Lake Stechlin indicates that pollen plays an important ecological role in nutrient cycling of temperate lakes. This requires further attention in aquatic ecology.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.