Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean 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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

A Realization of Bias Correction Method in the GMAO Coupled System

Over the past several decades, a tremendous effort has been made to improve model performance in the simulation of the climate system. The cold or warm sea surface temperature (SST) bias in the tropics is still a problem common to most coupled ocean atmosphere general circulation models (CGCMs). The precipitation biases in CGCMs are also accompanied by SST and surface wind biases. The deficiencies and biases over the equatorial oceans through their influence on the Walker circulation likely contribute the precipitation biases over land surfaces. In this study, we introduce an approach in the CGCM modeling to correct model biases. This approach utilizes the history of the model's short-term forecasting errors and their seasonal dependence to modify model's tendency term and to minimize its climate drift. The study shows that such an approach removes most of model climate biases. A number of other aspects of the model simulation (e.g. extratropical transient activities) are also improved considerably due to the imposed pre-processed initial 3-hour model drift corrections. Because many regional biases in the GEOS-5 CGCM are common amongst other current models, our approaches and findings are applicable to these other models as well.

Chang, Yehui↗
DOE OSTI2020

The Impact of Parameterized Lateral Mixing on the Antarctic Circumpolar Current in a Coupled Climate Model

This study examines the impact of changing the lateral diffusion coefficient ARedi on the transport of the Antarctic Circumpolar Current (ACC). The lateral diffusion coefficient ARedi is poorly constrained, with values ranging across an order of magnitude in climate models. The ACC is difficult to accurately simulate, and there is a large spread in eastward transport in the Southern Ocean (SO) in these models. This paper examines how much of that spread can be attributed to different eddy parameterization coefficients. A coarse-resolution, fully coupled model suite was run with A Redi = 400, 800, 1200, and 2400 m 2 s -1 . Additionally, two simulations were run with two-dimensional representations of the mixing coefficient based on satellite altimetry. Relative to the 400 m 2 s -1 case, the 2400 m 2 s -1 case exhibits 1) an 11% decrease in average wind stress from 50° to 65°S, 2) a 20% decrease in zonally averaged eastward transport in the SO, and 3) a 14% weaker transport through the Drake Passage. The decrease in transport is well explained by changes in the thermal current shear, largely due to increases in ocean density occurring on the northern side of the ACC. In intermediate waters these increases are associated with changes in the formation of intermediate waters in the North Pacific. We hypothesize that the deep increases are associated with changes in the wind stress curl allowing Antarctic Bottom Water to escape and flow northward.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Tracking Improvement in Simulated Marine Biogeochemistry Between CMIP5 and CMIP6

Purpose of Review The changes or updates in ocean biogeochemistry component have been mapped between CMIP5 and CMIP6 model versions, and an assessment made of how far these have led to improvements in the simulated mean state of marine biogeochemical models within the current generation of Earth system models (ESMs). Recent Findings The representation of marine biogeochemistry has progressed within the current generation of Earth system models. However, it remains difficult to identify which model updates are responsible for a given improvement. In addition, the full potential of marine biogeochemistry in terms of Earth system interactions and climate feedback remains poorly examined in the current generation of Earth system models. Summary Increasing availability of ocean biogeochemical data, as well as an improved understanding of the underlying processes, allows advances in the marine biogeochemical components of the current generation of ESMs. The present study scrutinizes the extent to which marine biogeochemistry components of ESMs have progressed between the 5th and the 6th phases of the Coupled Model Intercomparison Project (CMIP).

Marine Biogeochemistry↗
NASA NTRS2013

Mid-Pliocene Atlantic Meridional Overturning Circulation Not Unlike Modern

In the Pliocene Model Intercomparison Project (PlioMIP), eight state-of-the-art coupled climate models have simulated the mid-Pliocene warm period (mPWP, 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), northward ocean heat transport and ocean stratification simulated with these models. None of the models participating in PlioMIP simulates a strong mid-Pliocene AMOC as suggested by earlier proxy studies. Rather, there is no consistent increase in AMOC maximum among the PlioMIP models. The only consistent change in AMOC is a shoaling of the overturning cell in the Atlantic, and a reduced influence of North Atlantic Deep Water (NADW) at depth in the basin. Furthermore, the simulated mid-Pliocene Atlantic northward heat transport is similar to the pre-industrial. These simulations demonstrate that the reconstructed high-latitude mid-Pliocene warming can not be explained as a direct response to an intensification of AMOC and concomitant increase in northward ocean heat transport by the Atlantic.

Heat transfer↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2022

Observations Indicate That Clouds Amplify Mechanisms of Southern Ocean Heat Uptake

Abstract The Southern Ocean has absorbed most of the excess heat associated with anthropogenic greenhouse gas emissions. Since Southern Ocean observations are sparse in certain regions and seasons, much of our knowledge of ocean heat uptake is based on climate model simulations. However, climate models still inadequately represent some properties of Southern Ocean clouds, and they have not identified the mechanisms by which clouds may affect Southern Ocean heat uptake (SOHU). Here, we use the ERA5 and JRA‐55 reanalyses to assess the influence of clouds and other atmospheric processes on SOHU from 1979 to 2020. We find that years with the highest SOHU between 45° and 65°S are dominated by ocean heat uptake anomalies during winter and spring, but not during summer or fall. Winter and spring cloud cover are up to 7% higher when SOHU is up to 5.5 W/m 2 higher than the climatological seasonal mean, with the largest increases in the South Pacific Ocean. Clouds also contain more liquid water. These changes in cloud properties increase downwelling longwave radiation, amplifying ocean heat uptake. Cloud changes are also concomitant with a more stable lower atmosphere, which suppresses turbulent heat fluxes out of the surface. Overall, we find that SOHU is likely not mediated by enhanced surface shortwave absorption over the observational time period. A better understanding of how atmospheric processes impact ocean heat uptake may help improve our understanding of ocean heat uptake mechanisms in the current generation of climate models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Impact Study of the Assimilation of Surface Sensitive Microwave Radiances in the GEOS

The aim of this study is to investigate the feasibility of assimilating low frequency microwave observations from different satellite Microwave radiometers such as the Advanced Microwave Sounding Unit-A (AMSU-A). These observations are relevant to the description of air temperature, humidity, and surface parameters such as ocean surface temperature. Their assimilation into Goddard Earth Observing System (GEOS) modeling and assimilation system helps better constrain models in regions where very few observations are assimilated. In recent years, Channels 1–4 and 15 have not been assimilated in GEOS because of their large sensitivities to uncertain surface parameters such as emissivity and skin temperature. Here, the impact of using different channels will be discussed. Also we investigate the validity of the emissivity model in these channels over ocean and sea ice.Radiative transfer calculations necessary for the assimilation of satellite radiances are performed using the Community Radiative Transfer Model (CRTM). The CRTM performs forward radiative transfer calculations of brightness temperature from atmospheric profiles (temperature, moisture, and ozone) and surface characteristics. Additionally, the CRTM calculates the Jacobians of brightness temperature with respect to the input . The Jacobians are for the inversion of these radiances within the Gridpoint Statistical Interpolation (GSI) analysis system. Here, For direct assimilation of the radiance observations, the interface between GSI and CRTM, is slightly modified. Finally, we validate the assimilation results, and discuss the impact of each channels in detail.

Ebrahimi, Hamideh↗
NASA NTRS2000

Modeling Studies of the Effects of Tropical Rainfall on Ocean-Atmosphere Interactions and Oceanic Hydrological Cycle

This research project is a joint effort of UMD/JCESS, NASA/GSFC, NOAA/PMEL and LDEO, with UMD/JCESS being the leader and LDEO being responsible for coupled modeling. Although the project as a whole is for three years, the LDEO part was terminated at the end of the second year because of a shift of focus. The detailed description of the whole project and the effort made at LDEO have been given in the two joint annual reports and the two annual reports from LDEO. Here we only provide a brief summary of what we have done here at LDEO and a list of publications that resulted partly from this project.

Chen, Dake↗
NASA NTRS2008

Incorporating Satellite Time-Series Data into Modeling

In situ time series observations have provided a multi-decadal view of long-term changes in ocean biology. These observations are sufficiently reliable to enable discernment of even relatively small changes, and provide continuous information on a host of variables. Their key drawback is their limited domain. Satellite observations from ocean color sensors do not suffer the drawback of domain, and simultaneously view the global oceans. This attribute lends credence to their use in global and regional model validation and data assimilation. We focus on these applications using the NASA Ocean Biogeochemical Model. The enhancement of the satellite data using data assimilation is featured and the limitation of tongterm satellite data sets is also discussed.

Gregg, Watson↗
DOE OSTI2022

Equation‐Free Surrogate Modeling of Geophysical Flows at the Intersection of Machine Learning and Data Assimilation

Abstract There is a growing interest in developing data‐driven reduced‐order models for atmospheric and oceanic flows that are trained on data obtained either from high‐resolution simulations or satellite observations. The data‐driven models are non‐intrusive in nature and offer significant computational savings compared to large‐scale numerical models. These low‐dimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equation‐free surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an end‐to‐end non‐intrusive reduced‐order modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long short‐term memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for long‐term forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.

Pawar, Suraj↗
NASA NTRS1999

Wind Forcing of the Pacific Ocean Using Scatterometer Wind Data

The long-term objective of this research was an understanding of the wind-forced ocean circulation, particularly for the Pacific Ocean. To determine the ocean's response to the winds, we first needed to generate accurate maps of wind stress. For the ocean's response to wind stress we examined the sea surface height (SSH) both from altimeters and from numerical models for the Pacific Ocean.

Kelly, Kathryn A.↗
NASA NTRS1999

Satellite Altimetry, Ocean Circulation, and Data Assimilation

Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.

Fu, Lee-Lueng↗
NASA NTRS1991

The phytoplankton bloom in the northwestern Arabian Sea during the southwest monsoon of 1979

The present study investigates the biological variability of the northwestern Arabian Sea during the 1979 southwest monsoon by the synthesis of satellite ocean color remote sensing with an analysis of in situ hydrographic and meteorological data sets and the results of wind-driven modeling of upper-ocean circulation. The phytoplankton bloom peaked during August-September, extended from the Oman coast to about 65 deg E, and lagged behind the development of open-sea upwelling by at least 1 mo. The pigment distributions, hydrographic data, and model results all suggest that the boom was driven by spatially distinct upward nutrient fluxes to the euphotic zone forced by the physical processes of coastal upwelling and offshore Ekman pumping. Coastal upwelling was evident from May through September, yielded the most extreme concentrations of phytoplankton biomass, and, along the Arabian coast, was limited to the continental shelf in the promotion of high concentrations of phytoplankton.

Brock, John C.↗
NASA NTRS2021

Development and Validation of an Empirical Ocean Color Algorithm with Uncertainties: A Case Study with the Particulate Backscattering Coefficient

We explored how algorithm (model) and in situ measurement (observation) uncertainties can effectively be incorporated into empirical ocean color model development and assessment. In this study we focused on methods for deriving the particulate backscattering coefficient at 555 nm, b(bp)(555)/(m). We developed a simple empirical algorithm for deriving b(bp)(555) as a function of a remote sensing reflectance line height (LH) metric. Model training was performed using a high-quality bio-optical dataset that contains coincident in situ measurements of the spectral remote sensing reflectances, R(rs)(λ)/(sr), and the spectral particulate backscattering coefficients, b(bp)(λ). The LH metric used is defined as the magnitude of Rrs(555) relative to a linear baseline drawn between R(rs)(490) and R(rs)(670). Using an independent validation dataset, we compared the skill of the LH-based model with two other models. We used contemporary validation metrics, including bias and mean absolute error (MAE), that were corrected for model and observation uncertainties. The results demonstrated that measurement uncertainties do indeed impact contemporary validation metrics such as mean bias and MAE. Zeta-scores and z-tests for overlapping confidence intervals were also explored as potential methods for assessing model skill.

ocean color↗
DOE OSTI2020

Scaling Global Warming Impacts on Ocean Ecosystems: Lessons From a Suite of Earth System Models

An important technique used by climate modelers to isolate the impacts of increasing greenhouse gasses on Earth System processes is to simulate the impact of an abrupt increase in carbon dioxide. The spatial pattern of change provides a “fingerprint” that is generally much larger than natural variability. Insofar as the response to radiative forcing is linear (the impact of quadrupling CO 2 is twice the impact of doubling CO 2 ) this fingerprint can then be used to estimate the impact of historical greenhouse gas forcing. However, the degree to which biogeochemical cycles respond linearly to radiative forcing has rarely been tested. In this paper, we evaluate which ocean biogeochemical fields are likely to respond linearly to changing radiative forcing, which ones do not, and where linearity breaks down. We also demonstrate that the representation of lateral mixing by mesoscale eddies, which varies significantly across climate models, plays an important role in modulating the breakdown of linearity. Globally integrated surface rates of biogeochemical cycling (primary productivity, particulate export) respond in a relatively linear fashion and are only moderately sensitive to mixing. By contrast, the habitability of the interior ocean (as determined by hypoxia and calcite supersaturation) behaves non-linearly and is very sensitive to mixing. This is because the deep ocean, as well as certain regions in the surface ocean, are very sensitive to the magnitude of deep wintertime convection. The cessation of convection under global warming is strongly modulated by the representation of eddy mixing.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Tidal Prediction

Tides are the most predictable of oceanographic phenomena, due both to the simplicity and predictability of the astronomical forcing and to the near linearity of the ocean's dynamical response. In the classical and simplest scenario, tidal prediction is based on harmonic analysis of past measurements at a fixed location. Limits to predictability arise because isolated astronomical spectral lines are broadened into "cusps" of incoherent energy, for example through interactions with non-tidal flows. Tidal prediction at locations without past measurements has historically been a major challenge, but, owing to near-global observations of modern satellite altimeters, the empirical harmonic approach now yields reasonably accurate predictions throughout most of the open ocean. Advances in numerical modeling and data assimilation allow these predictions to be refined (especially in shallow seas where observations remain insufficient to directly constrain tidal wave structure) and extended to include tidal currents. We review recent progress in the development of global and regional-scale tidal prediction capabilities, summarize accuracy of available charts, and briefly consider outstanding issues. Satellite altimetry has also helped unravel the global tidal energy budget and has clarified the role of internal tides as a sink of tidal energy. We summarize these results, and then turn to the challenging problem of predicting internal tides. To the extent that low modes remain coherent with the surface tide, elevations can be directly mapped. We review current efforts in this direction, which are already producing charts with predictive capability. However, internal tides (especially higher modes) interact much more strongly with lower-frequency ocean flows, so a significant fraction of this tidal signal is intermittent and incoherent. We close with a brief review of ongoing efforts to model global tides in combination with wind-forced ocean motions.

Harmonic Constants↗
DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Onset of magma ocean solidification on Mars inferred from Mn-Cr chronometry

The mantle of Mars probably differentiated through the crystallization of a magma ocean during the first tens of million years (Ma) of Solar System evolution. However, the exact timescale of large-scale silicate differentiation of the martian mantle is debated, and in particular, it remains unclear when differentiation commenced. In this paper, we applied the short-lived 53 Mn- 53 Cr system to martian meteorites in order to date the onset of large-scale mantle differentiation on Mars. The new Cr isotope data demonstrate that martian meteorites exhibit no resolvable radiogenic 53 Cr variations, and instead have a uniform +20.2±1.2 (95% conf.) parts-per-million excess in 53 Cr/ 52 Cr relative to the terrestrial mantle. The investigated groups of martian meteorites are lithologically varied and derive from diverse mantle sources that probably had variable Mn/Cr. Hence, the lack of 53 Cr variability among martian meteorites demonstrates that silicate differentiation on Mars occurred after the extinction of 53 Mn. Provided that the sources of the martian meteorites have Mn/Cr variations that are typical of the terrestrial planets, this result implies that the onset of large-scale silicate differentiation must have occurred later than 20±5 Ma after Solar System formation. The onset of silicate differentiation on Mars inferred here is significantly later than time estimates for segregation of the martian core which conservatively occurred within <10 Ma after Solar System formation. Thus, the new Mn-Cr data imply that there was a small, but resolvable, time gap of at least 5 Ma between core formation and magma ocean solidification on Mars. If the age of core segregation is taken at face value, our results imply that the martian magma ocean remained mostly molten over several Ma. This inferred longevity of the magma ocean is inconsistent with thermal models predicting rapid (<1 Ma) solidification of the martian magma ocean. Although there is currently no unique solution to this conundrum, our results can potentially be explained by a protracted history of impact bombardment that delayed differentiation in a shallow magma ocean on Mars, or perhaps more readily, by the presence of an early and dense atmosphere that acted as an insulator and prevented the magma ocean from cooling quickly.

58 GEOSCIENCES↗