Role of reservoir regulation and groundwater feedback in a simulated ground-soil-vegetation continuum: A long-term regional scale analysis

The regional terrestrial water cycle is strongly altered by human activities. Among them, reservoir regulation is a way to spatially and temporally allocate water resources in a basin for multi-purposes. However, it is still not sufficiently understood how reservoir regulation modifies the regional terrestrial- and subsequently, the atmospheric water cycle. To address this question, the representation of reservoir regulation into the terrestrial component of fully coupled regional Earth system models is required. In this study, an existing process-based reservoir network module is implemented into NOAH-HMS, that is, the terrestrial component of an atmospheric–hydrologic modelling system, namely, the WRF-HMS. It allows to quantitatively differentiate role of reservoir regulation and of groundwater feedback in a simulated ground-soil-vegetation continuum. Our study focuses on the Poyang Lake basin, where the largest freshwater lake of China and reservoirs of different sizes are located. As compared to streamflow observations, the newly extended NOAH-HMS slightly improves the streamflow and streamflow duration curves simulation for the Poyang Lake basin for the period 1979–1986. The inclusion of reservoir regulation leads to major changes in the simulated groundwater recharges and evaporation from reservoirs at local scale, but has minor effects on the simulated soil moisture and surface runoff at basin scale. The performed groundwater feedback sensitivity analysis shows that the strength of the groundwater feedback is not altered by the consideration of reservoir regulation. Furthermore, both reservoir regulation and groundwater feedback modify the partitioning of the simulated evapotranspiration, thus affecting the atmospheric water cycle in the Poyang Lake region. This finding motivates future research with our extended fully coupled atmospheric–hydrologic modelling system by the community.     

Relating millimeter-scale turbulence to meter-scale subtidal erosion and accretion across the fringe of a coastal mangrove forest

Within a wave-exposed mangrove forest, novel field observations are presented, comparing millimeter-scale turbulent water velocity fluctuations with contemporaneous subtidal bed elevation changes. High-resolution velocity and bed level measurements were collected from the unvegetated mudflat, at the mangrove forest fringe, and within the forest interior over multiple tidal cycles (flood–ebb) during a 2-week period. Measurements demonstrated that the spatial variability in vegetation density is a control on sediment transport at sub-meter scales. Scour around single and dense clusters of pneumatophores was predicted by a standard hydraulic engineering equation for wave-induced scour around regular cylinders, when the cylinder diameter in the equations was replaced with the representative diameter of the dense pneumatophore clusters. Waves were dissipated as they propagated into the forest, but dissipation at infragravity periods (> 30 s) was observed to be less than dissipation at shorter periods (< 30 s), consistent with the predictions of a simple model. Cross-wavelet analysis revealed that infragravity-frequency fluctuations in the bed level were occasionally coherent with velocity, possibly indicating scour upstream of dense pneumatophore patches when infragravity waves reinforced tidal currents. Consequently, infragravity waves were a likely driver of sediment transport within the mangrove forest. Near-bed turbulent kinetic energy, estimated from the turbulent dissipation rate, was also correlated with bed level changes. Specifically, within the mangrove forest and over the unvegetated mudflat, high-energy events were associated with erosion or near-zero bed level change, whereas low-energy events were associated with accretion. In contrast, no single relationship between bed level changes and mean current velocity was applicable across both vegetated and unvegetated regions. These observations support the theory that sediment mobilization scales with turbulent energy, rather than mean velocity, a distinction that becomes important when vegetation controls the development of turbulence.     

Weighting of Observed Heads for the Inverse Problem

In using least-squares parameter estimation techniques to solve for hydrogeologic parameters, one may use a weighting function to reflect differing reliabilities of head measurements. In studies published to date, the weighting function has been used in an ad boc manner or not at all. The inverse square of the observed hydraulic gradient, adjusted to reflect the modeler's perception of geologic heterogeneity and data reliability, is typically an appropriate weighting function.     

Pits,rifts and slumps: the summit structure of Piton de la Fournaise

A clear model of structures and associated stress fields of a volcano can provide a framework in which to study and monitor activity. We propose a volcano-tectonic model for the dynamics of the summit of Piton de la Fournaise (La Reunion Island, Indian Ocean). The summit contains two main pit crater structures (Dolomieu and Bory), two active rift zones, and a slumping eastern sector, all of which contribute to the actual fracture system. Dolomieu has developed over 100 years by sudden large collapse events and subsequent smaller drops that include terrace formation. Small intra-pit collapse scars and eruptive fissures are located along the southern floor of Dolomieu. The western pit wall of Dolomieu has a superficial inward dipping normal fault boundary connected to a deeper ring fault system. Outside Dolomieu, an oval extension zone containing sub-parallel pit-related fractures extends to a maximum distance of 225 m from the pit. At the summit the main trend for eruptive fissures is N80°, normal to the north–south rift zone. The terraced structure of Dolomieu has been reproduced by analogue models with a roof to width ratio of approximately 1, suggesting an original magma chamber depth of about 1 km. Such a chamber may continue to act as a storage location today. The east flank has a convex–concave profile and is bounded by strike-slip fractures that define a gravity slump. This zone is bound to the north by strike-slip fractures that may delineate a shear zone. The southern reciprocal shear zone is probably marked by an alignment of large scoria cones and is hidden by recent aa lavas. The slump head intersects Dolomieu pit and may slide on a hydrothermally altered layer known to be located at a depth of around 300 m. Our model has the summit activity controlled by the pit crater collapse structure, not the rifts. The rifts become important on the mid-flanks of the cone, away from pit-related fractures. On the east flank the superficial structures are controlled by the slump. We suggest that during pit subsidence intra-pit eruptions may occur. During tumescence, however, the pit system may become blocked and a flank eruption is more likely. Intrusions along the rift may cause deformation that subsequently increases the slump’s potential to deform. Conversely, slumping may influence the east flank stress distribution and locally control intrusion direction. These predictions can be tested with monitoring data to validate the model and, eventually, improve monitoring.     

Capturing the residence time boundary layer—application to the Scheldt Estuary

At high Peclet number, the residence time exhibits a boundary layer adjacent to incoming open boundaries. In a Eulerian model, not resolving this boundary layer can generate spurious oscillations that can propagate into the area of interest. However, resolving this boundary layer would require an unacceptably high spatial resolution. Therefore, alternative methods are needed in which no grid refinement is required to capture the key aspects of the physics of the residence time boundary layer. An extended finite element method representation and a boundary layer parameterisation are presented and tested herein. It is also explained how to preserve local consistency in reversed time simulations so as to avoid the generation of spurious residence time extrema. Finally, the boundary layer parameterisation is applied to the computation of the residence time in the Scheldt Estuary (Belgium/The Netherlands). This timescale is simulated by means of a depth-integrated, finite element, unstructured mesh model, with a high space–time resolution. It is seen that the residence time temporal variations are mainly affected by the semi-diurnal tides. However, the spring–neap variability also impacts the residence time, particularly in the sandbank and shallow areas. Seasonal variability is also observed, which is induced by the fluctuations over the year of the upstream flows. In general, the residence time is an increasing function of the distance to the mouth of the estuary. However, smaller-scale fluctuations are also present: they are caused by local bathymetric features and their impact on the hydrodynamics.     

Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

Homogeneous nucleation of amorphous solid water particles in the upper mesosphere

Condensed water particles are known to exist in the high latitude upper mesosphere during the summer months. However, the mechanism or mechanisms through which they nucleate remains uncertain. It is postulated here that particles of amorphous solid water (ASW, condensed water with a non-crystalline structure) may nucleate homogeneously in the summer mesosphere. Using classical nucleation theory and a one-dimensional model, it is shown that more than 10⁵ cm^?3 amorphous solid water particles can nucleate homogeneously under mesopause conditions. Furthermore, it is shown that homogeneous nucleation competes with heterogeneous nucleation on meteoric smoke particles when the cooling rate is >0.5 K/h. The homogeneous nucleation of amorphous solid water could provide an explanation for the high density of ice particles (many thousands per cm³) thought to be required for electron depletions in the upper mesosphere. A parameterisation for homogeneous nucleation is presented which can be used in other mesospheric cloud models.     

Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Pit crater structure and processes governing persistent activity at Masaya Volcano, Nicaragua

 Persistent activity at Masaya Volcano, Nicaragua, is characterised by cycles of intense degassing, lava lake development and pit crater formation. It provides a useful site to study the processes which govern such activity, because of its easy accessibility and relatively short cycles (years to decades). An understanding of the present activity is important because Masaya is visited by large numbers of tourists, is located close to major cities and has produced voluminous lavas, plinian eruptions and ignimbrites in the recent past. We provide structural and geophysical data that characterise the "normal" present state of activity. These indicate that the ongoing degassing phase (1993 to present) was not caused by fresh magma intrusion. It was associated with shallow density changes within the active Santiago pit crater. The activity appears to be associated predominantly with shallow changes in the pit crater structure. More hazardous activity will occur only if there are significant departures from the present gravity, deformation and seismic signatures. Received: 16 May 1997 / Accepted: 29 October 1997     

Comparative hydrogeology – reference analysis of groundwater dynamics from neighbouring observation wells

ABSTRACT

Groundwater level fluctuations are caused by spatial and temporal superposition of processes within and outside the aquifer system. Most of the subsurface processes are usually observed on a small scale. Upscaling to the regional scale, as required for future climate change scenarios, is difficult due to data scarcity and increasing complexity. In contrast to the limited availability of system characteristics, high-resolution data records of groundwater hydrographs are more generally available. Exploiting the information contained in these records should thus be a priority for analysis of the chronical lack of data describing groundwater system characteristics. This study analyses the applicability of 63 indices derived from daily hydrographs to quantify different dynamics of groundwater levels in unconfined gravel aquifers from three groundwater regions (Bavaria, Germany). Based on the results of two different skill tests, the study aids index selection for different dynamic components of groundwater hydrographs.