Search NASASearch

SEARCH · Search NASA

Results for “proxy data”

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 145 records · Page 8

Data-Based Modeling of MMS/FPI Microchannel Plate Detector Resistance Using Temperature Proxy

The Magnetospheric MultiScale (MMS) space physics mission studies magnetic reconnection in the Earth’s magnetosphere. This work presents the results of an analysis of the resistance of the microchannel plates (MCP) included in the particle spectrometers, which are components of the Fast Plasma Investigation (FPI) aboard each MMS spacecraft. We found that a proxy for the resistance of the MCPs can be estimated from the MCP bias voltage and MCP high voltage power supply (HVPS) temperature. This estimation of MCP resistance allows for an accurate determination of nominal MCP current.

Joseph Patton

Empirical Indicators of Transmission Value in the Southeast United States

Concurrent differences in energy price between different parts of the electric grid are a key indicator of the value of additional transmission. In areas without a wholesale electricity market, such as the Southeast, an alternative indicator to price is the Federal Energy Regulatory Commission’s (FERC) system lambda data. This economic metric represents the minimized marginal production costs of thermal generators, including fuel and other variable operation and maintenance expenses. Balancing Authorities report a single system lambda for their entire balancing area. Most Southeastern lambdas exhibit sufficient price variation to support a transmission valuation analysis, although incomplete accounting of congestion costs or scarcity rents during peak load hours may underestimate the true value of transmission capacity. With transmission value defined as the annual average hourly absolute price difference between two regions and FERC’s system lambda data used as a price proxy, we find the following results in the Southeast region during 2012-2023 (reported in $\$2024$/MWh): Intra‐regional findings: Annual averages historically span $\$2$–$\$28$/MWh and average $\$12$/MWh in SERTP and span $\$4$–$\$19$/MWh and average $\$9$/MWh in FRCC, disregarding transmission value driven by anomalous data. The ranges of transmission value reported here are large, spanning an order of magnitude in some cases. Much of this variation is driven by year-to-year changes, with 2022 having a particularly high intra-regional transmission value due to elevated natural gas prices. Inter‐regional corridors: Annual average transmission values across three broader regions range from $\$6$ to $\$28$/MWh with a long-term average of $\$11$/MWh. Much of the transmission value is concentrated in a small portion of hours. Across all regions, severe weather—particularly polar vortex events in January 2018, February 2021, and December 2022—drives the largest price spreads. Seasonal patterns also emerge, with summer afternoons and fall mornings contributing consistently to transmission value, as for example between MISO and SOCO in 2023.

24 POWER TRANSMISSION AND DISTRIBUTION

Reevaluation of Stratospheric Ozone Trends From SAGE II Data Using a Simultaneous Temporal and Spatial Analysis

This paper details a new method of regression for sparsely sampled data sets for use with time-series analysis, in particular the Stratospheric Aerosol and Gas Experiment (SAGE) II ozone data set. Non-uniform spatial, temporal, and diurnal sampling present in the data set result in biased values for the long-term trend if not accounted for. This new method is performed close to the native resolution of measurements and is a simultaneous temporal and spatial analysis that accounts for potential diurnal ozone variation. Results show biases, introduced by the way data is prepared for use with traditional methods, can be as high as 10%. Derived long-term changes show declines in ozone similar to other studies but very different trends in the presumed recovery period, with differences up to 2% per decade. The regression model allows for a variable turnaround time and reveals a hemispheric asymmetry in derived trends in the middle to upper stratosphere. Similar methodology is also applied to SAGE II aerosol optical depth data to create a new volcanic proxy that covers the SAGE II mission period. Ultimately this technique may be extensible towards the inclusion of multiple data sets without the need for homogenization.

Damadeo, R. P.

Continuous snow depth and temperature measurements from dense network of above-ground distributed temperature profiling systems from 2021-09-23 to 2024-08-23, Seward Peninsula, Alaska

The dataset contains temperature measurements from distributed temperature profiling (DTP) systems (Dafflon et al., 2022; Wielandt et al., 2022; Wang et al., 2024a; Fiolleau et al., 2024) deployed vertically above the ground surface at a large number of locations from 2021 to 2024. The research is designed to improve understanding of the local heterogeneity in snow depth and snow thermal insulation dynamics, as well as their interactions in a discontinuous permafrost region (Wang et al., 2025). The DTP systems were deployed at 96 locations in a watershed along the Nome-Teller road at mile marker 27 (T27) and at 54 locations on a hillslope along the Kougarok road at mile marker 64 (K64) in the Seward Peninsula, Alaska. The probe location information is stored in Probe_locations_T27.csv and Probe_locations_K64.csv. Temperature measurements were recorded at 15-minute intervals using high-precision digital sensors (accuracy: ±0.1°C, resolution: 0.0078°C). The temperature probes, either 1.4 m or 1.6 m long, contain sensors spaced every 5 cm or 10 cm along their length. The temperature data are stored in compressed files following the format: DTP_snow_air_temperature_(site)_(start)_(end).zip, where site is either T27 or K64, and start and end represent the time series period. Within each ZIP file, individual CSV files are named by probe ID and contain temperature records at different heights above the ground surface.This dataset also includes derived snow depth time series over three snow seasons, estimated from temperature measurements. Snow depth was estimated by identifying the consecutive sensor pair that exhibited the largest drop in high-frequency temperature fluctuations (detailed in the methods). These data are stored in: Snow_depths_flags_(site)_(start)_(end).csv, which includes snow depth time series and corresponding quality flags (defined in the methods) from different probes. Additionally, the dataset includes derived metrics and supporting measurements at selected locations over two snow seasons, contributing to the manuscript of Wang et al., 2025. These locations were chosen based on the availability of high-quality snow depth time series during both seasons. The additional data include: (1) Air temperature proxies measured from the top sensors on the pole when they were not buried by snow, stored in Air_temperature_proxies_(site)_(start)_(end).csv (2) Ground interface temperature, recorded at 3 cm above the ground, stored in Ground_interface_temperature_(site)_(start)_(end).csv (3) Site characteristics, including vegetation height, elevation, and the topographic position index (TPI) within a 50 m radius, stored in Selected_probe_locations_gps_vegheight_tpi_elevation_(site).csv. These metrics were derived from 1 m resolution summer LiDAR-based digital elevation models and digital surface models from Singhania et al., 2023, DOI:10.5440/1832016. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv.This dataset is an updated version of a previous archive (Wang et al., 2024b, DOI: 10.15485/2475020), incorporating multiple seasons and improved snow depth estimation. Please note that due to large amount of information present in this dataset, many specificities associated with the acquisition of snow temperature, air temperature proxy and estimation of snow depth, and the future archiving of additional datasets on the soil temperature, thaw depth and soil characteristics at these locations, the author would welcome being contacted by people planning to use this dataset.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES

Preparations for Integrating Space-Based Total Lightning Observations into Forecast Operations

NASA's Short-term Prediction Research and Transition (SPoRT) Center has been a leader in collaborating with the United States National Weather Service (NWS) offices to integrate ground-based total lightning (intra-cloud and cloud-to-ground) observations into the real-time operational environment. For much of these collaborations, the emphasis has been on training, dissemination of data to the NWS AWIPS system, and focusing on the utility of these data in the warning decision support process. A shift away from this paradigm has occurred more recently for several reasons. For one, SPoRT's collaborations have expanded to new partners, including emergency managers and the aviation community. Additionally, and most importantly, is the impending launch of the GOES-R Geostationary Lightning Mapper (GLM). This has led to collaborative efforts to focus on additional forecast needs, new data displays, develop training for GLM uses based on the lessons learned from ground-based lightning mapping arrays, and ways to better relate total lightning data to other meteorological parameters. This presentation will focus on these efforts to prepare the operational end user community for GLM with an eye towards sharing lessons learned as EUMETSAT prepares for the Meteosat Third Generation Lightning Imager. This will focus on both software and training needs. In particular, SPoRT has worked closely with the Meteorological Development Laboratory to create the total lightning tracking tool. This software allows for NWS forecasters to manually track storms of interest and display a time series trend of observations. This tool also has been expanded to work on any gridded data set allowing for easy visual comparisons of multiple parameters in addition to total lightning. A new web display has been developed for the ground-based observations that can be easily extended to satellite observations. This paves the way for new collaborations outside of the NWS, both domestically and internationally, as the web display will be functional on PCs and mobile devices. Furthermore, SPoRT has helped developed the software plug-in to visualize GLM data. Examples using the official GLM proxy product will be used to provide a glimpse as to what real-time GLM and likely MTG-LI data will be in the near future.

GLM

Does The Nile Reflect Solar Variability?

Historical records of the Nile water level provide a unique opportunity to investigate the possibility that solar variability influences the Earth's climate. Particularly important are the annual records of the water level, which are uninterrupted for the years 622-1470 A.D. These records are non-stationary, so that standard spectral analyses cannot adequately characterize them. Here the Empirical Mode Decomposition technique, which is designed to deal with nonstationary, nonlinear time series, becomes useful. It allows the identification of two characteristic time scales in the water level data that can be linked to solar variability: the 88 year period and a time scale of about 200 years. These time scales are also present in the concurrent aurora data. Auroras are driven by coronal mass ejections and the rate of auroras is an excellent proxy for solar variability. Analysis of auroral data contemporaneous with the Nile data shows peaks at 88 years and about 200 years. This suggests a physical link between solar variability and the low frequency variations of the Nile water level. The link involves the influence of solar variability on the North Annual Mode of atmospheric variability and its North Atlantic and Indian Oceans patterns that affect rainfall over Eastern Equatorial Africa where the Nile originates.

statistical methods

Estimating Return on Investment for Energy Technical Assistance Programs

The U.S. Department of Energy's Office of State and Community Energy Programs engaged the National Laboratory of the Rockies to assess the return on investment (ROI) of technical assistance (TA) programs that support state, local, and Tribal energy planning. Although TA delivers value through capacity building, stakeholder engagement, and knowledge transfer, these benefits are often intangible and challenging to monetize. This study reviews existing ROI frameworks and synthesizes the most relevant elements into a hybrid approach tailored to energy TA programs. The proposed framework integrates monetary and non-monetary outcomes through early logic model development, baseline data collection, and the use of proxies for intangible benefits. As a case study, this paper applies this approach to the Communities Local Energy Action Program (Communities LEAP), demonstrating how ROI can inform program design, data strategy, and performance assessment. Findings underscore that ROI should be applied selectively and planned from the outset to ensure data alignment and attribution accuracy. The framework offers TA practitioners a structured approach that can be leveraged for future programs to evaluate and communicate the multifaceted value of TA investments.

29 ENERGY PLANNING, POLICY, AND ECONOMY

Synergistic Use of Hyperspectral UV-Visible OMI and Broadband Meteorological Imager MODIS Data for a Merged Aerosol Product

The retrieval of optimal aerosol datasets by the synergistic use of hyperspectral ultraviolet(UV)–visible and broadband meteorological imager (MI) techniques was investigated. The Aura Ozone Monitoring Instrument (OMI) Level 1B (L1B) was used as a proxy for hyperspectral UV–visible instrument data to which the Geostationary Environment Monitoring Spectrometer (GEMS) aerosol algorithm was applied. Moderate-Resolution Imaging Spectroradiometer (MODIS) L1B and dark target aerosol Level 2 (L2) data were used with a broadband MI to take advantage of the consistent time gap between the MODIS and the OMI. First, the use of cloud mask information from the MI infrared (IR) channel was tested for synergy. High-spatial-resolution and IR channels of the MI helped mask cirrus and sub-pixel cloud contamination of GEMS aerosol, as clearly seen in aerosol optical depth (AOD) validation with Aerosol Robotic Network (AERONET) data. Second, dust aerosols were distinguished in the GEMS aerosol-type classification algorithm by calculating the total dust confidence index (TDCI) from MODIS L1B IR channels. Statistical analysis indicates that the Probability of Correct Detection (POCD) between the forward and inversion aerosol dust models (DS) was increased from 72% to 94% by use of the TDCI for GEMS aerosol-type classification, and updated aerosol types were then applied to the GEMS algorithm. Use of the TDCI for DS type classification in the GEMS retrieval procedure gave improved single-scattering albedo (SSA) values for absorbing fine pollution particles (BC) and DS aerosols. Aerosol layer height (ALH) retrieved from GEMS was compared with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data, which provides high-resolution vertical aerosol profile information. The CALIOP ALH was calculated from total attenuated backscatter data at 1064 nm, which is identical to the definition of GEMS ALH. Application of the TDCI value reduced the median bias of GEMS ALH data slightly. The GEMS ALH bias approximates zero, especially for GEMS AOD values of>~0.4 and GEMS SSA values of<~0.95.Finally, the AOD products from the GEMS algorithm and MI were used in aerosol merging with the maximum-likelihood estimation method, based on a weighting factor derived from the standard deviation of the original AOD products. With the advantage of the UV–visible channel in retrieving aerosol properties over bright surfaces, the combined AOD products demonstrated better spatial data availability than the original AOD products, with comparable accuracy. Furthermore, pixel-level error analysis of GEMS AOD data indicates improvement through MI synergy.

aerosol

The diatom record from beneath the West Antarctic Ice Sheet and the global proxy perspective

Recent glaciological evaluation and modeling of the marine-based West Antarctic Ice Sheet (WAIS) support the possibility that the WAIS disintegrated during one or more Pleistocene interglacial period(s). The magnitude of sea level and oxygen isotope variation during certain late-Pleistocene interglacial periods is also consistent with the possibility of major retreat of the WAIS. Although oxygen isotopes from deep-sea sediments provide the best available proxy record for global ice volume (despite the ambiguities in the record), the source of ice volume changes must be hypothesized. Based on the intensity of interglacial isotopic shifts recorded in Southern Ocean marine sedimentary records, stage 11 (400,000 years ago) is the strongest candidate for WAIS collapse, but the records for stages 9, 7, and 5.5 are all consistent with the possibility of multiple late-Pleistocene collapses. Seismic reflection studies through the WAIS have revealed thick successions of strata with seismic characteristics comparable to upper Tertiary marine sediments. Small samples of glacial diamictons from beneath the ice sheet have been collected via hot-water drilled access holes. These sediments include mixed diatom assemblages of varying ages. Late-Miocene diatoms dominate many samples, probably reflecting marine deposition in West Antarctic basins prior to development of a dominantly glacial phase in West Antarctica. In addition to late-Miocene diatoms, samples from Upstream B (1988/89) contain rare post-Miocene diatoms, many of which imply deposition in the West Antarctic interior during one or more Pleistocene deglaciation periods. Age-diagnostic fossils in glacial sediments beneath ice sheets provide relatively coarse chronostratigraphic control, but they do contain direct evidence of regional deglaciation. Thus, sub-glacial till samples provide the evidence regarding the source of ice sheet variability seen in well-dated proxy records. Combined, these independent data sets can provide a more comprehensive and less speculative interpretation of the history of past glacial minima in currently glaciated polar regions.

Scherer, Reed P.

Investigation of several proxies to estimate sulfuric acid concentration in volcanic plume conditions

Sulfuric acid (H2SO4) is commonly accepted as a key precursor for atmospheric new particle formation (NPF). However, direct measurements of [H2SO4] remain challenging, thus preventing the determination of this important quantity, and, consequently, a complete understanding of its contribution to the NPF process. Several proxies have been developed to bridge the gaps, but their ability to predict [H2SO4] in very specific conditions such as those encountered in volcanic plumes (including in particular high sulphur dioxide mixing ratios) has not been evaluated so far. In this context, the main objective of the present study was to develop new proxies for daytime [H2SO4] in volcanic plume conditions and compare their performance to that of the proxies available in the literature. In specific, the data collected at Maïdo during the OCTAVE 2018 campaign, in the volcanic eruption plume of the Piton de la Fournaise, were first used to derive seven proxies based on the knowledge of sulphur dioxide (SO2) mixing ratio, global radiation, condensation sink (CS) and relative humidity (RH). In three of the seven proxies (F1–F3), all variables were given equal weight in the prediction of [H2SO4], while adjusted powers were allowed for the different variables in the other four proxies (A1–A4). Proxies A1–A4 were overall found to perform better compared to F1–F3, with, in specific, improved predictive ability for [H2SO4] > 2 × 108 per cu.cm. The CS was observed to play an important role in regulating [H2SO4], while, in contrast, the inclusion of RH did not improve the predictions. A last expression accounting for an additional sink term related to cluster formation, S1, was also tested and showed a very good predictive ability over the whole range of measured [H2SO4]. The newly developed proxies were in a second step further evaluated using airborne measurements performed in the passive degassing plume of Etna during the STRAP 2016 campaign. Increased correlations between observed and predicted [H2SO4] were obtained when the dependence of predicted [H2SO4] over CS was the lowest, and when the dependence over [SO2] was concurrently the highest. The best predictions were finally retrieved by the simple formulation of F2 (in which [SO2] and radiation alone were assumed to explain the variations of [H2SO4] with equal contributions), with a pre factor adapted to the STRAP data. All in all, our results illustrate the fairly good capacity of the proxy available in the literature to describe [H2SO4] in volcanic plume conditions, but highlight at the same time the benefit of the newly developed proxies for the prediction of the highest concentrations ([H2SO4] > 2–3 × 108 per cu.cm). Also, the contrasting behaviours of the new proxies in the two investigated datasets indicate that in volcanic plumes like in other environments, the relevance of a proxy can be affected by changes in environmental conditions, and that location specific coefficients do logically improve the predictions.

new particle formation

X-ray bright points and He I lambda 10830 dark points

Using near-simultaneous full disk Solar X-ray images and He I 10830 lambda, spectroheliograms from three recent rocket flights, dark points identified on the He I maps were compared with X-ray bright points identified on the X-ray images. It was found that for the largest and most obvious features there is a strong correlation: most He I dark points correspond to X-ray bright points. However, about 2/3 of the X-ray bright points were not identified on the basis of the helium data alone. Once an X-ray feature is identified it is almost always possible to find an underlying dark patch of enhanced He I absorption which, however, would not a priori have been selected as a dark point. Therefore, the He I dark points, using current selection criteria, cannot be used as a one-to-one proxy for the X-ray data. He I dark points do, however, identify the locations of the stronger X-ray bright points.

Golub, L.

X-ray bright points and He I lambda 10830 dark points

Using near-simultaneous full disk Solar X-ray images and He I 10830 lambda, spectroheliograms from three recent rocket flights, dark points identified on the He I maps were compared with x-ray bright points identified on the X-ray images. It was found that for the largest and most obvious features there is a strong correlation: most He I dark points correspond to X-ray bright points. However, about 2/3 of the X-ray bright points were not identified on the basis of the helium data alone. Once an X-ray feature is identified it is almost always possible to find an underlying dark patch of enhanced He I absorption which, however, would not a priori have been selected as a dark point. Therefore, the He I dark points, using current selection criteria, cannot be used as a one-to-one proxy for the X-ray data. He I dark points do, however, identify the locations of the stronger X-ray bright points.

Golub, L.

Updates in Developing a Prototype Science Pipeline and Full-Volume, Global Hyperspectral Synthetic Data Sets for NASA’s Earth System Observatory’s Upcoming Surface, Biology and Geology Mission

The Surface Biology and Geology (SBG) mission recently passed mission confirmation review and has entered phase A – design and development. SBG will acquire high resolution solar-reflected spectroscopy and thermal infrared observations at a data rate of ~2.5 TB/day and generate products at ~40 TB/day. Given that the per-day volume is greater than NASA’s total extant airborne hyperspectral data collection, collecting, processing, disseminating, and exploiting the SBG data present new challenges. To meet these challenges, we have developed a prototype science pipeline and a full-volume global hyperspectral synthetic data set to help prepare for SBG’s flight (see poster GC42D-0730). Our science pipeline is based on the science processing technology developed for NASA’s Kepler and TESS planet-hunting missions. The pipeline infrastructure, Ziggy, provides a scalable architecture for robust, repeatable, and replicable science and application products that can be run on a range of systems from a laptop to the cloud or a supercomputer. Ziggy is compliant with NASA Procedural Requirement (NPR) 7150.2C, is at a technical readiness level (TRL) of 7 and has been released to github.com/nasa/ziggy. We integrated Ziggy with EO-1/Hyperion workflows to build a prototype pipeline and ingested the 17-year mission archive that provides globally sampled visible through shortwave infrared spectra that are representative of SBG data types and volumes. We fully implemented the first stage and processed the entire 55 TB Hyperion data set from the raw data (Level 0) to top-of-the-atmosphere radiance (Level 1R). We are currently evaluating the ISOFIT atmospheric correction module to convert the L1R data to surface reflectance (Level 2) before reprocessing the full data set to L2. Crosschecks are being performed with RadCalNet as well as with coincident observations by AVIRIS. We are also investigating modern methods for georectifying the Hyperion scenes. Finally, we describe an analysis of the cost to conduct forward processing and reprocessing campaigns for SBG on HECC with dedicated compute and storage resources using the resurrected Hyperion pipeline as a proxy for full-volume SBG data. The analysis demonstrates that SBG L0 data can be processed to L2 on HECC with full reprocessing campaigns every two years for ~$2.6M over a 7-year lifespan. Moreover, 69% of the system capacity would be available for other activities, possibly enabling future open-source science activities, including algorithm development, L3+ processing, .etc.

ESD

Modeling Quiet Solar Luminosity Variability from TSI Satellite Measurements and Proxy Models during 1980-2018

A continuous record of direct total solar irradiance (TSI) observations began with a series of satellite experiments in 1978. This record requires comparisons of overlapping satellite observations with adequate relative precisions to provide useful long term TSI trend information. Herein we briefly review the active cavity radiometer irradiance monitor physikalisch-meteorologisches observatorium davos (ACRIM-PMOD) TSI composite controversy regarding how the total solar irradiance (TSI) has evolved since 1978 and about whether TSI significantly increased or slightly decreased from 1980 to 2000. The main question is whether TSI increased or decreased during the so-called ACRIM-gap period from 1989 to 1992. There is significant discrepancy between TSI proxy models and observations before and after the gap, which requires a careful revisit of the data analysis and modeling performed during the ACRIM-gap period. In this study, we use three recently proposed TSI proxy models that do not present any TSI increase during the ACRIM-gap, and show that they agree with the TSI data only from 1996 to 2016. However, these same models significantly diverge from the observations from 1981 and 1996. Thus, the scaling errors must be different between the two periods, which suggests errors in these models. By adjusting the TSI proxy models to agree with the data patterns before and after the ACRIM-gap, we found that these models miss a slowly varying TSI component. The adjusted models suggest that the quiet solar luminosity increased from the 1986 to the 1996 TSI minimum by about 0.45 W/m2 reaching a peak near 2000 and decreased by about 0.15 W/m2 from the 1996 to the 2008 TSI cycle minimum. This pattern is found to be compatible with the ACRIM TSI composite and confirms the ACRIM TSI increasing trend from 1980 to 2000, followed by a long-term decreasing trend since.

total solar irradiance

Assessing Uncertainties in Gridded Emissions: A Case Study for Fossil Fuel Carbon Dioxide (FFCO2) Emission Data

Fossil fuel carbon dioxide (CO2) emissions (FFCO2) are the largest input to the global carbon cycle on a decadal time scale. Because total emissions are assumed to be reasonably well constrained by fuel statistics, FFCO2 often serves as a reference in order to deduce carbon uptake by poorly understood terrestrial and ocean sinks. Conventional atmospheric CO2 flux inversions solve for spatially explicit regional sources and sinks and estimate land and ocean fluxes by subtracting FFCO2. Thus, errors in FFCO2 can propagate into the final inferred flux estimates. Gridded emissions are often based on disaggregation of emissions estimated at national or regional level. Although national and regional total FFCO2 are well known, gridded emission fields are subject to additional uncertainties due to the emission disaggregation. Assessing such uncertainties is often challenging because of the lack of physical measurements for evaluation. We first review difficulties in assessing uncertainties associated with gridded FFCO2 emission data and present several approaches for evaluation of such uncertainties at multiple scales. Given known limitations, inter-emission data differences are often used as a proxy for the uncertainty. The popular approach allows us to characterize differences in emissions, but does not allow us to fully quantify emission disaggregation biases. Our work aims to vicariously evaluate FFCO2 emission data using atmospheric models and measurements. We show a global simulation experiment where uncertainty estimates are propagated as an atmospheric tracer (uncertainty tracer) alongside CO2 in NASA's GEOS model and discuss implications of FFCO2 uncertainties in the context of flux inversions. We also demonstrate the use of high resolution urban CO2 simulations as a tool for objectively evaluating FFCO2 data over intense emission regions. Though this study focuses on FFCO2 emission data, the outcome of this study could also help improve the knowledge of similar gridded emissions data for non-CO2 compounds with similar emission characteristics.

FFCO

Historic Global Biomass Burning Emissions for CMIP6 (BB4CMIP) Based on Merging Satellite Observations with Proxies and Fire Models (1750-2015)

Fires have influenced atmospheric composition and climate since the rise of vascular plants, and satellite data has shown the overall global extent of fires. Our knowledge of historic fire emissions has progressively improved over the past decades due mostly to the development of new proxies and the improvement of fire models. Currently there is a suite of proxies including sedimentary charcoal records, measurements of fire-emitted trace gases and black carbon stored in ice and firn, and visibility observations. These proxies provide opportunities to extrapolate emissions estimates based on satellite data starting in 1997, but each proxy has strengths and weaknesses regarding, for example, the spatial and temporal extents over which they are representative. We developed a new historic biomass burning emissions dataset starting in 1750 that merges the satellite record with several existing proxies, and uses the average of six models from the Fire Model Intercomparison Project (FireMIP) protocol to estimate emissions when the available proxies had limited coverage. According to our approach, global biomass burning emissions were relatively constant with 10-year averages varying between 1.8 and 2.3 petagrams of carbon per year. Carbon emissions increased only slightly over the full time period and peaked during the 1990's after which they decreased gradually. There is substantial uncertainty in these estimates and patterns varied depending on choices regarding data representation, especially on regional scales. The observed pattern in fire carbon emissions is for a large part driven by African fires, which accounted for 58 percent of global fire carbon emissions. African fire emissions declined since about 1950 due to conversion of savanna to cropland, and this decrease is partially compensated for by increasing emissions in deforestation zones of South America and Asia. These global fire emissions estimates are mostly suited for global analyses and will be used in the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations.

Biomass burning; Satellite observations; computer

Coincident 1.3-year Periodicities in the ap Geomagnetic Index and the Solar Wind

Recent observations show an approximately 1.3-year period in the speed of the solar wind detected by the IMP 8 and Voyager 2 spacecraft. A similar period is also seen in the north-south (GSE) component of the magnetic field observed by IMP 8. Since both parameters are commonly used as input to models of geomagnetic activity, the 'ap' index (a measure of geomagnetic disturbance) is examined to look for this periodicity. The Lomb-Scargle periodogram method is used on the ap, plasma, and magnetic field data during the 1973-1994 time range. A dynamic FFT periodogram method is also used to analyze the ap data during this time, as well as to look for periods present between 1932 and 1972. A clear 1.3-year periodicity is present in the post-1986 data when the same period is observed in the plasma and field data. The V(2)B(zsm) and V(2)B(s) proxies for geomagnetic activity also show this periodicity. However, the southward (GSM) component of the magnetic field does not have a 1.3-year period, and neither do solar wind or ap data from 1973-1985. This demonstrates that the ap geomagnetic index can act as a proxy for solar wind periodicities at this time scale. Historic ap data are examined, and show that a similar periodicity in ap exists around 1942. Since auroral data show a 1.4-year periodicity, all these similar periods may result from a common underlying solar mechanism.

Paularena, K. I.

A Testbed for Geomagnetic Data Assimilation

Geomagnetic data assimilation merges past and present-day observations of the Earth’s mag-netic field with numerical geodynamo models and the results are used to initialize forecasts.We present a new ‘proxy model’ that can be used to test, or rapidly prototype, numerical tech-niques for geomagnetic data assimilation. The basic idea for constructing a proxy is to capturethe conceptual difficulties one encounters when assimilating observations into high-resolution,3-D geodynamo simulations, but at a much lower computational cost. The framework of usingproxy models as ‘gate-keepers’ for numerical methods that could/should be considered formore extensive testing on operational models has proven useful in numerical weather predic-tion, where advances in data assimilation and, hence, improved forecast skill, are at least inpart enabled by the common use of a wide range of proxy models. We also present a largeset of systematic data assimilation experiments with the proxy to reveal the importance oflocalization and inflation in geomagnetic data assimilation.

K Gwirtz