Geomorphological control of subsurface hydrology in the creekbank zone of tidal marshes

A combination of field and numerical modeling methods was used to assess porewater movement in a narrow (20 m) Spartina marsh which was flooded regularly by tidal waters. Soil composition and soil hydraulic properties did not vary across the marsh or with depth. Hydraulic head was monitored on a transect perpendicular to the creekbank. During exposure of the marsh surface, hydraulic gradients were predominantly horizontal; vertical gradients were small or zero. Subsurface flow was directed from the marsh interior toward the creekbank. Approximately 141 of pore water were discharged laterally to the adjacent tidal creek per meter of creekbank over a complete tidal cycle.A numerical hydrological model was modified to simulate subsurface hydraulics in the creekbank vicinity of regularly flooded tidal marshes. The model was parameterized to represent soil conditions, tidal fluctuations and topography at the field site. Observed changes in hydraulic head over complete tidal cycles were accurately predicted by the model. Model simulations identified the vertical infiltration of creek water into the marsh surface at the onset of tidal flooding as the primary source (66%) for the replacement of water drained at the creekbank. Significant replacement (31%) also occurred as discharge from the interior marsh. Horizontal recharge at the creekbank was minimal (3%).A sensitivity analysis was conducted with the model to assess the relative importance of geomorphological factors and soil properties in controlling pore water export at the creekbank of tidal marsh soils. Each parameter was varied systematically over a realistic range for field conditions. Changes in marsh elevation exerted greater control over creekbank discharge than changes in soil hydraulic properties. More rapid turnover of pore water near creekbanks of higher elevation marshes is hypothesized.     

Impact of lithology and soil properties on abandoned dryland terraces during the early stages of soil erosion by water in south‐east Spain

Soil erosion by water in abandoned dry terraces is one of the most important environmental problems in semiarid areas, enhancing biological degradation and reducing possible resources that can be obtained. However, little is known about the effects of the types of lithology and soil properties on the early stages of soil erosion. Therefore, the main aim of this research was to assess the effect of different lithologies (marls, limestones, and metamorphic—phyllites, schists, and greywackes—materials) and soil properties on the early stages of soil erosion by water in abandoned dry terraces, compared with similar terraces still in agricultural use. Soil analyses (texture, aggregate stability, and bulk density) and 22 rainfall simulations were carried out under dry conditions. During the experiments, local inclination, vegetation and stone cover, total organic matter, and antecedent soil moisture were also quantified. The results showed that the highest soil loss (41.41 g/m² in cultivated plots and 17.05 g/m² in the abandoned plots) and runoff (3.79 L/m² in the abandoned plot) occurred on marl substrata. Marls also showed the shallowest infiltration front (9 cm) and lowest infiltration rate (4.3 cm/min). Limestones and, especially, metamorphic areas, showed a lower degree of soil erosion, higher infiltration rates, and deeper infiltration fronts.     

A discrete element approach to elastic properties of layered geomaterials

Discrete element methods (DEMs) are used for layered geomaterials to investigate the dependency of traditional engineering constants on material properties and loading conditions. Shear deformations and compression tests parallel and perpendicular to layering are conducted on samples of varying kerogen volume fractions, confining pressures, porosities, and layer geometries. The goal of this article is to develop a method to better characterize oil shale (a transversely isotropic layered geomaterial) while eliminating high experimental costs. The DEM simulations conducted in this study demonstrate strong dependencies of Young's modulus, Poisson's ratio, and shear modulus on kerogen volume fraction and porosity. Furthermore, a rule of thumb for layer thickness and particle resolution is proposed for simulation design. Results agree well with robust effective medium theories, solidify the ability of DEM to model the mechanical properties of layered heterogenous materials, and encourage the use of DEM to study more complicated layered media and material failure. Copyright © 2011 John Wiley & Sons, Ltd.     

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.     

Application of the exact constitutive relationship of modified Cam clay to the undrained expansion of a spherical cavity

The modified Cam clay (MCC) model is used to study the response of virgin‐compressed clay subjected to undrained triaxial compression. The MCC constitutive relationship is obtained in a closed form. Both elastic and plastic deviatoric strains are considered in the analysis. The solution allows to obtain total and effective stress paths followed by the clay in undrained spherical expansion. Pore water pressures are determined from the difference between total and effective mean stresses. For illustration purposes, the analysis is also applied to the well‐known reconstituted normally consolidated London clay and the results are compared with the recently published data obtained by a numerical approach. In addition, the Almansi large strains are used in the analysis, as these allow to obtain limit expansion and pore pressures, whereas both small‐strain and logarithmic‐strain approaches do not permit to determine them. Copyright © 2010 John Wiley & Sons, Ltd.     

Using estimated thermodynamic properties to model accessory phases: the case of tourmaline

Accessory phases and minor components in minerals are commonly ignored in thermodynamic modelling. Such an approach seems unwarranted, as accessory phases can represent a significant element reservoir and minor components can substantially change their host mineral's stability field. However, a lack of thermodynamic data prohibits assessment of these effects. In this contribution, the polyhedron method is used to estimate the thermodynamic properties of tourmaline, a common and widespread accessory phase, stable over a range of P–T–X conditions. The polyhedron method allows Δ H , S , V , C _P and V _m ( T , P ) properties to be estimated from a linear stoichiometric summation over the fractional properties of its polyhedron constituents. To allow for estimates of tourmaline, fractional thermodynamic properties for B^III and B^IV polyhedra were derived. Mixing contributions to molar volume were evaluated and symmetrical mixing parameters derived for Al-Mg, Al-Fe and Al-Li interaction on tourmaline's Y-site and T-site Al-Si interaction. Evaluation of the estimated properties using experimental and natural equilibria between tourmaline and melts, minerals and hydrothermal fluids, shows that reliable semi-quantitative results are obtained. The boron contents in fluids coexisting with tourmaline are calculated to within an order of magnitude of measured content, and where anchor-points are available, agreement improves to within a factor of 2. Including tourmaline in petrogenetic modelling of metamorphic rocks indicates that its presence leads to disappearance of staurolite and garnet, among others, and modifies the X _Mg of coexisting phases, in line with observations on natural rocks.     

Elastic settlements under linear surface pressures on rectangular areas

This paper presents the practical closed‐form solutions for elastic settlement under the linear full‐contact pressures on rectangular areas resting upon an elastic mass. The linear pressure distribution is mostly assumed in practice for the rectangular footing having biaxial bending. The equations presented here are determined by evaluating the integration of Steinbrenner's strain approximation based on the stress solutions from the Boussinesq equations. The presented formulae are validated to be used for the elastic settlement under any point of linear full‐contact loading, triangular and embankment types of pressure cases. In this respect, it represents a significant step forward in the understanding of elastic settlement and rotation under the practical loading cases. The presented solutions are concise and easy to use. The numerical examples are demonstrated and the resulting influence factors in graphical forms are presented for their practical use especially in the elastic settlement calculations of rectangular footing subjected to the biaxial bending. Copyright © 2008 John Wiley & Sons, Ltd.     

An alternative interpretation for the formation of doubly plunging folds in sandstone terrains

Folds are marvellous features of mountain terrains, but despite extensive research, many fundamental problems have still not been solved. In terrains of sandstone, fold hinges are rarely straight lines but curved, forming a pattern characterized by doubly plunging, elliptical dome‐and‐basin structures. Such structures are an obvious manifestation of coeval or successive shortening deformation in two orthogonal principal directions in the horizontal plane. Based on an anatomic investigation of the fold pattern of sandstone beds at the rocky beaches of Saint‐Jean‐Port‐Joli (Quebec, Canada), we propose that the doubly plunging folds may result from a transition from a plane deformation to a constrictional deformation due to auxetic effects of quartz‐rich rocks. The sandstone beds possessed potentially such negative values of Poisson's ratio that, when placed under compression in one direction, they become contracted in the transverse direction, producing a series of doubly plunging folds. Further work is needed to approve or disapprove the interpretation.     

Controls on evapotranspiration from jack pine forests in the Boreal Plains Ecozone

The exchanges of water, energy and carbon between the land surface and the atmosphere are tightly coupled, so that errors in simulating evapotranspiration lead to errors in simulating both the water and carbon balances. Areas with seasonally frozen soils present a particular challenge due to the snowmelt-dominated hydrology and the impact of soil freezing on the soil hydraulic properties and plant root water uptake. Land surface schemes that have been applied in high latitudes often have reported problems with simulating the snowpack and runoff. Models applied at the Boreal Ecosystem Research and Monitoring Sites in central Saskatchewan have consistently over-predicted evapotranspiration as compared with flux tower estimates. We assessed the performance of two Canadian land surface schemes (CLASS and CLASS-CTEM) for simulating point-scale evapotranspiration at an instrumented jack pine sandy upland site in the southern edge of the boreal forest in Saskatchewan, Canada. Consistent with past reported results, these models over-predicted evapotranspiration, as compared with flux tower observations, but only in the spring period. Looking systematically at soil properties and vegetation characteristics, we found that the dominant control on evapotranspiration within these models was the canopy conductance. However, the problem of excessive spring ET could not be solved satisfactorily by changing the soil or vegetation parameters. The model overestimation of spring ET coincided with the overestimation of spring soil liquid water content. Improved algorithms for the infiltration of snowmelt into frozen soils and plant-water uptake during the snowmelt and soil thaw periods may be key to addressing the biases in spring ET.     

The impact of flow discharge on the hydraulic characteristics of headcut erosion processes in the gully region of the Loess Plateau

Headcut erosion is associated with major hydraulic changes induced by the gully head of concentrated flow. However, the variation in the hydraulic characteristics of the headcut erosion process is still not clear in the gully region of the Loess Plateau. A series of rainfall combined scouring experiments (flow discharges ranging from 3.6 to 7.2 m³ hr⁻¹, with 0.8 mm min⁻¹ rainfall intensity) were conducted on experimental plots to clarify the variation in the hydraulic parameters induced by gully head and erosion processes under different flow discharges. The results showed that concentrated flows in the catchment area and gully bed were turbulent (Reynolds number ranging from 1,876 to 6,693) and transformed between supercritical and subcritical (Froude number ranging from 0.96 to 3.73). The hydraulic parameters, such as the flow velocity, Reynolds number, shear stress, stream power, Darcy–Weisbach friction factor, and unit stream power in the catchment area were 0.45–0.59 m s⁻¹, 2086–6693, 1.96–5.33 Pa, 0.89–2.86 W m⁻², 0.08–0.16, and 0.023–0.031 m s⁻¹, respectively. When the concentrated flows dropped from the gully head, the hydraulic parameters in the gully bed decreased by 3.39–26.07%, 1.49–29.99%, 65.19–67.14%, 67.25–74.96%, 28.53–61.31%, and 67.82–77.14%, respectively, which contributed to the flow energy consumption at the gully head. As flow discharge increased, Reynolds number, shear stress, and stream power increased, while flow velocity, Froude number, unit stream power, and Darcy–Weisbach friction factor did not. The flow energy consumption at the gully head was 9.66–10.13, 13.25–13.74, 15.68–16.41, and 19.28–20.25 J s⁻¹, respectively, under different flow discharges and accounted for 60.58–68.50% of the flow energy consumption of the experimental plots. Generally, the sediment discharges increased rapidly at the initial stage, then increased slowly, and finally reached a steady state condition, which showed a significant declining logarithmic trend with experimental duration (P<.01) and increased with increasing flow discharge. Accordingly, the flow energy consumption was significantly correlated with the sediment yield. These findings could improve our understanding of the hydraulic properties and flow energy characteristics of headcut erosion.