Search NASA⌕ Search

Engineering topics

Leonardo Calle

Publications and source records attributed to Leonardo Calle.

6 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗
NASA NTRS2019

State of the science in reconciling top‐down and bottom‐up approaches for terrestrial CO2 budget

Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: “top‐down” atmospheric inversions and “bottom‐up” biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere–biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink–source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top‐down and bottom‐up approaches and yield more robust knowledge of regional CO2 budgets.

atmospheric inversion↗
NASA NTRS2019

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: an Application to Xco2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

segmentation algorithm↗
NASA NTRS2019

Role of forest regrowth in global carbon sink dynamics

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981–2010. For 2001–2010 we find a carbon sink of 0.85 (0.66–0.96) Pg/year located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03–1.96) Pg/year located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44–131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981–2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

forest↗
NASA NTRSDate not supplied

Emissions Characterization and Smoke Transport of A Prescribed Fire

Under the background of climate change, vast regions of the world will face a hotter and probably drier future that favors ignition and spread of wildfires. This poses additional challenges to the communities who are fighting the ever-larger wildfires. Fire and smoke models are valuable tools to fire managers and decision-makers to mitigate the impact of fires. Existing fire modeling systems bear large uncertainties due to the difficulty in collecting data to characterize fire behavior and smoke transport. More observations are urgently needed for improvement of fire modeling and reduction in model uncertainty. The coordinated prescribed burn experiment, such as the Fire and Smoke Model and Measurement Evaluation Experiment (FSMMEE), will allow the concurrent collection of critical measurements of fuel, fire behavior, smoke, and meteorology to better understand and model fires. In this preliminary investigation, we employed the data from the prescribed burn experiment at Langdon Mountain, Utah, in November 2019, to reconstruct the fire emissions and characterize the smoke transport. The NASA Unified Weather Research and Forecasting Model (NU-WRF) was utilized to assist in identifying the potential pre-fire remote sensing capabilities, such as fuel load, moisture, and fire radiative power, which are helpful to model development and advancement. Along the way, two sets of NU-WRF experiments would be done by applying either a default fire emissions inventory (i.e., NASA’s Quick Fire Emissions Dataset, QFED) or reconstructed fire emissions using the data collected from the prescribed fire. The results would help answer the questions such as “How do reconstructed emissions compare to QFED emissions and their impact on plume transport”.

NU-WRF↗