Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatiotemporal model”

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 37 records · Page 2

Optimal displacement in apparent motion and quadrature models of motion sensing

A grating appears to move if it is displaced by some amount between two brief presentations, or between multiple successive presentations. A number of recent experiments have examined the influence of displacement size upon either the sensitivity to motion, or upon the induced motion aftereffect. Several recent motion models are based upon quadrature filters that respond in opposite quadrants in the spatiotemporal frequency plane. Predictions of the quadrature model are derived for both two-frame and multiframe displays. Quadrature models generally predict an optimal displacement of 1/4 cycle for two-frame displays, but in the multiframe case the prediction depends entirely on the frame rate.

Watson, Andrew B.↗

A Spatiotemporal Indexing Approach for Efficient Processing of Big Array-Based Climate Data with MapReduce

Climate observations and model simulations are producing vast amounts of array-based spatiotemporal data. Efficient processing of these data is essential for assessing global challenges such as climate change, natural disasters, and diseases. This is challenging not only because of the large data volume, but also because of the intrinsic high-dimensional nature of geoscience data. To tackle this challenge, we propose a spatiotemporal indexing approach to efficiently manage and process big climate data with MapReduce in a highly scalable environment. Using this approach, big climate data are directly stored in a Hadoop Distributed File System in its original, native file format. A spatiotemporal index is built to bridge the logical array-based data model and the physical data layout, which enables fast data retrieval when performing spatiotemporal queries. Based on the index, a data-partitioning algorithm is applied to enable MapReduce to achieve high data locality, as well as balancing the workload. The proposed indexing approach is evaluated using the National Aeronautics and Space Administration (NASA) Modern-Era Retrospective Analysis for Research and Applications (MERRA) climate reanalysis dataset. The experimental results show that the index can significantly accelerate querying and processing (10 speedup compared to the baseline test using the same computing cluster), while keeping the index-to-data ratio small (0.0328). The applicability of the indexing approach is demonstrated by a climate anomaly detection deployed on a NASA Hadoop cluster. This approach is also able to support efficient processing of general array-based spatiotemporal data in various geoscience domains without special configuration on a Hadoop cluster.

big data↗

Probabilistic Machine Learning Estimation of Ocean Mixed Layer Depth from Dense Satellite and Sparse In-Situ Observations

The ocean mixed layer plays an important role in the coupling between the upper ocean and atmosphere across a wide range of time scales. Estimation of the variability of the ocean mixed layer is therefore important for atmosphere-ocean prediction and analysis. The increasing coverage of in situ Argo profile data allows for an increasingly accurate analysis of the mixed layer depth (MLD) variability associated with deviations from the seasonal climatology. However, sampling rates are not sufficient to fully resolve subseasonal (<90 day) MLD variability. Yet, many multivariate observations-based analyses include implicit modeled subseasonal MLD variability. One analysis method is optimal interpolation of in situ data, but the interior analysis can be improved by leveraging surface data with regression or variational approaches. Here, we demonstrate how machine learning methods and satellite sea surface temperature, salinity, and height facilitate MLD estimation in a pilot study of two regions: the mid-latitude southern Indian and the eastern equatorial Pacific Oceans. We construct multiple machine learning architectures to produce weekly 1/2° gridded MLD anomaly fields (relative to a monthly climatology) with uncertainty estimates. We test multiple traditional and probabilistic machine learning techniques to compare both accuracy and probabilistic calibration. We validate our methodology by applying it to ocean model simulations. We find that incorporating sea surface data through a machine learning model improves the performance of spatiotemporal MLD variability estimation compared to optimal interpolation of Argo observations alone. These preliminary results are a promising first step for the application of machine learning to MLD prediction.

Machine Learning↗

A GeoNEX-Based High-Spatiotemporal-Resolution Product of Land Surface Downward Shortwave Radiation and Photosynthetically Active Radiation

Surface downward shortwave radiation (DSR) and photosynthetically active radiation (PAR) play critical roles in the Earth’s surface processes. As the main inputs of various ecological, hydrological, carbon, and solar photovoltaic models, increasing requirements for high-spatiotemporal-resolution DSR and PAR es- timation with high accuracy have been observed in recent years. However, few existing products satisfy all of these requirements. This study employed a well-established physical-based lookup table (LUT) approach to the GeoNEX gridded top-of-atmosphere bidirectional reflectance factor data acquired by the Advanced Hi- mawari Imager (AHI) and Advanced Baseline Imager (ABI) sensors. It produced a data product of DSR and PAR over both AHI and ABI coverage at an hourly temporal step with a 1 km spatial resolution. GeoNEX DSR data were validated over 63 stations, and GeoNEX PAR data were validated over 27 stations. The vali- dation showed that the new GeoNEX DSR and PAR products have accuracy higher than other existing prod- ucts, with root mean square error (RMSE) of hourly GeoNEX DSR achieving 74.3Wm−2 (18.0%), daily DSR estimation achieving 18.0 W m−2 (9.2 %), hourly GeoNEX PAR achieving 34.9 W m−2 (19.6 %), and daily PAR achieving 9.5 W m−2 (10.5 %). The study also demonstrated the application of the high-spatiotemporal- resolution GeoNEX DSR product in investigating the spatial heterogeneity and temporal variability of surface solar radiation. The data product can be freely accessed through the NASA Advanced Supercomputing Division GeoNEX data portal: https://data.nas.nasa.gov/geonex/geonexdata/GOES16/GEONEX-L2/DSR-PAR/ (last ac- cess: 12 March 2023) and https://data.nas.nasa.gov/geonex/geonexdata/HIMAWARI8/GEONEX-L2/DSR-PAR/ (last access: 12 March 2023) (https://doi.org/10.5281/zenodo.7023863; Wang and Li, 2022).

Surface downward shortwave radiation (DSR)↗

Going With the Flow: Using a Particle Dispersion Model to Reveal Links Between Ocean Ecosystems and Marine Aerosols

Marine biogenic particle contributions to atmospheric aerosol concentrations are not well understood though they are important for determining cloud optical and cloud nucleating properties. Here we examine the relationship between marine aerosol measurements with satellite and model fields of ocean biology and meteorological variables during the North Atlantic Aerosols and Marine Ecosystems Study (NAAMES). NAAMES consisted of four field campaigns between November 2015 and April 2018 that aligned with the four major phases of the annual phytoplankton bloom cycle. The FLEXPART Lagrangian particle dispersion model is used to connect these variables spatiotemporally to ship-based aerosol. We find that correlations between some aerosol measurements with satellite measured and modelled variables increase with increasing trajectory length. This indicates biological and meteorological processes over the air mass history are influential to measured particle properties. For example, results show organic aerosol mass is driven by biogenic volatile organic compound (VOC) emissions from photosynthesis by phytoplankton stocks during advection into the region. In contrast, primary marine aerosol (PMA) concentrations showed better correlations to variables weighted with shorter air mass residence times, which reflects their localized origin as primary emissions. Evidence also suggest variability in PMA concentrations likely play an important role in VOC emission removal. The removal of VOCs lessen their impact on particle number and mass concentrations. The sensitivity of marine aerosols to the main ecosystem-atmosphere drivers and its dependence on air mass history are to be further discussed.

Kevin J Sanchez↗

Spatiotemporal Characteristics of the Association between AOD and PM over the California Central Valley

Many air pollution health effects studies rely on exposure estimates of particulate matter (PM) concentrations derived from remote sensing observations of aerosol optical depth (AOD). Simple but robust calibration models between AOD and PM are therefore important for generating reliable PM exposures. We conduct an in-depth examination of the spatial and temporal characteristics of the AOD-PM2.5 relationship by leveraging data from the Distributed Regional Aerosol Gridded Observation Networks (DRAGON) field campaign where eight NASA Aerosol Robotic Network (AERONET) sites were co-located with EPA Air Quality System (AQS) monitoring sites in California’s Central Valley from November 2012 to April 2013. With this spatiotemporally rich data we found that linear calibration models (R(exp 2) = 0.35, RMSE = 10.38 µg/cu.m) were significantly improved when spatial (R(exp 2) = 0.45, RMSE = 9.54 µg/cu.m), temporal (R(exp 2) = 0.62, RMSE = 8.30 µg/cu.m), and spatiotemporal (R(exp 2) = 0.65, RMSE = 7.58 µg/cu.m) functions were included. As a use-case we applied the best spatiotemporal model to convert space-borne MultiAngle Imaging Spectroradiometer (MISR) AOD observations to predict PM2.5 over the region (R(exp 2) = 0.60, RMSE = 8.42 µg/cu.m). Our results imply that simple AERONET AOD-PM2.5 calibrations are robust and can be reliably applied to space-borne AOD observations, resulting in PM2.5 prediction surfaces for use in downstream applications.

Meytar Sorek-Hamer↗

The Hourly GHI, DHI and DNI: Intercomparison of the CERES-Based Data and the NSRDB Data through Comparison with the Ground-Based BSRN Data

The NASA Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg(Ed4.1) provides hourly global horizontal irradiance (GHI), diffuse horizontal irradiance (DHI) and direct horizontal irradiance (DirHI) at 1-degree latitude by 1-degree longitude resolution, and the time span is from March 2000 to near present. While the GHI agrees well with the ground-based Baseline Surface Radiation Network (BSRN) data, the DHI and DirHI are, respectively, positively and negatively biased against the BSRN data, and so is the direct normal irradiance (DNI) derived simply by dividing the DirHI by the cosine of the hourly solar zenith angle (SZA). We found that the biases tend to be well-defined functions of cos(SZA) and the cloud fraction (CLFR). We thus performed a bias-correction on the hourly DHI and DNI in the latitude-cos(SZA)-CLFR phase space, and therefrom, we derived the hourly global tilted irradiance (GTI) at a number of tilt angles using the isotropic diffuse irradiance model, and the results agree well with the GTI derived from the original BSRN near-instantaneous records. Meanwhile, the National Renewable Energy Laboratory (NREL) has produced the National Solar Radiation Database (NSRDB) with finer spatiotemporal resolution using its Physical Solar Model (PSM) and surface-based measured data, model results and satellite-based data as inputs. Now the NSRDB data covers not only the United States and its neighboring regions, but Europe, Africa, Asia and Oceania. The available temporal resolutions are 5, 10, 15, 30 and 60 minutes and the spatial resolutions are 2, 4 and 10 km, depending on the selected region. In this presentation, we compare both the CERES-based and NSRDB GHI, DHI and DNI with their BSRN counterparts and show how finer spatial resolution may give us an advantage in data accuracy and usability.

Taiping Zhang↗

Using remotely sensed information to support landslide hazard and exposure assessment throughout the disaster lifecycle

The global coverage and temporal frequency that satellites provide offers a unique opportunity to estimate landslide hazard and exposure throughout the disaster lifecycle, from pre-event planning and forecasting to post-event mapping and impact assessment, and finally to recovery and mitigation. The relevance of satellite-derived data and model products is largely contingent on the spatiotemporal sampling, the hazard characteristics, and the needs from the research or applications community. This work presents an advanced Landslide Hazard Assessment for Situational Awareness (LHASA) framework that brings together satellite and model products with new machine learning techniques and global inventory data to better model landslide hazard and exposure. We present several new ways to map, model, and assess landslide hazard using a range of satellite data. Two new thrusts of this work are to better account for the exacerbating impacts of fires and to provide a multi-day forecast of potential hazard. Together, these additional components blend information from a suite of satellite and model sources to improve early warning of potentially hazardous areas, identify landslide occurrence and impacts in near real-time, and better characterize the spatiotemporal patterns of landslide hazard and exposure more broadly for future awareness and planning. This suite of tools and products is open to the public and provides information to better assess the potential occurrence and impacts of landslides within different regions of the world. This presentation explores both the architecture behind this framework and examples of how the model components and products have been used by different stakeholders around the world.

Thomas Stanley↗

Ion Dynamics Model for Collisionless Radio Frequency Sheaths

Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.

Bose, Deepak↗

Land Surface Process and Air Quality Research and Applications at MSFC

This viewgraph presentation provides an overview of land surface process and air quality research at MSFC including atmospheric modeling and ongoing research whose objective is to undertake a comprehensive spatiotemporal analysis of the effects of accurate land surface characterization on atmospheric modeling results, and public health applications. Land use maps as well as 10 meter air temperature, surface wind, PBL mean difference heights, NOx, ozone, and O3+NO2 plots as well as spatial growth model outputs are included. Emissions and general air quality modeling are also discussed.

Quattrochi, Dale↗

Projection-Based Reduced Order Modeling for Spacecraft Thermal Analysis

This paper presents a mathematically rigorous, subspace projection-based reduced order modeling (ROM) methodology and an integrated framework to automatically generate reduced order models for spacecraft thermal analysis. Two key steps in the reduced order modeling procedure are described: (1) the acquisition of a full-scale spacecraft model in the ordinary differential equation (ODE) and differential algebraic equation (DAE) form to resolve its dynamic thermal behavior; and (2) the ROM to markedly reduce the dimension of the full-scale model. Specifically, proper orthogonal decomposition (POD) in conjunction with discrete empirical interpolation method (DEIM) and trajectory piece-wise linear (TPWL) methods are developed to address the strong nonlinear thermal effects due to coupled conductive and radiative heat transfer in the spacecraft environment. Case studies using NASA-relevant satellite models are undertaken to verify the capability and to assess the computational performance of the ROM technique in terms of speed-up and error relative to the full-scale model. ROM exhibits excellent agreement in spatiotemporal thermal profiles (<0.5% relative error in pertinent time scales) along with salient computational acceleration (up to two orders of magnitude speed-up) over the full-scale analysis. These findings establish the feasibility of ROM to perform rational and computationally affordable thermal analysis, develop reliable thermal control strategies for spacecraft, and greatly reduce the development cycle times and costs.

0000↗

A Framework for Deep Learning Emulation of Numerical Models With a Case Study in Satellite Remote Sensing

Numerical models based on physics represent the state of the art in Earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning (DL), across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in Earth systems. A difficult test for DL-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A DL emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here, we examine, with a case study in satellite-based remote sensing, the hypothesis that DL approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the DL emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of DL and the growing desire for higher resolution simulations in the Earth sciences.

Bayesian Deep Learning↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

Spatiotemporal Analysis of High-Resolution Pressure-Sensitive Paint Measurements on a Launch Vehicle Wind Tunnel Model

An accurate understanding of unsteady aerodynamic loads is critical for launch vehicle structural design and for development of the vibro-acoustic environment for payload components. Analysts commonly use point-wise surface pressure measurements from wind tunnel testing to build statistical models of the in-flight turbulent boundary layer flow structures, often extrapolating over large swaths of the surface due to measurement sparsity. We present novel results from a transonic wind tunnel test using fast-response pressure-sensitive paint (PSP) and show how it can be used to visualize and quantify localized turbulent boundary layer wall pressure statistics with unprecedented spatial resolution. The wind tunnel model is a scaled model of the Space Launch System (SLS) rocket Block 1B Crew configuration and analysis focused on the compression corner between the Service Module fairing and the Universal Stage Adapter.

Corcos Model↗

Quantitative Measures of Immersion in Cloud and the Biogeography of Cloud Forests

Sites described as tropical montane cloud forests differ greatly, in part because observers tend to differ in their opinion as to what constitutes frequent and prolonged immersion in cloud. This definitional difficulty interferes with hydrologic analyses, assessments of environmental impacts on ecosystems, and biogeographical analyses of cloud forest communities and species. Quantitative measurements of cloud immersion can be obtained on site, but the observations are necessarily spatially limited, although well-placed observers can examine 10 50 km of a mountain range under rainless conditions. Regional analyses, however, require observations at a broader scale. This chapter discusses remote sensing and modeling approaches that can provide quantitative measures of the spatiotemporal patterns of cloud cover and cloud immersion in tropical mountain ranges. These approaches integrate remote sensing tools of various spatial resolutions and frequencies of observation, digital elevation models, regional atmospheric models, and ground-based observations to provide measures of cloud cover, cloud base height, and the intersection of cloud and terrain. This combined approach was applied to the Monteverde region of northern Costa Rica to illustrate how the proportion of time the forest is immersed in cloud may vary spatially and temporally. The observed spatial variation was largely due to patterns of airflow over the mountains. The temporal variation reflected the diurnal rise and fall of the orographic cloud base, which was influenced in turn by synoptic weather conditions, the seasonal movement of the Intertropical Convergence Zone and the north-easterly trade winds. Knowledge of the proportion of the time that sites are immersed in clouds should facilitate ecological comparisons and biogeographical analyses, as well as land use planning and hydrologic assessments in areas where intensive on-site work is not feasible.

Lawton, R. O.↗

Soil Moisture Estimation in South Asia via Assimilation of SMAP Retrievals

A soil moisture retrieval assimilation framework is implemented across South Asia in an attempt to improve regional soil moisture estimation as well as to provide a consistent regional soil moisture dataset. This study aims to improve the spatiotemporal variability of soil moisture estimates by assimilating Soil Moisture Active Passive (SMAP) near-surface soil moisture retrievals into a land surface model. The Noah-MP (v4.0.1) land surface model is run within the NASA Land Information System software framework to model regional land surface processes. NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA2) and Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals (IMERG) provide the meteorological boundary conditions to the land surface model. Assimilation is carried out using both cumulative distribution function (CDF)-corrected (DA-CDF) and uncorrected SMAP retrievals (DA-NoCDF). CDF matching is applied to correct the statistical moments of the SMAP soil moisture retrieval relative to the land surface model. Comparison of assimilated and model-only soil moisture estimates with publicly available in situ measurements highlights the relative improvement in soil moisture estimates by assimilating SMAP retrievals. Across the Tibetan Plateau, DA-NoCDF reduced the mean bias and RMSE by 8.4 % and 9.4 %, even though assimilation only occurred during less than 10 % of the study period due to frozen (or partially frozen) soil conditions. The best goodness-of-fit statistics were achieved for the IMERG DA-NoCDF soil moisture experiment. The general lack of publicly available in situ measurements across irrigated areas limited a domain-wide direct model validation. However, comparison with regional irrigation patterns suggested correction of biases associated with an unmodeled hydrologic phenomenon (i.e., anthropogenic influence via irrigation) as a result of SMAP soil moisture retrieval assimilation. The greatest sensitivity to assimilation was observed in cropland areas. Improvements in soil moisture potentially translate into improved spatiotemporal patterns of modeled evapotranspiration, although limited influence from soil moisture assimilation was observed on modeled processes within the carbon cycle such as gross primary production. Improvement in fine-scale modeled estimates by assimilating coarse-scale retrievals highlights the potential of this approach for soil moisture estimation over data-scarce regions.

Jawairia Ahmad↗

Modeling Landslide Susceptibility Over Large Regions with Fuzzy Overlay

Landslide susceptibility mapping is most effective if detailed surface and subsurface information can be combined with authoritative landslide catalogs or a deep understanding of local conditions. However, these types of homogeneous input data and catalogs are frequently not available over large areas. In this study, we model landslide susceptibility in Central America and the Caribbean islands by combining three globally available datasets and one regional dataset with fuzzy overlay. This primarily heuristic model provides the flexibility to test a range of different contributing variables and the capability to compare landslide inventories within the model framework that vary greatly in their size, spatiotemporal scope, and collection methods. We create a regional susceptibility map and evaluate its performance using receiver operating characteristics for both continuous and binned susceptibility values. This susceptibility map forms the basis for a near-real-time landslide hazard assessment system that couples susceptibility with rainfall and soil moisture triggers to estimate potential landslide activity at a regional scale. The application of this susceptibility model at the regional scale provides a foundation for transferring the methodology to other geographic areas.

Landslide↗

Characterizing Martian Volcanic Provinces’ Magmatic Evolution and Chemistry through Equations of State Modeling Initial Study

Here we discuss a novel interdisciplinary approach to investigating igneous compositions of large volcanic provinces on Mars using remote sensing data sets (Mars Odyssey Gamma Ray and neutron Spectrometer suite--GRS, and gravity) to inform petrologic and thermoelastic modeling. Martian volcanic provinces, starting from Noachian to Amazonian age, including Elysium (EVP), are locations of great geologic interest which have been active over long time scales from Hesperian to Amazonian. Regional scale change in eruptive processes are poorly understood. Compared to large igneous provinces on Earth, the martian volcanic activity has persisted for orders of magnitude longer. Therefore, changes in mantle chemistry, pressure, and temperature during lithospheric cooling are expected to produce significant changes in the conditions of magma production, storage and ascent, affecting the degree of fractional crystallization and crustal contamination. In this study, we perform a detailed modeling to test the hypothesis that compositional variability within volcanic provinces resulted from spatiotemporal changes in the depth of magma formation and present initial results. Our methods include constraining the pressure and temperature conditions of EVP as a case study of geologically recent magmatic evolution on Mars using GRS informed surface chemistry constraints and pMELTS modeling. Second, we place constraints on the density and seismic velocities of the EVP melt through thermoelastic modeling with a local gravity analysis. This analysis is also extended to Noachian aged volcanic provinces on Mars initial results estimates mantle pressure for each sub-region is 16 kbar pressure, whereas, degree of partial melting is low and varies between 10 to 12. This study aims to validate our theoretical model and develop a perspective view of spatiotemporal changes at the interior of Mars throughout the time.

Martian science, Mars, geology, geophysics↗