Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temporal analysis of products”

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 19 records

Global Validation of Columnar Water Vapor Derived from EOS MODIS-MAIAC Algorithm Against the Ground-Based AERONET Observations

The water vapor is a relevant greenhouse gas in the Earth's climate system, and satellite products become one of the most effective way to characterize and monitor the columnar water vapor (CWV) content at global scale. Recently, a new product (MCD19) was released as part of MODIS (Moderate Resolution Imaging Spectroradiometer) Collection 6 (C6). This operational product from the Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm includes a high 1-kilometer resolution CWV retrievals. This study presents the first global validation of MAIAC C6 CWV obtained from MODIS MCD19A2 product. This evaluation was performed using Aerosol Robotic Network (AERONET) observations at 265 sites (2000-2017). Overall, the results show a good agreement between MAIAC/AERONET CWV retrievals, with correlation coefficient higher than 0.95 and RMS (Root Mean Square) error lower than 0.250 centimeters. The binned error analysis revealed an underestimation (approximately 10 percent) of Aqua CWV retrievals with negative bias for CWV higher than 3.0 centimeters. In contrast, Terra CWV retrievals show a slope of regression close to unity and a low mean bias of 0.075 centimeters. While the accuracy is relatively similar between 1.0 and 5.0 centimeters for both sensor products, Terra dataset is more reliable for applications in humid tropical areas (less than 5.0 centimeters). The expected error was defined as plus or minus 15 percent, with less than 68 percent of retrievals falling within this envelope. However, the accuracy is regionally dependent, and lower error should be expected in some regions, such as South America and Oceania. Since MODIS instruments have exceeded their design lifetime, time series analysis was also presented for both sensor products. The temporal analysis revealed a systematic offset of global average between Terra and Aqua CWV records. We also found an upward trend (approximately 0.2 centimeters per decade) in Terra CWV retrievals, while Aqua CWV retrievals remain stable over time. The sensor degradation influences the ability to detect climate signals, and this study indicates the need for revisiting calibration of the MODIS bands 17-19, mainly for Terra instrument, to assure the quality of the MODIS water vapor product. Finally, this study presents a comprehensive validation analysis of MAIAC CWV over land, raising the understanding of its overall quality.

MAIAC↗

Studies of the net surface radiative flux from satellite radiances during FIFE

Studies of the net surface radiative flux from satellite radiances during First ISLSCP Field Experiment (FIFE) are presented. Topics covered include: radiative transfer model validation; calibration of VISSR and AVHRR solar channels; development and refinement of algorithms to estimate downward solar and terrestrial irradiances at the surface, including photosynthetically available radiation (PAR) and surface albedo; verification of these algorithms using in situ measurements; production of maps of shortwave irradiance, surface albedo, and related products; analysis of the temporal variability of shortwave irradiance over the FIFE site; development of a spectroscopy technique to estimate atmospheric total water vapor amount; and study of optimum linear combinations of visible and near-infrared reflectances for estimating the fraction of PAR absorbed by plants.

Frouin, Robert↗

A Spatio-Temporal Approach for Global Validation and Analysis of MODIS Aerosol Products

With the launch of the MODIS sensor on the Terra spacecraft, new data sets of the global distribution and properties of aerosol are being retrieved, and need to be validated and analyzed. A system has been put in place to generate spatial statistics (mean, standard deviation, direction and rate of spatial variation, and spatial correlation coefficient) of the MODIS aerosol parameters over more than 100 validation sites spread around the globe. Corresponding statistics are also computed from temporal subsets of AERONET-derived aerosol data. The means and standard deviations of identical parameters from MOMS and AERONET are compared. Although, their means compare favorably, their standard deviations reveal some influence of surface effects on the MODIS aerosol retrievals over land, especially at low aerosol loading. The direction and rate of spatial variation from MODIS are used to study the spatial distribution of aerosols at various locations either individually or comparatively. This paper introduces the methodology for generating and analyzing the data sets used by the two MODIS aerosol validation papers in this issue.

Ichoku, Charles↗

Making Better Use of Satellite Data: The Satellite Needs Working Group

The U.S. Group on Earth Observations (USGEO) initiated in 2016 the Satellite Needs Working Group (SNWG) to identify and communicate the Earth observation needs of U.S. federal agencies. The SNWG identifies such needs through a biennial survey followed by interviews and follow-up discussions by the satellite Earth data providers of the U.S. Government: the National Aeronautics and Space Administration (NASA), the National Oceanic and Atmospheric Administration (NOAA), and the U.S. Geological Survey (USGS). Solutions and services are identified that leverage current or upcoming satellite missions to meet the identified needs; implementation of services that are estimated to significantly increase the level of satisfaction of multiple U.S. agencies are funded by NASA. The SNWG process has resulted in the implementation of numerous services that have impacted operations of not only U.S. agencies but of academic and international institutions as well. Notable examples include the Harmonized Landsat Sentinel-2 project that leverages European Space Agency (ESA) and NASA satellite assets to generate a global, analysis-ready, surface reflectance product with a temporal resolution of two days; the Airborne Data Management Group that curates and provides access to relevant resources, information and data from existing and past NASA field campaigns; and the Dynamic Surface Water Extent product which consists of harmonized but independent water extents derived from both from optical and radar data. These and other products are hosted at the NASA Distributed Active Archive Centers (DAACs) for free and open access. The SNWG Management Office at NASA’s Interagency Implementation and Advanced Concepts Team (IMPACT) manages the implementation of selected solutions and, importantly, NASA’s response to the needs of federal agencies through a Stakeholder Engagement Program. The impact of implemented services and solutions is hampered without efforts to build capacity around the use of the services. The Stakeholder Engagement Program ensures relevant training and outreach to SNWG agencies in collaboration with each solution implementation team.

remote sensing↗

Land Surface Temperature Measurements from EOS MODIS Data

This report summarizes the accomplishments made by the MODIS LST (Land-Surface Temperature) group at University of California, Santa Barbara, under NASA Contract. Version 1 of the MODIS Land-Surface Temperature Algorithm Theoretical Basis Document (ATBD) was reviewed in June 1994, version 2 reviewed in November 1994, version 3.1 in August 1996, and version 3.3 updated in April 1999. Based on the ATBD, two LST algorithms were developed, one is the generalized split-window algorithm and another is the physics-based day/night LST algorithm. These two LST algorithms were implemented into the production generation executive code (PGE 16) for the daily standard MODIS LST products at level-2 (MODII-L2) and level-3 (MODIIA1 at 1 km resolution and MODIIB1 at 5km resolution). PGE codes for 8-day 1 km LST product (MODIIA2) and the daily, 8-day and monthly LST products at 0.05 degree latitude/longitude climate model grids (CMG) were also delivered. Four to six field campaigns were conducted each year since 2000 to validate the daily LST products generated by PGE16 and the calibration accuracies of the MODIS TIR bands used for the LST/emissivity retrieval from versions 2-4 of Terra MODIS data and versions 3-4 of Aqua MODIS data. Validation results from temperature-based and radiance-based methods indicate that the MODIS LST accuracy is better than 1 C in most clear-sky cases in the range from -10 to 58 C. One of the major lessons learn from multi- year temporal analysis of the consistent V4 daily Terra MODIS LST products in 2000-2003 over some selected target areas including lakes, snow/ice fields, and semi-arid sites is that there are variable numbers of cloud-contaminated LSTs in the MODIS LST products depending on surface elevation, land cover types, and atmospheric conditions. A cloud-screen scheme with constraints on spatial and temporal variations in LSTs was developed to remove cloud-contaminated LSTs. The 5km LST product was indirectly validated through comparisons to the 1 km LST product. Twenty three papers related to the LST research work were published in journals over the last decade.

Wan, Zheng-Ming↗

Weather analysis and interpretation procedures developed for the US/Canada wheat and barley exploratory experiment

Procedures and techniques for providing analyses of meteorological conditions at segments during the growing season were developed for the U.S./Canada Wheat and Barley Exploratory Experiment. The main product and analysis tool is the segment-level climagraph which depicts temporally meteorological variables for the current year compared with climatological normals. The variable values for the segment are estimates derived through objective analysis of values obtained at first-order station in the region. The procedures and products documented represent a baseline for future Foreign Commodity Production Forecasting experiments.

Trenchard, M. H.↗

Northern Brazil Agriculture: Measuring Soybean Yields in Northern Brazil During El Niño-Southern Oscillation Conditions to Evaluate Trends in Agricultural Production and Support Crop Forecasting, 1984 – 2023

As one of the largest agricultural exporters in the world, Brazil’s crop productivity is highly influential to the world’s food supply. Crop productivity across Brazil is heavily influenced by climatic factors, such as the El Niño-Southern Oscillation (ENSO), which drives spatially heterogenous effects on growing conditions. To understand the relationship between ENSO conditions and staple crop productivity, the team used 1984 – 2023 climatic data and vegetation indices for four Brazilian states: Bahia, Mato Grosso, Pará, and Tocantins. The team focused on soybean growing areas in each state derived from annual land use/land cover classifications. Monthly mean Normalized Difference Vegetation Index (NDVI) per state were calculated using multispectral imagery from the Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), and Landsat 8 Operational Land Imagery (OLI) sensors. In addition, the team used crop production indices, ENSO anomalies, temperature, and precipitation data in this study. The team found that although monthly ENSO anomaly did have a positive relationship between temperature and a slight negative relationship with precipitation, ENSO and monthly soy NDVI were not significantly correlated. Results find no association between ENSO conditions and NDVI in the study area at the state level spatial resolution during the study period. Furthermore, the team found that there was no correlation between cumulative growing season soy NDVI and detrended soy production yield for this study region at the spatial and temporal scale of the analysis. Findings from this project will be used by the United States Department of Agriculture’s Foreign Agriculture Service to inform crop productivity forecasting.

remote sensing↗

SatCORPS Global Cloud Composite (GCC): the Design and Delivery of A High Quality, High Resolution, Global Cloud Product Available in Near-Real Time

The NASA Satellite ClOud and Radiation Property retrieval System (SatCORPS) supports the development of an analysis ready and cloud-optimized data transformation pipeline and geospatial service enablement of a global cloud composite (GCC) product derived from global geostationary satellite imagery. This geospatial service will be available at high temporal and spatial resolution via the SatCORPS web mapping application for visualization and analysis as well as direct ingestion to common geospatial software and custom programming. The resulting global cloud composite products from the processing pipeline can then be geospatially-service enabled as ArcGIS Image Services and Open Geospatial Consortium (OGC) Web Mapping/Coverage Services for visualization and analysis via a web mapping application and common geospatial software. Near real time global observations are created through the composition of five geostationary satellites that provides modelling and forecasting communities with the capability to provide high quality and timely information to start the projection process. The Global Cloud Composite product combines information from geostationary satellites, GOES-16, GOES-17, Himawari-8, Meteosat-11 and Meteosat-9 to create a single global composite netcdf file and images using the different products within the netcdf file. The SatCORPS team, though our Global Cloud Composite (GCC) product and web-based visualization tools including Geographical Information System (GIS) services provide near real time global cloud product information to both automated processes and traditional web users that is timely and high quality derived from geostationary satellites. The Global Cloud Composite product takes advantage of the scalable processing resources provided by the AWS batch service to provide new composites every thirty minutes. Because information from each of the low earth orbiting satellites is available on schedules tuned to the specific satellite, the processing algorithm temporally composites the final dataset as each satellite’s information becomes available. The SatCORPS team has leveraged our experience using Amazon Web Services (AWS) to build a low latency high availability tool that allows end users both human and automated to acquire high quality and high-resolution Geostationary Earth Orbiting (GEO) information at zero cost to the end user. This presentation will describe how we architected and implemented the service as well as lessons learned based on our experiences both developing and operating the system. The lessons learned include how we integrated multiple services including Amazon Batch, Amazon S3 and Amazon Lambda service to create a low cost but high-performance processing system that is capable of identifying and processing the most appropriate satellite overpass information into global cloud composites. We will also describe our web-based tools including our Geographic Information System that can be used for visualization and analysis. The products from the processing can be geospatially-service enabled as ArcGIS Image Services and Open Geospatial Consortium (OGC) Web Mapping/Coverage Services for visualization and analysis via a web mapping application and common geospatial software. The SatCORPS Global Composite Cloud product provides sophisticated global composited cloud research products with very low latency that we see that as filling a rapidly growing need in the research and modelling community with no up-front nor ongoing costs associated with downloading or using the information.

AWS AMCE SMCE GCC SATCORPS GLOBAL CLOUD COMPOSITE ↗

Comparisons of Temporal and Spatial Trends in the Spatially Complete Global Spectral Surface Albedos Products

Five years of spatially complete snow-free land surface albedo data have been prepared using high quality white-sky and black-sky land surface albedo observations (MOD43B3) from the MODIS instrument aboard NASA's Terra satellite platform. The data were generated using an updated ecosystem-dependent temporal interpolation technique. In this paper, we describe the refinements in the technique and the creation of a spatially complete snow-free five-year aggregate climatology product. This paper also describes an error analysis of the interpolation technique. The filled albedo products are examined through comparisons of temporal and spatial trends for pixels that have been filled versus pixels that have been retained in the original MOD43B3 values. The variability in the trends showcase how the filling technique maintains the pixel-level spatial, spectral, and temporal integrity of the MOD43B3 data. These comparisons are made for both a single-year of filled data, year 2002, and for the five-year aggregate climatology product.

Moody, Eric G.↗

Developing A Continuous Ozone Record Through the SAGE and Aura Missions With NASA Reanalysis Products

During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the loss and the subsequent recovery of the stratospheric ozone layer. Due to the employed solar occultation and self-calibration method, the SAGE monitors have produced stable data throughout the lifetime of each instrument. However, over ten years passed between the end of the SAGE II and SAGE III/M3M missions in 2005 and the launch of SAGE III/ISS instrument in 2017, leaving a gap in the data that much be bridged in order to assess the trends in the ozone record. Reanalysis products, such as the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), are attractive candidates for trend analysis due to the statistically optimized combination of multiple observing systems and the regular temporal and spatial coverage. In this study, we explore using the SAGE records to develop a stable reanalysis data product, suitable for trend analysis, from the start of the SAGE II record in 1984 through the present. Changes in the assimilated observation systems can introduce discontinuities within the MERRA-2 ozone record, such as in 2004 when the MERRA-2 system shifted from assimilating ozone retrievals collected by SBUV instruments to those collected by instruments onboard the Aura satellite. We follow the radiative transfer procedure outlined by Wargan et al. (2018) to address discontinuities in the MERRA-2 ozone dataset at the 2004 transition and during the Aura record. SAGE II ozone profiles are used to address discontinuities in upper stratospheric ozone associated with changes in the MERRA-2 meteorological observing system in 1998 and 1995. Lastly, we will use the resulting bias-corrected MERRA-2 ozone fields to assess the relative performance of the data from different SAGE sensors.

SAGE↗

Barometric Tides from ECMWF Operational Analyses

The solar diurnal and semidiurnal tidal oscillations in surface pressure are extracted from the the operational analysis product of the European Centre for Medium Range Weather Forecasting (ECMWF). For the semidiurnal tide this involves a special temporal interpolation, following Van den Dool and colleagues. The resulting tides are compared with a ground truth tide dataset, a compilation of well-determined tide estimates deduced from long time series of station barometer measurements. These comparisons show that the ECMWF tides are significantly more accurate than the tides deduced from two other widely available reanalysis products. Spectral analysis of ECMWF pressure series shows that the tides consist of sharp central peaks with modulating sidelines at integer multiples of 1 cycle/year, superimposed on a broad cusp of stochastic energy. The integrated energy in the cusp dominates that of the sidelines. This complicates development of a simple model that can characterize the full temporal variability of the tides.

Ray, R. D.↗

Calibrated Noise Measurements with Induced Receiver Gain Fluctuations

The lack of well-developed techniques for modeling changing statistical moments in our observations has stymied the application of stochastic process theory in science and engineering. These limitations were encountered when modeling the performance of radiometer calibration architectures and algorithms in the presence of non stationary receiver fluctuations. Analyses of measured signals have traditionally been limited to a single measurement series. Whereas in a radiometer that samples a set of noise references, the data collection can be treated as an ensemble set of measurements of the receiver state. Noise Assisted Data Analysis is a growing field of study with significant potential for aiding the understanding and modeling of non stationary processes. Typically, NADA entails adding noise to a signal to produce an ensemble set on which statistical analysis is performed. Alternatively as in radiometric measurements, mixing a signal with calibrated noise provides, through the calibration process, the means to detect deviations from the stationary assumption and thereby a measurement tool to characterize the signal's non stationary properties. Data sets comprised of calibrated noise measurements have been limited to those collected with naturally occurring fluctuations in the radiometer receiver. To examine the application of NADA using calibrated noise, a Receiver Gain Modulation Circuit (RGMC) was designed and built to modulate the gain of a radiometer receiver using an external signal. In 2010, an RGMC was installed and operated at the National Institute of Standards and Techniques (NIST) using their Noise Figure Radiometer (NFRad) and national standard noise references. The data collected is the first known set of calibrated noise measurements from a receiver with an externally modulated gain. As an initial step, sinusoidal and step-function signals were used to modulate the receiver gain, to evaluate the circuit characteristics and to study the performance of a variety of calibration algorithms. The receiver noise temperature and time-bandwidth product of the NFRad are calculated from the data. Statistical analysis using temporal-dependent calibration algorithms reveals that the natural occurring fluctuations in the receiver are stationary over long intervals (100s of seconds); however the receiver exhibits local non stationarity over the interval over which one set of reference measurements are collected. A variety of calibration algorithms have been applied to the data to assess algorithms' performance with the gain fluctuation signals. This presentation will describe the RGMC, experiment design and a comparative analysis of calibration algorithms.

Racette, Paul↗

Thick-target bremsstrahlung interpretation of short time-scale solar hard X-ray features

Steady-state analyses of bremsstrahlung hard X-ray production in solar flares are appropriate only if the lifetime of the high energy electrons in the X-ray source is much shorter than the duration of the observed X-ray burst. For a thick-target nonthermal model, this implies that a full time-dependent analysis is required when the duration of the burst is comparable to the collisional lifetime of the injected electrons, in turn set by the lengths and densities of the flaring region. In this paper we present the results of such a time-dependent analysis, and we point out that the intrinsic temporal signature of the thick-target production mechanism, caused by the finite travel time of the electrons through the target, may indeed rule out such a mechanism for extremely short duration hard X-ray events.

Emslie, A. G.↗

Weekly Gridded Aquarius L-band Radiometer-Scatterometer Observations and Salinity Retrievals over the Polar Regions - Part 2: Initial Product Analysis

Following the development and availability of Aquarius weekly polar-gridded products, this study presents the spatial and temporal radiometer and scatterometer observations at L band (frequency1.4 GHz) over the cryosphere including the Greenland and Antarctic ice sheets, sea ice in both hemispheres, and over sub-Arctic land for monitoring the soil freeze-thaw state. We provide multiple examples of scientific applications for the L-band data over the cryosphere. For example, we show that over the Greenland Ice Sheet, the unusual 2012 melt event lead to an L-band brightness temperature (TB) sustained decrease of 5 K at horizontal polarization. Over the Antarctic ice sheet, normalized radar cross section (NRCS) observations recorded during ascending and descending orbits are significantly different, highlighting the anisotropy of the ice cover. Over sub-Arctic land, both passive and active observations show distinct values depending on the soil physical state (freeze-thaw). Aquarius sea surface salinity (SSS) retrievals in the polar waters are also presented. SSS variations could serve as an indicator of fresh water input to the ocean from the cryosphere, however the presence of sea ice often contaminates the SSS retrievals, hindering the analysis. The weekly grided Aquarius L-band products used a redistributed by the US Snow and Ice Data Center at http:nsidc.orgdataaquariusindex.html, and show potential for cryospheric studies.

Aquarius↗

Midwest Water Resources II: Evaluating Evapotranspiration with NASA Earth Observations and In Situ Observations to Understand Water Balance in Midwest Agriculture

Seasonal water variability in the midwestern United States extensively affects the agricultural community, as it impacts irrigation schedules, growing seasons, and overall ecosystem function. Evapotranspiration (ET) is a critical climatic variable in the water cycle and is used to evaluate spatiotemporal trends in drought and flood conditions. The NASA DEVELOP team partnered with the United States Department of Agriculture (USDA) Midwest Climate Hub, the Minnesota Department of Agriculture, the Illinois State Water Survey, and Michigan State University to compare remotely sensed ET products with in situ observations from January 2001 through December 2020. Remotely sensed actual ET (aET) data were sourced from NASA’s Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and reference ET (refET) data were derived from the Gridded Surface Meteorological (gridMET) dataset. For in situ comparison, aET data were downloaded from the AmeriFlux database while refET data were collected from the Illinois Climate Network and Michigan State University’s Enviro-weather database. For a holistic assessment of ET, this project generated comparisons between remotely sensed and in situ observations, calculated descriptive statistics for validation between refET datasets, and spatially produced statistical validation maps regarding in situ sites. The temporal and spatial gaps of AmeriFlux data limited aET analysis. This comparative assessment of ET products across the Midwest can be used by project partners to assess regional water trends and guide future land management decisions.

Addison Pletcher↗

The Version 2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979-Present)

The Global Precipitation Climatology Project (GPCP) Version 2 Monthly Precipitation Analysis is described. This globally complete, monthly analysis of surface precipitation at 2.5 degrees x 2.5 degrees latitude-longitude resolution is available from January 1979 to the present. It is a merged analysis that incorporates precipitation estimates from low-orbit-satellite microwave data, geosynchronous-orbit-satellite infrared data, and rain gauge observations. The merging approach utilizes the higher accuracy of the low-orbit microwave observations to calibrate, or adjust, the more frequent geosynchronous infrared observations. The data set is extended back into the premicrowave era (before 1987) by using infrared-only observations calibrated to the microwave-based analysis of the later years. The combined satellite-based product is adjusted by the raingauge analysis. This monthly analysis is the foundation for the GPCP suite of products including those at finer temporal resolution, satellite estimate, and error estimates for each field. The 23-year GPCP climatology is characterized, along with time and space variations of precipitation.

Adler, Robert F.↗

PBL Height from AIRS, GPS RO, and MERRA-2 Products in NASA GES DISC and Their 10 Year Seasonal Mean Intercomparison

Within the planetary boundary layer (PBL), surface forcing response, drag, turbulence, and vertical mixing are important processes and play a more critical role here than in the overlying “free atmosphere”. The PBL Height (PBLH) is an important parameter in climate models, weather forecasts, and air quality prediction. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provides data processing, archiving, and distribution services for numerous Earth science products. PBLH is a parameter in three products served by GES DISC, which are from the Atmospheric Infrared Sounder (AIRS), the Global Positioning System (GPS) radio occultation (RO) experiment, and the NASA reanalysis product Modern-Era Retrospective analysis for Research and Applications – 2 (MERRA-2). These products have different spatial and temporal resolutions and coverages, and their PBLH definitions are also different. To better serve the PBL research community, we have summarized the specifications of these products. A ten-year seasonal mean intercomparison is also conducted to provide further guidance to users. The intercomparison results show that MERRA-2 has a much shallower PBL than AIRS and GPS RO. An experimental study indicates the different PBLH definition in MERRA-2 caused smaller values of PBLH. The improvement of the water vapor retrieval in AIRS version 7 over version 6 results in the version 7 PBLH agreeing better with GPS RO and MERRA-2 than version 6, especially near the equator and low latitudes.

Feng Ding↗

The connection between forest degradation and urban energy demand in sub-Saharan Africa: a characterization based on high-resolution remote sensing data

Charcoal is a key energy source for urban households in sub-Saharan Africa and charcoal production is the main cause of forest degradation across the region. We used multitemporal high-resolution remote sensing optical imagery to quantify the extent and intensity of forest degradation associated with charcoal production and its impact on forest carbon stocks for the main supplying area of an African capital. This analysis documents the advance of forest degradation and quantifies its aboveground biomass removals over a seven-year period, registering that, between 2013 and 2016, the average annual area under charcoal production was 103 sq.km and the annual aboveground biomass removals reached 1081 000 (SD = 2461) Mg. Kiln densities in the study area rose to 2 kilns-ha, with an average of 90.7 Mg ha^(−1) of extracted aboveground biomass. Charcoal production was responsible for the degradation of 55.5% of the mopane woodlands in the study area between 2013 and 2019. We estimated post-disturbance recovery times using an ecosystem model calibrated for the study area. The simulations showed that recovery times could require up to 150 years for current aboveground biomass extraction rates. The results of the remote sensing analysis and the simulations of the ecosystem model corroborate the unsustainability of the present patterns of charcoal production. The detailed characterization of the spatial and temporal patterns of charcoal production was combined with household survey information to quantify the impact of the urban energy demand of the Maputo urban area on forest carbon stocks. The analysis shows that Maputo charcoal demand was responsible for the annual degradation of up to 175.3 sq.km and that the contribution of the study area to this demand fluctuated between 75% and 33% over the study period. The extent, advance pace and distance from urban centers documented in this study support the idea that forest degradation from charcoal production cannot merely be considered a peri-urban process. The intensity of the aboveground biomass (AGB) removals and its contribution to forest carbon stocks changes is significant at the national and regional levels.

F Sedano↗