A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.     

A mesocosm experiment investigating the effects of substratum quality and wave exposure on the survival of fish eggs

In a mesocosm experiment, the attachment of bream (Abramis brama) eggs to spawning substrata with and without periphytic biofilm coverage and their subsequent survival with and without low-intensity wave exposure were investigated. Egg attachment was reduced by 73% on spawning substrata with a natural periphytic biofilm, compared to clean substrata. Overall, this initial difference in egg numbers persisted until hatching. The difference in egg numbers was even increased in the wave treatment, while it was reduced in the no-wave control treatment. Exposure to a low-intensity wave regime affected egg development between the two biofilm treatments differently. Waves enhanced egg survival on substrata without a biofilm but reduced the survival of eggs on substrata with biofilm coverage. In the treatment combining biofilm-covered substrata and waves, no attached eggs survived until hatching. In all treatments, more than 75% of the eggs became detached from the spawning substrata during the egg incubation period, and <1% of these detached eggs survived within the substratum interstices. Hence, detached eggs contributed little to the reproductive success we observed. The implications of these results on the spawning success of fish species that use freshly inundated spawning substrata in shallow waters are discussed.     

Long‐term variability of proglacial groundwater‐fed hydrological systems in an area of glacier retreat,Skeiðarársandur,Iceland

Proglacial groundwater‐fed features, such as seeps, substantially impact proglacial geomorphology, hydrology, and ecology. However, there is a paucity of research on the impacts of climate change and glacier retreat on the extent of these important features. This paper aims to investigate the impact of glacier retreat on proglacial groundwater levels and on the extent of groundwater‐fed seeps. Research has taken place in western Skeiðarársandur, the large proglacial outwash plain of Skeiðarárjökull, a retreating temperate glacier in southeast Iceland. Changes in the extent of proglacial groundwater seeps were mapped using historical aerial photographs from 1986, 1997, and 2012. Proglacial groundwater levels were monitored in shallow boreholes between 2000 and 2012. The western margin of Skeiðarárjökull has retreated approximately 1 km beyond its position in 1986. However, this retreat was punctuated by short periods of readvance. The geomorphology and groundwater systems at the site were substantially impacted by the November 1996 jökulhlaup, whose deposits altered approximately 18% of the area of groundwater seeps. The surface areas of groundwater seeps and lakes in the study area have declined by ~97% between 1986 and 2012. Most of the decline took place after 1997, when the mean annual rate of retreat increased three‐fold. Groundwater levels also declined substantially between 2000 and 2012, although this trend varies spatially. The paper provides a conceptual model of the controls on proglacial shallow groundwater systems. Direct impacts of glacier retreat are suggested as the main cause for the declines in proglacial groundwater levels and in the extent of groundwater seeps. These declines are expected to adversely impact sandur ecology. Copyright © 2014 John Wiley & Sons, Ltd.     

Minor element- and carbonaceous material thermometry of high-grade metapelites from the Sauwald Zone, Southern Bohemian Massif (Upper Austria)

The Sauwald area is located at the southern rim of the Bohemian Massif and contains migmatites and high-grade metapelitic and granitic gneisses. These rocks were metamorphosed during the post-collisional high-T/low-P stage of the Variscan metamorphic event (~330 Ma). Metapelitic samples were taken from two localities near Kössldorf and Pyret in Upper Austria and the investigated samples contain the mineral assemblage garnet + cordierite + spinel + sillimanite + K-feldspar + quartz + biotite + muscovite + magnetite + graphite. An important aspect of this study is the evaluation of previously published P-T estimates for these high-grade metapelites (Knop et al. 1995; Tropper et al. 2006) involving Ti-in-biotite, Na-in-cordierite thermometry and the micro-Raman thermometer based upon the degree of crystallization of carbonaceous material. In the two samples studied three texturally and chemically different biotites are distinguished. Biotite inclusions in garnet have the highest Ti contents of 5–6 wt.% TiO₂. Matrix biotites contain 2–4 wt.% TiO₂ and biotites from late-stage muscovite-biotite symplectites contain <2 wt.% TiO₂. This corresponds to temperatures of 730–760°C (stage 1), 600–700°C (stage 2), and 550–610°C (stage 3). Since the Ti-in-biotite thermometer strongly depends on X _Mg of biotite, which is susceptible to changes during retrogression the calculated temperatures for stage 1 are interpreted as minimum temperatures of the peak metamorphic stage. The Na contents of the studied cordierites vary from 0.1 to 0.2 wt.% Na₂O. Application of the Na-in-cordierite thermometer yields temperatures in the range of 770–900°C; they are strongly dependent on the bulk Na₂O content of the samples. The micro-Raman geothermometer of graphite was applied to carbonaceous material, which occurs as inclusions in garnet and cordierite. It yielded a maximum temperature >650°C, i.e. above the calibration limit of this method. This study shows that the obtained temperature estimates agree well with the P-T estimates based on phase equilibrium thermobarometry (Knop et al. 1995; Tropper et al. 2006), thus illustrating the validity of these thermometers. Nevertheless in order to more precisely constrain the metamorphic evolution of these high-grade rocks, better constrained experimental calibrations of, for instance the Na-in-cordierite thermometer, are clearly needed.     

Earthquakes in Switzerland and surrounding regions during 2008

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2008. During this period, 451 earthquakes and 75 quarry blasts were detected and located in the region under consideration. The three strongest events occurred in the Valais, near Lac des Toules (M_L 3.6), and in Graubünden, near Ilanz (M_L 3.7) and Paspels (M_L 4.0). Although felt by the population, they were not reported to have caused any damage. However, with a total of only 15 events with M_L ≥ 2.5, the seismic activity in the year 2008 was far below the average over the previous 33 years.