Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Indexes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78

Sources of Subseasonal-To-Seasonal Predictability of Atmospheric Rivers and Precipitation in the Western United States

Atmospheric rivers (ARs) account for a large portion of winter precipitation in the western US. To evaluate the sources of AR and precipitation predictability at subseasonal-to-seasonal timescales, we examine the relationships between two climate modes, the El Niño/Southern Oscillation (ENSO) and Madden-Julian oscillation (MJO), and winter hydroclimate in the western US. Our analysis uses a large ensemble of the Weather Research and Forecast (WRF) model simulations from 1981 to 2017, facilitating the assessment of uncertainty in climate mode-AR relationship due to climate variability and short length of observations. Over the North Pacific basin, we find ENSO-related latitudinal shifts of ARs, but ENSO has little effect on basin-wide averaged AR frequency. Over the western US, there is some uncertainty in the ENSO-AR connection due to the impacts of ENSO indices and datasets. However, extreme El Niño events defined by the ENSO longitude index are consistently linked to increased landfalling AR activity. The MJO can lead to significantly enhanced or suppressed landfalling AR activity in the large-ensembles, depending on the phase of MJO and time lag, while observations are too short in length to robustly show this signal. ENSO substantially modulates the MJO-AR relationship, triggering variegated responses of landfalling ARs and AR precipitation in La Niña and El Niño years. Our findings highlight the need to evaluate concurrent effects of different climate modes on ARs and precipitation, and may shed light on a path toward more accurate subseasonal-to-seasonal prediction of ARs and precipitation over the western US.

54 ENVIRONMENTAL SCIENCES↗

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗

Tower-Based Remote Sensing Reveals Mechanisms Behind a Two-phased Spring Transition in a Mixed-Species Boreal Forest

The boreal forest is a major contributor to the global climate system, therefore, reducing uncertainties in how the forest will respond to a changing climate is critical. One source of uncertainty is the timing and drivers of the spring transition. Remote sensing can provide important information on this transition, but persistent foliage greenness, seasonal snow cover, and a high prevalence of mixed forest stands (both deciduous and evergreen species) complicate interpretation of these signals. We collected tower-based remotely sensed data (reflectance-based vegetation indices and Solar-Induced Chlorophyll Fluorescence [SIF]), stem radius measurements, gross primary productivity, and environmental conditions in a boreal mixed forest stand. Evaluation of this data set shows a two-phased spring transition. The first phase is the reactivation of photosynthesis and transpiration in evergreens, marked by an increase in relative SIF, and is triggered by thawed stems, warm air temperatures, and increased available soil moisture. The second phase is a reduction in bulk photoprotective pigments in evergreens, marked by an increase in the Chlorophyll-Carotenoid Index. Deciduous leaf-out occurs during this phase, marked by an increase in all remotely sensed metrics. The second phase is controlled by soil thaw. Our results demonstrate that remote sensing metrics can be used to detect specific physiological changes in boreal tree species during the spring transition. The two-phased transition explains inconsistencies in remote sensing estimates of the timing and drivers of spring recovery. Our results imply that satellite-based observations will improve by using a combination of vegetation indices and SIF, along with species distribution information.

54 ENVIRONMENTAL SCIENCES↗

Exploring Hydrologic Model Process Connectivity at the Continental Scale Through an Information Theory Approach

Exploring water fluxes between hydrological model (HM) components is essential to assess and improve model realism. Many classical metrics for HM diagnosis rely solely on streamflow and hence provide limited insights into model performance across processes. This study applies an information theory measure known as “transfer entropy” (TE) to systematically quantify the transfer of information among major HM components. To test and demonstrate the benefits of TE, we use the Framework for Understanding Structural Errors (FUSE) model to mimic and compare four commonly used HM structures, VIC, PRMS, SACRAMENTO, and TOPMODEL, across 671 catchments spanning a variety of hydrologic regimes in the conterminous United States. We explore connections between HM components and catchment landscape characteristics (e.g., climate, topography, soil, and vegetation) and characterize their nonlinear associations using distance correlation and Spearman correlation coefficients. Our results indicate that while the information transferred from precipitation to runoff is similar across model structures (likely as a result of calibration), the information transferred among other components can vary significantly from a FUSE structure to another. We find that aridity, precipitation duration and frequency, snow fraction, mean elevation, forest area, and leaf area index are often significantly associated with TE between the main HM components. Here, we propose that the presence of meaningful nonlinear associations can be used to diagnose process representation in HMs. Our results highlight the necessity to enhance the conventional streamflow-only calibration approach for a more realistic representation of water dynamics in the models.

large-sample hydrology↗

Impact of vegetation physiology and phenology on watershed hydrology in a semiarid watershed in the Pacific Northwest in a changing climate

Concurrent with changes in climate conditions, increases in atmospheric composition and transport can affect land surface processes by regulating physiological and phenological responses of plants, and consequently perturb water, carbon, and nitrogen (N) dynamics at multiple scales. Howerver, the relationships bwtween plants physiological and phenological and land surface processes in a changing climate are less well known due to the limitation of representing palnts properties in traditional hydrological models. In this study, we use the Community Land Model version 5 (CLM5) to investigate how projected changes in large scale atmospheric conditions can affect hydrologic and biogeochemical dynamics in the Upper Columbia-Priest Rapids (UCPR) watershed, a typical semi-arid watershed located in the north-western United States dominated by croplands. Six numerical experiments considering climate change, CO2 concentration, and N deposition consistent with the Representative Concentration Pathways (RCPs) 4.5 and 8.5 are performed at 1-km resolution over the period of 2015-2099. Our results show that, compared with historical results (1985-2014), the elevated CO2 directly increases plants N uptake rates that provide additional leaf N to enhance photosynthesis by the end of 21st century (2070-2099). The higher GPP leads to changes in stomatal resistance and increased leaf area index (LAI) that affect evaporatranspiration (ET) and its components. As a result of stomatal responses to a higher CO2 concentration, ET over corn fields declines by 24.6 mm yr-1 under RCP4.5, and by 64.0 mm yr-1 under RCP8.5, respectively by the end of the century from its historical mean of 428 mm yr-1. Meanwhile, ET over soybean fields increases by 27.7 mm yr-1 under RCP4.5 and 9.98 mm yr-1 under RCP8.5 from its historical mean of 507 mm yr-1 as a result of CO2 fertilization as soybeans are able to uptake more N to increase LAI significantly. Furtherfore, the changes in physiological and phenological in corn and soybean can strongly modulate water budget terms over the watershed in response to expected changes in greenhouse gas concentrations and atmospheric composition. Additionally, increases in CO2 concentration and N deposition are beneficial for improving water use efficiency of corn and soybean, therefore lead to less water consumption per unit gram of carbon produced from these crops. Our study elucidates the need to jointly consider effects of surface climate and atmospheric composition and transport in understanding and predicting watershed hydrologic and biogeochemical responses, especially for an agricultural watershed.

Zhu, Bowen↗

Vegetation Response to Rising CO 2 Amplifies Contrasts in Water Resources Between Global Wet and Dry Land Areas

Rising atmospheric CO 2 impacts on vegetation physiological processes can alter land feedbacks on precipitation and water resources, but understanding of regional differences in these changes is uncertain. We investigate the impact of rising CO 2 on land water resources for different wetness levels using four Earth system models. We find an overall tendency of runoff to increase across all wetness levels. However, runoff increases in wet regions are much larger than those in dry regions, especially in wet seasons. This substantial amplification of contrasts between wet and dry regions increases at 3% per 100 ppm increase in CO 2 relative to the historical period, reaching 18% for a quadrupling of CO 2 , quantified by a new wetting contrast index (WCI). Physiological effects suppress evapotranspiration more in wet than dry regions, which has a larger contribution than radiative forcing to the amplification of runoff contrast, reshaping the spatial distribution of future land water resources.

58 GEOSCIENCES↗

Warm Phase of AMV Damps ENSO Through Weakened Thermocline Feedback

Abstract Interactions between ocean basins affect El Niño–Southern Oscillation (ENSO), altering its impacts on society. Here, we explore the effect of Atlantic Multidecadal Variability (AMV) on ENSO behavior using idealized experiments performed with the NCAR‐CESM1 model. Comparing warm (AMV+) to cold (AMV−) AMV conditions, we find that ENSO sea surface temperature (SST) anomalies are reduced by ∼10% and ENSO precipitation anomalies are shifted to the west during El Niño and east during La Niña. Using the Bjerknes stability index, we attribute the reduction in ENSO variability to a weakened thermocline feedback in boreal autumn. In AMV+, the Walker circulation and trade winds strengthen over the tropical Pacific, increasing the background zonal SST gradient. The background changes shift ENSO anomalies westwards, with wind stress anomalies more confined to the west. We suggest the changes in ENSO‐wind stress decrease the strength of the thermocline feedback in the east, eventually reducing ENSO growth rate.

Trascasa‐Castro, Paloma↗

Seasonality of Interbasin SST Contributions to Atlantic Tropical Cyclone Activity

Abstract Recent studies have demonstrated that the difference in sea surface temperature anomalies (SSTAs) between the tropical Atlantic main development region (MDR) and the tropical Pacific (Niño 3) modulates Atlantic tropical cyclone activity. This study further explores the seasonality of Pacific and Atlantic contributions to Atlantic hurricane activity. Our analysis shows that while MDR and Niño 3 SSTAs are equally important for late‐season (September–November) activity, early‐season (June‐August) activity is largely modulated by MDR SSTAs. This reflects the increased (reduced) variance of MDR (Niño 3) SSTAs in the early‐season due to their phase locking to the seasonal cycle. Further analysis yields skillful forecasts using an MDR‐Niño 3 interbasin index derived from hindcasts of the North American Multi‐Model Ensemble with May initial conditions. However, the prediction skill for MDR SSTAs is lower than that of Niño 3 SSTAs, suggesting that increasing the prediction skill for MDR SSTAs is key to improving seasonal outlooks.

West, Robert↗

Earth System Model Overestimation of Cropland Temperatures Scales With Agricultural Intensity

Intensive crop growth can modify regional climate by partitioning energy to latent heating through transpiration, cooling growing season temperatures. Recent work shows that cooling associated with agriculture can dampen anthropogenic warming over breadbasket regions. However, it is unknown whether climate models reproduce crop influences on regional climate, and thus the future risk of extreme climate events over global breadbasket regions. Here, we show that models overestimate growing season temperatures and underestimate evapotranspiration (ET) over global croplands, and that these differences increase with cropped area. We trace this warm and dry difference through each model's representation of the surface energy budget, showing that model differences in transpiration, leaf area index, and the ratio of transpiration to total ET drive the overall effect. While the implications of these model deficiencies for future projections are uncertain, they point to the importance of improving representations of crop-climate processes to better assess breadbasket vulnerability to climate change.

54 ENVIRONMENTAL SCIENCES↗

Examination of Radiation Belt Dynamics During Substorm Clusters: Activity Drivers and Dependencies of Trapped Flux Enhancements

Here, dynamical variations of radiation belt trapped electron fluxes are examined to better understand the variability of enhancements linked to substorm clusters. Analysis is undertaken using the Substorm Onsets and Phases from Indices of the Electrojet substorm cluster algorithm for event detection. Observations from low earth orbit are complemented by additional measurements from medium earth orbit to allow a major expansion in the energy range considered, from medium energy energetic electrons up to ultra-relativistic electrons. The number of substorms identified inside a cluster does not depend strongly on solar wind drivers or geomagnetic indices either before, during, or after the cluster start time. Clusters of substorms linked to moderate (100 nT < AE ≤ 300 nT) or strong AE (AE ≥ 300 nT) disturbances are associated with radiation belt flux enhancements, including up to ultra-relativistic energies by the strongest substorms (as measured by strong southward Bz and high AE). These clusters reliably occur during times of high speed solar winds streams with associated increased magnetospheric convection. However, substorm clusters associated with quiet AE disturbances (AE ≤ 100 nT) lead to no significant chorus whistler mode intensity enhancements, or increases in energetic, relativistic, or ultra-relativistic electron flux in the outer radiation belts. In these cases the solar wind speed is low, and the geomagnetic Kp index indicates a lack of magnetospheric convection. Our study clearly indicates that clusters of substorms occurring outside of high speed wind streams are not by themselves sufficient to drive acceleration, which may be due to the lack of pre-cluster convection.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep Learning Estimation of Daily Ground–Level NO 2 Concentrations from Remote Sensing Data

The limited number of nitrogen dioxide (NO 2 ) surface measurements calls for the development of highly accurate approaches to estimating surface NO 2 concentrations. In this study, we leverage a new satellite instrument, the TROPOspheric Monitoring Instrument (TROPOMI), along with other predictor variables, to estimate daily surface NO 2 concentrations over Texas in 2019. We use the deep convolutional neural network (Deep-CNN), an advanced deep learning algorithm, to obtain estimates and achieve a correlation coefficient (R) of 0.91, an index of agreement (IOA) of 0.95, and a mean absolute bias (MAB) of 1.75 ppb in surface NO 2 estimation. Additionally, we leverage a novel approach, SHapley Additive exPlanations (SHAP), to describe how Deep-CNN understands each predictor variable. The SHAP results show that the Deep-CNN model has an advanced understanding of the dataset, revealing that TROPOMI closely captures levels of NO 2 . In addition, we show the superiority of our Deep-CNN model at estimating surface NO 2 over other well-known machine learning and regression models in the field, including the support vector machines (SVM), random forest (RF), and multiple linear regression (MLR). Although SVM and RF show strong capabilities at estimating surface NO 2 concentrations, their accuracy is inferior to that of the Deep-CNN model, ranking second and third in model accuracy in this study. The MLR, however, shows a poor ability at NO 2 estimation and ranks last among all models. Furthermore, testing the impact of sample size on model performance, we also show that, compared to other models, Deep-CNN needs more samples to trigger its strength at surface NO 2 estimation.

54 ENVIRONMENTAL SCIENCES↗

Spatiotemporal Variations of Evapotranspiration in Amazonia Using the Wavelet Phase Difference Analysis

The relationships and seasonal-to-annual variations among evapotranspiration (ET), precipitation (P), terrestrial water storage anomalies (TWSA), radiation (downward shortwave radiation, DSR), and phenology (leaf area index, LAI) are complex across the Amazon basin. Here, to analyze how ET is controlled by these influencing factors, we used wavelet phase difference (WPD) to investigate the effects of P, TWSA, DSR, and LAI on ET at different spatiotemporal scales. The Amazon-scale averaged ET has strong correlations with these factors at the annual and multi-year periodicities. The patterns of WPDs have south-north and west-east divides due to the significant variation in climatic conditions. The results demonstrate that ET is mainly affected by water and energy availability while vegetation regulates both processes. The deep soil moisture/groundwater can provide strong subsidies to ET during the meteorological dry season in the water-limited area of Amazon. The WPD can well reflect the responses of ET to the variations of P, TWSA, DSR, and LAI, and the process of vegetation sustaining ET in the dry years in the water-limited area of the Amazon.

54 ENVIRONMENTAL SCIENCES↗

The Behaviors of Intraseasonal Cloud Organization During DYNAMO/AMIE

This study investigates the organization of tropical convection associated with the Madden-Julian Oscillation (MJO) during the Dynamics of the MJO/Atmospheric Radiation Measurement MJO Investigation Experiment field campaign. While it is known that tropical clouds can organize and impact the large-scale environment, how this occurs, and its underlying mechanism are not fully understood. Application of several existing cloud organization indices showed inconsistent evolutions in the measured degree of organization with the MJO. The inconsistency arises from the varying definitions and assumptions of cloud organization behaviors that are applied to each index. While these indices often combine different properties of clouds, such as their number density, size, and distance between them, the analysis of these properties separately provided further understanding of how clouds organize with the MJO. Using the rainfall clusters identified from the S-Polka radar, we find that deep convective rainfall clusters begin to increase their number density before the arrival of MJO enhanced convective center, which is accompanied by increased proximity (shorter distance to each other) and followed by growth in their size. However, the nonrandomness in the spatial distribution of rainfall clusters maximizes as MJO convection decays. Deep convective clusters become the least randomly distributed as the clusters decay because of the suppression and decay of isolated deep convective cells, while clustered deep convective cells exist longer. This evolution of cloud organization is analogous to mesoscale convective systems, indicating that the duration and frequency of their organization stages are altered by the large-scale environmental perturbations associated with the MJO.

54 ENVIRONMENTAL SCIENCES↗

Impact of Changing Winds on the Mauna Loa CO 2 Seasonal Cycle in Relation to the Pacific Decadal Oscillation

Abstract Long‐term measurements at the Mauna Loa Observatory (MLO) show that the CO 2 seasonal cycle amplitude (SCA) increased from 1959 to 2019 at an overall rate of 0.22 0.034 ppm decade −1 while also varying on interannual to decadal time scales. These SCA changes are a signature of changes in land ecological CO 2 fluxes as well as shifting winds. Simulations with the TM3 tracer transport model and CO 2 fluxes from the Jena CarboScope CO 2 Inversion suggest that shifting winds alone have contributed to a decrease in SCA of −0.10 0.022 ppm decade −1 from 1959 to 2019, partly offsetting the observed long‐term SCA increase associated with enhanced ecosystem net primary production. According to these simulations and MIROC‐ACTM simulations, the shorter‐term variability of MLO SCA is nearly equally driven by varying ecological CO 2 fluxes (49%) and varying winds (51%). We also show that the MLO SCA is strongly correlated with the Pacific Decadal Oscillation (PDO) due to varying winds, as well as with a closely related wind index (U‐PDO). Since 1980, 44% of the wind‐driven SCA decrease has been tied to a secular trend in the U‐PDO, which is associated with a progressive weakening of westerly winds at 700 mbar over the central Pacific from 20°N to 40°N. Similar impacts of varying winds on the SCA are seen in simulations at other low‐latitude Pacific stations, illustrating the difficulty of constraining trend and variability of land CO 2 fluxes using observations from low latitudes due to the complexity of circulation changes.

54 ENVIRONMENTAL SCIENCES↗

Cloud Morphology Evolution in Arctic Cold‐Air Outbreak: Two Cases During COMBLE Period

Abstract Cloud feedbacks play an important role in Arctic warming. Cloud morphology, for example, cloud size and spatial distributions, is among key factors that directly impact their radiative effects. In this work, we use two cases observed during the Cold‐air Outbreak (CAO) in the Marine Boundary Layer Experiment (COMBLE) to study the evolution of cloud size distributions as an air mass is advected from the Arctic over a comparatively warm ocean and cloud mesoscale organization changes from rolls to cells. Cloud objects are identified from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance images through an object segmentation procedure and roll breakup is identified by homogeneities in cloud water path (CWP). Roll breakup is found to be accompanied by a local minimum in wind shear and local maxima in cloud size and marine cold air outbreak index. The mean cloud horizontal aspect ratio has weak fetch dependency and is around 2 in roll, transition, and cell regimes. Regardless of distance from the ice edge, smaller clouds (<10 km 2 ) dominate the population number but not cloud cover. Cloud size distributions show bimodality in transition and cell regimes. For clouds with comparable sizes, mean nearest neighbor distances normalized by equivalent cloud radius converge to a single value for all regimes and for all but the smallest clouds, suggesting that clouds of comparable sizes in CAOs are separated by distance proportional to their sizes. The presented statistical results pave the way to evaluating model simulated cloud organizations during CAO events.

54 ENVIRONMENTAL SCIENCES↗

Strong Local Evaporative Cooling Over Land Due to Atmospheric Aerosols

Abstract Aerosols can enhance terrestrial productivity through increased absorption of solar radiation by the shaded portion of the plant canopy—the diffuse radiation fertilization effect. Although this process can, in principle, alter surface evaporation due to the coupling between plant water loss and carbon uptake, with the potential to change the surface temperature, aerosol‐climate interactions have been traditionally viewed in light of the radiative effects within the atmosphere. Here, we develop a modeling framework that combines global atmosphere and land model simulations with a conceptual diagnostic tool to investigate these interactions from a surface energy budget perspective. Aerosols increase the terrestrial evaporative fraction, or the portion of net incoming energy consumed by evaporation, by over 4% globally and as much as ∼40% regionally. The main mechanism for this is the increase in energy allocation from sensible to latent heat due to global dimming (reduction in global shortwave radiation) and slightly augmented by diffuse radiation fertilization. In regions with moderately dense vegetation (leaf area index >2), the local surface cooling response to aerosols is dominated by this evaporative pathway, not the reduction in incident radiation. Diffuse radiation fertilization alone has a stronger impact on gross primary productivity (+2.18 Pg C y −1 or +1.8%) than on land evaporation (+0.18 W m −2 or +0.48%) and surface temperature (−0.01 K). Our results suggest that it is important for land surface models to distinguish between quantity (change in total magnitude) and quality (change in diffuse fraction) of radiative forcing for properly simulating surface climate.

Chakraborty, TC↗

Representation of Leaf-to-Canopy Radiative Transfer Processes Improves Simulation of Far-Red Solar-Induced Chlorophyll Fluorescence in the Community Land Model Version 5

Recent advances in satellite observations of solar-induced chlorophyll fluorescence (SIF) provide a new opportunity to evaluate and constrain the simulation of terrestrial gross primary productivity (GPP). Accurate representation of the processes driving SIF emission and the radiative transfer of SIF to remote sensing sensors is an essential prerequisite for the evaluation and data assimilation. Recently, SIF simulations have been incorporated into several land surface models, but the scaling of SIF from leaf-level to canopy level is usually not well-represented. In this work, we incorporate the simulation of far-red SIF observed at nadir into the Community Land Model version 5 (CLM5). An efficient and accurate method based on escape probability is developed to scale SIF from leaf-level to top-of-canopy while taking clumping and the radiative transfer processes into account. SIF simulated by CLM5 and a canopy-level model agreed well at sites except one in needle leaf forest (R2>0.91, root-mean-square error < 0.19W m -2 sr -1 um -1 ), and captured the day-to-day variation of tower-measured SIF at temperate forest sites (R 2 >0.68). At the global scale, simulated SIF generally captured the spatial and seasonal pat37 terns of satellite-observed SIF (R2 > 0.76 except for tropical forest). Factors including the fluorescence emission model, clumping, bidirectional effect, and canopy properties (leaf optical properties and leaf area index) had considerable impacts on SIF simulation, and the discrepancies between simulated and observed SIF varied with plant functional type. By improving the representation of radiative transfer for SIF simulation, our model allows better comparisons between simulated and observed SIF towards constraining and evaluating GPP simulations.

59 BASIC BIOLOGICAL SCIENCES↗

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗