AmeriFlux Sites Research Articles

Evaluation of the Uncertainty in Satellite-Based Crop State Variable Retrievals Due to Site and Growth Stage Specific Factors and Their Potential in Coupling with Crop Growth Models.

Coupling crop growth models and remote sensing provides the potential to improve our understanding of the genotype x environment x management (G × E × M) variability of crop growth on a global scale. Unfortunately, the uncertainty in the relationship between the satellite measurements and the crop state variables across different sites and growth stages makes it difficult to perform the coupling. In this study, we evaluate the effects of this uncertainty with MODIS data at the Mead, Nebraska Ameriflux sites (US-Ne1, US-Ne2, and US-Ne3) and accurate, collocated Hybrid-Maize (HM) simulations of leaf area index (LAI) and canopy light use efficiency (LUECanopy). The simulations are used to both explore the sensitivity of the satellite-estimated genotype × management (G × M) parameters to the satellite retrieval regression coefficients and to quantify the amount of uncertainty attributable to site and growth stage specific factors. Additional ground-truth datasets of LAI and LUECanopy are used to validate the analysis. The results show that uncertainty in the LAI/satellite measurement regression coefficients lead to large uncertainty in the G × M parameters retrievable from satellites. In addition to traditional leave-one-site-out regression analysis, the regression coefficient uncertainty is assessed by evaluating the retrieval performance of the temporal change in LAI and LUECanopy. The weekly change in LAI is shown to be retrievable with a correlation coefficient absolute value (|r|) of 0.70 and root-mean square error (RMSE) value of 0.4, which is significantly better than the performance expected if the uncertainty was caused by random error rather than secondary effects caused by site and growth stage specific factors (an expected |r| value of 0.36 and RMSE value of 1.46 assuming random error). As a result, this study highlights the importance of accounting for site and growth stage specific factors in remote sensing retrievals for future work developing methods coupling remote sensing with crop growth models.

Lessons Learned From Modeling Irrigation From Field to Regional Scales

AbstractCorrectly calculating the timing and amount of crop irrigation is crucial for capturing irrigation effects on surface water and energy budgets and land‐atmosphere interactions. This study incorporated a dynamic irrigation scheme into the Noah with multiparameterization land surface model and investigated three methods of determining crop growing season length by agriculture management data. The irrigation scheme was assessed at field scales using observations from two contrasting (irrigated and rainfed) AmeriFlux sites near Mead, Nebraska. Results show that crop‐specific growing‐season length helped capture the first application timing and total irrigation amount, especially for soybeans. With a calibrated soil‐moisture triggering threshold (IRR_CRI), using planting and harvesting dates alone could reasonably predict the first application for maize. For soybeans, additional constraints on growing season were required to correct an early bias in the first modeled application. Realistic leaf area index input was essential for identifying the leaf area index‐based growing season. When transitioning from field to regional scales, the county‐level calibrated IRR_CRI helped mitigate overestimated (underestimated) total irrigation amount in southeastern Nebraska (lower Mississippi River Basin). In these two heavily irrigated regions, irrigation produced a cooling effect of 0.8–1.4 K, a moistening effect of 1.2–2.4 g/kg, a reduction in sensible heat flux by 60–105 W/m2, and an increase in latent heat flux by 75–120 W/m2. Most of irrigation water was used to increase soil moisture and evaporation, rather than runoff. Lacking regional‐scale irrigation timing and crop‐specific parameters makes transferring the evaluation and parameter‐constraint methods from field to regional scales difficult.

The Arctic-Boreal vulnerability experiment model benchmarking system

NASA’s Arctic-Boreal Vulnerability Experiment (ABoVE) integrates field and airborne data into modeling and synthesis activities for understanding Arctic and Boreal ecosystem dynamics. The ABoVE Benchmarking System (ABS) is an operational software package to evaluate terrestrial biosphere models against key indicators of Arctic and Boreal ecosystem dynamics, i.e.: carbon biogeochemistry, vegetation, permafrost, hydrology, and disturbance. The ABS utilizes satellite remote sensing data, airborne data, and field data from ABoVE as well as collaborating research networks in the region, e.g.: the Permafrost Carbon Network, the International Soil Carbon Network, the Northern Circumpolar Soil Carbon Database, AmeriFlux sites, the Moderate Resolution Imaging Spectroradiometer, the Orbiting Carbon Observatory 2, and the Soil Moisture Active Passive mission. The ABS is designed to be interactive for researchers interested in having their models accurately represent observations of key Arctic indicators: a user submits model results to the system, the system evaluates the model results against a set of Arctic-Boreal benchmarks outlined in the ABoVE Concise Experiment Plan, and the user then receives a quantitative scoring of model strengths and deficiencies through a web interface. This interactivity allows model developers to iteratively improve their model for the Arctic-Boreal Region by evaluating results from successive model versions. We show here, for illustration, the improvement of the Lund–Potsdam–Jena-Wald Schnee und Landschaft (LPJwsl) version model through the ABoVE ABS as a new permafrost module is coupled to the existing model framework. The ABS will continue to incorporate new benchmarks that address indicators of Arctic-Boreal ecosystem dynamics as they become available.

Investigating old-growth ponderosa pine physiology using tree-rings, δ13 C, δ18 O, and a process-based model.

In dealing with predicted changes in environmental conditions outside those experienced today, forest managers and researchers rely on process‐based models to inform physiological processes and predict future forest growth responses. The carbon and oxygen isotope ratios of tree‐ring cellulose (δ13Ccell, δ18Ocell) reveal long‐term, integrated physiological responses to environmental conditions. We incorporated a submodel of δ18Ocell into the widely used Physiological Principles in Predicting Growth (3‐PG) model for the first time, to complement a recently added δ13Ccell submodel. We parameterized the model using previously reported stand characteristics and long‐term trajectories of tree‐ring growth, δ13Ccell, and δ18Ocell collected from the Metolius AmeriFlux site in central Oregon (upland trees). We then applied the parameterized model to a nearby set of riparian trees to investigate the physiological drivers of differences in observed basal area increment (BAI) and δ13Ccell trajectories between upland and riparian trees. The model showed that greater available soil water and maximum canopy conductance likely explain the greater observed BAI and lower δ13Ccell of riparian trees. Unexpectedly, both observed and simulated δ18Ocell trajectories did not differ between the upland and riparian trees, likely due to similar δ18O of source water isotope composition. The δ18Ocell submodel with a Peclet effect improved model estimates of δ18Ocell because its calculation utilizes 3‐PG growth and allocation processes. Because simulated stand‐level transpiration (E) is used in the δ18O submodel, aspects of leaf‐level anatomy such as the effective path length for transport of water from the xylem to the sites of evaporation could be estimated.

LandMOD ET mapper: A new matlab-based graphical user interface (GUI) for automated implementation of SEBAL and METRIC models in thermal imagery

Remote sensing based surface energy balance (SEB) models, such as SEB Algorithm for Land (SEBAL) and Mapping evapotranspiration (ET) at high Resolution with Internalized Calibration (METRIC) are widely used to monitor ET. However, their operational use is challenged by model complexity and the tedious process of selecting hot and cold pixels. We developed a new Matlab-based ET mapping toolbox named “LandMOD ET Mapper” for automated implementation of these models to facilitate their widespread applications among new users. The toolbox is demonstrated using one Landsat and one Moderate Resolution Imaging Spectroradiometer (MODIS) Terra image covering three AmeriFlux sites in Nebraska, U.S. Daily ET estimates were comparable with the observed ET for five days when near cloud-free Landsat and MODIS Terra images were available in the 2006 maize growing season. The toolbox provides a viable ET mapping option to new and inexperienced users with reduced time and labor demands.

Impact of the Revisit of Thermal Infrared Remote Sensing Observations on Evapotranspiration Uncertainty—A Sensitivity Study Using AmeriFlux Data

Thermal infrared remote sensing observations have been widely used to provide useful information on surface energy and water stress for estimating evapotranspiration (ET). However, the revisit time of current high spatial resolution (<100 m) thermal infrared remote sensing systems, sixteen days for Landsat for example, can be insufficient to reliably derive ET information for water resources management. We used in situ ET measurements from multiple Ameriflux sites to (1) evaluate different scaling methods that are commonly used to derive daytime ET estimates from time-of-day observations; and (2) quantify the impact of different revisit times on ET estimates at monthly and seasonal time scales. The scaling method based on a constant evaporative ratio between ET and the top-of-atmosphere solar radiation provided slightly better results than methods using the available energy, the surface solar radiation or the potential ET as scaling reference fluxes. On average, revisit time periods of 2, 4, 8 and 16 days resulted in ET uncertainties of 0.37, 0.55, 0.73 and 0.90 mm per day in summer, which represented 13%, 19%, 23% and 31% of the monthly average ET calculated using the one-day revisit dataset. The capability of a system to capture rapid changes in ET was significantly reduced for return periods higher than eight days. The impact of the revisit on ET depended mainly on the land cover type and seasonal climate, and was higher over areas with high ET. We did not observe significant and systematic differences between the impacts of the revisit on monthly ET estimates that are based on morning or afternoon observations. We found that four-day revisit scenarios provided a significant improvement in temporal sampling to monitor surface ET reducing by around 40% the uncertainty of ET products derived from a 16-day revisit system, such as Landsat for instance.

Evaluating the Spatial Representativeness of the MODerate Resolution Image Spectroradiometer Albedo Product (MCD43) at AmeriFlux Sites

Land surface albedo is a key parameter in regulating surface radiation budgets. The gridded remote sensing albedo product often represents information concerning an area larger than the nominal spatial resolution because of the large viewing angles of the observations. It is essential to quantify the spatial representativeness of remote sensing products to better guide the sampling strategy in field experiments and match products from different sources. This study quantifies the spatial representativeness of the MODerate Resolution Image Spectroradiometer (MODIS) (collection V006) 500 m daily albedo product (MCD43A3) using the high-resolution product as intermediate data for different land cover types. A total of 1820 paired high-resolution Landsat Thematic Mapper (TM) and coarse-resolution (MODIS) albedo data from five land cover types were used. The TM albedo data was used as the spatial-complete high resolution data to evaluate the spatial representativeness of the MODIS albedo product. Semivarioagrams were estimated from 30 m Landsat data at different spatial scales. Surface heterogeneity was evaluated with sill value and relative coefficient of variation. The 30 m Landsat albedo data was aggregated to 450 m–1800 m using two different methods and compared with MODIS albedo product. The spatial representativeness of MODIS albedo product was determined according to the surface heterogeneity and the consistency of MODIS data and the aggregated TM value. Results indicated that for evergreen broadleaf forests, deciduous broadleaf forests, open shrub lands, woody savannas and grasslands, the MODIS 500 m daily albedo product represents a spatial scale of approximately 630 m. For mixed forests and croplands, the representative spatial scale was approximately 690 m. The difference obtained was primarily because of the complexity of the landscape structure. For mixed forests and croplands, the structure of the landscape was relatively complex due to the presence of different forest and plant types in the pixel area, whereas the other landscape structures were considerably simpler.

Bayesian Analyses of 17 Winters of Water Vapor Fluxes Show Bark Beetles Reduce Sublimation

AbstractSublimation is an important hydrological flux in cold, snow‐dominated ecosystems. In high‐elevation spruce‐fir forests of western North America, spruce beetle outbreaks have killed trees, reduced the canopy, and altered processes that control sublimation. We evaluated two hypotheses related to effects of disturbance on sublimation in this ecosystem: (1) the dominant source for sublimation is canopy intercepted snow and (2) the loss of canopy following a beetle disturbance leads to less total sublimation. To incorporate uncertainty hierarchically across multiple data sources and address phenomenological parsimony, Bayesian statistics were used to analyze 17 years (2000–2016) of winter eddy covariance flux data at the Glacier Lakes Ecosystem Experiments Sites AmeriFlux sites where a spruce beetle outbreak caused 75–85% basal area mortality. Our analysis revealed that resistances to sublimate snow from the canopy were an order of magnitude less than from the snowpack, and the maximum snow loading capacity in disturbed canopies was reduced to 34% of its pre‐outbreak value. Total sublimation has decreased since 2010, 2 years after the main outbreak, declining 24% (with a 95% credible interval, C.I., between 18% and 38%) during 2014–2016 due to a 32% decrease in canopy sublimation. Snowpack sublimation only increased 3% over this period. With less total sublimation, the forest retained 6.1% (4.5–12.3% C.I.) more snowpack mass or equivalently 4.4% (3.2–8.8 C.I.) of the annual precipitation. Considering tree growth and ecological succession are slow in spruce‐fir forests, this decrease in sublimation should persist as an increased snowpack for decades, with substantial impacts on catchment hydrologic processes and potentially streamflow.

Evaluating a new algorithm for satellite-based evapotranspiration for North American ecosystems: Model development and validation

We introduce a different operational approach to estimate 8-day average daily evapotranspiration (ET) using both routinely available data and the Penman-Monteith (P-M) equation for canopy transpiration and evaporation of intercepted water and Priestley and Taylor for soil evaporation. Our algorithm considered the environmental constraints on canopy resistance and ET by (1) including vapor pressure deficit (VPD), incoming solar radiation, soil moisture, and temperature constraints on stomatal conductance; (2) using leaf area index (LAI) to scale from the leaf to canopy conductance; and (3) calculating canopy resistance as a function of environmental variables such as net radiation and VPD. Remote sensing data from the Moderate Resolution Spectroradiometer (MODIS) and satellite soil moisture data were used to derive the ET model. The algorithm was calibrated and evaluated using measured ET data from 20 AmeriFlux Eddy covariance flux sites for the period of 2003–2012. We found good agreements between our 8-day ET estimates and observations with mean absolute error (MAE) ranges from 0.17 mm/day to 0.94 mm/day compared with MAE ranging from 0.28 mm/day to 1.50 mm/day for MODIS ET. Compared to MODIS ET, our proposed algorithm has higher correlations and higher Willmott’s index of agreement with observations for the majority of the Ameriflux sites. The strong relationship between the model estimated ET and the flux tower observations implies that our model has the potential to be applied to different ecosystems and at different temporal scales.

Advancing Terrestrial Ecosystem Science With a Novel Automated Measurement System for Sun‐Induced Chlorophyll Fluorescence for Integration With Eddy Covariance Flux Networks

AbstractSun‐induced chlorophyll fluorescence (SIF) provides critical information on the dynamics of gross primary productivity, a unique role not readily achievable using other methods. Long‐term continuous SIF observations have the potential to advance terrestrial ecosystem science. Realizing this potential, however, requires synergistic implementation of SIF measurements within eddy covariance (EC) flux networks. There is a need for SIF systems that can integrate seamlessly with EC instrumentation to maximize synergistic use of obtained data. Here we introduce the Fluorescence Auto‐Measurement Equipment (FAME) and protocol that fulfill such a purpose. FAME is designed specifically for plug‐and‐play integration with existing EC data acquisition systems. Its innovative hardware and software designs provide versatility, extensibility, autonomous operation, and ease of maintenance for acquiring SIF data of high quality and quantity. A major novel feature of FAME is its synchronized sampling of spectral irradiance and environmental variables, allowing for more precise interpretation of the SIF signal. FAME has been deployed since September 2016 at the Missouri Ozark AmeriFlux site, providing high‐quality measurements even when air temperatures approached 40 °C. Results reveal that canopy SIF saturated or even slightly decreased at high light, similar to leaf‐level photosynthesis. Clear diurnal hysteresis was observed: For the same light, morning SIF was higher than afternoon. Dynamic energy dissipation processes and stress‐induced movements of chloroplasts and leaves may explain the observed pattern. The technology and measurement protocol introduced here advances the coordinated observation of SIF and EC fluxes and represents a step change in observational ecosystem and carbon cycle science research.

Development of the DayCent-Photo model and integration of variable photosynthetic capacity

We integrated a photosynthetic sub-model into the daily Century model (DayCent) to improve the estimations of seasonal changes in carbon fluxes at the Niwot Ridge LTER site and the Harvard forest LTER site (DayCent-Photo). The photosynthetic sub-model, adapted from the SIPNET/PnET family of models, includes solar radiation and vapor pressure deficit controls on production, as well as temperature and water stress terms. A key feature we added to the base photosynthetic equations is the addition of a variable maximum net photosynthetic rate (Amax). We optimized the parameters controlling photosynthesis using a variation of the Metropolis-Hastings algorithm along with data-assimilation techniques. The model was optimized and validated against observed net ecosystem exchange (NEE) and estimated gross primary production (GPP) and ecosystem respiration (RESP) values for AmeriFlux sites at Niwot Ridge and Harvard forest. The inclusion of a variable Amax rate greatly improved model performance (NEE RMSE = 0.63 gC·m–2, AIC = 2099) versus a version with a single Amax parameter (NEE RMSE = 0.74 gC·m–2, AIC = 3724). DayCent-Photo was able to capture the inter-annual and seasonal flux patterns for NEE, GPP, ecosystem respiration (RESP), and daily actual evapotranspiration (AET), but tended to overestimate yearly NEE uptake. The DayCent-Photo model has been successfully set up to simulate daily NEE, GPP, RESP, and AET for deciduous forest, conifer forests, and grassland systems in the US using AmeriFlux data sets and has recently been improved to include the impact of UV radiation surface litter decay (DayCent-UV). The simulated influence of a variable Amax rate suggests a need for further studies on the process controls affecting the seasonal photosynthetic rates. The results for all of the forest and grassland sites show that maximum Amax values occurs early during the growing period and taper off toward the end of the growing season.

Influences of Root Hydraulic Redistribution on N2O Emissions at AmeriFlux Sites

AbstractIt has long been suspected that root hydraulic redistribution (HR) affects the carbon and nitrogen cycles. Nitrous oxide (N2O) is an important greenhouse gas and is the primary stratospheric ozone‐depleting substance. To our knowledge, the influences of HR on N2O emissions have not been investigated. Here we use the HR schemes of Ryel et al. and Amenu and Kumar incorporated into CLM4.5 to examine N2O emissions at five AmeriFlux sites. The results show that HR reduced N2O emissions by 28–92% in the four natural ecosystems experiencing a dry season, whereas it had a very limited effect on the Corn Belt site that has strong emissions but with no distinct dry season. We hypothesize that N2O emissions in ecosystems with a distinct dry season are likely overestimated by CENTURY‐based Earth system models.

Hydraulic redistribution affects modeled carbon cycling via soil microbial activity and suppressed fire.

Hydraulic redistribution (HR) of water from moist to drier soils, through plant roots, occurs world-wide in seasonally dry ecosystems. Although the influence of HR on landscape hydrology and plant water use has been amply demonstrated, HR's effects on microbe-controlled processes sensitive to soil moisture, including carbon and nutrient cycling at ecosystem scales, remain difficult to observe in the field and have not been integrated into a predictive framework. We incorporated a representation of HR into the Community Land Model (CLM4.5) and found the new model improved predictions of water, energy, and system-scale carbon fluxes observed by eddy covariance at four seasonally dry yet ecologically diverse temperate and tropical AmeriFlux sites. Modeled plant productivity and microbial activities were differentially stimulated by upward HR, resulting at times in increased plant demand outstripping increased nutrient supply. Modeled plant productivity and microbial activities were diminished by downward HR. Overall, inclusion of HR tended to increase modeled annual ecosystem uptake of CO2 (or reduce annual CO2 release to the atmosphere). Moreover, engagement of CLM4.5's ground-truthed fire module indicated that though HR increased modeled fuel load at all four sites, upward HR also moistened surface soil and hydrated vegetation sufficiently to limit the modeled spread of dry season fire and concomitant very large CO2 emissions to the atmosphere. Historically, fire has been a dominant ecological force in many seasonally dry ecosystems, and intensification of soil drought and altered precipitation regimes are expected for seasonally dry ecosystems in the future. HR may play an increasingly important role mitigating development of extreme soil water potential gradients and associated limitations on plant and soil microbial activities, and may inhibit the spread of fire in seasonally dry ecosystems.

A method to calculate a locally relevant water year for ecohydrological studies using eddy covariance data

AbstractIncreased rainfall variability due to climate change significantly impacts carbon and water cycling in ecosystems, but these impacts may be masked when using arbitrary annual reporting periods such as the calendar year, which may not have any relevance to natural annual ecosystem processes. A variety of alternative annual integration periods have been described for specific purposes or locations, but are of limited general applicability. Here, we present an eddy covariance data‐driven empirical method to determine a locally relevant annual time period. The method selects a start date for a locally relevant water year (LRWY) that maximizes correlation between annual precipitation (AP), and annual evapotranspiration (AET) and annual gross primary production. The method was tested using data from 2004 to 2013 for 2 Ameriflux sites (woodland and grassland) in Central Texas. The timeframe included periods of unusually high rainfall and periods of extreme drought. The highest correlation between AP, and AET and annual gross primary production was obtained with an LRWY beginning in mid‐September. Use of the LRWY better captured the impact of soil water recharge in the autumn on photosynthesis the following spring than did calendar years. Use of the LRWY also identified more annual periods in which AET exceeded AP, which more accurately reflected the impact of drought on ecosystems processes than did analysis based on calendar years.

Regional evapotranspiration from an image-based implementation of the Surface Temperature Initiated Closure (STIC1.2) model and its validation across an aridity gradient in the conterminous US

Abstract. Recent studies have highlighted the need for improved characterizations of aerodynamic conductance and temperature (gA and T0) in thermal remote-sensing-based surface energy balance (SEB) models to reduce uncertainties in regional-scale evapotranspiration (ET) mapping. By integrating radiometric surface temperature (TR) into the Penman–Monteith (PM) equation and finding analytical solutions of gA and T0, this need was recently addressed by the Surface Temperature Initiated Closure (STIC) model. However, previous implementations of STIC were confined to the ecosystem-scale using flux tower observations of infrared temperature. This study demonstrates the first regional-scale implementation of the most recent version of the STIC model (STIC1.2) that integrates the Moderate Resolution Imaging Spectroradiometer (MODIS) derived TR and ancillary land surface variables in conjunction with NLDAS (North American Land Data Assimilation System) atmospheric variables into a combined structure of the PM and Shuttleworth–Wallace (SW) framework for estimating ET at 1 km × 1 km spatial resolution. Evaluation of STIC1.2 at 13 core AmeriFlux sites covering a broad spectrum of climates and biomes across an aridity gradient in the conterminous US suggests that STIC1.2 can provide spatially explicit ET maps with reliable accuracies from dry to wet extremes. When observed ET from one wet, one dry, and one normal precipitation year from all sites were combined, STIC1.2 explained 66 % of the variability in observed 8-day cumulative ET with a root mean square error (RMSE) of 7.4 mm/8-day, mean absolute error (MAE) of 5 mm/8-day, and percent bias (PBIAS) of −4 %. These error statistics showed relatively better accuracies than a widely used but previous version of the SEB-based Surface Energy Balance System (SEBS) model, which utilized a simple NDVI-based parameterization of surface roughness (zOM), and the PM-based MOD16 ET. SEBS was found to overestimate (PBIAS = 28 %) and MOD16 was found to underestimate ET (PBIAS = −26 %). The performance of STIC1.2 was better in forest and grassland ecosystems as compared to cropland (20 % underestimation) and woody savanna (40 % overestimation). Model inter-comparison suggested that ET differences between the models are robustly correlated with gA and associated roughness length estimation uncertainties which are intrinsically connected to TR uncertainties, vapor pressure deficit (DA), and vegetation cover. A consistent performance of STIC1.2 in a broad range of hydrological and biome categories, as well as the capacity to capture spatio-temporal ET signatures across an aridity gradient, points to the potential for this simplified analytical model for near-real-time ET mapping from regional to continental scales.

Modeling canopy-induced turbulence in the Earth system: a unified parameterization of turbulent exchange within plant canopies and the roughness sublayer (CLM-ml v0)

Abstract. Land surface models used in climate models neglect the roughness sublayer and parameterize within-canopy turbulence in an ad hoc manner. We implemented a roughness sublayer turbulence parameterization in a multilayer canopy model (CLM-ml v0) to test if this theory provides a tractable parameterization extending from the ground through the canopy and the roughness sublayer. We compared the canopy model with the Community Land Model (CLM4.5) at seven forest, two grassland, and three cropland AmeriFlux sites over a range of canopy heights, leaf area indexes, and climates. CLM4.5 has pronounced biases during summer months at forest sites in midday latent heat flux, sensible heat flux, gross primary production, nighttime friction velocity, and the radiative temperature diurnal range. The new canopy model reduces these biases by introducing new physics. Advances in modeling stomatal conductance and canopy physiology beyond what is in CLM4.5 substantially improve model performance at the forest sites. The signature of the roughness sublayer is most evident in nighttime friction velocity and the diurnal cycle of radiative temperature, but is also seen in sensible heat flux. Within-canopy temperature profiles are markedly different compared with profiles obtained using Monin–Obukhov similarity theory, and the roughness sublayer produces cooler daytime and warmer nighttime temperatures. The herbaceous sites also show model improvements, but the improvements are related less systematically to the roughness sublayer parameterization in these canopies. The multilayer canopy with the roughness sublayer turbulence improves simulations compared with CLM4.5 while also advancing the theoretical basis for surface flux parameterizations.

Estimation of Turbulent Heat Fluxes by Assimilation of Land Surface Temperature Observations From GOES Satellites Into an Ensemble Kalman Smoother Framework

AbstractIn different studies, land surface temperature (LST) observations have been assimilated into the variational data assimilation (VDA) approaches to estimate turbulent heat fluxes. The VDA methods yield accurate turbulent heat fluxes, but they need an adjoint model, which is difficult to derive and code. They also cannot directly calculate the uncertainty of their estimates. To overcome the abovementioned drawbacks, this study assimilates LST data from Geostationary Operational Environmental Satellite into the ensemble Kalman smoother (EnKS) data assimilation system to estimate turbulent heat fluxes. EnKS does not need to derive the adjoint term and directly generates statistical information on the accuracy of its predictions. It uses the heat diffusion equation to simulate LST. EnKS with the state augmentation approach finds the optimal values for the unknown parameters (i.e., evaporative fraction and neutral bulk heat transfer coefficient, CHN) by minimizing the misfit between LST observations from Geostationary Operational Environmental Satellite and LST estimations from the heat diffusion equation. The augmented EnKS scheme is tested over six Ameriflux sites with a wide range of hydrological and vegetative conditions. The results show that EnKS can predict not only the model parameters and turbulent heat fluxes but also their uncertainties over a variety of land surface conditions. Compared to the variational method, EnKS yields suboptimal turbulent heat fluxes. However, suboptimality of EnKS is small, and its results are comparable to those of the VDA method. Overall, EnKS is a feasible and reliable method for estimation of turbulent heat fluxes.

Merging a mechanistic enzymatic model of soil heterotrophic respiration into an ecosystem model in two AmeriFlux sites of northeastern USA

Heterotrophic respiration (Rh), microbial processing of soil organic matter to carbon dioxide (CO2), is a major, yet highly uncertain, carbon (C) flux from terrestrial systems to the atmosphere. Temperature sensitivity of Rh is often represented with a simple Q10 function in ecosystem models and earth system models (ESMs), sometimes accompanied by an empirical soil moisture modifier. More explicit representation of the effects of soil moisture, substrate supply, and their interactions with temperature has been proposed as a way to disentangle the confounding factors of apparent temperature sensitivity of Rh and improve the performance of ecosystem models and ESMs. The objective of this work was to insert into an ecosystem model a more mechanistic, but still parsimonious, model of environmental factors controlling Rh and evaluate the model performance in terms of soil and ecosystem respiration. The Dual Arrhenius and Michaelis-Menten (DAMM) model simulates Rh using Michaelis-Menten, Arrhenius, and diffusion functions. Soil moisture affects Rh and its apparent temperature sensitivity in DAMM by regulating the diffusion of oxygen, soluble C substrates, and extracellular enzymes to the enzymatic reaction site. Here, we merged the DAMM soil flux model with a parsimonious ecosystem flux model, FöBAAR (Forest Biomass, Assimilation, Allocation and Respiration). We used high-frequency soil flux data from automated soil chambers and landscape-scale ecosystem fluxes from eddy covariance towers at two AmeriFlux sites (Harvard Forest, MA and Howland Forest, ME) in the northeastern USA to estimate parameters, validate the merged model, and to quantify the uncertainties in a multiple constraints approach. The optimized DAMM-FöBAAR model better captured the seasonal and inter-annual dynamics of soil respiration (Soil R) compared to the FöBAAR-only model for the Harvard Forest, where higher frequency and duration of drying events significantly regulate substrate supply to heterotrophs. However, DAMM-FöBAAR showed improvement over FöBAAR-only at the boreal transition Howland Forest only in unusually dry years. The frequency of synoptic-scale dry periods is lower at Howland, resulting in only brief water limitation of Rh in some years. At both sites, the declining trend of soil R during drying events was captured by the DAMM-FöBAAR model; however, model performance was also contingent on site conditions, climate, and the temporal scale of interest. While the DAMM functions require a few more parameters than a simple Q10 function, we have demonstrated that they can be included in an ecosystem model and reduce the model-data mismatch. Moreover, the mechanistic structure of the soil moisture effects using DAMM functions should be more generalizable than the wide variety of empirical functions that are commonly used, and these DAMM functions could be readily incorporated into other ecosystem models and ESMs.

Evapotranspiration partitioning for three agro-ecosystems with contrasting moisture conditions: a comparison of an isotope method and a two-source model calculation

Quantification of the contribution of transpiration (T) to evapotranspiration (ET) is a requirement for understanding changes in carbon assimilation and water cycling in a changing environment. So far, few studies have examined seasonal variability of T/ET and compared different ET partitioning methods under natural conditions across diverse agro-ecosystems. In this study, we apply a two-source model to partition ET for three agro-ecosystems (rice, wheat and corn). The model is coupled with a plant physiology scheme for the canopy conductance. The model-estimated T/ET ranges from 0 to 1, with a near continuous increase over time in the early growing season when leaf area index (LAI) is less than 2.5 and then convergence towards a stable value beyond LAI of 2.5. The seasonal change in T/ET can be described well as a function of LAI, implying that LAI is a first-order factor affecting ET partitioning. Application of the model to seven other Ameriflux sites reveals that soil moisture and canopy conductance also influence the ET partitioning. The two-source model results show that the growing-season (May–September for rice, April–June for wheat and June–September for corn) T/ET is 0.50, 0.84 and 0.64, while an isotopic approach shows that T/ET is 0.74, 0.93 and 0.81 for rice, wheat and maize, respectively. The two-source model results are supported by soil lysimeter and eddy covariance measurements made during the same time period for wheat (0.87). Uncertainty analysis suggests that further improvements to the Craig-Gordon model prediction of the evaporation isotope composition and to measurement of the isotopic composition of ET are necessary to achieve accurate flux partitioning at the ecosystem scale using water isotopes as tracers.The program code for the two-source model is available at the open-source platform https://github.com/zhongwangwei/SiLSM_v3.

Evaluation of ozone deposition models over a subalpine forest in Niwot Ridge, Colorado

In this study, we evaluated three conceptually similar ozone gas deposition models. These dry deposition models are frequently used with chemical transport models for calculations over large spatial domains. However, large scale applications of surface-atmosphere exchange of reactive gases require modeling results as accurate as possible to avoid nonlinear accumulation of errors in the spatially representative results. In this paper, model evaluation and comparison against measured data over a coniferous forest at Niwot Ridge AmeriFlux site (Colorado, USA) is carried out. At this site, no previous model calibration took place for any of the models, therefore, we can test and compare their performances under similar conditions as they would perform in a spatial application. Our results show systematic model errors in all the three cases, model performance varies with time of the day, and the errors show a pronounced seasonal pattern as well. The introduction of soil moisture content stress in the model improved model performance regarding the magnitude of fluxes, but the correlation between measured and modeled ozone deposition values remains low. Our results suggest that ozone dry deposition model results should be interpreted carefully in large scale applications, where the accuracy can vary with land cover sometimes are biased. © 2018, Hungarian Meteorological Service. All rights reserved.