Search NASA⌕ Search

Engineering topics

Reto Ruedy

Publications and source records attributed to Reto Ruedy.

The Sensitivity of the Equatorial Pacific ODZ to Particulate Organic Matter Remineralization in A Climate Model Under Pre-Industrial Conditions

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ O2 ol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Oxygen-deficient zone↗

Future Climate Change Under SSP Emission Scenarios with GISS-E2.1

This paper presents the response to anthropogenic forcing in the GISS-E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non-interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One-Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1-G) and HYCOM (E2.1-H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high-mitigation scenarios SSP1-1.9 and SSP1-2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1-H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5-8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1-G (−1.21 × 106 km2/°C) and E2.1-H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5-8.5 in the E2.1-G model, comparable to the reduction in the corresponding CMIP5 GISS-E2 RCP8.5 simulation. Both low-end scenarios SSP1-1.9 and SSP1-2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).

climate models↗

Effects of Spin-Orbit Resonances and Tidal Heating on the Inner Edge of the Habitable Zone

Much attention has been given to the climate dynamics and habitable boundaries of synchronously rotating planets around low mass stars. However, other rotational states are possible, including spin–orbit resonant configurations, particularly when higher eccentricity orbits can be maintained in a system. Additionally, the oscillating strain as a planet moves from periastron to apoastron results in friction and tidal heating, which can be an important energy source. Here, we simulate the climate of ocean-covered planets near the inner edge of the habitable zone around M to solar stars with the NASA GISS ROCKE-3D general circulation model, and leverage the planetary evolution software package, VPLanet, to calculate tidal heating rates for Earth-sized planets orbiting 2600 and 3000 K stars. This study is the first to use a 3D general circulation model that implements tidal heating to investigate habitability for multiple resonant states. We find that for reference experiments without tidal heating, the resonant state has little impact on the radial position of the inner edge because for a given stellar flux, higher-order states tend to be warmer than synchronous rotators, but for a given temperature, have drier upper atmospheres. However, when strong tidal heating is present, the rotational component implies a strong dependence of habitable conditions on the system evolution and rotational state. Since tidal and stellar heating both decrease rapidly with orbital distance, this results in a compact orbital width separating temperate and uninhabitable climates. We summarize these results and also compare ROCKE-3D to previously published simulations of the inner edge.

Spin-Orbit Resonances↗

Drivers of Air-Sea CO2 Flux Seasonality and its Long-Term Changes in the NASA-GISS Model CMIP6 Submission.

Climate change will affect both the mean state and seasonality of marine physical and biogeochemical properties, with important implications for the oceanic sink of atmospheric CO2. Here, we investigate the seasonal cycle of the air-sea exchange of CO2 and pCO2,sw (surface seawater pCO2) and their long term changes using the CMIP6 submission of the NASA-GISS modelE (GISS-E2.1-G). In comparison to the CMIP5 submission (GISS-E2-R), we find that on the global scale, the seasonal cycles of the CO2 flux and NPP have improved, while the seasonal cycles of DIC, alkalinity, and macronutrients have deteriorated. Moreover, for all ocean biogeochemietry fields, changes in skill between E2.1-G and E2-R display large regional variability. For E2.1-G, we find similar modeled and observed CO2 flux seasonal cycles in the subtropical gyres, where seasonal anomalies of pCO2,sw and the flux are temperature-driven, and the Southern Ocean, where anomalies are DIC-driven. Biases in these seasonal cycles are largest in the subpolar and equatorial regions, driven by a combination of biases in temperature, DIC, alkalinity, and wind speed. When comparing the historical simulation to a simulation with an idealized increase in atmospheric pCO2, we find that the seasonal amplitudes of the CO2 flux and pCO2,sw generally increase. These changes are produced by increases in the sensitivity of pCO2,sw to its respective drivers. These findings are consistent with the notion that the sensitivity of pCO2,sw is expected to increase due to the increase of atmospheric pCO2, with changes in the seasonality of temperature, DIC, and alkalinity having secondary influences.

Air-Sea CO2↗

CMIP6 Historical Simulations (1850‐2014) with GISS‐E2.1

Simulations of the CMIP6 historical period 1850‐‐2014, characterized by the emergence of anthropogenic climate drivers like greenhouse gases, are presented for different configurations of the NASA Goddard Institute for Space Studies (GISS) Earth System ModelE2.1. The GISS‐E2.1 ensembles are more sensitive to forcing than their CMIP5 predecessors (GISS‐E2), but warm less during recent decades due to reduced total forcing. This forcing reduction is attributed to an increase in longwave opacity in pre‐industrial simulations, resulting in an atmosphere less sensitive to further increases in opacity that result from rising greenhouse gas concentrations. This demonstrates the importance of the base climatology to forcing and forced climate trends. Most model versions match observed temperature trends since 1979 from the ocean to the stratosphere. The choice of ocean model is important to the transient climate response, as found in CMIP5 GISS‐E2: the model that more efficiently exports heat to the deep ocean shows a smaller rise in tropospheric temperature. Model sea level rise over the historical period is traced to excessive drawdown of aquifers to meet irrigation demand with a smaller contribution from thermal expansion. This shows how fully coupled models can provide indirect observational constraints upon forcing, in this case, constraining irrigation rates with observed sea level changes. The overall agreement of GISS‐E2.1 with observed trends is familiar from evaluation of its predecessors, as is the conclusion that these trends are almost entirely anthropogenic in origin.

NASA Goddard Institute for Space Studies (GISS) Ea↗

Global Carbon Cycle and Climate Feedbacks in the NASA GISS ModelE2.1

We present results from the NASA GISS ModelE2.1‐G‐CC Earth System Model with coupled climate‐carbon cycle simulations that were submitted to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) Coupled Climate‐Carbon Cycle MIP (C4MIP). Atmospheric CO2 concentration and carbon budgets for the land and ocean in the historical simulations were generally consistent with observations. Low simulated atmospheric CO2 concentrations during 1850‐1950 were due to excess uptake from prescribed land cover change, which erroneously replaced arid shrublands with higher biomass crops, and assumed high 2004 LAI values in vegetated lands throughout the historical simulation. At the end of the historical period, slightly higher simulated CO2 than observed resulted from the land being an insufficient net carbon sink, despite the net effect of CO2 fertilization and warming‐induced increases to leaf photosynthetic capacity. The global ocean carbon uptake agreed well with the observations with the largest discrepancies in the low latitudes. Future climate projection at 2091‐2100 agreed with CMIP5 models in the northward shift, of temperate deciduous forest climate and expansion across Eurasia along 60 °N latitude, and dramatic regional biome shifts from drying and warming in continental Europe. Carbon feedback parameters were largely similar to the CMIP5 model ensemble. For our model, the variation of land feedback parameters within the uncertainty arises from the fertilization feedback being less sensitive due to lack of increased vegetation growth, and the comparably more negative ocean carbon‐climate feedback is due to the large slowdown of the Atlantic overturning circulation.

Carbon cycle↗

GISS‐E2.1: Configurations and Climatology

This paper describes the GISS‐E2.1 contribution to the Coupled Model Intercomparison Project, Phase 6 (CMIP6). This model version differs from the predecessor model (GISS‐E2) chiefly due to parameterization improvements to the atmospheric and ocean model components, while keeping atmospheric resolution the same. Model skill when compared to modern era climatologies is significantly higher than in previous versions. Additionally, updates in forcings have a material impact on the results. In particular, there have been specific improvements in representations of modes of variability (such as the Madden‐Julian Oscillation and other modes in the Pacific) and significant improvements in the simulation of the climate of the Southern Oceans, including sea ice. The effective climate sensitivity to 2xCO2 is slightly higher than previously at 2.7‐‐3.1°C (depending on version), and is a result of lower CO2 radiative forcing and stronger positive feedbacks.

General Circulation Model↗

Influence of Heat Flow from the Planetary Interior on Surface Habitability

an primordial and/or radiogenic heat flow from a rocky planet’s interior help maintain a surface environment amenable to life?On rocky Earth-like planets with shallow oceans (≤5 km depth) that are on the outer edge of the habitable zone, the additional energy provided might make the difference between a world supporting regional habitat space and a world uninhabitable by any known life-form. To test this hypothesis, we introduce heat flux from the planetary interior as a forcing to the NASA/GISSROCKE-3D exoplanet GCM, and explore whether it can enhance the habitability of marginal terrestrial rocky planets (e.g.,paleo Earth in a “snowball” climate state).GCM simulations of modern Earth have not included geothermal heat flux as a climate forcing, as it is more than an order of magnitude less than the 2.9 W/m heating caused by anthropogenic greenhouse gases [Flanner, GRL 2009]. Modern Earth measurements in the oceans [Downes et al., GRL 2019] and results from ocean GCMs [Barnes et al., Ocean Model. 2017]suggest that localized higher heat flux can have measurable thermodynamic and dynamic impacts, such as weakening deep ocean stratification and increased poleward ocean heat transports. A similar result with ROCKE-3D would be important, since our snowball Earth simulations without heat flux produce features like ocean stratification that should have hampered the survival of life. We have created a set of spatially variable geothermal heat fluxes at varying horizontal resolutions based on modern Earth heat flow patterns [Davies, Geochem. Geophys. Geosys. 2013] to test the sensitivity of the ROCKE-3D GCM initially in a modern Earth context. Preliminary results show that even modest heat fluxes have a measurable impact on deep ocean temperatures (upto +0.9°C locally) and that this heat can be transported at depth along ocean circulation paths. Estimates of heat flux impacts on vertical ocean mixing and surface ocean conditions require long simulations at the higher spatial resolutions needed to express heat flow at the scale of most ocean ridges. Ongoing simulations are focused on testing the impacts of horizontal resolution and flux scaling on surface environments, especially for “snowball” climate scenarios.

Heat flow↗

A NASA GISTEMP Observational Uncertainty Ensemble: Regional and Monthly Uncertainty

The historical global temperature record is an essential data product for quantifying the variability and change of the Earth system. In recent years, better characterization of observational uncertainty in global and hemispheric trends has become available, but the methodologies are not necessarily applicable to analyses at smaller regional areas, or monthly means, where station sparsity and other systematic issues contribute to greater uncertainty. This work details a gridded uncertainty ensemble of historical temperature anomalies from the Goddard Institute for Space Studies (GISS) Surface Temperature product (GISTEMP) product. This ensemble characterizes the complex spatial and temporal correlation structure of uncertainty in gridded historical temperature, enabling proper uncertainty propagation for climate and social science at regional and monthly scales. This work details the methodology for generating the uncertainty ensemble, key statistics of the uncertainty evolution over space and time, and provides best practices for using the uncertainty ensemble in future studies. Summary statistics from the uncertainty ensemble are in good agreement with production GISTEMP. Two applications of the uncertainty ensemble are also presented. First, the warmest year on record is shown to most likely be 2016 with a 53.2% chance and 2020 as the second most likely with a 44.4% chance. Second, it is shown that the arctic is warming 2.5 - 5 times faster than the globe, significantly faster than the regularly quoted twice as fast.

GISTEMP↗