Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

Cohesive‐frictional model for bond and splitting action of prestressing wire

This paper presents a numerical procedure for bond between indented wires and concrete, and the coupled splitting process of the surrounding concrete. The bond model is an interface, non‐associative, plasticity model. It is coupled with a cohesive fracture model for concrete to take into account the splitting of such concrete. Bond between steel and concrete is fundamental for the transmission of stresses between both materials in precast prestressed concrete. Indented wires are used to improve the bond in these structural elements. The radial component of the prestressing force, increased by Poisson's effect, may split the surrounding concrete, decreasing the wire confinement and diminishing the bonding. The combined action of the bond and the splitting is studied with the proposed model. The results of the numerical model are compared with the results of a series of tests, such as those which showed splitting induced by the bond between wire and concrete. Tests with different steel indentation depths were performed. The numerical procedure accurately reproduces the experimental records and improves knowledge of this complex process. Copyright © 2010 John Wiley & Sons, Ltd.     

Evolution of scouring process downstream of grade‐control structures under steady and unsteady flows

For many incised channels, one of the most common strategies is to install some hard structures, such as grade‐control structures (GCSs), in the riverbed to resist further incision. In this study, a series of experiments, including both steady and unsteady flow conditions, were conducted to investigate the scouring process downstream of a GCS. Three distinct phases, including the initial, developing and equilibrium phases, during the evolution of scour holes were identified. In addition, a semi‐empirical method was proposed to predict the equilibrium scour‐hole profile for the scour countermeasure design. In general, the comparisons between the experimental and simulated results are reasonably consistent. As the studies on temporal variation of the scour depth at GCSs caused by floods are limited, the effect of flood hydrograph shapes on the scour downstream of GCSs without upstream sediment supply was also investigated experimentally in this study. Based on the dimensional analysis and the concept of superposition, a methodology is proposed to simulate the time evolution of the maximum scour depth downstream of a GCS for steady flows. Moreover, the proposed scheme predicts reasonably well the temporal variations of the maximum scour depth for unsteady flows with both single and multiple peak. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil saturation index of salt marshes subjected to spring‐neap tides: a new variable for describing marsh soil aeration condition

Organized spatial distribution of plants (plant zonation) in salt marshes has been linked to the soil aeration condition in the rhizosphere through simplistic tidal inundation parameters. Here, a soil saturation index (ratio of saturation period to tidal period at a soil depth) is introduced to describe the soil aeration condition. This new index captures the effects of not only the tidal inundation period and frequency but also the flow dynamics of groundwater in the marsh soil. One‐dimensional numerical models based on saturated flow with the Boussinesq approximations and a two‐dimensional variably saturated flow model were developed to explore the behaviour of this new soil aeration variable under the influence of spring‐neap tides. Simulations revealed two characteristic zones of soil aeration across the salt marsh: a relatively well aerated near‐creek zone and a poorly aerated interior zone. In the near‐creek zone, soils undergo periodic wetting and drying as the groundwater table fluctuates throughout the spring‐neap cycle. In the interior, the soil remains largely water saturated except for neap tide periods when limited drainage occurs. Although such a change of soil aeration condition has been observed in previous numerical simulations, the soil saturation index provides a clear delineation of the zones that are separated by an ‘inflexion point’ on the averaged index curve. The results show how the saturation index represents the effects of soil properties, tidal parameters and marsh platform elevation on marsh soil aeration. Simulations of these combined effects have not been possible with traditional tidal inundation parameters. The saturation index can be easily derived using relatively simple models based on five non‐dimensional variables. Copyright © 2010 John Wiley & Sons, Ltd.     

Study on low‐strain integrity testing of pipe‐pile using the elastodynamic finite integration technique

This study focuses on the three‐dimensional (3‐D) characteristics of wave propagation in pipe‐pile using elastodynamic finite integration technique. First, a real 3‐D pile‐soil model in cylindrical coordinate system is presented. Then, the governing equations are established. With the boundary and initial conditions, the numerical solution is obtained. The accuracy and feasibility of the self‐written code are further verified via comparing with the measured data. Velocity histories at different angles of pile top and pile tip are illustrated, and the snapshots reflecting the 3‐D characteristics of wave propagation are given. It shows that the interferences of Rayleigh waves can confuse the result interpretation for pile integrity testing. The increase of hammer contact time can effectively mitigate the interferences, and the interferences of Rayleigh waves are weakest at an angle of 90° from where hammer hits. Besides, surrounding soil can partly mitigate the wave interferences. Copyright © 2012 John Wiley & Sons, Ltd.     

Effective block diagonal preconditioners for Biot's consolidation equations in piled‐raft foundations

The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd.     

Copper‐gold mineralisation in New Guinea: Numerical modelling of collision,fluid flow and intrusion‐related hydrothermal systems

Two‐ and three‐dimensional numerical modelling techniques, constrained by key geodynamic data, provide insights into the controls on development of porphyry‐related Cu–Au mineralisation in the Tertiary collision zone of New Guinea. Modelling shows that the creation of local dilation to facilitate magma emplacement can be caused by reactivation of arc‐normal transfer faults, where they cut the weakened fold belt. Additionally, dilation occurs where fluid overpressuring is caused by collision‐related, south‐directed fluid flow being localised into the more permeable units of the Mesozoic passive‐margin sedimentary succession. Rapid uplift and erosion, which may be a mechanism for magmatic fluid release in these systems, is shown to be greatest in the west of West Papua, where the stronger Australian crust acts as a buttress. Within the Papuan Fold Belt, uplift is greatest near the margins, where the weaker fold belt abuts the stronger crust and/or major faults have been reactivated. Increased orographically induced precipitation and erosion exposes the lower parts of the stratigraphy within or on the margins of these uplifted zones. On a smaller scale, 2–D coupled fluid‐flow ‐ thermal‐chemical modelling uses a scenario of fluid mixing to calculate metal precipitation distribution and magnitude around an individual intrusive complex. Modelling highlights the interdependence of the spatial permeability structure, the regional temperature gradient, and the geometry of the convection cells and how this impacts on the distribution of metal precipitation.     

A mass‐balance approach to estimate in‐stream processes in a large river

A mass‐balance approach was used to estimate in‐stream processes related to inorganic nitrogen species (NH₄⁺, NO₂^? and NO₃^?) in a large river characterized by highly variable hydrological conditions, the Garonne River (south‐west France). Studies were conducted in two consecutive reaches of 30 km located downstream of the Toulouse agglomeration (population 760 000, seventh order), impacted by modification of discharge regime and high nitrogen concentrations. The mass‐balance was calculated by two methods: the first is based on a variable residence time (VRT) simulated by a one‐dimensional (1‐D) hydraulic model; the second is a based on a calculation using constant residence time (CRT) evaluated according to hydrographic peaks. In the context of the study, removal of dissolved inorganic nitrogen (DIN) for a reach of 30 km is underestimated by 11% with the CRT method. In sub‐reaches, the discrepancy between the two methods led to a 50% overestimation of DIN removal in the upper reach (13 km) and a 43% underestimation in the lower reach (17 km) using the CRT method. The study highlights the importance of residence time determination when using modelling approaches in the assessment of whole stream processes in short‐duration mass‐balance for a large river under variable hydrological conditions. Copyright © 2007 John Wiley & Sons, Ltd.