SEARCH · Search NASA
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.
Ocean modeling plays an important role in both understanding the current climatic conditions and predicting future climate change. However, modeling the ocean circulation at various spatial and temporal scales is a very challenging computational task.
The Accelerated Climate Model for Energy (ACME), a new initiative by the U.S. Department of Energy, includes unstructured-mesh ocean, land-ice, and sea-ice components using the Model for Prediction Across Scales (MPAS) framework. The ability to run coupled high-resolution global simulations efficiently on large, high-performance computers is a priority for ACME. Sub-ice shelf ocean cavities are a significant new capability in ACME, and will be used to better understand how changing ocean temperature and currents influence glacial melting and retreat. These simulations take advantage of the horizontal variable-resolution mesh and adaptive vertical coordinate in MPAS-Ocean, in order to place high resolution below ice shelves and near grounding lines.
With the growing power and shrinking cost of personal computers. the availability of fast ethernet interconnections, and public domain software packages, it is now possible to combine them to build desktop parallel computers (named Beowulf or PC clusters) at a fraction of what it would cost to buy systems of comparable power front supercomputer companies. This led as to build and assemble our own sys tem. specifically for climate ocean modeling. In this article, we present our experience with such a system, discuss its network performance, and provide some performance comparison data with both HP SPP2000 and Cray T3E for an ocean Model used in present-day oceanographic research.
The focus of this paper is on the representation of Antarctic coastal polynyas in global ice-ocean general circulation models (OGCMs), in particular their local, regional, and high-frequency behavior. This is verified with the aid of daily ice concentration derived from satellite passive microwave data using the NASATeam 2 (NT2) and the bootstrap (BS) algorithms. Large systematic regional and temporal discrepancies arise, some of which are related to the type of convection parameterization used in the model. An attempt is made to improve the fresh-water flux associated with melting and freezing in Antarctic coastal polynyas by ingesting (assimilating) satellite ice concentration where it comes to determining the thermodynamics of the open-water fraction of a model grid cell. Since the NT2 coastal open-water fraction (polynyas) tends to be less extensive than the simulated one in the decisive season and region, assimilating NT2 coastal ice concentration yields overall reduced net freezing rates, smaller formation rates of Antarctic Bottom Water, and a stronger southward flow of North Atlantic Deep Water across 30 S. Enhanced net freezing rates occur regionally when NT2 coastal ice concentration is assimilated, concomitant with a more realistic ice thickness distribution and accumulation of High-Salinity Shelf Water. Assimilating BS rather than NT2 coastal ice concentration, the differences to the non-assimilated simulation are generally smaller and of opposite sign. This suggests that the model reproduces coastal ice concentration in closer agreement with the BS data than with the NT2 data, while more realistic features emerge when NT2 data are assimilated.
With the growing power and shrinking cost of personal computers, the availability of fast ethernet interconnections, and public domain software packages, it is now possible to combine them to build desktop parallel computers at a fraction of what it would cost to buy systems of comparable power from supercomputer companies.
A prequel study showed that dynamic downscaling using a regional climate model (RCM) over Africa improved the Goddard Institute for Space Studies Atmosphere-Ocean Global Climate Model (GISS AOGCM: ModelE) simulation of June–September rainfall patterns over Africa. The current study applies bias corrections to the lateral and lower boundary data from the AOGCM driving the RCM, based on the comparison of a 30-year simulation to the actual climate. The analysis examines the horizontal pattern of June–September total accumulated precipitation, the time versus latitude evolution of zonal mean West Africa (WA) precipitation (showing monsoon onset timing), and the latitude versus altitude cross-section of zonal winds over WA (showing the African Easterly Jet and the Tropical Easterly Jet). The study shows that correcting for excessively warm AOGCM Atlantic sea-surface temperatures (SSTs) improves the simulation of key features, whereas applying 30-year mean bias corrections to atmospheric variables driving the RCM at the lateral boundaries does not improve the RCM simulations. We suggest that AOGCM climate projections for Africa should benefit from downscaling by nesting an RCM that has demonstrated skill in simulating African climate, driven with bias-corrected SST.
Climate models exhibit major radiative biases over the Southern Ocean owing to a poor representation of mixed-phase clouds. In this study, we used the remote-sensing dataset from the Measurements of Aerosols, Radiation and Clouds over the Southern Ocean (MARCUS) campaign to assess the ability of the Weather Research and Forecasting (WRF) model to reproduce frontal clouds off Antarctica. It focuses on the modeling of thin mid-level supercooled liquid water layers which precipitate ice. The standard version of WRF produces almost fully glaciated clouds and cannot reproduce cloud top turbulence. Our work demonstrates the importance of adapting the ice nucleation parameterization to the pristine austral atmosphere to reproduce the supercooled liquid layers. Once simulated, droplets significantly impact the cloud radiative effect by increasing downwelling longwave fluxes and decreasing downwelling shortwave fluxes at the surface. The net radiative effect is a warming of snow and ice covered surfaces and a cooling of the ocean. Despite improvements in our simulations, the local turbulent circulation related to cloud-top radiative cooling is not properly reproduced, advocating for the need to develop a parameterization for top-down convection to capture the turbulence-microphysics interplay at cloud top.
We introduce a new climate model (GISS E2.2) that has been specially optimized for the Middle Atmosphere, and whose output is being contributed to the CMIP6 archive. The top of the model is at a geopotential altitude of 89 km, and parameterizations of moist convection and various forms of gravity wave drag based on tropospheric processes are chosen specifically for this optimization. We first evaluate the model in its configuration as a coupled atmosphere-chemistry model with respect to its simulation of the mean state of the Middle Atmosphere, from the mesosphere down through the Upper Troposphere/Lower Stratosphere. Then we assess its use as a coupled atmosphere-ocean climate model by exploring its mean ocean climatology. To evaluate its variability, we report on its simulation of the primary modes in the troposphere, stratosphere and ocean. Two climate change simulations are presented, the responses to instantaneous increases of 2xCO2 and 4xCO2, run with two different ocean models. Sensitivity studies are performed to illustrate the effect of parameterizations on the model results. We compare these results to the lower vertical resolution/top GISS Model E2.1, whose output has also been submitted to CMIP6. The different choices made for these models are explored. It is shown that important improvements in the circulation above and below the tropopause can be obtained when attention is paid to representation of Middle Atmosphere processes in climate model development.
The Ocean Heat-Carbon Nexus, linking ocean heat and carbon uptake, is crucial for understanding climate responses to cumulative carbon dioxide (CO 2 ) emissions and to net zero CO 2 emissions. It results from a suite of processes involving the exchange of heat and carbon across the sea-air interface as well as their storage below the mixed layer and redistribution by the ocean large-scale circulation. The Ocean Heat and Carbon Nexus is assumed to be consistently represented across two modelling platforms used in the latest IPCC assessments: the Earth System Models (ESMs) and the Simple Climate Models (SCMs). However, our research shows significant deficiencies in state-of-the-art SCMs in replicating the ocean heat-carbon nexus of ESMs due to a crude treatment of the ocean thermal and carbon cycle coupling. With one SCM, we show that a more realistic heat-to-carbon uptake ratio exacerbates the projected warming by 0.1 °C in low overshoot scenarios and up to 0.2 °C in high overshoot scenarios. It is therefore critical to explore how SCMs’ physical inconsistencies, such as the representation of the ocean heat-carbon nexus, can affect future warming projections used in climate assessments, not just by SCMs in Working Group 3 but also by ESMs in Working Group 1 via SCM-driven emission-to-concentration translation.
Convective transition statistics, which describe the relation between column-integrated water vapor (CWV) and precipitation, are compiled over tropical oceans using satellite and ARM site measurements to quantify the temperature and resolution dependence of the precipitation–CWV relation at fast time scales relevant to convection. At these time scales, and for precipitation especially, uncertainties associated with observational systems must be addressed by examining features with a variety of instrumentation and identifying robust behaviors versus instrument sensitivity at high rain rates. Here the sharp pickup in precipitation as CWV exceeds a certain critical threshold is found to be insensitive to spatial resolution, with convective onset occurring at higher CWV but at lower column relative humidity as bulk tropospheric temperature increases. Mean tropospheric temperature profiles conditioned on precipitation show vertically coherent structure across a wide range of temperature, reaffirming the use of a bulk temperature measure in defining the convective transition statistics. The joint probability distribution of CWV and precipitation develops a peak probability at low precipitation for CWV above critical, with rapidly decreasing probability of high precipitation below and near critical, and exhibits systematic changes under spatial averaging. The precipitation pickup with CWV is reasonably insensitive to time averaging up to several hours but is smoothed at daily time scales. This work demonstrates that CWV relative to critical serves as an effective predictor of precipitation with only minor geographic variations in the tropics, quantifies precipitation-related statistics subject to different spatial–temporal resolution, and provides a baseline for model comparison to apply these statistics as observational constraints on precipitation processes.
Process models of phytoplankton production are discussed with respect to their suitability for incorporation into global-scale numerical ocean circulation models. Exact solutions are given for integrals over the mixed layer and the day of analytic, wavelength-independent models of primary production. Within this class of model, the bias incurred by using a triangular approximation (rather than a sinusoidal one) to the variation of surface irradiance through the day is computed. Efficient computation algorithms are given for the nonspectral models. More exact calculations require a spectrally sensitive treatment. Such models exist but must be integrated numerically over depth and time. For these integrations, resolution in wavelength, depth, and time are considered and recommendations made for efficient computation. The extrapolation of the one-(spatial)-dimension treatment to large horizontal scale is discussed.
The Southern Ocean surface has freshened in recent decades, increasing water column stability and reducing upwelling of warmer subsurface waters. The majority of CMIP5 models underestimate or fail to capture this historical surface freshening, yet little is known about the impact of this model bias on regional ocean circulation and hydrography. Here experiments are performed using a global coupled climate model with additional freshwater applied to the Southern Ocean to assess the influence of recent surface freshening. The simulations explore the impact of persistent and long-term broad-scale freshening as a result of processes including precipitation minus evaporation changes. Thus, unlike previous studies, the freshening is applied as far north as 55°S, beyond the Antarctic ice margin. It is found that imposing a large-scale surface freshening causes a surface cooling and sea ice increase under preindustrial conditions, because of a reduction in ocean convection and weakened entrainment of warm subsurface waters into the surface ocean. This is consistent with intermodel relationships between CMIP5 models and the simulations, suggesting that models with larger surface freshening also exhibit stronger surface cooling and increased sea ice. Additional experiments are conducted with surface salinity restoration applied to capture observed regional salinity trends. Remarkably, without any mechanical wind trend forcing, these simulations accurately represent the spatial pattern of observed surface temperature and sea ice trends around Antarctica. This study highlights the importance of accurately simulating changes in Southern Ocean salinity to capture changes in ocean circulation, sea surface temperature, and sea ice.
Hydrological mass transport in the geophysical fluids of the atmosphere-hydrosphere-solid Earth surface system can excite Earth's rotational variations in both length-of-day and polar motion. These effects can be computed in terms of the hydrological angular momentum by proper integration of global meteorological data. We do so using the 40-year NCEP data and the 18-year NASA GEOS-1 data, where the precipitation and evapotranspiration budgets are computed via the water mass balance of the atmosphere based on Oki et al.'s (1995) algorithm. This hydrological mass redistribution will also cause geocenter motion and changes in Earth's gravitational field, which are similarly computed using the same data sets. Corresponding geodynamic effects due to the oceanic mass transports (i.e. oceanic angular momentum and ocean-induced geocenter/gravity changes) have also been computed in a similar manner. We here compare two independent sets of the result from: (1) non-steric ocean surface topography observations based on Topex/Poseidon, and (2) the model output of the mass field by the Parallel Ocean Climate Model. Finally, the hydrological and the oceanic time series are combined in an effort to better explain the observed non-atmospheric effects. The latter are obtained by subtracting the atmospheric angular momentum from Earth rotation observations, and the atmosphere- induced geocenter/gravity effects from corresponding geodetic observations, both using the above-mentioned atmospheric data sets.
To assess deep convective parameterizations in a variety of GCMs and examine the fast-time-scale convective transition, a set of statistics characterizing the pickup of precipitation as a function of column water vapor (CWV), PDFs and joint PDFs of CWV and precipitation, and the dependence of the moisture–precipitation relation on tropospheric temperature is evaluated using the hourly output of two versions of the GFDL Atmospheric Model, version 4 (AM4), NCAR CAM5 and superparameterized CAM (SPCAM). The 6-hourly output from the MJO Task Force (MJOTF)/GEWEX Atmospheric System Study (GASS) project is also analyzed. Contrasting statistics produced from individual models that primarily differ in representations of moist convection suggest that convective transition statistics can substantially distinguish differences in convective representation and its interaction with the large-scale flow, while models that differ only in spatial–temporal resolution, microphysics, or ocean–atmosphere coupling result in similar statistics. Most of the models simulate some version of the observed sharp increase in precipitation as CWV exceeds a critical value, as well as that convective onset occurs at higher CWV but at lower column RH as temperature increases. While some models quantitatively capture these observed features and associated probability distributions, considerable intermodel spread and departures from observations in various aspects of the precipitation–CWV relationship are noted. For instance, in many of the models, the transition from the low-CWV, nonprecipitating regime to the moist regime for CWV around and above critical is less abrupt than in observations. Additionally, some models overproduce drizzle at low CWV, and some require CWV higher than observed for strong precipitation. For many of the models, it is particularly challenging to simulate the probability distributions of CWV at high temperature.
Explore the source record for details and available documents.
Improvements in time-stepping schemes can provide substantial improvements in the throughput of global ocean-climate models, but is challenging due to large variations in time scales in the dynamics and spatial scales in the global, variable resolution mesh. Two funded projects aim to speed up the Model for Prediction Across Scales-Ocean (MPAS-Ocean). The first will add a variable-timestepping formulation to MPAS-Ocean, which is a critical step towards running efficiently on global variable-resolution meshes with regional and coastal refinement. The second optimizes the barotropic (depth-integrated) mode, which is currently a performance bottleneck at high core counts. The proposed compute time will be used for extensive model validation and parameter studies for both methods.
This study compares the impacts of Arctic sea ice decline on the Atlantic meridional overturning circulation (AMOC) in two configurations of the Community Earth System Model with different horizontal resolution. In a suite of model experiments, we impose radiative imbalance at the ice surface, replicating a loss of sea ice cover comparable to that observed during 1979–2014, and we find dramatic differences in the AMOC response between the two models. In the lower-resolution configuration, the AMOC weakens by about one-third over the first 100 years, approaching a new quasi-equilibrium. By contrast, in the higher-resolution configuration, the AMOC weakens by ~10% during the first 20–30 years followed by a full recovery driven by invigorated deep water formation in the Labrador Sea and adjacent regions. We investigate these differences using a diagnostic AMOC stability indicator, which reflects the AMOC freshwater transport in and out of the basin and hence the strength of the basin-scale salt-advection feedback. This indicator suggests that the AMOC in the lower-resolution model is less stable and more sensitive to surface perturbations, as confirmed by hosing experiments mimicking Arctic freshening due to sea ice decline. Differences between the models’ mean states, including the Atlantic Ocean mean surface freshwater fluxes, control the differences in AMOC stability. Our results demonstrate that the AMOC stability indicator is indeed useful for evaluating AMOC sensitivity to perturbations. We emphasize that, despite the differences in the long-term adjustment, both models simulate a multidecadal AMOC weakening caused by Arctic sea ice decline, relevant to climate change.