Estimation of groundwater evaporation and salt flux from Owens Lake, California, USA

Groundwater evaporation and subsequent precipitation of soluble salts at Owens Lake in eastern California have created one of the single largest sources of airborne dust in the USA, yet the evaporation and salt flux have not been fully quantified. In this study, we compare eddy correlation, microlysimeters and solute profiling methods to determine their validity and sensitivity in playa environments. These techniques are often used to estimate evaporative losses, yet have not been critically compared at one field site to judge their relative effectiveness and accuracy. Results suggest that eddy correlation methods are the most widely applicable for the variety of conditions found on large playa lakes. Chloride profiling is shown to be highly sensitive to thermal and density-driven fluxes in the near surface and, as a result, appears to underestimate yearly groundwater evaporation. Yearly mean groundwater evaporation from the playa surface estimated from the three study areas was found to range from 88 to 104 mm year⁻¹, whereas mean evaporation from the brine-covered areas was 872 mm year⁻¹. Uncertainties on these mean rates were estimated to be ±25%, based on comparisons between eddy correlation and lysimeter estimates. On a yearly basis, evaporation accounts for approximately 47 × 10⁶ m³ of water loss from the playa surface and open-water areas of the lake. Over the playa area, as much as 7.5 × 10⁸ kg (7.5 × 10⁵ t) of salt are annually concentrated by evaporation at or near the playa surface, much of which appears to be lost during dust storms in area.     

Local isotropy and anisotropy in the sheared and heated atmospheric surface layer

Longitudinal velocity and temperature measurements above a uniform dry lakebed were used to investigate sources of eddy-motion anisotropy within the inertial subrange. Rather than simply test the adequacy of locally isotropic relations, we investigated directly the sources of anisotropy. These sources, in a daytime desert-like climate, include: (1) direct interaction between the large-scale and small-scale eddy motion, and (2) thermal effects on the small-scale eddy motion. In order to explore these two anisotropy sources, we developed statistical measures that are sensitive to such interactions. It was found that the large-scale/small-scale interaction was significant in the inertial subrange up to 3 decades below the production scale, thus reducing the validity of the local isotropy assumption. The anisotropy generated by thermal effects was also significant and comparable in magnitude to the former anisotropy source. However, this thermal anisotropy was opposite in sign and tended to counteract the anisotropy generated by the large-scale/smallscale interaction. The thermal anisotropy was attributed to organized ramp-like patterns in the temperature measurements. The impact of this anisotropy cancellation on the dynamics of inertial subrange eddy motion was also considered. For that purpose, the Kolmogorov-Obukhov structure function equation, as derived from the Navier-Stokes equations for locally isotropic turbulence, was employed. The Kolmogorov-Obukhov structure function equation in conjunction with Obukhov's constant skewness closure hypothesis reproduced the measured second- and third-order structure functions. Obukhov's constant skewness closure scheme, which is also based on the local isotropy assumption, was verified and was found to be in good agreement with the measurements. The accepted 0.4 constant skewness value derived from grid turbulence experiments overestimated our measurements. A suggested 0.26 constant skewness value, which we derived from Kolmogorov's constant, was found to be adequate.     

A Theoretical and Experimental Investigation of Energy-Containing Scales in the Dynamic Sublayer of Boundary-Layer Flows

The existence of universal power laws at low wavenumbers (K) in the energy spectrum (E_u) of the turbulent longitudinal velocity (u) is examined theoretically and experimentally for the near-neutral atmospheric surface layer. Newly derived power-law solutions to Tchen's approximate integral spectral budget equation are tested for strong- and weak-interaction cases between the mean flow and turbulent vorticity fields. To verify whether these solutions reproduce the measured E_u at low wavenumbers, velocity measurements were collected in the dynamic sublayer of the atmosphere at three sites and in the inner region of a laboratory open channel. The atmospheric surface layer measurements were carried out using triaxial sonic anemometers over tall corn, short grass, and smooth desert-like sandy soil. The open channel measurements were performed using a two-dimensional boundary-layer probe above a smooth stainless steel bed. Comparisons between the proposed analytical solution for E_u, the dimensional analysis by Kader and Yaglom, and the measured Haar wavelet E_u spectra are presented. It is shown that when strong interaction between the mean flow and turbulent vorticity field occurs, wavelet spectra measurements, predictions by the analytical solution, and predictions by the dimensional analysis of Kader-Yaglom (KY) are all in good agreement and confirm the existence of a -1 power law in E_u(= C_uuu² _* K^-1, where C_uu is a constant and u_* is the friction velocity). The normalized upper wavenumber limit of the -1 power law (Kz = 1, where z is the height above the zero-plane displacement) is estimated using two separate approaches and compared to the open channel and atmospheric surface-layer measurements. It is demonstrated that the measured upper wavenumber limit is consistent with Tchen's budget but not with the KY assumptions. The constraints as to whether the mean flow and turbulent vorticity strongly interact are considered using a proposed analysis by Panchev. It is demonstrated that the arguments by Panchev cannot be consistent with surface-layer turbulence. Using dimensional analysis and Heisenberg's turbulent viscosity model, new constraints are proposed. The new constraints agree with the open channel and atmospheric surface-layer measurements, Townsend's inactive eddy motion hypothesis, and the Perry et al. analysis.     

Flume Experiments on Turbulent Flows Across Gaps of Permeable and Impermeable Boundaries

Laser Doppler anemometery and laser-induced fluorescence techniques were used to explore the spatial structure of the flow within and above finite cavities created within porous and solid media. The cavities within these two configurations were identical in size and were intended to mimic flow disturbances created by finite gaps and forest clearing. Because flows over permeable boundaries differ from their solid counterparts, the study here addresses how these differences in boundary conditions produce differences in, (i) bulk flow properties including the mean vorticity within and adjacent to the gaps, (ii) second-order statistics such as the standard deviations and turbulent stresses, (iii) the relative importance of advective to turbulent stress terms across various regions within and above the gaps, and (iv) the local imbalance between ejections and sweeps and momentum transport efficiencies of updrafts and downdrafts. Both configurations exhibited a primary recirculation zone of comparable dimensions inside the gap. The mean vorticity spawned at the upstream corner of the gap was more intense for the solid configuration when compared to its porous counterpart. The free-shear layer spawned from the upstream corner-edge deeper into the gap for the porous configuration. The momentum flux at the interface within and above the gap was enhanced by a factor of 1.5–2.0 over its upstream value, and this enhancement zone was much broader in size for the porous configuration. For the turbulent transport terms in the longitudinal and vertical mean momentum balances, these transport terms were significant inside the gap for both boundary configurations when compared to their upstream counterpart. The effectiveness of using incomplete cumulant expansion methods to describe the momentum transport efficiencies, and the relative contributions of ejections and sweeps to turbulent stresses, especially in this zone, were also demonstrated. The flatness factor for both velocity components, often used as a measure of intermittency, was highest in the vicinity of the upstream corner in both configurations. However, immediately following the downstream corner, the flatness factor remained large for the porous configuration, in contrast to its solid configuration counterpart.     

Scaling Properties of Biologically Active Scalar Concentration Fluctuations in the Atmospheric Surface Layer over a Managed Peatland

The higher-order scalar concentration fluctuation properties are examined in the context of Monin–Obukhov similarity theory for a variety of greenhouse gases that have distinct and separate source/sink locations along an otherwise ideal micrometeorological field site. Air temperature and concentrations of water vapour, carbon dioxide and methane were measured at high frequency (10 Hz) above a flat and extensive peat-land soil in the San Joaquin–Sacramento Delta (California, USA) area, subjected to year-round grazing by beef cattle. Because of the heterogeneous distribution of the sources and sinks of CO₂ and especially CH₄ emitted by cattle, the scaling behaviour of the higher-order statistical properties diverged from predictions based on a balance between their production and dissipation rate terms, which can obtained for temperature and H₂O during stationary conditions. We identify and label these departures as ‘exogenous’ because they depend on heterogeneities and non-stationarities induced by boundary conditions on the flow. Spectral analysis revealed that the exogenous effects show their signatures in regions with frequencies lower than those associated with scalar vertical transport by turbulence, though the two regions may partially overlap in some cases. Cospectra of vertical fluxes appear less influenced by these exogenous effects because of the modulating role of the vertical velocity at low frequencies. Finally, under certain conditions, the presence of such exogenous factors in higher-order scalar fluctuation statistics may be ‘fingerprinted’ by a large storage term in the mean scalar budget.     

Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Turbulent Intensities and Velocity Spectra for Bare and Forested Gentle Hills: Flume Experiments

To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.     

A Note on the Flux-Variance Similarity Relationships for Heat and Water Vapour in the Unstable Atmospheric Surface Layer

Atmospheric surface layer (ASL) experiments over the past 10 years demonstrate that the flux-variance similarity functions for water vapour are consistently larger in magnitude than their temperature counterpart. In addition, latent heat flux calculations using the flux-variance method do not compare as favorably to eddy- correlation measurements when compared to their sensible heat counterpart. These two findings, in concert with measured heat to water vapour transport efficiencies in excess of unity, are commonly used as evidence of dissimilarity between heat and water vapour transport in the unstable atmospheric surface layer. In this note, it is demonstrated that even if near equality in flux-profile similarity functions for heat and water vapour is satisfied, the flux-variance similarity functions for water vapour are larger in magnitude than temperature for a planar, homogeneous, unstably-stratified turbulent boundary-layer flow.     

Landslide Research at the British Geological Survey:Capture,Storage and Interpretation on a National and Site-Specific Scale

Abstract: Landslide research at the British Geological Survey (BGS) is carried out through a number of activities, including surveying, database development and real-time monitoring of landslides. Landslide mapping across the UK has been carried out since BGS started geological mapping in 1835. Today, BGS geologists use a combination of remote sensing and ground-based investigations to survey landslides. The development of waterproof tablet computers (BGS·SIGMAmobile), with inbuilt GPS and GIS for field data capture provides an accurate and rapid mapping methodology for field surveys. Regional and national mapping of landslides is carried out in conjunction with site-specific monitoring, using terrestrial LiDAR and differential GPS technologies, which BGS has successfully developed for this application. In addition to surface monitoring, BGS is currently developing geophysical ground-imaging systems for landslide monitoring, which provide real-time information on subsurface changes prior to failure events. BGS’s mapping and monitoring activities directly feed into the BGS National Landslide Database, the most extensive source of information on landslides in Great Britain. It currently holds over 14?000 records of landslide events. By combining BGS’s corporate datasets with expert knowledge, BGS has developed a landslide hazard assessment tool, GeoSure, which provides information on the relative landslide hazard susceptibility at national scale.