Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean Modeling”

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,027 records · Page 57

The Chlorine Isotopic Composition Of Lunar UrKREEP

Since the long standing paradigm of an anhydrous Moon was challenged there has been a renewed focus on investigating volatiles in a variety of lunar samples. Numerous studies have examined the abundances and isotopic compositions of volatiles in lunar apatite, Ca5(PO4)3(F,Cl,OH). In particular, apatite has been used as a tool for assessing the sources of H2O in the lunar interior. However, current models for the Moon's formation have yet to fully account for its thermal evolution in the presence of H2O and other volatiles. For ex-ample, in the context of the lunar magma ocean (LMO) model, it is anticipated that chlorine (and other volatiles) should have been concentrated in the late-stage LMO residual melts (i.e., the dregs enriched in incompatible elements such as K, REEs (Rare Earth Elements), and P, collectively called KREEP, and in its primitive form - urKREEP, given its incompatibility in mafic minerals like olivine and pyroxene, which were the dominant phases that crystallized early in the cumulate pile of the LMO. When compared to chondritic meteorites and terrestrial rocks, lunar samples have exotic chlorine isotope compositions, which are difficult to explain in light of the abundance and isotopic composition of other volatile species, especially H, and the current estimates for chlorine and H2O in the bulk silicate Moon (BSM). In order to better understand the processes involved in giving rise to the heavy chlorine isotope compositions of lunar samples, we have performed a comprehensive in situ high precision study of chlorine isotopes in lunar apatite from a suite of Apollo samples covering a range of geochemical characteristics and petrologic types.

Barnes, J. J.↗

Timing of Crystallisation of the Lunar Magma Ocean Constrained by the Oldest Zircon

The presently favoured concept for the early evolution of the Moon involves consolidation of debris from a giant impact of a Mars sized body with Earth forming a primitive Moon with a thick global layer of melt referred to as the Lunar Magma Ocean1 . It is widely accepted that many significant features observed on the Moon today are the result of crystallisation of this magma ocean. However, controversy exists over the precise timing and duration of the crystallisation process. Resolution of this problem depends on the establishment of precise and robust key crystallisation time points. We report a 4417 6 Myr old zircon in lunar breccia sample 72215,195, which provides a precisely determined younger limit for the solidification of the Lunar Magma Ocean. A model based on these data, together with the age of the Moon forming giant impact, defines an exponential time frame for crystallisation and suggests formation of anorthositic crust after about 80-85% of the magma ocean was solidified. In combination with other zircon ages the 4417 +/- 6 Myr age also suggests that the very small (less than a few per cent) residual portion of the magma ocean continued to solidify during the following 300-500 m.y.

Nemchin, A.↗

Sensitivity of the Atmospheric Response to Warm Pool El Nino Events to Modeled SSTs and Future Climate Forcings

Warm pool El Nino (WPEN) events are characterized by positive sea surface temperature (SST) anomalies in the central equatorial Pacific. Under present-day climate conditions, WPEN events generate poleward propagating wavetrains and enhance midlatitude planetary wave activity, weakening the stratospheric polar vortices. The late 21st century extratropical atmospheric response to WPEN events is investigated using the Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM), version 2. GEOSCCM simulations are forced by projected late 21st century concentrations of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) and by SSTs and sea ice concentrations from an existing ocean-atmosphere simulation. Despite known ocean-atmosphere model biases, the prescribed SST fields represent a best estimate of the structure of late 21st century WPEN events. The future Arctic vortex response is qualitatively similar to that observed in recent decades but is weaker in late winter. This response reflects the weaker SST forcing in the Nino 3.4 region and subsequently weaker Northern Hemisphere tropospheric teleconnections. The Antarctic stratosphere does not respond to WPEN events in a future climate, reflecting a change in tropospheric teleconnections: The meridional wavetrain weakens while a more zonal wavetrain originates near Australia. Sensitivity simulations show that a strong poleward wavetrain response to WPEN requires a strengthening and southeastward extension of the South Pacific Convergence Zone; this feature is not captured by the late 21st century modeled SSTs. Expected future increases in GHGs and decreases in ODSs do not affect the polar stratospheric responses to WPEN.

stratosphere↗

The active role of the ocean in the temporal evolution of climate sensitivity

Transient climate sensitivity has been shown to be influenced by the changing pattern of SST and Ocean Heat Uptake (OHU), which in turn have been attributed to ocean circulation changes. A set of novel experiments are performed to isolate the active role of the ocean by comparing a fully-coupled CO2-quadrupling CESM simulation against a partially-coupled one, where the effect of the ocean circulation change and its impact on surface fluxes are disabled. The active OHU is the main factor responsible for the reduced transient climate sensitivity and weaker surface warming response in the fully-coupled simulation. The passive OHU excites qualitatively similar feedbacks to CO2 forcing in a slab ocean model configuration due to the similar SST spatial pattern response in both experiments. Additionally, the non-unitary forcing efficacy of the active OHU (1.7) explains the very different net feedback parameters in the fully-and partially-coupled responses.

climate, ocean heat uptake, OHU, ocean, SST↗

Southern Ocean Sulfate Aerosol Sources Quantified From Sulfur Isotopes in Antarctic Ice Cores

The Southern Ocean has emerged as a key region for constraining aerosol-climate interactions due to its relatively low anthropogenic influence. Sulfate is an important aerosol over the Southern Ocean, and models suggest dimethyl sulfide (DMS) is the largest source of sulfate during summer. However, sulfur isotopes of sulfate (δ 34 S(SO 4 2− )) in Antarctic ice cores suggest a significant contribution from a previously unexplained non-DMS source. Here we show that the fractional contribution from passive volcanic degassing (f volc ) explains observed δ 34 S(SO 4 2− ), and that a global chemical transport model underestimates f volc across Antarctica. Underestimated f volc implies that the model mischaracterizes sulfate sources in this important region. The discrepancy between observed and modeled sulfate sources can be reconciled by increasing passive volcanic sulfur degassing emissions and decreasing DMS emissions. Our results imply that current biases in emissions inventories could bias assessments of aerosol-cloud interactions in the Southern Ocean region and globally.

Jongebloed, Ursula A. [Univ. of Washington, Seattl↗

Lagrangian Timescales of Southern Ocean Upwelling in a Hierarchy of Model Resolutions

Here, we study upwelling pathways and timescales of Circumpolar Deep Water (CDW) in a hierarchy of models using a Lagrangian particle tracking method. Lagrangian timescales of CDW upwelling decrease from 87 years to 31 years to 17 years as the ocean resolution is refined from 1° to 0.25° to 0.1°. We attribute some of the differences in timescale to the strength of the eddy fields, as demonstrated by temporally degrading high-resolution model velocity fields. Consistent with the timescale dependence, we find that an average Lagrangian particle completes 3.2 circumpolar loops in the 1° model in comparison to 0.9 loops in the 0.1° model. In conclusion, these differences suggest that advective timescales and thus interbasin merging of upwelling CDW may be overestimated by coarse-resolution models, potentially affecting the skill of centennial scale climate change projections.

58 GEOSCIENCES↗

Soil Moisture Active/Passive (SMAP) L-Band Microwave Radiometer Post-Launch Calibration Upgrade

The Soil Moisture Active/Passive (SMAP) microwave radiometer is a fully polarimetric L-band radiometer flown on theSMAP satellite in a 6 AM/6 PM sun-synchronous orbit at 685 km altitude. After the SMAP L1B_TB data product version 3 was released in 2016, the radiometer has been undergoing further calibration and validation with the goal of reducing both the bias in the cold-sky measurements and calibration drift in the global ocean measurements experienced during eclipse seasons in data product version 3. The post-launch calibration algorithm has been upgraded by using new estimates of the reflector emissivity as well as using multiple scenes to calibrate the radiometer internal reference sources and antenna gain simultaneously. In addition, a correction offset is applied to the ocean roughness model for horizontal polarization based on nadir observations. Test and validation results show that the goal is achieved (e.g., biases are removed and the calibration stability achieved for data release version 4 is 0.1 K(rms) over both the global ocean and CS).

Calibratio↗

Mixed-phase cloud physics and Southern Ocean cloud feedback in climate models

Increasing optical depth poleward of 45° is a robust response to warming in global climate models. Much of this cloud optical depth increase has been hypothesized to be due to transitions from ice-dominated to liquid-dominated mixed-phase cloud. In this study, the importance of liquid-ice partitioning for the optical depth feedback is quantified for 19 Coupled Model Intercomparison Project Phase 5 models. All models show a monotonic partitioning of ice and liquid as a function of temperature, but the temperature at which ice and liquid are equally mixed (the glaciation temperature) varies by as much as 40 K across models. Models that have a higher glaciation temperature are found to have a smaller climatological liquid water path (LWP) and condensed water path and experience a larger increase in LWP as the climate warms. The ice-liquid partitioning curve of each model may be used to calculate the response of LWP to warming. It is found that the repartitioning between ice and liquid in a warming climate contributes at least 20% to 80% of the increase in LWP as the climate warms, depending on model. Intermodel differences in the climatological partitioning between ice and liquid are estimated to contribute at least 20% to the intermodel spread in the high-latitude LWP response in the mixed-phase region poleward of 45°S. As a result, it is hypothesized that a more thorough evaluation and constraint of global climate model mixed-phase cloud parameterizations and validation of the total condensate and ice-liquid apportionment against observations will yield a substantial reduction in model uncertainty in the high-latitude cloud response to warming.

58 GEOSCIENCES↗

Evaluation of Preindustrial to Present-day Black Carbon and its Albedo Forcing from Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

As a part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against the observations including 12 ice core records, a long-term surface mass concentrations and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using the NCAR Community Land and Sea-Ice model 4 with prescribed meteorology from 1996-2000, which includes the SNICAR BC-snow model. We evaluated the vertical profile of BC snow concentrations from these offline simulations to using recent BC snowpack measurements. Despite using the same BC emissions, global BC burden differs by approximately a factor of 3 among models due to the differences in aerosol removal parameterizations and simulated meteorology among models; 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However,models agree well on 2.5~3 times increase in the global BC burden from preindustrial to present-day, which matches with the 2.5 times increase in BC emissions. We find a large model diversity at both NH and SH high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC mass concentrations well in Europe and North America except at Jungfrauch and Ispra. However, the models fail to capture the Arctic BC seasonality due tosevere underestimations during winter and spring. Compared to recent snowpack measurements, the simulated vertically resolved BC snow concentrations are, on average, within a factor of 2-3 of observations except for Greenland and Arctic Ocean. However, model and observation differ widely due to missing interannual variations in emissions and possibly due to the choice of the prescribed meteorology period (i.e., 1996-2000).

ACCMIP↗

On isostatic geoid anomalies

In regions of slowly varying lateral density changes, the gravity and geoid anomalies may be expressed as power series expansions in topography. Geoid anomalies in isostatically compensated regions can be directly related to the local dipole moment of the density-depth distribution. This relationship is used to obtain theoretical geoid anomalies for different models of isostatic compensation. The classical Pratt and Airy models give geoid height-elevation relationships differing in functional form but predicting geoid anomalies of comparable magnitude. The thermal cooling model explaining ocean floor subsidence away from mid-ocean ridges predicts a linear age-geoid height relationship of 0.16 m/m.y. Geos 3 altimetry profiles were examined to test these theoretical relationships. A profile over the mid-Atlantic ridge is closely matched by the geoid curve derived from the thermal cooling model. The observed geoid anomaly over the Atlantic margin of North America can be explained by Airy compensation. The relation between geoid anomaly and bathymetry across the Bermuda Swell is consistent with Pratt compensation with a 100-km depth of compensation.

Haxby, W. F.↗

A Multi-Season Study of the Effects of MODIS Sea-Surface Temperatures on Operational WRF Forecasts at NWS Miami, FL

Studies at the Short-term Prediction Research and Transition (SPORT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) sea-surface temperature (SST) composites in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. Recent work by LaCasse et al (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPORT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The project's goal is to determine whether more accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running WRF in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software. Twenty-seven hour forecasts are run dally initialized at 0300, 0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at 1/12deg resolution (approx.9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPORT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water, The MODIS SST composites for initializing the SPORT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST data into the SPORT WRF runs is staggered such that SSTs are updated with a new composite every six hours in each of the WRF runs. From mid-February to July 2007, over 500 parallel WRF simulations have been collected for analysis and verification. This paper will present verification results comparing the NWS MIA operational WRF runs to the SPORT experimental runs, and highlight any substantial differences noted in the predicted mesoscale phenomena for specific cases.

Case, Jonathan L.↗

Update on GOCART Model Development and Applications

Recent results from the GOCART and GMI models are reported. They include: Updated emission inventories for anthropogenic and volcano sources, satellite-derived vegetation index for seasonal variations of dust emission, MODIS-derived smoke AOT for assessing uncertainties of biomass-burning emissions, long-range transport of aerosol across the Pacific Ocean, and model studies on the multi-decadal trend of regional and global aerosol distributions from 1980 to 2010, volcanic aerosols, and nitrate aerosols. The document was presented at the 2013 AEROCENTER Annual Meeting held at the GSFC Visitors Center, May 31, 2013. The Organizers of the meeting are posting the talks to the public Aerocentr website, after the meeting.

emission inventories↗

Data from: Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

This data collection represents all UH8to2 model data used for the publication "Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model." This includes ocean and sea ice characteristics from the 2017-2020 UH8to2 model run. Further details about the model setup, forcing, and our interpretation of these data can be found in the publication.

58 GEOSCIENCES↗

SST Variation Due to Interactive Convective-Radiative Processes

The recent linking of Cloud-Resolving Models (CRMs) to Ocean-Mixed Layer (OML) models has provided a powerful new means of quantifying the role of cloud systems in ocean-atmosphere coupling. This is due to the fact that the CRM can better resolve clouds and cloud systems and allow for explicit cloud-radiation interaction. For example, Anderson (1997) applied an atmospheric forcing associated with a CRM simulated squall line to a 3-D OML model (one way or passive interaction). His results suggested that the spatial variability resulting from the squall forcing can last at least 24 hours when forced with otherwise spatially uniform fluxes. In addition, the sea surface salinity (SSS) variability continuously decreased following the forcing, while some of the SST variability remained when a diurnal mixed layer capped off the surface structure. The forcing used in the OML model, however, focused on shorter time (8 h) and smaller spatial scales (100-120 km). In this study, the 3-D Goddard Cumulus Ensemble Model (GCE; 512 x 512 x 23 cu km, 2-km horizontal resolution) is used to simulate convective active episodes occurring in the Western Pacific warm pool and Eastern Atlantic regions. The model is integrated for seven days, and the simulated results are coupled to an OML model to better understand the impact of precipitation and changes in the planetary boundary layer upon SST variation. We will specifically examine and compare the results of linking the OML model with various spatially-averaged outputs from GCE simulations (i.e., 2 km vs. 10-50 km horizontal resolutions), in order to help understand the SST sensitivity to multi-scale influences. This will allow us to assess the importance of explicitly simulated deep and shallow clouds, as well as the subgrid-scale effects (in coarse-model runs) upon SST variation. Results using both 1-D and 2-D OML models will be evaluated to assess the effects of horizontal advection.

Tao, W.-K.↗

Upper Atmosphere Heating From Ocean-Generated Acoustic Wave Energy

We present that colliding sea surface waves generate the ocean microbarom, an acoustic signal that may transmit significant energy to the upper atmosphere. Previous estimates of acoustic energy flux from the ocean microbarom and mountain-wind interactions are on the order of 0.01 to 1 mW/m 2 , heating the thermosphere by tens of Kelvins per day. We captured upgoing ocean microbarom waves with a balloon-borne infrasound microphone; the maximum acoustic energy flux was approximately 0.05 mW/m 2 . This is about half the average value reported in previous ground-based microbarom observations spanning 8 years. The acoustic flux from the microbarom episode described here may have heated the thermosphere by several Kelvins per day while the source persisted. Lastly, we suggest that ocean wave models could be used to parameterize acoustically generated heating of the upper atmosphere based on sea state.

58 GEOSCIENCES↗

Prediction and Predictability of the Madden Julian Oscillation in the NASA GEOS-5 Seasonal-to-Subseasonal System

In this study, we examine the prediction skill and predictability of the Madden Julian Oscillation (MJO) in a recent version of the NASA GEOS-5 atmosphere-ocean coupled model run at at 1/2 degree horizontal resolution. The results are based on a suite of hindcasts produced as part of the NOAA SubX project, consisting of seven ensemble members initialized every 5 days for the period 1999-2015. The atmospheric initial conditions were taken from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and the ocean and the sea ice were taken from a GMAO ocean analysis. The land states were initialized from the MERRA-2 land output, which is based on observation-corrected precipitation fields. We investigated the MJO prediction skill in terms of the bivariate correlation coefficient for the real-time multivariate MJO (RMM) indices. The correlation coefficient stays at or above 0.5 out to forecast lead times of 26-36 days, with a pronounced increase in skill for forecasts initialized from phase 3, when the MJO convective anomaly is located in the central tropical Indian Ocean. A corresponding estimate of the upper limit of the predictability is calculated by considering a single ensemble member as the truth and verifying the ensemble mean of the remaining members against that. The predictability estimates fall between 35-37 days (taken as forecast lead when the correlation reaches 0.5) and are rather insensitive to the initial MJO phase. The model shows slightly higher skill when the initial conditions contain strong MJO events compared to weak events, although the difference in skill is evident only from lead 1 to 20. Similar to other models, the RMM-index-based skill arises mostly from the circulation components of the index. The skill of the convective component of the index drops to 0.5 by day 20 as opposed to day 30 for circulation fields. The propagation of the MJO anomalies over the Maritime Continent does not appear problematic in the GEOS-5 hindcasts implying that the Maritime Continent predictability barrier may not be a major concern in this model. Finally, the MJO prediction skill in this version of GEOS-5 is superior to that of the current seasonal prediction system at the GMAO; this could be partly attributed to a slightly better representation of the MJO in the free running version of this model and partly to the improved atmospheric initialization from MERRA-2.

Achuthavarier, Deepthi↗

Dynamics of coupled ice-ocean system in the marginal ice zone: Study of the mesoscale processes and of constitutive equations for sea ice

This study is aimed at the modelling of mesoscale processed such as up/downwelling and ice edge eddies in the marginal ice zones. A 2-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model (f-plane) through interfacial stresses. The constitutive equations of the sea ice are formulated on the basis of the Reiner-Rivlin theory. The internal ice stresses are important only at high ice concentrations (90-100%), otherwise the ice motion is essentially free drift, where the air-ice stress is balanced by the ice-water stress. The model was tested by studying the upwelling dynamics. Winds parallel to the ice edge with the ice on the right produce upwilling because the air-ice momentum flux is much greater that air-ocean momentum flux, and thus the Ekman transport is bigger under the ice than in the open water. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4 day period. This forcing resembles successive cyclone passings. In the model with a thin oceanic upper layer, ice bands were formed.

Hakkinen, S.↗

MASK3 for Advanced WEC Dynamics and Controls

This report outlines the "MASK3" wave tank test within the Advanced WEC Dynamics and Controls (AWDC) project. This test represents the final test in the AWDC project. The focus of the MASK3 test was to consider coordinated 3-degree-of-freedom (3DOF) control of a WEC in a realistic ocean environment. A key aspect of this test was the inclusion of a "self-tuning" mechanism which uses an optimization algorithm to update controller gains based on a changing sea state. The successful implementation of the self-tuning mechanism is the last crucial step required for such a controller to be implemented in real ocean environments.

16 TIDAL AND WAVE POWER↗