Validity limits for the van Genuchten–Mualem model and implications for parameter estimation and numerical simulation

The calculation of the relative hydraulic conductivity function based on water retention data is an attractive and widely used approach, since direct measurements of unsaturated conductivities are difficult. We show theoretically under which conditions an air-entry value for water retention data is definitely required when using the statistical approach of Mualem. Moreover we rigorously specify the conditions for which the classical van Genuchten–Mualem model leads to wrong predictions of relative hydraulic conductivity and, hence, an alternative formulation including an air-entry value should be used. Significant consequences are demonstrated for the inverse parameter estimation based on multistep outflow experiments. Furthermore it is shown that the use of a physically correct formulation of the water retention curve including an air-entry value and the derived hydraulic conductivity function influences not only the stability of numerical simulations but also their final results. This is especially grave as simulations with van Genuchten–Mualem parameters are frequently used to compare experiments and simulations and to draw conclusions on the correctness of Richards’ equation.     

A review of: “Turbulence and random processes in fluid mechanics”

Abstract

by M. T. Landahl and E. Mollo-Christensen, Cambridge University Press, vi+154 pp., ?34.50 (ISBN (0-521-26306-9) 1986.     

Techniques for assessing water infrastructure for nonstationary extreme events: a review

Statistical and physically-based methods have been used for designing and assessing water infrastructure such as spillways and stormwater drainage systems. Traditional approaches assume that hydrological processes evolve in an environment where the hydrological cycle is stationary over time. However, in recent years, it has become increasingly evident that in many areas of the world the foregoing assumption may no longer apply, due to the effect of anthropogenic and climatic induced stressors that cause nonstationary conditions. This has attracted the attention of national and international agencies, research institutions, academia, and practicing water specialists, which has led to developing new techniques that may be useful in those cases where there is good evidence and attribution of nonstationarity. We review the various techniques proposed in the field and point out some of the challenges ahead in future developments and applications. Our review emphasizes hydrological design to protect against extreme events such as floods and low flows.     

Spatial Variability of Turbulent Fluxes in the Roughness Sublayer of an Even-Aged Pine Forest

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (_u, _w) and temperature (_T) are more planar homogeneous than their vertical flux of momentum (u_* ²) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO₂ and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the canonical structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the moving-equilibrium hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u_*, especially when u_* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (I_w) and the arrival frequency of organized structures (/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.     

The climate of Namaqualand in the nineteenth century

Southern African climatic change research is hampered by a lack of long-term historical data sets. This paper aims to extend the historical climate record for southern Africa to the semi-arid area of Namaqualand in the Northern Cape province of South Africa. This is achieved through extensive archival research, making use of historical documentary sources such as missionary journals and letters, traveller’s writings and government reports and letters. References to precipitation and other climatic conditions have been extracted and categorised, providing a proxy precipitation data set for Namaqualand for the nineteenth century. Notwithstanding problems of data accuracy and interpretation the reconstruction enables the detection of severe and extreme periods. Measured meteorological data, available from the late 1870s, was compared to the data set derived from documentary sources in order to ascertain the accuracy of the data set and monthly rainfall data has been used to identify seasonal anomalies. Confidence ratings on derived dry and wet periods, where appropriate, have been assigned to each year. The study extends the geographical area of existing research and extracts the major periods of drought and climatic stress, from the growing body of historical climate research. The most widespread drought periods affecting the southern and eastern Cape, Namaqualand and the Kalahari were 1820–1821; 1825–1827; 1834; 1861–1862; 1874–1875; 1880–1883 and 1894–1896. Finally, a possible correspondence is suggested between some of the widespread droughts and the El Nino Southern Oscillation (ENSO).     

Extraction of crystal-poor rhyolite from a hornblende-bearing intermediate mush: a case study of the caldera-forming Matahina eruption,Okataina volcanic complex

The Matahina Ignimbrite (~160 km³ rhyolite magma, 330 ka) was deposited during a caldera-forming eruption from the Okataina Volcanic Centre, Taupo Volcanic Zone (TVZ), New Zealand. Juvenile clasts are divided into three groups: Group (1) the dominant crystal-poor rhyolite type, Group (2) a minor coarse-grained, mingled/mixed intermediate type, and Group (3) a rare fine-grained basalt. The ignimbrite consists of the Group 1 type and is divided into three members: a lower and middle member, which is high-silica, crystal-poor (<10 vol.%) rhyolite, and the upper member, which is low-silica and slightly more crystal-rich (up to 21 vol.%). Cognate, crystal-rich (up to 50 vol.%) basalt to intermediate pumice occurs on top of lag breccias and within lithic-rich pyroclastic density current deposits along the caldera margin (Groups 2 and 3). Several lines of evidence indicate that the intermediate clasts represent the cumulate complement to the melt-rich rhyolite: (1) continuity in the compositions of plagioclase, orthopyroxene, hornblende, and oxides and normal zoning of individual phenocrysts; (2) the silicic glass from the intermediate magma (interstitial melt) overlaps compositionally with the bulk rock rhyolite and glass; (3) high Zr and a slight positive Eu anomaly in the intermediate magma relative to quenched enclaves from other intermediate TVZ eruptions indicates zircon and plagioclase accumulation, respectively; (4) an increase in the Cl contents in glass from the least evolved to most evolved is consistent with the concentration of volatiles during magma evolution. Most of the compositional variations in the low- to high-silica rhyolites can be accounted for by continued Rayleigh fractionation (up to 15%), following melt extraction from the underlying mush, under varying fO₂–fH₂O conditions to form a slightly compositionally zoned rhyolitic cap. This link to the varying fO₂–fH₂O conditions is evidenced by the strong correlation between key geochemical parameters (e.g. Dy, Y), that qualitatively reflect fH₂O conditions (presence or absence of hornblende/biotite), and fO₂ estimated from Fe–Ti oxide equilibrium. Magma mingling/mixing between the basalt–andesite and the main slightly compositionally zoned rhyolitic magma occurred during caldera-collapse, modifying the least-evolved rhyolite at the lower portion of the reservoir and effectively destroying any pre-eruptive gradients.     

The great sulfur depletion of Earth’s mantle is not a signature of mantle–core equilibration

The extreme depletion of the Earth’s mantle in sulfur is commonly seen as a signature of metal segregation from Earth’s mantle to Earth’s core. However, in addition to S, the mantle contains other elements as volatile as S that are hardly depleted relative to the lithophile volatility trend although they are potentially as siderophile as sulfur. We report experiments in metal-sulfide–silicate systems to show that the CI normalized abundances of S, Pb, and Sn in Earth’s mantle cannot be reproduced by element partitioning in Fe ± S–silicate systems, neither at low nor at high pressure. Much of the volatile inventory of the Earth’s mantle must have been added late in the accretion history, when metal melt segregation to the core had become largely inactive. The great depletion in S is attributed to the selective segregation of a late sulfide matte from an oxidized and largely crystalline mantle. Apparently, the volatile abundances of Earth’s mantle are not in redox equilibrium with Earth’s core.     

Long‐Term Sediment Export Estimates from Northern Jordan using Roman Cisterns as Sediment Traps

Roman cisterns served as rainwater storage devices for centuries and are densely distributed in parts of northern Jordan. A major earthquake hit the region ca. A.D. 750 and in a short time many settlements were abandoned. As a consequence, most cisterns were not maintained, and they filled with sediments that today provide a postabandonment depositional record. In two field surveys, we mapped the locations of more than 100 cisterns in the Wadi Al‐Arab basin and selected two for detailed stratigraphic analysis that included ¹⁴C and optically stimulated luminescence dating. Catchment basin area for each cistern was determined by differential GPS. Both cisterns filled with sediments after the great earthquake and consequent abandonment of the region. Calculated sediment volumes are translated to long‐term average sediment export rates of 2.6–6.6 t ha⁻¹a⁻¹, which are comparable to erosion and sediment yield rates from other studies within the Mediterranean region. Our pilot study suggests that this approach can be applied elsewhere to calculate long‐term sediment export rates on hill slopes containing relict cisterns.     

Mechanisms and Model Diversity of Trade-Wind Shallow Cumulus Cloud Feedbacks: A Review

Shallow cumulus clouds in the trade-wind regions are at the heart of the long standing uncertainty in climate sensitivity estimates. In current climate models, cloud feedbacks are strongly influenced by cloud-base cloud amount in the trades. Therefore, understanding the key factors controlling cloudiness near cloud-base in shallow convective regimes has emerged as an important topic of investigation. We review physical understanding of these key controlling factors and discuss the value of the different approaches that have been developed so far, based on global and high-resolution model experimentations and process-oriented analyses across a range of models and for observations. The trade-wind cloud feedbacks appear to depend on two important aspects: (1) how cloudiness near cloud-base is controlled by the local interplay between turbulent, convective and radiative processes; (2) how these processes interact with their surrounding environment and are influenced by mesoscale organization. Our synthesis of studies that have explored these aspects suggests that the large diversity of model responses is related to fundamental differences in how the processes controlling trade cumulus operate in models, notably, whether they are parameterized or resolved. In models with parameterized convection, cloudiness near cloud-base is very sensitive to the vigor of convective mixing in response to changes in environmental conditions. This is in contrast with results from high-resolution models, which suggest that cloudiness near cloud-base is nearly invariant with warming and independent of large-scale environmental changes. Uncertainties are difficult to narrow using current observations, as the trade cumulus variability and its relation to large-scale environmental factors strongly depend on the time and/or spatial scales at which the mechanisms are evaluated. New opportunities for testing physical understanding of the factors controlling shallow cumulus cloud responses using observations and high-resolution modeling on large domains are discussed.     

Feldspar geothermometry of the Hell Canyon Pluton,Boulder Batholith,Montana

The bulk compositions of the groundmass alkali feldspar from the Hell Canyon Pluton is 0.146mole% albite. The composition of the outermost zone of the oscillatory zoned plagioclase is 0.686 mole% albite, whereas the most calcic cores have a composition of 0.43 mole% albite. The structural state of the alkali feldspar is near orthoclase. Both composition of coexisting feldspars and structural state of the alkali feldspar are nearly constant throughout the pluton.Exsolved albite in the alkali feldspar have a composition of 0.965 mole% albite and the orthoclase host has a composition of 0.032 mole%. Singe crystal X-ray studies indicate that the albite intergrowths are coherent with the host.Equilibrium temperatures derived from the coexisting feldspar average 554 ° C; about 150 ° C, too low for the minimum solidus temperatures for reasonable emplacement pressures (2 kb). If this minimum solidus temperature is assumed, then the alkali feldspar has lost about 0.15 mole% albite. This loss was most likely caused by hydrothermal solutions associated with the crystallizing magma and equilibrated at about 550 ° C. However, based on the coherent albite intergrowths and the orthoclase structure state it can be inferred that the system was relatively free of volatiles below 500 ° C. Final equilibirium between orthoclase host and albite intergrowths occurred at about 311 ° C.