Modeling plant transpiration under limited soil water: comparison of different plant and soil hydraulic parameterizations and preliminary implications for their use in land surface models

Accurate estimates of how soil water stress affects plant transpiration are crucial for reliable land surface model (LSM) predictions. Current LSMs generally use a water stress factor, β, dependent on soil moisture content, θ, that ranges linearly between β = 1 for unstressed vegetation and β = 0 when wilting point is reached. This paper explores the feasibility of replacing the current approach with equations that use soil water potential as their independent variable, or with a set of equations that involve hydraulic and chemical signaling, thereby ensuring feedbacks between the entire soil–root–xylem–leaf system. A comparison with the original linear θ-based water stress parameterization, and with its improved curvi-linear version, was conducted. Assessment of model suitability was focused on their ability to simulate the correct (as derived from experimental data) curve shape of relative transpiration versus fraction of transpirable soil water. We used model sensitivity analyses under progressive soil drying conditions, employing two commonly used approaches to calculate water retention and hydraulic conductivity curves. Furthermore, for each of these hydraulic parameterizations we used two different parameter sets, for 3 soil texture types; a total of 12 soil hydraulic permutations. Results showed that the resulting transpiration reduction functions (TRFs) varied considerably among the models. The fact that soil hydraulic conductivity played a major role in the model that involved hydraulic and chemical signaling led to unrealistic values of β, and hence TRF, for many soil hydraulic parameter sets. However, this model is much better equipped to simulate the behavior of different plant species. Based on these findings, we only recommend implementation of this approach into LSMs if great care with choice of soil hydraulic parameters is taken

The effect of thermal conductivity of RIM moulds in kinetics cure

Reaction Injection Moulding (RIM) is a moulding technology used for the production of large size and complex plastic parts. The RIM process is characterized essentially by the injection of a highly reactive chemical system (usually polyurethane) and fast cure, in a mould properly closed and thermally controlled. Several studies show that rapid manufacturing moulds obtained in epoxy resins for Thermoplastic Injection Moulding (TIM) affect the moulding process and the final properties of parts. The cycle time and mechanical properties of final parts are reduced, due to a low thermal conductivity of epoxy materials. In contrast, the low conductivity of materials usually applied for the rapid manufacturing of RIM moulds, increase the mechanical properties of final injected parts and reduce the cycle time. This study shows the effect of the rapid manufacturing moulds material during the RIM process. Several materials have been tested for rapid manufacturing of RIM moulds and the analysis of both, temperature profile of moulded parts during injection and the cure data experimentally obtained in a mixing and reaction cell, allow to determine and model the real effect of the mould material on the RIM process.

Benchmarking 2D hydraulic models for urban flooding

This paper describes benchmark testing of six two-dimensional (2D) hydraulic models (DIVAST, DIVASTTVD, TUFLOW, JFLOW, TRENT and LISFLOOD-FP) in terms of their ability to simulate surface flows in a densely urbanised area. The models are applied to a 1·0 km × 0·4 km urban catchment within the city of Glasgow, Scotland, UK, and are used to simulate a flood event that occurred at this site on 30 July 2002. An identical numerical grid describing the underlying topography is constructed for each model, using a combination of airborne laser altimetry (LiDAR) fused with digital map data, and used to run a benchmark simulation. Two numerical experiments were then conducted to test the response of each model to topographic error and uncertainty over friction parameterisation. While all the models tested produce plausible results, subtle differences between particular groups of codes give considerable insight into both the practice and science of urban hydraulic modelling. In particular, the results show that the terrain data available from modern LiDAR systems are sufficiently accurate and resolved for simulating urban flows, but such data need to be fused with digital map data of building topology and land use to gain maximum benefit from the information contained therein. When such terrain data are available, uncertainty in friction parameters becomes a more dominant factor than topographic error for typical problems. The simulations also show that flows in urban environments are characterised by numerous transitions to supercritical flow and numerical shocks. However, the effects of these are localised and they do not appear to affect overall wave propagation. In contrast, inertia terms are shown to be important in this particular case, but the specific characteristics of the test site may mean that this does not hold more generally.

An overview of Earth's global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

A double-Gerdien instrument for simultaneous bipolar air conductivity measurements on balloon platforms

A bipolar air conductivity instrument is described for use with a standard disposable meteorological radiosonde package. It is intended to provide electrical measurements at cloud boundaries, where the ratio of the bipolar air conductivities is affected by the presence of charged particles. The sensors are two identical Gerdien-type electrodes, which, through a voltage decay method, measure positive and negative air conductivities simultaneously. Voltage decay provides a thermally stable approach and a novel low current leakage electrometer switch is described which initiates the decay sequence. The radiosonde supplies power and telemetry, as well as measuring simultaneous meteorological data. A test flight using a tethered balloon determined positive (σ+) and negative (σ−) conductivities of σ+ = 2.77±0.2 fS m−1 and σ− = 2.82±0.2 fS m−1, respectively, at 400 m aloft, with σ+/σ− = 0.98±0.04.

Exploring the spatial relations between soil physical properties and apparent electrical conductivity

Relations between the apparent electrical conductivity of the soil (ECa) and top- and sub-soil physical properties were examined for two arable fields in southern England (Crowmarsh Battle Farms and the Yattendon Estate). The spatial variation of ECa and the soil properties was explored geostatistically. The variogram ranges showed that ECa varied on a similar spatial scale to many of the soil physical properties in both fields. Several features in the map of kriged predictions of ECa were also evident in maps of the soil properties. In addition, the correlation coefficients showed a strong relation between ECa and several soil properties. A moving correlation analysis enabled differences in the relations between ECa and the soil properties to be examined within the fields. The results indicated that relations were inconsistent; they were stronger in some areas than others. A regression of ECa on the principal component scores of the leading components for both fields showed that the first two components accounted for a large proportion of the variance in ECa, whereas the others accounted for little or none. For Crowmarsh topsoil sand and clay, loss on ignition and volumetric water measured in the autumn had large correlations on the first component, and for Yattendon they were large for topsoil sand and clay, and autumn and spring volumetric water. The cross-variograms suggested strong coregionalization between ECa and several soil physical properties; in particular subsoil sand and silt at Crowmarsh, and subsoil sand and clay at Yattendon. The structural correlations from the linear model of coregionalization confirmed the strength of the relations between ECa and the subsoil properties. Nevertheless, no one property was consistently important for both fields. Although a map of ECa can indicate the general patterns of spatial variation in the soil, it is not a substitute for information on soil properties obtained by sampling and analysing the soil. Nevertheless, it could be used to guide further sampling. (c) 2005 Elsevier B.V. All rights reserved.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

Phosphorus dynamics in the ditch system of a restored peat wetland

The restoration of wetlands as bird habitats often involves the maintenance of a fluctuating water regime by careful, localised ditch water management using pumps and sluices. However, there is evidence in the literature to Suggest that alternate flood/drainage cycles can accelerate nutrient cycling and transport within the soil and, therefore, pose a threat to water quality through the process of eutrophication. This study focused on the dynamics and losses of soil P in a recently re-wetted, eutrophic fen peat developed on alluvium ill South west England. During the 2-year Study (2001 and 2002), soil water tensiometry revealed that the field water table (fluctuating annually between +20 and 60 cm relative to ground level) was extensively influenced across most of the 8.4 ha field site by the management of the adjacent ditch water levels. This conservation-led, prescribed water balance was facilitated by the high hydraulic conductivity (1.1 x 10(-s) ms(-1)) of the lower (70-140 cm), degraded layer of peat. However, only during a 7-day period of water table drawdown by intermittent pump drainage, approximately 45 g ha(-1) of dissolved reactive P (DRP) entered the pumped ditch from the field via this degraded layer. Summer rainfall events >35 mm d(-1) also coincided with significant peaks ill ditch water P concentration (up to 200 mu g L-1 DRP). Even larger peaks (Up to 700 mu g L-1 DRP) Occurred With the annual onset of autumn reflooding. These episodic P loss events pose a serious potential threat to biological water quality. (C) 2009 Elsevier B.V. All rights reserved.

Wetland restoration within agricultural watersheds: Balancing water quality protection with habitat conservation

Peat wetlands that have been restored from agricultural Land have the potential to act as Long term sources of phosphorus (P) and, therefore have to potenital to accelerate freshwater eutrophication. During a two-year study the water table in a eutrophic fen peat that was managed by pump drainage fluctuated annually between +20 cm and -60 cm relative to ground Level. This precise management was facilitated by the high hydraulic conductivity (K) of the humified peat (1.1 x 10(-5) m s(-1)) below around 60 cm depth. However, during one week of intermittent pumping, as much as 50 g ha(-1) dissolved P entered the pumped ditch. Summer. rainfall events and autumn reflooding also triggered P losses. The P Losses were attributed to the low P sorption capacity (217 mg kg(-1)) of the saturated peat below 60 cm, combined with its high K and the reductive dissolution of Fe bound P.

Hydraulic modeling of runoff over a rough surface under partial inundation

This study presents a numerical method to derive the Darcy- Weisbach friction coefficient for overland flow under partial inundation of surface roughness. To better account for the variable influence of roughness with varying levels of emergence, we model the flow over a network which evolves as the free surface rises. This network is constructed using a height numerical map, based on surface roughness data, and a discrete geometry skeletonization algorithm. By applying a hydraulic model to the flows through this network, local heads, velocities, and Froude and Reynolds numbers over the surface can be estimated. These quantities enable us to analyze the flow and ultimately to derive a bulk friction factor for flow over the entire surface which takes into account local variations in flow quantities. Results demonstrate that although the flow is laminar, head losses are chiefly inertial because of local flow disturbances. The results also emphasize that for conditions of partial inundation, flow resistance varies nonmonotonically but does generally increase with progressive roughness inundation.