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.     

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.     

Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.     

Thermal history of the Hodgkinson Province and Laura Basin,Queensland: Multiple cooling episodes identified from apatite fission track analysis and vitrinite reflectance data

Apatite fission track analysis and vitrinite reflectance data from outcrop and well samples in the Hodgkinson Province and Laura Basin reveal regional Cretaceous cooling. Apatite fission track analysis appears to define two discrete cooling episodes, in the mid‐Cretaceous (110–100 Ma) and Late Cretaceous (80–70 Ma), although in most samples data allow only definition of a single episode. Rocks now at outcrop cooled from Cretaceous palaeotemperatures generally between 50 and 130°C in the south of the region, and from >100°C in the north. Some samples from the Hodgkinson Province also show evidence for an Early Jurassic cooling episode, characterised by maximum palaeotemperatures varying from at least 95°C (from apatite fission track analysis) to ～200–220°C (from vitrinite reflectance), with cooling beginning at around 200 Ma. Apatite fission track analysis data do not reveal the earlier event in the Laura Basin, but on the basis of vitrinite reflectance data from Permian? units this event is also inferred to have affected the pre‐Jurassic basin units in this region. The regional extent of the Cretaceous cooling episode in the Hodgkinson Province suggests that the elevated palaeotemperatures in this region were most likely due to greater depth of burial, with subsequent cooling due to kilometre‐scale denudation. For a palaeogeothermal gradient of 30°C/km and a palaeosurface temperature of 25°C the total degree of Cretaceous cooling experienced by the samples corresponds to removal of between ～0.8 and >3.0 km of Triassic and younger section removed by denudation, beginning some time between ca 110 and 80 Ma. Higher palaeogradients would require correspondingly lower amounts of removed section. The geology of the Laura Basin suggests that an explanation of the observed Cretaceous palaeotemperatures in this region in terms of deeper burial may be untenable. Heating due to hot fluid flow may be a more realistic mechanism for producing the observed Cretaceous palaeothermal effects in the Laura Basin.     

Structural kinematics,metamorphic P–T profiles and zircon geochronology across the Greater Himalayan Crystalline Complex in south‐central Tibet: implication for a revised channel flow

A specific question about the Himalayas is whether the orogeny grew by distributed extrusion or discrete thrusting. To place firm constraints on tectonic models for the orogeny, kinematic, thermobarometric and geochronological investigations have been undertaken across the Greater Himalayan Crystalline Complex (GHC) in the Nyalam region, south‐central Tibet. The GHC in this section is divided into the lower, upper and uppermost GHC, with kinematically top‐to‐the‐south, alternating with top‐to‐the‐north shear senses. A new thrust named the Nyalam thrust is recognized between the lower and upper GHC, with a 3 kbar pressure reversion, top‐to‐the‐south thrust sense, and was active after the exhumation of the GHC. Peak temperature reached ～749 °C in the cordierite zone, and decreased southwards to 633–667 °C in the kyanite and sillimanite‐muscovite zones, and northwards to greenschist facies at the top of the South Tibetan Detachment System (STDS). Pressure at peak temperature reached a maximum value in the kyanite zone of 9.0–12.6 kbar and decreased northwards to ～4.1 kbar in the cordierite zone. Zircon U‐Pb ages of a sillimanite migmatite and an undeformed leucogranite dyke cutting the mylonitized rocks in the STDS reveal a long‐lived partial melting of the GHC, which initiated at 39.7–34 Ma and ceased at 14.1 Ma. Synthesizing the obtained and collected results, a revised channel flow model is proposed by considering the effect of heat advection and convection by melt and magma migration. Our new model suggests that distributed processes like channel flow dominated during the growth of the Himalayan orogen, while discrete thrusting occurred in a later period as a secondary process.     

One-dimensional dispersion in unsteady flow in an adsorbing porous medium: An analytical solution

Analytical solutions describing the concentration distribution along one-dimensional unsteady seepage flow through adsorbing saturated finite porous medium have been obtained. An exponential function concentration is enforced at the source of the dispersion, while the change in the concentration is zero at the other boundary. A new time variable has been introduced to solve the unsteady flow problem and the solution is illustrated.     

State of the art of modelling two-phase flow in fractured rock

 Two-phase flow in fractured rock is an important phenomenon for a range of practical problems, not the least of which is non-aqueous phase liquid (NAPL) contamination of groundwater. Although multiphase systems have long been studied in the petroleum field, in the hydrogeological field progress has only just reaches the point where models are being developed. Scale effect is one of the main issues of concern. Although models presented in this paper have the potential to provide useful predictions, they can only be used to investigate a variety of possible scenarios with parameters being specified in the form of distribution of values. The calibration and validation of all but the simplest of these models poses a formidable task, with great demands on hydrogeologists and geophysicists to provide adequate data. Received: 20 May 1996 · Accepted: 19 August 1996     

The effect of terrace geology on ground‐water movement and on the interaction of ground water and surface water on a mountainside near Mirror Lake,New Hampshire,USA

The west watershed of Mirror Lake in the White Mountains of New Hampshire contains several terraces that are at different altitudes and have different geologic compositions. The lowest terrace (FSE) has 5 m of sand overlying 9 m of till. The two next successively higher terraces (FS2 and FS1) consist entirely of sand and have maximum thicknesses of about 7 m. A fourth, and highest, terrace (FS3) lies in the north‐west watershed directly adjacent to the west watershed. This highest terrace has 2 m of sand overlying 8 m of till. All terraces overlie fractured crystalline bedrock. Numerical models of hypothetical settings simulating ground‐water flow in a mountainside indicated that the presence of a terrace can cause local ground‐water flow cells to develop, and that the flow patterns differ based on the geologic composition of the terrace. For example, more ground water moves from the bedrock to the glacial deposits beneath terraces consisting completely of sand than beneath terraces that have sand underlain by till. Field data from Mirror Lake watersheds corroborate the numerical experiments. The geology of the terraces also affects how the stream draining the west watershed interacts with ground water. The stream turns part way down the mountainside and passes between the two sand terraces, essentially transecting the movement of ground water down the valley side. Transects of water‐table wells were installed across the stream's riparian zone above, between, and below the sand terraces. Head data from these wells indicated that the stream gains ground water on both sides above and below the sand terraces. However, where it flows between the sand terraces the stream gains ground water on its uphill side and loses water on its downhill side. Biogeochemical processes in the riparian zone of the flow‐through reach have resulted in anoxic ground water beneath the lower sand terrace. Results of this study indicate that it is useful to understand patterns of ground‐water flow in order to fully understand the flow and chemical characteristics of both ground water and surface water in mountainous terrain. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.