Semi-arid forests are found to sustain a massive sensible heat flux in spite of having a low surface to air temperature difference by lowering the aerodynamic resistance to heat transfer (rH) – a ...property called the canopy convector effect (CCE). In this work large-eddy simulations are used to demonstrate that the CCE appears more generally in canopy turbulence. It is indeed a generic feature of canopy turbulence: rH of a canopy is found to reduce with increasing unstable stratification, which effectively increases the aerodynamic roughness for the same physical roughness of the canopy. This relation offers a sufficient condition to construct a general description of the CCE. In addition, we review existing parameterizations for rH from the evapotranspiration literature and test to what extent they are able to capture the CCE, thereby exploring the possibility of an improved parameterization.
Rising temperatures and changes in snow cover, as can be expected under a warmer global climate, may have large impacts on mountain grassland productivity limited by cold and long winters. Here, we ...combined two existing models, the multi-layer atmosphere-SOiL-VEGetation model (SOLVEG) and the BASic GRAssland model (BASGRA), which accounts for snow, freeze–thaw events, grass growth, and soil carbon balance. The model was applied to simulate the responses of managed grasslands to anomalously warm winter conditions. The grass growth module considered key ecological processes under a cold environment, such as leaf formation, elongation and death, tillering, carbon allocation, and cold acclimation, in terms of photosynthetic activity. Input parameters were derived for two pre-Alpine grassland sites in Germany, for which the model was run using 3 years of data that included a winter with an exceptionally small amount of snow. The model reproduced the temporal variability of observed daily mean heat fluxes, soil temperatures, and snow depth throughout the study period. High physiological activity levels during the extremely warm winter led to a simulated CO2 uptake of 100 gC m−2, which was mainly allocated into the belowground biomass and only to a minor extent used for additional plant growth during early spring. If these temporary dynamics are representative of long-term changes, this process, which is so far largely unaccounted for in scenario analysis using global terrestrial biosphere models, may lead to carbon accumulation in the soil and/or carbon loss from the soil as a response to global warming.
Overview of the PALM model system 6.0 Maronga, Björn; Banzhaf, Sabine; Burmeister, Cornelia ...
Geoscientific Model Development,
03/2020, Letnik:
13, Številka:
3
Journal Article
Recenzirano
Odprti dostop
In this paper, we describe the PALM model system 6.0. PALM (formerly an abbreviation for Parallelized Large-eddy Simulation Model and now an independent name) is a Fortran-based code and has been ...applied for studying a variety of atmospheric and oceanic boundary layers for about 20 years. The model is optimized for use on massively parallel computer architectures. This is a follow-up paper to the PALM 4.0 model description in Maronga et al. (2015). During the last years, PALM has been significantly improved and now offers a variety of new components. In particular, much effort was made to enhance the model with components needed for applications in urban environments, like fully interactive land surface and radiation schemes, chemistry, and an indoor model. This paper serves as an overview paper of the PALM 6.0 model system and we describe its current model core. The individual components for urban applications, case studies, validation runs, and issues with suitable input data are presented and discussed in a series of companion papers in this special issue.
It is known from arid and semi-arid ecosystems that atmospheric water vapor can directly be adsorbed by the soil matrix. Soil water vapor adsorption was typically neglected and only recently received ...attention because of improvements in measurement techniques. One technique rarely explored for the measurement of soil water vapor adsorption is eddy covariance (EC). Soil water vapor adsorption may be detectable as downwardly directed (i.e., negative) EC latent heat (λE) flux measurements under dry conditions, but a systematic assessment of the use of negative λE fluxes from EC flux stations to characterize adsorption is missing. We propose a classification method to characterize soil water vapor adsorption, excluding conditions of dew and fog when λE derived from EC is not trustworthy due to stable atmospheric conditions. We compare downwardly directed λE fluxes from EC with measurements from weighing lysimeters for 4 years in a Mediterranean savanna ecosystem and 3 years in a temperate agricultural site. Our aim is to assess if overnight water inputs from soil water vapor adsorption differ between ecosystems and how well they are detectable by EC.
This study tests the sensitivity of a Lagrangian Stochastic footprint model to the turbulence statistics describing the flow field, with a focus on the within canopy processes. Representative ...profiles of the input velocity statistics are taken from a long-term dataset of turbulence measurements within and above a tall spruce canopy. Based on a wavelet tool, which allows a detailed analysis of coherent structures along the vertical profile, we characterize several typical states of coupling and decoupling between surface, canopy and atmosphere. For each coupling regime, three flux footprints using different sources for turbulence statistics are compared: the first based on conditionally-averaged measurements, the second on a simple numerical solution and the third on measurements taken from literature. The effects of profile smoothing and connecting measured canopy data to parametrized atmospheric surface layer profiles are considered. Significant differences between footprints based on modelled and measured profiles were found for exchange regimes with the lower section of the profiles decoupled from the atmospheric surface layer. As such cases are likely to occur for tall canopies with moderate density, our results suggest that the accuracy of Lagrangian Stochastic footprint modelling could be improved by using better turbulence profiles for different exchange regimes.
The eddy-covariance method provides the most direct estimates for fluxes between ecosystems and the atmosphere. However, dispersive fluxes can occur in the presence of secondary circulations, which ...can inherently not be captured by such single-tower measurements. In this study, we present options to correct local flux measurements for such large-scale transport based on a non-local parametric model that has been developed from a set of idealized large-eddy simulations. This method is tested for three real-world sites (DK-Sor, DE-Fen, and DE-Gwg), representing typical conditions in the mid-latitudes with different measurement heights, different terrain complexities, and different landscape-scale heterogeneities. Two ways to determine the boundary-layer height, which is a necessary input variable for modelling the dispersive fluxes, are applied, which are either based on operational radio soundings and local in situ measurements for the flat sites or from backscatter-intensity profiles obtained from co-located ceilometers for the two sites in complex terrain. The adjusted total fluxes are evaluated by assessing the improvement in energy balance closure and by comparing the resulting latent heat fluxes with evapotranspiration rates from nearby lysimeters. The results show that not only the accuracy of the flux estimates is improved but also the precision, which is indicated by RMSE values that are reduced by approximately 50 %. Nevertheless, it needs to be clear that this method is intended to correct for a bias in eddy-covariance measurements due to the presence of large-scale dispersive fluxes. Other reasons potentially causing a systematic underestimated or overestimation, such as low-pass filtering effects and missing storage terms, still need to be considered and minimized as much as possible. Moreover, additional transport induced by surface heterogeneities is not considered.
Display omitted
•Aircraft-based CH4 measurements yield high resolution flux estimates.•Comparable CH4 flux measurements using EC, wavelet covariance and REA techniques.•Wavelet analysis separates CH4 ...emissions from managed and unmanaged CH4 sources.•Larger top-down CH4 flux than agricultural bottom-up CH4 emissions.•Wetland area and surface temperature influence the top-down flux estimates.
Agriculture is estimated to produce more than 40% of anthropogenic methane (CH4) emissions, contributing to global climate change. Bottom-up, IPCC based methodologies are typically used to estimate the agriculture sector’s contribution, but these estimates are rarely verified beyond the farm gate, due to the challenge of separating interspersed sources. We present flux measurements of CH4, using eddy covariance (EC), relaxed eddy accumulation (REA) and wavelet covariance obtained using an aircraft-based measurement platform and compare these top-down estimates with bottom-up footprint adjusted inventory estimates of CH4 emissions for an agricultural region in eastern Ontario, Canada. Top-down CH4 fluxes agree well (mean±1 standard error: EC=17±4mg CH4m−2d−1; REA=19±3mg CH4m−2d−1, wavelet covariance=16±3mgCH4m−2d−1), and are not statistically different, but significantly exceed bottom-up inventory estimates of CH4 emissions based on animal husbandry (8±1mgCH4m−2d−1). The discrepancy between top-down and bottom-up estimates was found to be related to both increasing fractional area of wetlands in the flux footprint, and increasing surface temperature. For the case when the wetland area in the flux footprint was less than 10% fractional coverage, the top-down and bottom-up estimates were within the measurement error. This result provides the first independent verification of agricultural methane emissions inventories at the regional scale. Wavelet analysis, which provides spatially resolved fluxes, was used to attempt to separate CH4 emissions from managed and unmanaged CH4 sources. Opportunities to minimize the challenges of verifying agricultural CH4 emissions inventories using aircraft flux measuring systems are discussed.
The role of secondary circulations has recently been studied in the context of well-defined surface heterogeneity in a semiarid ecosystem where it was found that energy balance closure over a ...desert–forest system and the structure of the boundary layer was impacted by advection and flux divergence. As a part of the CliFF (“Climate feedbacks and benefits of semi-arid forests”, a collaboration between KIT, Germany, and the Weizmann Institute, Israel) campaign, we studied the boundary layer dynamics and turbulent transport of energy corresponding to this effect in Yatir Forest situated in the Negev Desert in Israel. The forest surrounded by small shrubs presents a distinct feature of surface heterogeneity, allowing us to study the differences between their interactions with the atmosphere above by conducting measurements with two eddy covariance (EC) stations and two Doppler lidars. As expected, the turbulence intensity and vertical fluxes of momentum and sensible heat are found to be higher above the forest compared to the shrubland. Turbulent statistics indicative of nonlocal motions are also found to differ over the forest and shrubland and also display a strong diurnal cycle. The production of turbulent kinetic energy (TKE) over the forest is strongly mechanical, while buoyancy effects generate most of the TKE over the shrubland. Overall TKE production is much higher above the forest compared to the shrubland. The forest is also found to be more efficient in dissipating TKE. The TKE budget appears to be balanced on average both for the forest and shrubland, although the imbalance of the TKE budget, which includes the role of TKE transport, is found to be quite different in terms of diurnal cycles for the forest and shrubland. The difference in turbulent quantities and the relationships between the components of TKE budget are used to infer the characteristics of the turbulent transport of energy between the desert and the forest.
Long-term tall-tower eddy-covariance (EC) measurements have been recently established in three European pilot cities as part of the ICOS-Cities project. We conducted a comparison of EC software to ...ensure a reliable generation of interoperable flux estimates, which is the prerequisite for avoiding methodological biases and improving the comparability of the results. We analyzed datasets covering 5 months collected from EC tall-tower installations located in urbanized areas of Munich, Zurich, and Paris. Fluxes of sensible heat, latent heat, and CO2 were calculated using three software packages (i.e., TK3, EddyPro, and eddy4R) to assess the uncertainty of flux estimations attributed to differences in implemented postprocessing schemes. A very good agreement on the mean values and standard deviations was found across all three sites, which can probably be attributed to a uniform instrumentation, data acquisition, and preprocessing. The overall comparison of final flux time series products showed a good but not yet perfect agreement among the three software packages. TK3 and EddyPro both calculated fluxes with low-frequency spectral correction, resulting in better agreement than between TK3 and the eddy4R workflow with disabled low-frequency spectral treatment. These observed flux discrepancies indicate the crucial role of treating low-frequency spectral loss in flux estimation for tall-tower EC systems.
Large-eddy simulation (LES) has become a well-established tool in the atmospheric boundary layer research community to study turbulence.
It allows three-dimensional realizations of the turbulent ...fields, which large-scale models and most experimental studies cannot yield.
To resolve the largest eddies in the mixed layer, a moderate grid resolution in the range of 10 to 100 m is often sufficient, and these simulations can be run on a computing cluster with a few hundred processors or even on a workstation for simple configurations.
The desired resolution is usually limited by the computational resources. However, to compare with tower measurements of turbulence and exchange fluxes in the surface layer, a much higher resolution is required.
In spite of the growth in computational power, a high-resolution LES of the surface layer is often not feasible: to fully resolve the energy-containing eddies near the surface, a grid spacing of O(1 m) is required.
One way to tackle this problem is to employ a vertical grid nesting technique, in which the surface is simulated at the necessary fine grid resolution, and it is coupled with a standard, coarse, LES that resolves the turbulence in the whole boundary layer.
We modified the LES model PALM (Parallelized Large-eddy simulation Model) and implemented a two-way nesting technique, with coupling in both directions between the coarse and the fine grid. The coupling algorithm has to ensure correct boundary conditions for the fine grid.
Our nesting algorithm is realized by modifying the standard third-order Runge–Kutta time stepping to allow communication of data between the two grids.
The two grids are concurrently advanced in time while ensuring that the sum of resolved and sub-grid-scale kinetic energy is conserved. We design a validation test and show that the temporally averaged profiles from the fine grid agree well compared to the reference simulation with high resolution in the entire domain.
The overall performance and scalability of the nesting algorithm is found to be satisfactory.
Our nesting results in more than 80 % savings in computational power for 5 times higher resolution in each direction in the surface layer.