Social Contagion and General Diffusion Models of Adolescent Religious Transitions: A Tutorial, and EMOSA Applications.

Epidemic Models of the Onset of Social Activities (EMOSA) describe behaviors that spread through social networks. Two social influence methods are represented, social contagion (one-to-one spread) and general diffusion (spread through cultural channels). Past models explain problem behaviors-smoking, drinking, sexuality, and delinquency. We provide review, and a tutorial (including examples). Following, we present new EMOSA models explaining changes in adolescent and young adult religious participation. We fit the model to 10 years of data from the 1997 U.S. National Longitudinal Survey of Youth. Innovations include a three-stage bi-directional model, Bayesian Markov Chain Monte Carlo (MCMC) estimation, graphical innovations, and empirical validation. General diffusion dominated rapid reduction in church attendance during adolescence; both diffusion and social contagion explained church attendance stability in early adulthood.

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming.

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79-86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. These findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

Tree rings provide no evidence of a CO2 fertilization effect in old-growth subalpine forests of western Canada.

Atmospheric CO2 concentrations are now 1.7 times higher than the preindustrial values. Although photosynthetic rates are hypothesized to increase in response to rising atmospheric CO2 concentrations, results from in situ experiments are inconsistent in supporting a CO2 fertilization effect of tree growth. Tree-ring data provide a historical record of tree-level productivity that can be used to evaluate long-term responses of tree growth. We use tree-ring data from old-growth, subalpine forests of western Canada that have not had a stand-replacing disturbance for hundreds of years to determine if growth has increased over 19th and 20th centuries. Our sample consisted of 5,858 trees belonging to five species distributed over two sites in the coastal zone and two in the continental climate of the interior. We calculated annual increments in tree basal area, adjusted these increments for tree size and age, and tested whether there was a detectable temporal trend in tree growth over the 19th and 20th centuries. We found a similar pattern in 20th century growth trends among all species at all sites. Growth during the 19th century was mostly stable or increasing, with the exception of one of the coastal sites, where tree growth was slightly decreasing; whereas growth during the 20th century consistently decreased. The unexpected decrease in growth during the 20th century indicates that there was no CO2 fertilization effect on photosynthesis. We compared the growth trends from our four sites to the trends simulated by seven Earth System Models, and saw that most of the models did not predict these growth declines. Overall, our results indicate that these old-growth forests are unlikely to increase their carbon storage capacity in response to rising atmospheric CO2 , and thus are unlikely to contribute substantially to offsetting future carbon emissions.

Microbial models with data-driven parameters predict stronger soil carbon responses to climate change.

Long-term carbon (C) cycle feedbacks to climate depend on the future dynamics of soil organic carbon (SOC). Current models show low predictive accuracy at simulating contemporary SOC pools, which can be improved through parameter estimation. However, major uncertainty remains in global soil responses to climate change, particularly uncertainty in how the activity of soil microbial communities will respond. To date, the role of microbes in SOC dynamics has been implicitly described by decay rate constants in most conventional global carbon cycle models. Explicitly including microbial biomass dynamics into C cycle model formulations has shown potential to improve model predictive performance when assessed against global SOC databases. This study aimed to data-constrained parameters of two soil microbial models, evaluate the improvements in performance of those calibrated models in predicting contemporary carbon stocks, and compare the SOC responses to climate change and their uncertainties between microbial and conventional models. Microbial models with calibrated parameters explained 51% of variability in the observed total SOC, whereas a calibrated conventional model explained 41%. The microbial models, when forced with climate and soil carbon input predictions from the 5th Coupled Model Intercomparison Project (CMIP5), produced stronger soil C responses to 95 years of climate change than any of the 11 CMIP5 models. The calibrated microbial models predicted between 8% (2-pool model) and 11% (4-pool model) soil C losses compared with CMIP5 model projections which ranged from a 7% loss to a 22.6% gain. Lastly, we observed unrealistic oscillatory SOC dynamics in the 2-pool microbial model. The 4-pool model also produced oscillations, but they were less prominent and could be avoided, depending on the parameter values.

Peat profile database from peatlands in Canada.

Peatlands cover approximately 12% of the Canadian landscape and play an important role in the carbon cycle through their centennial- to millennial-scale storage of carbon under waterlogged and anoxic conditions. In recognizing the potential of these ecosystems as natural climate solutions and therefore the need to include them in national greenhouse gas inventories, the Canadian Model for Peatlands module (CaMP v. 2.0) was developed by the Canadian Forest Service. Model parameterization included compiling peat profiles across Canada to calibrate peat decomposition rates from different peatland types, to define typical bulk density profiles, and to describe the hydrological (i.e., water table) response of peatlands to climatic changes. A total of 1217 sites were included in the dataset from published and unpublished sources. The CORESITES table contains site location and summary data for each profile, as well as an estimate of total carbon mass per unit area (in megagrams of C per hectare). Total carbon mass per unit area at each location was calculated using bulk density and carbon content through each profile. The PROFILES table contains data for depth (in centimeters), bulk density (in grams per cubic meter), ash and carbon content (in percentage), and material descriptions for contiguous samples through each peat profile. Data gaps for bulk density and C content were filled using interpolation, regression trees, and assigned values based on material description and/or soil classification to allow for the estimation of total carbon mass per unit area. A subset of the sites (N = 374) also have pH and pore water trace-elemental geochemistry data and are found in the WATER table. The REFERENCES table contains the full citation of each source of the data and is linked to each core location through the SOURCEDATA table. The LOOKUP table defines codes in the database that required more space that what was sufficient in the metadata tables. The data can be accessed on Open Government Canada and will be useful for future work on carbon stock mapping and ecosystem modeling. All metadata and data are provided © Her Majesty the Queen in Right of Canada, 2023 and information contained in this publication may be reproduced for personal or public noncommercial purposes with attribution, whereas commercial reproduction and distribution are prohibited except with written permission from NRCan; complete details are noted in the Supporting Information file Metadata S1 (see Class III.B.3: Copyright restrictions).

Traceable components of terrestrial carbon storage capacity in biogeochemical models.

Biogeochemical models have been developed to account for more and more processes, making their complex structures difficult to be understood and evaluated. Here, we introduce a framework to decompose a complex land model into traceable components based on mutually independent properties of modeled biogeochemical processes. The framework traces modeled ecosystem carbon storage capacity (Xss ) to (i) a product of net primary productivity (NPP) and ecosystem residence time (τE ). The latter τE can be further traced to (ii) baseline carbon residence times (τ'E ), which are usually preset in a model according to vegetation characteristics and soil types, (iii) environmental scalars (ξ), including temperature and water scalars, and (iv) environmental forcings. We applied the framework to the Australian Community Atmosphere Biosphere Land Exchange (CABLE) model to help understand differences in modeled carbon processes among biomes and as influenced by nitrogen processes. With the climate forcings of 1990, modeled evergreen broadleaf forest had the highest NPP among the nine biomes and moderate residence times, leading to a relatively high carbon storage capacity (31.5 kg cm(-2) ). Deciduous needle leaf forest had the longest residence time (163.3 years) and low NPP, leading to moderate carbon storage (18.3 kg cm(-2) ). The longest τE in deciduous needle leaf forest was ascribed to its longest τ'E (43.6 years) and small ξ (0.14 on litter/soil carbon decay rates). Incorporation of nitrogen processes into the CABLE model decreased Xss in all biomes via reduced NPP (e.g., -12.1% in shrub land) or decreased τE or both. The decreases in τE resulted from nitrogen-induced changes in τ'E (e.g., -26.7% in C3 grassland) through carbon allocation among plant pools and transfers from plant to litter and soil pools. Our framework can be used to facilitate data model comparisons and model intercomparisons via tracking a few traceable components for all terrestrial carbon cycle models. Nevertheless, more research is needed to develop tools to decompose NPP and transient dynamics of the modeled carbon cycle into traceable components for structural analysis of land models.