Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatio-Temporal”

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 289 records · Page 16

Negative extreme events in gross primary productivity and their drivers in China during the past three decades

Climate extremes have remarkable impacts on ecosystems and are expected to increase with future global warming. However, only few studies have focused on the ecological extreme events and their drivers in China. In this study, we carried out an analysis of negative extreme events in gross primary productivity (GPP) in China and the sub-regions during 1982–2015, using monthly GPP simulated by 12 process-based models (TRENDYv6) and an observation-based model (Yao-GPP). Extremes were defined as the negative 5th percentile of GPP anomalies, which were further merged into individual extreme events using a three-dimensional contiguous algorithm. Spatio-temporal patterns of negative GPP anomalies were analyzed by taking the 1000 largest extreme events into consideration. Results showed that the effects of extreme events decreased annual GPP by 2.8% (i.e. 208 TgC/year) in TRENDY models and 2.3% (i.e. 151 TgC/year) in Yao-GPP. Hotspots of extreme GPP deficits were mainly observed in North China (−53 gC/sq.m year) in TRENDY models and Northeast China (−42 gC/sq.m year) in Yao-GPP. For China as a whole, attribution analyses suggested that extreme low precipitation was associated with 40%–50% of extreme negative GPP events. Most events in northern and western China could be explained by meteorological droughts (i.e. low precipitation) while GPP extreme events in southern China were more associated with temperature extremes, in particular with cold spells. GPP was revealed to be much more sensitive to heat/drought than to cold/wet extreme events. Combined with projected changes in climate extremes in China, GPP negative anomalies caused by drought events in northern China and by temperature extremes in southern China might be more prominent in the future.

Climate change↗

Mars Atmospheric Chemistry Simulations with the GEM-Mars General Circulation Model

General Circulation Models with interactive physical and chemistry processes are the state-of-the-art tools for an integrated view and understanding of the Martian atmosphere and climate system. The GEM-Mars model currently includes 16 tracers for chemical composition and applies a fully online, interactive calculation of the photo- and gas phase chemistry of carbon dioxide (CO2) and water vapor (H2O). These species largely control the chemical composition of the neutral Mars atmosphere through their photolysis products and their subsequent interactions. Water vapor undergoes a complex cycle on Mars as it is transported and interacts with ice reservoirs both at the planet's surface and at water ice clouds, which in turn provide radiative feedbacks. In the photochemical cycles involving CO2 and H2O, the abundances of 5 species have been reported by previous investigations with significant spatio-temporal coverage: CO2, H2O, CO, O3 and H2O2. This paper presents the current status of the atmospheric chemistry simulations in GEM-Mars by comparing them to a selection of these observational datasets as well as to oxygen dayglow emission from O2(a1∆g). The results are consistent with previous model-data comparisons and illustrate that the water cycle and the photochemistry are well implemented in the model. In particular, the simulation of the key reservoir species H2O2 provides a good match to the available data. Model-data biases for ozone columns and oxygen airglow are related to the simulated water vapor vertical profile, as these species have important column contributions from vertical layers at the top of the hygropause.

F Daerden↗

Global analysis and forecasts of carbon monoxide on Mars

Spatial and temporal variations in the Martian carbon monoxide (CO) cycle have been investigated through combining Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) retrievals of carbon monoxide mixing ratio and Mars Climate Sounder (MCS) temperature profiles with a Martian global circulation model (GCM) to produce the first global reanalysis of the Martian CO cycle. The reanalysis reduces the root mean square error between the forecast and CRISM CO retrievals by a factor of 2–4, dependent on the time of year. Forecasts initiated from the reanalysis show an improved match to standalone CRISM CO retrievals from northern winter solstice to northern summer solstice, indicating the benefit of data assimilation in simulating the CO abundance. The northern summer solstice CO minimum between 10°S – 50°S in the CRISM CO retrievals is found to be caused by a suppression of CO-enriched air leaking from the Hellas and Argyre basins, and is also strongly influenced by the amount of carbon dioxide condensing at the time. The GCM is the first reported to simulate the local CO minimum however it is weaker in strength, as CO-enriched air is also released northward in the area of Argyre basin in the GCM because of a shifted boundary between the southern polar vortex and mid-latitudes. The reanalysis, as a result of the assimilation of MCS temperature profiles, indicates that the polar vortex boundary is northward of Argyre basin and hence no northerly transport of CO-enriched air should be present over this location. Differences in local CO abundance between the reanalysis and GCM are also evident in spatio-temporal regions where no nearby CRISM CO retrievals are available as a result of alterations in local circulation and the sublimation/condensation of carbon dioxide through the combined assimilation of CRISM CO retrievals alongside MCS temperature profiles, illustrating how constraints can be imposed indirectly on the CO cycle through the powerful technique of data assimilation.

James A. Holmes↗

EOSDIS STAC Briefing

The Spatio-Temporal Asset Catalog (STAC) specification provides a common language to describe a range of geospatial information, so it can more easily be indexed and discovered. A 'spatiotemporal asset' is any file that represents information about the earth captured in a certain space and time. STAC provides a standards-based, web and cloud friendly cataloging specification that has seen a large amount of adoption in the cloud geospatial information system space. This presentation will demonstrate how NASA EOSDIS leverages STAC technology to provide value-added features to both our data discovery and transformation services. Finally, we will propose a way to improve our federated discovery capabilities using STAC.

NASA EOSDIS STAC Cloud catalog metadata↗

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We will provide a brief summary of the the progress made in the field to understand these systems theoretically and what new directions groups are pursuing. We will also report on our collaboration's progress to simulate these systems using general relativistic MHD and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure how the binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power. The novel computational methods enabling inclusion of the black holes in the domain, including multi-patch methods, will be described.

Scott C Noble↗

Large Eddy Simulations of Upper-Ocean Turbulence Across the Equatorial Pacific Cold Tongue

Turbulence in the equatorial Pacific cold tongue plays multiple roles in the sea-surface temperature budget. However, observations of turbulence are limited to a handful of field experiments lasting from a few days to one month at 0 N , 140 W, and the key turbulent fluxes can only be inferred. Although explicit simulations of turbulence resolve these fluxes, the required large eddy simulations (LES) are costly and thus also limited to a few weeks. A two-fold hypothesis motivates this contribution to the PUMP pre-field modelling experiments in support of TPOS: 1) that there exist useful relationships between the surface wind, wave, and buoyancy forcing, oceanic shear and stratification, and turbulent mixing, but 2) these relationships differ between regions of equatorial dynamics (f~0) and higher latitudes, such that the state space is insufficiently sampled to understand these relationships and represent them in ocean general circulation models (OGCM). To address the lack of observations, we simulated 35 continuous days of upper ocean turbulence using LES during October-November 1985 at sites in two regions (70 days in total) of the cold tongue: 0 N, 140 W and 3 N, 140 W. The model reproduces many features of the observed turbulence at 0,140W, including the “deep-cycle” observed by Gregg et al. (1985) Nature and the turbulent pulses every few hours during the night and in the early morning. Furthermore, the results at 3 N support the hypothesis advanced by Cherian et al. (2021) JPO that strong, but intermittent deep-cycle turbulence is observed to the north of the equator in the presence of strong shear associated with Tropical Instability Waves. Comparisons with the K-profile parameterization (KPP) scheme show that additional tuning of the KPP parameters can reduce differences between the turbulent fluxes produced by KPP and LES at both sites. The results will suggest modifications to OGCMs and already emphasize the need for a more extensive measurements of vertical profiles of shear, stratification and turbulent mixing to reveal the spatio-temporal structure of turbulence and how it relates to forcing, shear and stratification across the equatorial Pacific Cold tongue.

Large Eddy Simulations↗

Prompt Response of the Dayside Magnetosphere to Discrete Structures within the Sheath Region of a Coronal Mass Ejection

A sequence of discrete solar wind structures within the sheath region of an interplanetary coronal mass ejection (ICME) on Nov 06, 2015 caused a series of compressions and releases of the dayside magnetosphere. Each compression resulted in a brief adiabatic enhancement of ions (electrons) driving bursts of electromagnetic ion cyclotron (whistler mode chorus) wave growth across the dayside magnetosphere. Fine-structured rising tones were observed in the EMIC wave bursts, resulting in nonlinear scattering of relativistic electrons in the outer radiation belt. Multipoint observations allow us to study the spatial structure and evolution of these sheath structures as they propagate Earthward from L1 as well as the spatio-temporal characteristics of the magnetospheric response. This event highlights the importance of fine-scale solar wind structure, in particular within complex sheath regions, in driving dayside phenomena within the inner magnetosphere.

L. W. Blum↗

Validation of Soil Moisture Data Products from the NASA SMAP Mission

The National Aeronautics and Space Administration (NASA)Soil Moisture Active Passive(SMAP) mission has been validatingits soil moisture (SM) products since the start of data production onMarch 31, 2015. Prior to launch, the mission defined a set of criteria for core validation sites (CVS) that enable the testing of the key mission SM accuracy requirement(unbiased root-mean-square error <0.04 m3/m3). Thevalidation approach also includes other (“sparse network”) in situSM measurements, satellite SM products, model-based SM products, and field experiments. Over the past six years, the SMAP SM products have been analyzed with respect to thesereference data,and the analysis approaches themselves have been scrutinizedin an effort to best understand the products’ performance. Validation of themost recent SMAP Level 2 and 3 SMretrievalproducts (R17000) shows that the L-band (1.4 GHz) radiometer-based SM record continues to meet mission requirements. The products aregenerallyconsistentwith SM retrievals from the European Space Agency (ESA)Soil Moisture Ocean Salinity mission, althoughthere aredifferencesin some regions. The high-resolution (3-km) SM retrieval product,generated by combining Copernicus Sentinel-1 data with SMAP observations,performswithin expectations. Currently, however,there is limited availability of3-kmCVSdatato support extensive validation at this spatial scale. The most recent (version 5)SMAP Level4 SMdata assimilation productprovidingsurface and root-zone SM with complete spatio-temporal coverageat 9-km resolution also meets performance requirements. The SMAP SM validation program will continue throughout the mission life; futureplans include expanding ittoforestedand high-latituderegions

SMAP↗

Sea Surface Height Anomalies of the Arctic Ocean from Icesat-2: A First Examination and Comparisons with Cryosat-2

Accurately resolving spatio-temporal variations in sea surface height across the polar oceans is key to improving our understanding of ocean circulation variability and change. Here, we examine the first two years (2018-2020) of Arctic Ocean sea surface height anomalies (SSHA) from the photon-counting laser altimeter onboard NASA’s ICE, Cloud, and Land Elevation Satellite-2 (ICESat-2). ICESat-2 SSHA estimates are compared to estimates from ESA’s CryoSat-2 mission, including semi-synchronous along-track measurements from the recent CRYO2ICE orbit alignment campaign. There are documented residual centimeter-scale range biases between the ICESat-2 beams (in release-003 data) and we opted for a single-beam approach in our comparisons. We find good agreement in the along-track estimates (spatial correlation coefficient >0.8 and mean differences <0.03 m) as well as in the gridded monthly SSHA estimates (temporal correlation coefficient of 0.76 and a mean difference of 0.01 m) from the two altimeters, suggesting ICESat-2 adds to the CryoSat-2 SSHA estimates.

ICESat-2↗

Multi-Flight Common Routes

Flights often experience large delays when they are routed around weather. Multi-flight common routes advisories provide delay recovery by suggesting time-saving re-routes for groups of flights whose current weather-avoidance routes have become outdated because the weather has dissipated and/or moved away. The multi-flight common routes tool provides time-saving route change advisories taking into account flight plans, wind fields, and the spatio-temporal evolution of predicted convective weather. A graphical user interface enables these advisories to be easily modified by a traffic manager for possible operational implementation.

ATD-3↗

Air Pollution Scenario over Pakistan: Characterization and Ranking of Extremely Polluted Cities using Long-Term Concentrations of Aerosols and Trace Gases

Pakistan ranks third in the world in terms of mortality attributable to air pollution, with aerosol mass concentrations (PM2.5) consistently well above WHO (World Health Organization) air quality guidelines (AQG). However, regulation is dependent on a sparse network of air quality monitoring stations and insufficient ground data. This study utilizes long-term observations of aerosols and trace gases to characterize and rank the air pollution scenarios and pollution characteristics of 80 selected cities in Pakistan. Datasets used include (1) the Aqua and Terra (AquaTerra) MODIS (Moderate Resolution Imaging Spectroradiometer) Level 2 Collection 6.1 merged Dark Target and Deep Blue (DTB) aerosol optical depth (AOD) retrieval products; (2) the CAMS (Copernicus Atmosphere Monitoring Service) reanalysis PM1, PM2.5, and PM10 data; (3) the MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2) reanalysis PM2.5 data, (4) the OMI (Ozone Monitoring Instrument) tropospheric vertical column density (TVCD) of nitrogen dioxide (NO2), and VCD of sulfur dioxide (SO2) in the Planetary Boundary Layer (PBL), (5) the VIIRS (Visible Infrared Imaging Radiometer Suite) Nighttime Lights data, (6) MODIS Collection 6 Version 2 global monthly fire location data (MCD14ML), (7) population density, (8) MODIS Level 3 Collection 6 land cover types, (9) AERONET (AErosol RObotic NETwork) Version 3 Level 2.0 data, and (10) ground-based PM2.5 concentrations from air quality monitoring stations. Potential Source Contribution Function (PSCF) analyses were performed by integrating with ground-based PM2.5 concentrations and the NOAA (National Oceanic and Atmospheric Administration) HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) air parcel back trajectories to identify potential pollution source areas which are responsible for extreme air pollution in Pakistan. Results show that the ranking of the top polluted cities depends on the type of pollutant considered and the metric used. For example, Jhang, Multan, and Vehari were characterized as the top three polluted cities in Pakistan when considering AquaTerra DTB AOD products; for PM1, PM2.5, and PM10 Lahore, Gujranwala, and Okara were the top three; for tropospheric NO2 VCD Lahore, Rawalpindi, and Islamabad and for PBL SO2 VCD Lahore, Mirpur, and Gujranwala. The results demonstrate that Pakistan’s entire population has been exposed to high PM2.5 concentrations for many years, with a mean annual value of 54.7 μg/cu. m, over all Pakistan from 2003 to 2020. This value exceeds Pakistan’s National Environmental Quality Standards (Pak-NEQS, i.e., <15 μg/cu. m annual mean) for ambient air defined by the Pakistan Environmental Protection Agency (Pak-EPA) as well as the WHO Interim Target-1 (i.e., mean annual PM2.5 <35 μg/cu. m). The spatial analyses of the concentrations of aerosols and trace gases in terms of population density, nighttime lights, land cover types, and fire location data, and the PSCF analysis indicate that Pakistan’s air quality is strongly affected by anthropogenic sources inside of Pakistan, with contributions from surrounding countries. Statistically significant positive (increasing) trends in PM1, PM2.5, PM10, tropospheric NO2 VCD, and SO2 VCD were observed in ~89%, ~67%, ~48%, 91%, and ~88% of the Pakistani cities (80 cities), respectively. This comprehensive analysis of aerosol and trace gas levels, their characteristics in spatio-temporal domains, and their trends over Pakistan, is the first of its kind. Results will be helpful to the Ministry of Climate Change (Government of Pakistan), Pak-EPA, SUPARCO (Pakistan Space and Upper Atmosphere Research Commission), policymakers, and the local research community to mitigate air pollution and its effects on human health.

Muhammad Bilal↗

Assessment of aerosol burden over Ghana

Although air pollution in Ghana is ranked number one in environmental health threats to public health and sixth to cause of deaths, routine monitoring is rare. This paper presents fourteen years (2005-2018) assessment of aerosol optical depth (AOD) at 3 km resolution from MODIS Aqua and Terra satellites to ascertain the Spatio-temporal and seasonal distribution of aerosols over Ghana and its major cities. The MODIS AOD at 3 km were validated against ground-based Aerosol Robotic Network (AERONET) AODs to ascertain the suitability of the MODIS 3 km data for air quality application in the region. The contribution of distant aerosols to city aerosol loadings was also assessed with Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) backscatter model. A moderate-high aerosol burden (AODs ~ 0.50) was observed over Ghana with a significant contribution from the pre-monsoon season. City centres of Takoradi and Kumasi showed higher aerosol loads (AODs ~ 0.80) than Accra and Tamale. The HYSPLIT model showed that distant or transported aerosol sources to the city centres were of both marine and land generated origins. Linear regression analysis between MODIS AOD and AERONET AOD showed a reasonably good correlation of ~ 0.60 for Aqua and Terra. From the validation analysis, both Aqua and Terra satellites can be used for air quality monitoring over Ghana; however, more ground research must be conducted to ascertain better aerosol model assumptions for the region.

Aerosols↗

TLife-LSTM: Forecasting Future COVID-19 Progression with Topological Signatures of Atmospheric Conditions

Understanding the impact of atmospheric conditions on SARS-CoV2 is critical to model COVID-19 dynamics and sheds a light on the future spread around the world. Furthermore, geographic distri- butions of expected clinical severity of COVID-19 may be closely linked to prior history of respiratory diseases and changes in humidity, tem- perature, and air quality. In this context, we postulate that by tracking topological features of atmospheric conditions over time, we can provide a quanti?able structural distribution of atmospheric changes that are likely to be related to COVID-19 dynamics. As such, we apply the machinery of persistence homology on time series of graphs to extract topological signatures and to follow geographical changes in relative humidity and temperature. We develop an integrative machine learning framework named Topological Lifespan LSTM (TLife-LSTM) and test its predictive capabilities on forecasting the dynamics of SARS-CoV2 cases. We validate our framework using the number of con?rmed cases and hospitalization rates recorded in the states of Washington and California in the USA. Our results demonstrate the predictive potential of TLife-LSTM in forecasting the dynamics of COVID-19 and modeling its complex spatio-temporal spread dynamics.

Gel, Yulia R.↗

Hybrid phenology matching model for robust crop phenological retrieval

Crop phenology regulates seasonal agroecosystem carbon, water, and energy exchanges, and is a key component in empirical and process-based crop models for simulating biogeochemical cycles of farmlands, assessing gross and net primary production, and forecasting the crop yield. The advances in phenology matching models provide a feasible means to monitor crop phenological progress using remote sensing observations, with a priori information of reference shapes and reference phenological transition dates. Yet the underlying geometrical scaling assumption of models, together with the challenge in defining phenological references, hinders the applicability of phenology matching in crop phenological studies. The objective of this study is to develop a novel hybrid phenology matching model to robustly retrieve a diverse spectrum of crop phenological stages using satellite time series. The devised hybrid model leverages the complementary strengths of phenometric extraction methods and phenology matching models. It relaxes the geometrical scaling assumption and can characterize key phenological stages of crop cycles, ranging from farming practice-relevant stages (e.g., planted and harvested) to crop development stages (e.g., emerged and mature). To systematically evaluate the influence of phenological references on phenology matching, four representative phenological reference scenarios under varying levels of phenological calibrations in terms of time and space are further designed with publicly accessible phenological information. The results indicate that the hybrid phenology matching model can achieve high accuracies for estimating corn and soybean phenological growth stages in Illinois, particularly with the year- and region-adjusted phenological reference (R-squared higher than 0.9 and RMSE less than 5 days for most phenological stages). The inter-annual and regional phenological patterns characterized by the hybrid model correspond well with those in the crop progress reports (CPRs) from the USDA National Agricultural Statistics Service (NASS). Compared to the benchmark phenology matching model, the hybrid model is more robust to the decreasing levels of phenological reference calibrations, and is particularly advantageous in retrieving crop early phenological stages (e.g., planted and emerged stages) when the phenological reference information is limited. This innovative hybrid phenology matching model, together with CPR-enabled phenological reference calibrations, holds 3 considerable promise in revealing spatio-temporal patterns of crop phenology over extended geographical regions.

Phenology↗

LUVOIR-ECLIPS Closed-Loop Adaptive Optics Performance and Contrast Predictions

One of the primary science goals of the Large UV/Optical/Infrared Surveyor (LUVOIR) mission concept is to detect and characterize Earth-like exoplanets orbiting nearby stars with direct imaging. The success of its coronagraph instrument ECLIPS (Extreme Coronagraph for Living Planetary Systems) depends on the ability to stabilize the wavefront from a large segmented mirror such that optical path differences are limited to tens of picometers RMS during an exposure time of a few hours. In order to relax the constraints on the mechanical stability, ECLIPS will be equipped with a wavefront sensing and control (WS\&C) architecture to correct wavefront errors up to temporal frequencies $\gtrsim$1 Hz. These errors may be dominated by spacecraft structural dynamics exciting vibrations at the segmented primary mirror. In this work, we present detailed simulations of the WS\&C system within the ECLIPS instrument and the resulting contrast performance. This study assumes wavefront aberrations based on a finite element model of a simulated telescope with spacecraft structural dynamics. Wavefront residuals are then computed according to a model of the adaptive optics system that includes numerical propagation to simulate a realistic wavefront sensor and an analytical model of the temporal performance. An end-to-end numerical propagation model of ECLIPS is then used to estimate the residual starlight intensity distribution at the science detector. We show that the contrast performance depends strongly on the target star magnitude and the spatio-temporal distribution of wavefront errors from the telescope. In cases with significant vibration, we advocate for the use of laser metrology to mitigate high temporal frequency wavefront errors and increase the mission yield.

Sirbu, Dan↗

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

Discovering Physically Meaningful Structures from Climate Extreme Data

The original proposal described an interdisciplinary team spanning UC San Diego (lead), Columbia University, and UC Irvine, with Columbia investigators including Pierre Gentine, Elias Bareinboim, and Marcus van Lier-Walqui. The proposal further specified a leadership structure in which Columbia co-investigators contributed across the three aims, with Co-PI Gentine serving as a point of contact with science teams and with responsibilities distributed across aims.

42 ENGINEERING↗

Space-Time Data fusion for Remote Sensing Applications

NASA has been collecting massive amounts of remote sensing data about Earth's systems for more than a decade. Missions are selected to be complementary in quantities measured, retrieval techniques, and sampling characteristics, so these datasets are highly synergistic. To fully exploit this, a rigorous methodology for combining data with heterogeneous sampling characteristics is required. For scientific purposes, the methodology must also provide quantitative measures of uncertainty that propagate input-data uncertainty appropriately. We view this as a statistical inference problem. The true but notdirectly- observed quantities form a vector-valued field continuous in space and time. Our goal is to infer those true values or some function of them, and provide to uncertainty quantification for those inferences. We use a spatiotemporal statistical model that relates the unobserved quantities of interest at point-level to the spatially aggregated, observed data. We describe and illustrate our method using CO2 data from two NASA data sets.

carbon dioxide↗