Chad William Higgins

Evapotranspiration: A process driving mass transport and energy exchange in the soil‐plant‐atmosphere‐climate system

[1] The role of evapotranspiration (ET) in the global, continental, regional, and local water cycles is reviewed. Elevated atmospheric CO2, air temperature, vapor pressure deficit (D), turbulent transport, radiative transfer, and reduced soil moisture all impact biotic and abiotic processes controlling ET that must be extrapolated to large scales. Suggesting a blueprint to achieve this link is the main compass of this review. Leaf-scale transpiration (fe) as governed by the plant biochemical demand for CO2 is first considered. When this biochemical demand is combined with mass transfer formulations, the problem remains mathematically intractable, requiring additional assumptions. A mathematical “closure” that assumes stomatal aperture is autonomously regulated so as to maximize the leaf carbon gain while minimizing water loss is proposed, which leads to analytical expressions for leaf-scale transpiration. This formulation predicts well the effects of elevated atmospheric CO2 and increases in D on fe. The case of soil moisture stress is then considered using extensive gas exchange measurements collected in drought studies. Upscaling the fe to the canopy is then discussed at multiple time scales. The impact of limited soil water availability within the rooting zone on the upscaled ET as well as some plant strategies to cope with prolonged soil moisture stress are briefly presented. Moving further up in direction and scale, the soil-plant system is then embedded within the atmospheric boundary layer, where the influence of soil moisture on rainfall is outlined. The review concludes by discussing outstanding challenges and how to tackle them by means of novel theoretical, numerical, and experimental approaches.

Albedo effect on radiative errors in air temperature measurements

Most standard air temperature measurements are subject to significant errors mainly due to sensor heating by solar radiation, even when the measurement principle is accurate and precise. We present various air temperature measurements together with other measurements of meteorological parameters using different sensor systems at a snow-covered and a vegetated site. Measurements from naturally ventilated air temperature sensors in multiplate shields are compared to temperatures measured using sonic anemometers which are unaffected by solar radiation. Over snow, 30 min mean temperature differences can be as large as 10°C. Unshielded thermocouples were also tested and are generally less affected by shortwave radiation. Temperature errors decrease with decreasing solar radiation and increasing wind speed but do not completely disappear at a given solar radiation even in the presence of effective ventilation. We show that temperature errors grow faster for reflected than for incident solar radiation, demonstrating the influence of the surface properties on radiative errors, and we detect the albedo as a variable with major influence on the magnitude of the error as well as a key quantity in possible error correction schemes. An extension is proposed for an existing similarity regression model to correct for radiative errors; thus, surface-reflected shortwave radiation is identified as a principal source of error and the key variable for obtaining a unique nondimensional scaling of radiative errors.

Alignment trends of velocity gradients and subgrid-scale fluxes in the turbulent atmospheric boundary layer

Field experimental data in the atmospheric surface layer are analyzed using toolsfrom statistical geometry. The data consist of velocity measurements from sonicanemometer arrays. In the context of large eddy simulations (LES), these arrayspermit the spatial filtering needed to separate large from small scales. Time seriesof various quantities relevant to LES are evaluated from the data. Results show thatthe preferred filtered fluid deformation is axisymmetric extension and the preferredsubgrid stress state is axisymmetric contraction. The filtered fluctuating vorticityshows preferred alignments with the mean vorticity, with the streamwise direction,and with the intermediate strain-rate eigenvector. The alignment between eigenvectorsof the subgrid-scale stress and filtered strain rate is used to test eddy viscosity andmixed model formulations. In qualitative agreement with prior laboratory measurements at much lower Reynolds numbers, a bimodal distribution is observed, which can be reduced to good alignment with eddy viscosity closure using the mixed model.

Field study of the dynamics and modelling of subgrid-scale turbulence in a stable atmospheric surface layer over a glacier

A field experiment – the Snow Horizontal Array Turbulence Study (SnoHATS) – has been performed over an extensive glacier in Switzerland in order to study small-scale turbulence in the stable atmospheric surface layer, and to investigate the role, dynamics and modelling of the subgrid scales (SGSs) in the context of large-eddy simulations. The a priori data analysis aims at comparing the role and behaviour of the SGSs under stable conditions with previous studies under neutral or unstable conditions. It is found that the SGSs in a stable surface layer remain an important sink of temperature variance and turbulent kinetic energy from the resolved scales and carry a significant portion of the fluxes when the filter scale is larger than the distance to the wall. The fraction of SGS fluxes (out of the total fluxes) is found to be independent of stability. In addition, the stress–strain alignment is similar to the alignment under neutral and unstable conditions. The model coefficients vary considerably with stability but in a manner consistent with previous findings, which also showed that scale-dependent dynamic models can capture this variation. Furthermore, the variation of the coefficients for both momentum and heat SGS fluxes can be shown to be better explained by stability parameters based on vertical gradients, rather than vertical fluxes. These findings suggest that small-scale turbulence dynamics and SGS modelling under stable conditions share many important properties with neutral and convective conditions, and that a unified approach is thus possible. This paper concludes with a discussion of some other challenges for stable boundary-layer simulations that are not encountered in the neutral or unstable cases.

The Effect of Scale on the Applicability of Taylor's Frozen Turbulence Hypothesis in the Atmospheric Boundary Layer

Taylor’s frozen turbulence hypothesis is the central assumption invoked in most experiments designed to investigate turbulence physics with time resolving sensors. It is also frequently used in theoretical discussions when linking Lagrangian to Eulerian flow formalisms. In this work we seek to quantify the effectiveness of Taylor’s hypothesis on the field scale using water vapour as a passive tracer. A horizontally orientated Raman lidar is used to capture the humidity field in space and time above an agricultural region in Switzerland. High resolution wind speed and direction measurements are conducted simultaneously allowing for a direct test of Taylor’s hypothesis at the field scale. Through a wavelet decomposition of the lidar humidity measurements we show that the scale of turbulent motions has a strong influence on the applicability of Taylor’s hypothesis. This dependency on scale is explained through the use of dimensional analysis. We identify a ‘persistency scale’ that can be used to quantify the effectiveness of Taylor’s hypothesis, and present the accuracy of the hypothesis as a function of this non-dimensional length scale. These results are further investigated and verified through the use of large-eddy simulations.

Flow over Hills: A Large-Eddy Simulation of the Bolund Case

Simulation of local atmospheric flows around complex topography is important for several applications in wind energy (short-term wind forecasting and turbine siting and control), local weather prediction in mountainous regions and avalanche risk assessment. However, atmospheric simulation around steep mountain topography remains challenging, and a number of different approaches are used to represent such topography in numerical models. The immersed boundary method (IBM) is particularly well-suited for efficient and numerically stable simulation of flow around steep terrain. It uses a homogenous grid and permits a fast meshing of the topography. Here, we use the IBM in conjunction with a large-eddy simulation (LES) and test it against two unique datasets. In the first comparison, the LES is used to reproduce experimental results from a wind-tunnel study of a smooth three-dimensional hill. In the second comparison, we simulate the wind field around the Bolund Hill, Denmark, and make direct comparisons with field measurements. Both cases show good agreement between the simulation results and the experimental data, with the largest disagreement observed near the surface. The source of error is investigated by performing additional simulations with a variety of spatial resolutions and surface roughness properties.

The Effect of Filter Dimension on the Subgrid-Scale Stress, Heat Flux, and Tensor Alignments in the Atmospheric Surface Layer

Abstract In field experiments designed to study subgrid-scale parameterizations for large eddy simulation, the flow field is often measured and then filtered in two-dimensional planes. This two-dimensional filtering serves as a surrogate for three-dimensional filtering. The question of whether this will yield accurate results in subgrid-scale (SGS) models is addressed by analyzing data from a field experiment in which 16 sonic anemometers were deployed in a four by four grid. The experiment was held in July 2002 at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility in the Utah West Desert. The full SGS stress tensor and its parameterizations using both two- and three-dimensional filterings are obtained. Comparisons are given between two- and three-dimensional filterings of the field measurements based on probability density functions (PDFs) and energy spectra of the SGS stress elements. The PDFs reveal that quantities calculated with two-dimensional filtering exhibit greater int...

Towards oscillation-free implementation of the immersed boundary method with spectral-like methods

It is known that, when the immersed boundary method (IBM) is implemented within spectral-like methods, the Gibbs oscillation seriously deteriorates the calculation of derivatives near the body surface. In this paper, a radial basis function (RBF) based smoothing technique is proposed with the intention of eliminating or efficiently reducing the Gibbs oscillation without affecting the flow field outside the body. Based on this technique, a combined IBM/spectral scheme is developed to solve the incompressible Navier-Stokes equations. Numerical simulations of flow through a periodic lattice of cylinders of various cross sections are performed. The results demonstrate that the proposed methodology is able to give accurate and nearly oscillation-free numerical solutions of incompressible viscous flows.

The heat flux and the temperature gradient in the lower atmosphere

] MostparameterizationsusedinLargeEddySimulationsof the atmospheric boundary layer are based on theassumption that subgrid-scale fluxes are aligned againstspatial gradients of transported quantities (down-gradientclosures). Based on field experiments, we determine thedistribution and most probable relative orientations of thesubgrid-scale (SGS) heat flux relative to parameterizationsbased on the temperature gradient. We show that, underneutral and unstable atmospheric stability, the SGS heatflux most likely lies within the geometric span of the socalledmixedtensoreddydiffusivitymodel.

Measured and Estimated Water Vapor Advection in the Atmospheric Surface Layer

The flux of water vapor due to advection is measured using high-resolution Raman lidar that was orientated horizontally across a land‐lake transition.At the same time, a fullsurface energy balance is performed toassess the impact of scalar advection on energy budget closure. The flux of water vapor due to advection is then estimated with analytical solutions to the humidity transport equation that show excellent agreement with the field measurements. Although the magnitude of the advection was not sufficient to account for the total energy deficit for this field site, the analytical approach is used to explore situations where advection would be the dominanttransportmechanism.Theauthorsfindthatadvectionisatmaximumwhenthemeasurementheightis 0.036 times the distance to a land surfacetransition. The framework proposed in this paper can beused to predict the potential impact of advection on surface flux measurements prior to field deployment and can be used as a data analysis algorithm to calculate the flux of water vapor due to advection from field measurements.