Search NASA⌕ Search

Engineering topics

Paul Lerner

Publications and source records attributed to Paul Lerner.

Aera-Mip: Emission Pathways, Remaining Budgets, and Carbon Cycle Dynamics Compatible With 1.5 and 2 °C Global Warming Stabilization

While international climate policies now focus on limiting global warming to well below 2 °C or pursuing a 1.5 °C level of global warming, the climate modelling community has not provided an experimental design in which all Earth system models (ESMs) converge and stabilize at the same prescribed global warming levels. This gap hampers accurate estimations based on comprehensive ESMs of the carbon emission pathways and budgets needed to meet such agreed warming levels and of the associated climate impacts under temperature stabilization. Here, we apply the Adaptive Emission Reduction Approach (AERA) with ESMs to provide such simulations in which all models converge at 1.5 and 2.0 °C warming levels by adjusting their emissions over time. These emission-driven simulations provide a wide range of emission pathways and resulting atmospheric CO 2 projections for a given warming level, uncovering uncertainty ranges that were previously missing in the traditional Coupled Model Intercomparison Project (CMIP) scenarios with prescribed greenhouse gas concentration pathways. Meeting the 1.5 °C warming level requires a 40 % (full model range: 7 % to 76 %) reduction in multi-model mean CO 2 -forcing-equivalent (CO 2 -fe) emissions from 2025 to 2030, a 98 % (57 % to 127 %) reduction from 2025 to 2050, and a stabilization at 1.0 (−1.7 to 2.9) PgC yr−1 from 2100 onward after the 1.5 °C global warming level is reached. Meeting the 2.0 °C warming level requires a 47 % (8 % to 92 %) reduction in multi-model mean CO 2 -fe emissions until 2050 and a stabilization at 1.7 (−1.5 to 2.7) PgC yr−1 from 2100 onward. The on-average positive emissions under stabilized global temperatures are the result of a decreasing transient climate response to cumulative CO 2 -fe emissions over time under stabilized global warming. This evolution is consistent with a slightly negative zero emissions commitment – initially assumed to be zero – and leads to an increase in the post-2025 CO 2 -fe emission budget by a factor of 2.2 (−0.8 to 6.9) by 2150 for the 1.5 °C warming level and a factor of 1.4 (0.9 to 2.4) for the 2.0 °C warming level compared to its first estimate in 2025. The median CO 2 -only carbon budget by 2150, relative to 2020, is 800 GtCO 2 for the 1.5 °C warming level and 2250 GtCO 2 for the 2.0 °C warming level. These median values exceed the median IPCC AR6 estimates by 60 % for the 1.5 °C warming level and 67 % for 2.0 °C. Some of the differences may be explained by the choice of the mitigation scenario for non-CO 2 radiative agents. Our simulations highlight shifts in carbon uptake dynamics under stabilized temperature, such as a cessation of the carbon sinks in the North Atlantic and in tropical forests. On the other hand, the Southern Ocean remains a carbon sink centuries after temperatures stabilize. Overall, this new type of warming-level-based emission-driven simulation offers a more coherent assessment across climate models and opens up a wide range of possibilities for studying both the carbon cycle and climate impacts, such as extreme events, under climate stabilization.

Adaptive Emission Reduction Approach↗

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↗

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↗

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↗

Influence of sediment resuspension on deep-ocean Pa-231 and Th-230 cycling: A PROGRESS REPORT

The naturally-occurring particle-reactive radionuclides protactinium-231 (231Pa) and thorium-230 (230Th) are used as tracers of a variety of oceanic processes, both at present and in the past. Most notably, the sediment 231Pa/230Th ratio has been used to infer changes in the Atlantic Meridional Overturning Circulation over the last (de)glaciation. However, recent measurements along the U.S. GEOTRACES North Atlantic transect (GA03) revealed two features which are at odds with current understanding about 231Pa and 230Th behaviour in the ocean: (i) a sharp decrease in dissolved 231Pa and 230Th activities with depth below 2000-4000 m and (ii) very large particulate 231Pa and 230Th activities near the bottom, at a number of stations between the New England continental shelf and Bermuda. Concomitant measurements of particulate matter concentration and potential temperature showed that both features are associated with the benthic nepheloid layer (BNL) and the bottom mixed layer (BML) that are present at these stations. Here we develop and apply a simplified model of the exchange of particles, 231Pa, and 230Th between the BNL and the upper sediment, to explore the extent to which the radionuclide anomalies observed near the bottom at a number of GA03 stations can be explained by local sediment resuspension. We find that the model can broadly reproduce the observed anomalies at two stations where samples for radionuclide analyses were collected near the seafloor. Sensitivity tests with the model show that the 231Pa/230Th ratio of particles in the BML and the sediment varies by a factor of 3 as the sediment resuspension rate fluctuates within a range consistent with observational estimates. The modelled variability is comparable to the spatial variability of 231Pa/230Th of suspended particles in the modern North Atlantic and to the variability of Atlantic sediment 231Pa/230Th records across the last (de)glacial period. Two factors are found to contribute to the modelled sensitivity of the sediment 231Pa/230Th to sediment resuspension rate: the vertical turbulent mixing in the BML and the differential scavenging intensity of Pa and Th due to variation in particle concentration. Overall, our study indicates that the exchange of material between the BNL and the upper sediment can affect the particulate 231Pa/230Th ratio in the bottom water and the sediment, which may complicate the use of sediment 231Pa/230Th as a palaeoceanographic tracer.

particle-reactive radionuclides↗

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↗

Development of Digital Twin Technologies for Climate Projections

Climate projections are increasingly needed for adaptation, climate resilience and related decision making. However, existing projections have systematic biases, are limited in scope, and are not readily available for most potential users. While the ideal of an observational data-driven ‘digital twin’ for climate is initially attractive, there is only a very limited set of climate data available with which to train such a tool. Nonetheless, we are confident that there is a role for ‘digital twin technologies’ in removing biases, increasing computational efficiency, expanding scenarios and data accessibility.

digital twins↗

Sampling the ENSO response to tropical volcanic eruptions in the GISS-E2-1 model: physical & biogeochemical ocean effects

Despite observational evidence that an El-Niño is likely to follow tropical volcanic eruptions (Predybaylo et al, 2017), modelling studies have so far reached no consensus on either the sign or physical mechanism of El Niño response to volcanism (Khodri et al, 2017) . A previous modelling study found that following the Pinatubo eruption, the tropical Pacific developed El Niño-like conditions that triggered a decrease in the upwelling of DIC and increased carbon uptake (Eddebbar et al., 2019). However, the carbon cycle response to volcanisms as it depends on ENSO phase remains unclear. Here we perform thirty ensemble control simulations using the GISS-E2-1-G model to sample the pre-industrial control ENSO cycle and thirty more to superimpose a Pinatubo-like volcanic eruption. We find that the model tends to favor the development of a weak El-Niño following the volcanic eruption only in the case of a volcano that erupts at the start of a neutral ENSO phase. In response to the weak El-Niño, some anomalies have the same sign as those that develop from a strong El-Niño, such as weaker winds, shallower mixed layer depths, less direct but more diffuse radiation at the surface of the ocean, less nitrate, and lower primary productivity. However, other anomalies show the opposite sign as those from a strong El-Niño, such as more precipitation, less surface DIC, and less surface ocean pCO 2 . In contrast to the neutral phase, the volcanic response during positive ENSO and negative ENSO phases are swamped by the model’s internal variability.

ENSO↗