Search

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

Leaf-level gas exchange data support the mechanistic understanding of plant fluxes of carbon and water. These fluxes inform our understanding of ecosystem function, are an important constraint on parameterization of terrestrial biosphere models, are necessary to understand the response of plants to global environmental change, and are integral to efforts to improve crop production. Collection of these data using gas analyzers can be both technically challenging and time consuming, and individual studies generally focus on a small range of species, restricted time periods, or limited geographic regions. The high value of these data is exemplified by the many publications that reuse and synthesize gas exchange data, however the lack of metadata and data reporting conventions make full and efficient use of these data difficult. Here we propose a reporting format for leaf-level gas exchange data and metadata to provide guidance to data contributors on how to store data in repositories to maximize their discoverability, facilitate their efficient reuse, and add value to individual datasets. For data users, the reporting format will better allow data repositories to optimize data search and extraction, and more readily integrate similar data into harmonized synthesis products. The reporting format specifies data table variable naming and unit conventions, as well as metadata characterizing experimental conditions and protocols. For common data types that were the focus of this initial version of the reporting format, i.e., survey measurements, dark respiration, carbon dioxide and light response curves, and parameters derived from those measurements, we took a further step of defining required additional data and metadata that would maximize the potential reuse of those data types. To aid data contributors and the development of data ingest tools by data repositories we provided a translation table comparing the outputs of common gas exchange instruments. Extensive consult

Leaf respiration in the dark (Rdark) and light (Rday) is poorly characterized in diverse tropical ecosystems, and little to no information exists on the degree of light suppression in common tree species within the Amazon basin, and their dependences upon plant functional traits and position within the canopy. We quantified Rdark and apparent Rday using the Kok method and measured key leaf traits in 26 tree individuals of different species distributed in three different canopy positions: canopy, lower canopy, and understory. To explore the relationships between the leaf traits we used the standardized major axis (SMA). We found that canopy trees had significantly higher rates of Rdark and Rday than trees in the understory. The difference between Rdark and Rday (the light suppression of respiration) was greatest in the understory (68 ± 9%, 95% CI) and lower canopy (49 ± 9%, 95% CI) when compared to the canopy (37 ± 10%, 95% CI). We also found that Rday was significantly and strongly correlated with Rdark (p < 0.001) for all the canopy positions. Also, leaf mass per area (LMA) and leaf Phosphorus concentration (P) had a significant relationship with Rdark (p < 0.001; p = 0.003), respectively. In addition, a significant relationship was found for LMA in the canopy and lower canopy positions (p = 0.009; p = 0.048) while P was only significant in the canopy (p = 0.044). Finally, no significant relationship was found between Rdark and nitrogen, sugars, and starch. Our results highlight the importance of including representation of the light suppression of leaf respiration in terrestrial biosphere models and also of accounting for vertical gradients within forest canopies and connections with functional traits.

Leaf respiration in the dark (Rdark) and light (Rday) is poorly characterized in diverse tropical ecosystems, and little to no information exists on the degree of light suppression in common tree species within the Amazon basin, and their dependences upon plant functional traits and position within the canopy. We quantified Rdark and apparent Rday using the Kok method and measured key leaf traits in 26 tree individuals of different species distributed in three different canopy positions: canopy, lower canopy, and understory. To explore the relationships between the leaf traits we used the standardized major axis (SMA). We found that canopy trees had significantly higher rates of Rdark and Rday than trees in the understory. The difference between Rdark and Rday (the light suppression of respiration) was greatest in the understory (68 ± 9%, 95% CI) and lower canopy (49 ± 9%, 95% CI) when compared to the canopy (37 ± 10%, 95% CI). We also found that Rday was significantly and strongly correlated with Rdark (p < 0.001) for all the canopy positions. Also, leaf mass per area (LMA) and leaf Phosphorus concentration (P) had a significant relationship with Rdark (p < 0.001; p = 0.003), respectively. In addition, a significant relationship was found for LMA in the canopy and lower canopy positions (p = 0.009; p = 0.048) while P was only significant in the canopy (p = 0.044). Finally, no significant relationship was found between Rdark and nitrogen, sugars, and starch. Our results highlight the importance of including representation of the light suppression of leaf respiration in terrestrial biosphere models and also of accounting for vertical gradients within forest canopies and connections with functional traits.

Leaf-level gas exchange data support the mechanistic understanding of plant fluxes of carbon and water. These fluxes inform our understanding of ecosystem function, are an important constraint on parameterization of terrestrial biosphere models, are necessary to understand the response of plants to global environmental change, and are integral to efforts to improve crop production. Collection of these data using gas analyzers can be both technically challenging and time consuming, and individual studies generally focus on a small range of species, restricted time periods, or limited geographic regions. The high value of these data is exemplified by the many publications that reuse and synthesize gas exchange data, however the lack of metadata and data reporting conventions make full and efficient use of these data difficult. Here we propose a reporting format for leaf-level gas exchange data and metadata to provide guidance to data contributors on how to store data in repositories to maximize their discoverability, facilitate their efficient reuse, and add value to individual datasets. For data users, the reporting format will better allow data repositories to optimize data search and extraction, and more readily integrate similar data into harmonized synthesis products. The reporting format specifies data table variable naming and unit conventions, as well as metadata characterizing experimental conditions and protocols. For common data types that were the focus of this initial version of the reporting format, i.e., survey measurements, dark respiration, carbon dioxide and light response curves, and parameters derived from those measurements, we took a further step of defining required additional data and metadata that would maximize the potential reuse of those data types. To aid data contributors and the development of data ingest tools by data repositories we provided a translation table comparing the outputs of common gas exchange instruments. Extensive consult

Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While singlePFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early-and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states a

Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While singlePFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early-and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states a

Vegetation plays an important role in regulating global carbon cycles and is a key component of the Earth system models (ESMs) that aim to project Earth's future climate. In the last decade, the vegetation component within ESMs has witnessed great progress from simple "big-leaf" approaches to demographically structured approaches, which have a better representation of plant size, canopy structure, and disturbances. These demographically structured vegetation models typically have a large number of input parameters, and sensitivity analysis is needed to quantify the impact of each parameter on the model outputs for a better understanding of model behavior. In this study, we conducted a comprehensive sensitivity analysis to diagnose the Community Land Model coupled to the Functionally Assembled Terrestrial Simulator, or CLM4.5(FATES). Specifically, we quantified the first- and second-order sensitivities of the model parameters to outputs that represent simulated growth and mortality as well as carbon fluxes and stocks for a tropical site with an extent of 1×1°. While the photosynthetic capacity parameter (Vc;max25) is found to be important for simulated carbon stocks and fluxes, we also show the importance of carbon storage and allometry parameters, which determine survival and growth strategies within the model. The parameter sensitivity changes with different sizes of trees and climate conditions. The results of this study highlight the importance of understanding the dynamics of the next generation of demographically enabled vegetation models within ESMs to improve model parameterization and structure for better model fidelity.

Global terrestrial models currently predict that the Amazon rainforest will continue to act as a carbon sink in the future, primarily owing to the rising atmospheric carbon dioxide (CO2) concentration. Soil phosphorus impoverishment in parts of the Amazon basin largely controls its functioning, but the role of phosphorus availability has not been considered in global model ensembles—for example, during the Fifth Climate Model Intercomparison Project. Here we simulate the planned free-air CO2 enrichment experiment AmazonFACE with an ensemble of 14 terrestrial ecosystem models. We show that phosphorus availability reduces the projected CO2-induced biomass carbon growth by about 50% to 79 ± 63 g C m−2 yr−1 over 15 years compared to estimates from carbon and carbon–nitrogen models. Our results suggest that the resilience of the region to climate change may be much less than previously assumed. Variation in the biomass carbon response among the phosphorus-enabled models is considerable, ranging from 5 to 140 g C m−2 yr−1, owing to the contrasting plant phosphorus use and acquisition strategies considered among the models. The Amazon forest response thus depends on the interactions and relative contributions of the phosphorus acquisition and use strategies across individuals, and to what extent these processes can be upregulated under elevated CO2.

Vegetation plays an important role in regulating global carbon cycles and is a key component of the Earth system models (ESMs) that aim to project Earth's future climate. In the last decade, the vegetation component within ESMs has witnessed great progress from simple "big-leaf" approaches to demographically structured approaches, which have a better representation of plant size, canopy structure, and disturbances. These demographically structured vegetation models typically have a large number of input parameters, and sensitivity analysis is needed to quantify the impact of each parameter on the model outputs for a better understanding of model behavior. In this study, we conducted a comprehensive sensitivity analysis to diagnose the Community Land Model coupled to the Functionally Assembled Terrestrial Simulator, or CLM4.5(FATES). Specifically, we quantified the first- and second-order sensitivities of the model parameters to outputs that represent simulated growth and mortality as well as carbon fluxes and stocks for a tropical site with an extent of 1×1°. While the photosynthetic capacity parameter (Vc;max25) is found to be important for simulated carbon stocks and fluxes, we also show the importance of carbon storage and allometry parameters, which determine survival and growth strategies within the model. The parameter sensitivity changes with different sizes of trees and climate conditions. The results of this study highlight the importance of understanding the dynamics of the next generation of demographically enabled vegetation models within ESMs to improve model parameterization and structure for better model fidelity.

Global terrestrial models currently predict that the Amazon rainforest will continue to act as a carbon sink in the future, primarily owing to the rising atmospheric carbon dioxide (CO2) concentration. Soil phosphorus impoverishment in parts of the Amazon basin largely controls its functioning, but the role of phosphorus availability has not been considered in global model ensembles—for example, during the Fifth Climate Model Intercomparison Project. Here we simulate the planned free-air CO2 enrichment experiment AmazonFACE with an ensemble of 14 terrestrial ecosystem models. We show that phosphorus availability reduces the projected CO2-induced biomass carbon growth by about 50% to 79 ± 63 g C m−2 yr−1 over 15 years compared to estimates from carbon and carbon–nitrogen models. Our results suggest that the resilience of the region to climate change may be much less than previously assumed. Variation in the biomass carbon response among the phosphorus-enabled models is considerable, ranging from 5 to 140 g C m−2 yr−1, owing to the contrasting plant phosphorus use and acquisition strategies considered among the models. The Amazon forest response thus depends on the interactions and relative contributions of the phosphorus acquisition and use strategies across individuals, and to what extent these processes can be upregulated under elevated CO2.

We compared initial screening transferrin saturation (TfSat) and serum ferritin (SF) phenotypes and HFE C282Y and H63D genotypes of 645 Native American and 43,453 white Hemochromatosis and Iron Overload Screening Study participants who did not report a previous diagnosis of hemochromatosis or iron overload. Elevated measurements were defined as TfSat >50% in men and >45% in women and SF >300 ng/ml in men and >200 ng/ml in women. Mean TfSat was 31% in Native American men and 32% in white men (p = 0.0337) and 25% in Native American women and 27% in white women (p < 0.0001). Mean SF was 153 microg/l in Native American and 151 microg/l in white men (p = 0.8256); mean SF was 55 microg/l in Native American women and 63 microg/l in white women (p = 0.0015). The C282Y allele frequency was 0.0340 in Native Americans and 0.0683 in whites (p < 0.0001). The H63D allele frequency was 0.1150 in Native Americans and 0.1532 in whites (p = 0.0001). We conclude that the screening TfSat and SF phenotypes of Native Americans are similar to those of whites. The allele frequencies of HFE C282Y and H63D are significantly lower in Native Americans than in whites.

We compared initial screening transferrin saturation (TfSat) and serum ferritin (SF) phenotypes and HFE C282Y and H63D genotypes of 645 Native American and 43,453 white Hemochromatosis and Iron Overload Screening Study participants who did not report a previous diagnosis of hemochromatosis or iron overload. Elevated measurements were defined as TfSat >50% in men and >45% in women and SF >300 ng/ml in men and >200 ng/ml in women. Mean TfSat was 31% in Native American men and 32% in white men (p = 0.0337) and 25% in Native American women and 27% in white women (p < 0.0001). Mean SF was 153 microg/l in Native American and 151 microg/l in white men (p = 0.8256); mean SF was 55 microg/l in Native American women and 63 microg/l in white women (p = 0.0015). The C282Y allele frequency was 0.0340 in Native Americans and 0.0683 in whites (p < 0.0001). The H63D allele frequency was 0.1150 in Native Americans and 0.1532 in whites (p = 0.0001). We conclude that the screening TfSat and SF phenotypes of Native Americans are similar to those of whites. The allele frequencies of HFE C282Y and H63D are significantly lower in Native Americans than in whites.

A human gene with strong homology to the MAGE gene family located in XqZ7-qter has been isolated by using exon-trapping of cosmids in the Xp21.3 region. We have mapped and sequenced cDNA and genomic clones corresponding to this gene, MAGE-Xp, and shown that the last exon contains the open reading frame and is present in a minimum of five copies in a 30-kb interval, MAGE-Xp, is expressed only in testis and, unlike the Xq27-qter MAGE genes, it is not expressed in any of 12 different tumor tissues tested. However, the gene and predicted protein structure are conserved, suggesting a similar function, MAGE-Xp is located in the 160-kb critical interval defined for the locus involved in sex determination within Xp21 and is 50 kb distal to the DAX-1 gene, which is responsible for X-chromosome-linked adrenal hypoplasia congenita.