Modelling non‐linear ground response of non‐liquefiable soils

A non‐linear finite element (FE) model is presented to account for soil column effects on strong ground motion. A three‐dimensional bounding surface plasticity model with a vanishing elastic region, appropriate for non‐liquefiable soils, is formulated to accommodate the effects of plastic deformation right at the onset of loading. The elasto‐plastic constitutive model is cast within the framework of a FE soil column model, and is used to re‐analyse the downhole motion recorded by an array at a Large‐Scale Seismic Test (LSST) site in Lotung, Taiwan, during the earthquake of 20 May 1986; as well as the ground motion recorded at Gilroy 2 reference site during the Loma Prieta earthquake of 17 October 1989. Results of the analysis show maximum permanent shearing strains experienced by the soil column in the order of 0.15 per cent for the Lotung event and 0.8 per cent for the Loma Prieta earthquake, which correspond to modulus reduction factors of about 30 and 10 per cent respectively, implying strong non‐linear response of the soil deposit at the two sites. Copyright © 2000 John Wiley & Sons, Ltd.     

Influence of river channel geometry in stream flow modelling and guidelines for field investigation

Fully physics‐based, process‐level, distributed fluid flow and reactive transport hydrological models are rarely used in practice until recent years. These models are useful tools to help understand the fundamental physical, chemical, and biological processes that take place in nature. In this study, sensitivity analyses based on a mountain area river basin modelling study are performed to investigate the effect of river channel geometric characteristics on downstream water flow. Numerical experiments show that reduction in the river channel geometric measurement interval may not significantly affect the downstream water stage simulation as long as measurement accuracy at special nodes is guaranteed. The special upstream nodes include but are not limited to 1) nodes located close to the observation station, 2) nodes near the borders of different land covers with considerable riverbed roughness changes, 3) nodes at entering points of tributaries causing discharge jump and 4) nodes with a narrow cross‐section width that may control the flow conditions. This information provides guidelines for field investigation to efficiently obtain necessary geometric data for physics‐based hydrological modelling. It is especially useful in alpine areas such as the Tibetan Plateau where field investigation capability is limited under severe topography and climate condition. Copyright © 2013 John Wiley & Sons, Ltd.     

Bank erosion and channel width change in a tropical catchment

Catchment sediment budget models are used to predict the location and rates of bank erosion in tropical catchments draining to the Great Barrier Reef lagoon, yet the reliability of these predictions has not been tested due to a lack of measured bank erosion data. This paper presents the results of a 3 year field study examining bank erosion and channel change on the Daintree River, Australia. Three different methods were employed: (1) erosion pins were used to assess the influence of riparian vegetation on bank erosion, (2) bench‐marked cross‐sections were used to evaluate annual changes in channel width and (3) historical aerial photos were used to place the short term data into a longer temporal perspective of channel change (1972–2000). The erosion pin data suggest that the mean erosion rate of banks with riparian vegetation is 6·5 times (or 85%) lower than that of banks without riparian vegetation. The changes measured from cross‐section surveys suggest that channel width has increased by an average of 0·74 (±0·47) m a^?1 over the study period (or ～0·8% yr^?1). The aerial photo results suggest that over the last 30 years the Daintree River has undergone channel contraction of the order of 0·25 m a^?1. The cross‐section data were compared against modelled SedNet bank erosion rates, and it was found that the model underestimated bank erosion and was unable to represent the variable erosion and accretion processes that were observed in the field data. The reach averaged bank erosion rates were improved by the inclusion of locally derived bed slope and discharge estimates; however, the results suggest that it will be difficult for catchment scale sediment budget models to ever accurately predict the location and rate of bank erosion due to the variation in bank erosion rates in both space and time. Copyright © 2008 John Wiley & Sons, Ltd.     

Quantification of peatland water storage capacity using the water table fluctuation method

Peat specific yield (S_Y) is an important parameter involved in many peatland hydrological functions such as flood attenuation, baseflow contribution to rivers, and maintaining groundwater levels in surficial aquifers. However, general knowledge on peatland water storage capacity is still very limited, due in part to the technical difficulties related to in situ measurements. The objectives of this study were to quantify vertical S_Y variations of water tables in peatlands using the water table fluctuation (WTF) method and to better understand the factors controlling peatland water storage capacity. The method was tested in five ombrotrophic peatlands located in the St. Lawrence Lowlands (southern Québec, Canada). In each peatland, water table wells were installed at three locations (up‐gradient, mid‐gradient, and down‐gradient). Near each well, a 1‐m long peat core (8 cm × 8 cm) was sampled, and subsamples were used to determine S_Y with standard gravitational drainage method. A larger peat sample (25 cm × 60 cm × 40 cm) was also collected in one peatland to estimate S_Y using a laboratory drainage method. In all sites, the mean water table depth ranged from 9 to 49 cm below the peat surface, with annual fluctuations varying between 15 and 29 cm for all locations. The WTF method produced similar results to the gravitational drainage experiments, with values ranging between 0.13 and 0.99 for the WTF method and between 0.01 and 0.95 for the gravitational drainage experiments. S_Y was found to rapidly decrease with depth within 20 cm, independently of the within‐site location and the mean annual water table depth. Dominant factors explaining S_Y variations were identified using analysis of variance. The most important factor was peatland site, followed by peat depth and seasonality. Variations in storage capacity considering site and seasonality followed regional effective growing degree days and evapotranspiration patterns. This work provides new data on spatial variations of peatland water storage capacity using an easily implemented method that requires only water table measurements and precipitation data.     

Influence of shrub canopy morphology and rainfall characteristics on stemflow within a revegetated sand dune in the Tengger Desert,NW China

Stemflow of xerophytic shrubs was monitored on event basis within a revegetated sand dune. Quantity of stemflow showed a clear species‐specific dependence in combination with the rainfall characteristics. Results obtained revealed that for ovate‐leaved C. korshinskii with an inverted cone‐shaped canopy and smooth bark, the quantity of stemflow in depth accounted for 7.2% of the individual gross rainfall, while it was 2.0% for needle‐leaved A. ordosica with a cone‐shaped canopy and coarse bark. There were significant positive linear relationships between stemflow and individual gross rainfall and rainfall intensity for the two shrubs. An individual gross rainfall of 1.4 and 1.8 mm was necessary for stemflow generation for C. korshinskii and A. ordosica, respectively. Multiple regression analysis showed that the abiotic and biotic variables including the individual gross rainfall, mean windspeed (WS), canopy height, branch length, and canopy volume have significant influence on stemflow for C. korshinskii, whereas for A. ordosica, the notable influencing variables were individual gross rainfall, stem diameter, and leaf area index. Generally, WS has less effect on stemflow than that of rainfall for A. ordosica. The correlation relationship between individual gross rainfall and funneling ratio showed that the funneling ratio attains its peak when the gross rainfall is 13 and 16 mm for C. korshinskii and A. ordosica, respectively, implying that the canopy morphology emerged as determining factors on funneling ratio decrease when the individual gross rainfall exceeds these values. In comparison, higher WS increased the funneling ratio remarkably for C. korshinskii than A. ordosica due partly to the greater branch length and canopy projection area in C. korshinskii. Funneling ratio can be used as an integrated variable for the effects of canopy morphology and rainfall characteristics on stemflow. The implication of stemflow on water balance and its contribution to sustain the shrubs and the revegetation efforts was discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Testing the genetic relationship between fluid alteration and brecciation in CM chondrites

Boriskino is a poorly studied CM chondrite with numerous millimeter‐ to centimeter‐scale clasts exhibiting sharp boundaries. Clast textures and mineralogies attest to diverse geological histories with various degrees of aqueous alteration. We conducted a petrographic, chemical, and isotopic study on each clast type of the breccia to investigate if there exists a genetic link between brecciation and aqueous alteration, and to determine the controlling parameter of the extent of alteration. Boriskino is dominated by CM2 clasts for which no specific petrographic type could be assigned based on the chemical compositions and modal abundances of constituents. One clast stands out and is identified as a CM1 lithology, owing to its lack of anhydrous silicates and its overall abundance of dolomite‐like carbonates and acicular iron sulfides. We observe that alteration phases near clast boundaries exhibit foliation features, suggesting that brecciation postdated aqueous alteration. We measured the O‐isotopic composition of Ca‐carbonates and dolomite‐like carbonates to determine their precipitation temperatures following the methodology of Verdier‐Paoletti et al. (2017). Both types of carbonates yield similar ranges of precipitation temperatures independent of clast lithology, ranging from ?13.9 ± 22.4 (2σ) to 166.5 ± 47.3 °C, precluding that temperature alone accounts for the differences between the CM1 and CM2 lithologies. Instead, we suggest that initial water/rock ratios of 0.75 and 0.61 for the CM1 and CM2 clasts, respectively, might control the extent of aqueous alteration. Based on these estimates, we suggest that Boriskino clasts originated from a single parent body with heterogeneous distribution of water either due to local differences in the material permeability or in the initial content of ice available. These conditions would have produced microenvironments with differing geochemical conditions thus leading to a range of degrees of aqueous alteration.     

Determination of synthetic wind functions for estimating open water evaporation with Computational Fluid Dynamics

A computational fluid dynamics (CFD)‐based methodology is proposed to derive convective mass‐transfer coefficients (wind functions) that are required for estimating evaporation of water bodies with the mass‐transfer method. Three‐dimensional CFD was applied to model heat transfer in two water bodies: a Class‐A tank evaporimeter and an on‐farm artificial pond. The standard k–? model assuming isotropic turbulence was adopted to describe turbulent heat transport, whereas the heat and mass transfer analogy was assumed to derive the wind functions. The CFD‐derived wind functions were very similar to those empirically derived from the experimental water bodies. The evaporation rates calculated with the synthetic wind functions were in good agreement with hourly and daily evaporation measurements for the tank and pond, respectively. The proposed CFD‐approach is generalisable and cost effective, because it has low input data requirements. Besides, it provides additional capability of modelling the spatial distribution of the evaporation rate over the water surface. Although the application of CFD to water bodies evaporation modelling is still in development, it looks very promising. Copyright © 2012 John Wiley & Sons, Ltd.     

The linkage between fault throw and footwall scarp erosion patterns: an example from the Bremstein Fault Complex,offshore Mid‐Norway

Studies of normal fault systems in modern extensional regimes (e.g. Basin and Range), and in exhumed, ancient rift basins (e.g. Gulf of Suez Rift) have shown a link between the evolution of fault‐related footwall topography and associated erosional drainage systems. In this study, we use 3D seismic reflection data to image the footwall crest of a gravity‐driven fault system developed during late Middle Jurassic to Early Cretaceous rifting on the Halten Terrace, offshore Mid‐Norway. This 22‐km‐long fault system lacks significant footwall uplift, with hangingwall subsidence accommodating throw accumulation on the fault system. Significant erosion has occurred along the length of the footwall crest and is defined by 96 catchments characterized by erosional channels. These erosional channels consist of small, linear systems up to 750 m long located along the front of the fault footwall. Larger, dendritic channel systems extend further back (up to 3 km normal to fault strike) into the footwall. These channels are up to 7 km long, up to 50 m deep and up to 1 km wide. Fault throw varies along strike, with greatest throw in the centre of the fault decreasing towards the fault tips; localized throw minima are interpreted to represent segment linkage points, which were breached as the fault grew. Comparison of the catchment location to the throw distribution shows that the largest catchments are in the centre of the fault and decrease in size to the fault tips. There is no link between the location of the breached segment linkage points and the location and size of the footwall catchments, suggesting that the first‐order control on footwall erosion patterns is the overall fault‐throw distribution.     

Discrete simulations of a triaxial compression test for sand by DEM

A quasi‐static homogeneous drained triaxial compression test on cohesionless sand under constant lateral pressure was simulated using a three‐dimensional discrete element method. Grains were modelled by means of particle clusters composed of rigid spheres or spheres with contact moments imitating irregular particle shapes. Attention was paid to the effect of initial void ratio and grain shape mixture on the shear strength, volume changes, force chains, kinetic, elastic and dissipated energies. In addition, the effect of the mean grain size, grain size distribution, grain size range, specimen size and roughness and stiffness of boundaries was numerically analysed in initially dense sand. Some numerical results were compared with available experimental results. Copyright © 2014 John Wiley & Sons, Ltd.     

Removal of Volatile Organic Compounds in Vertical Flow Filters: Predictions from Reactive Transport Modeling

Vertical flow filters are containers filled with porous medium that are recharged from top and drained at the bottom, and are operated at partly saturated conditions. They have recently been suggested as treatment technology for groundwater containing volatile organic compounds (VOCs). Numerical reactive transport simulations were performed to investigate the relevance of different filter operation modes on biodegradation and/or volatilization of the contaminants and to evaluate the potential limitation of such remediation mean due to volatile emissions. On the basis of the data from a pilot‐scale vertical flow filter intermittently fed with domestic waste water, model predictions on the system’s performance for the treatment of contaminated groundwater were derived. These simulations considered the transport and aerobic degradation of ammonium and two VOCs, benzene and methyl tertiary butyl ether (MTBE). In addition, the advective‐diffusive gas‐phase transport of volatile compounds as well as oxygen was simulated. Model predictions addressed the influence of depth and frequency of the intermittent groundwater injection, degradation rate kinetics, and the composition of the filter material. Simulation results show that for unfavorable operation conditions significant VOC emissions have to be considered and that operation modes limiting VOC emissions may limit aerobic biodegradation. However, a suitable combination of injection depth and composition of the filter material does facilitate high biodegradation rates while only little VOC emissions take place. Using such optimized operation modes would allow using vertical flow filter systems as remediation technology suitable for groundwater contaminated with volatile compounds.