Engineering PapersSearch

Engineering topics

LeGrande, Allegra N.

Publications and source records attributed to LeGrande, Allegra N..

A Double Bootstrap Approach to Superposed Epoch Analysis to Evaluate Response Uncertainty

The association between climate variability and episodic events, such as the antecedent moisture conditions prior to wildfire or the cooling following volcanic eruptions, is commonly assessed using Superposed Epoch Analysis (SEA). In SEA the epochal response is typically calculated as the average climate conditions prior to and following all event years or their deviation from climatology. However, the magnitude and significance of the inferred climate association may be sensitive to the selection or omission of individual key years, potentially resulting in a biased assessment of the relationship between these events and climate. Here we describe and test a modified double-bootstrap SEA that generates multiple unique draws of the key years and evaluates the sign, magnitude, and significance of event-climate relationships within a probabilistic framework. This multiple resampling helps quantify multiple uncertainties inherent in conventional applications of SEA within dendrochronology and paleoclimatology. We demonstrate our modified SEA by evaluating the volcanic cooling signal in a Northern Hemisphere tree-ring temperature reconstruction and the link between drought and wildfire events in the western United States. Finally, we make our Matlab and R code available to be adapted for future SEA applications.

Dendroclimatology

Assessing the Impact of Large Volcanic Eruptions of the Last Millennium (850-1850 CE) on Australian Rainfall Regimes

Explosive volcanism is an important natural climate forcing, impacting global surface temperatures and regional precipitation. Although previous studies have investigated aspects of the impact of tropical volcanism on various ocean–atmosphere systems and regional climate regimes, volcanic eruptions remain a poorly understood climate forcing and climatic responses are not well constrained. In this study, volcanic eruptions are explored in particular reference to Australian precipitation, and both the Indian Ocean Dipole (IOD) and El Niño–Southern Oscillation (ENSO). Using nine realizations of the last millennium (LM) (850–1850CE) with different time-evolving forcing combinations, from the NASA GISS ModelE2-R, the impact of the six largest tropical volcanic eruptions of this period are investigated. Overall, we find that volcanic aerosol forcing increased the likelihood of El Niño and positive IOD conditions for up to four years following an eruption, and resulted in positive precipitation anomalies over north-west (NW) and south-east (SE) Australia. Larger atmospheric sulfate loading during larger volcanic eruptions coincided with more persistent positive IOD and El Niño conditions, enhanced positive precipitation anomalies over NW Australia, and dampened precipitation anomalies over SE Australia.

Blake, Stephanie A. P.

Multi-model Comparison of the Volcanic Sulfate Deposition from the 1815 Eruption of Mt. Tambora

The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyze both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-core derived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Marshall, Lauren

The PMIP4 Contribution to CMIP6-Part 4: Scientific Objectives and Experimental Design of the PMIP4-CMIP6 Last Glacial Maximum Experiments and PMIP4 Sensitivity Experiments

The Last Glacial Maximum (LGM, 21,000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Kageyama, Masa

The DeepMIP Contribution to PMIP4: Experimental Design for Model Simulations of the EECO, PETM, and pre-PETM (version 1.0)

Past warm periods provide an opportunity to evaluate climate models under extreme forcing scenarios, in particular high ( greater than 800 ppmv) atmospheric CO2 concentrations. Although a post hoc intercomparison of Eocene (approximately 50 Ma) climate model simulations and geological data has been carried out previously, models of past high-CO2 periods have never been evaluated in a consistent framework. Here, we present an experimental design for climate model simulations of three warm periods within the early Eocene and the latest Paleocene (the EECO, PETM, and pre-PETM). Together with the CMIP6 pre-industrial control and abrupt 4(times) CO2 simulations, and additional sensitivity studies, these form the first phase of DeepMIP - the Deep-time Model Intercomparison Project, itself a group within the wider Paleoclimate Modeling Intercomparison Project (PMIP). The experimental design specifies and provides guidance on boundary conditions associated with palaeogeography, greenhouse gases, astronomical configuration, solar constant, land surface processes, and aerosols. Initial conditions, simulation length, and output variables are also specified. Finally, we explain how the geological data sets, which will be used to evaluate the simulations, will be developed.

climate models

The Model Intercomparison Project on the Climatic Response to Volcanic Forcing (VolMIP): Experimental Design and Forcing Input Data for CMIP6

The enhancement of the stratospheric aerosol layer by volcanic eruptions induces a complex set of responses causing global and regional climate effects on a broad range of timescales. Uncertainties exist regarding the climatic response to strong volcanic forcing identified in coupled climate simulations that contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). In order to better understand the sources of these model diversities, the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) has defined a coordinated set of idealized volcanic perturbation experiments to be carried out in alignment with the CMIP6 protocol. VolMIP provides a common stratospheric aerosol data set for each experiment to minimize differences in the applied volcanic forcing. It defines a set of initial conditions to assess how internal climate variability contributes to determining the response. VolMIP will assess to what extent volcanically forced responses of the coupled ocean-atmosphere system are robustly simulated by state-of-the-art coupled climate models and identify the causes that limit robust simulated behavior, especially differences in the treatment of physical processes. This paper illustrates the design of the idealized volcanic perturbation experiments in the VolMIP protocol and describes the common aerosol forcing input data sets to be used.

Volcanoes

Supplement of: The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El NioSouthern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium.An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records.Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the amount effect. During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher

The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El Niño-Southern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850 CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium. An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records. Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the "amount effect". During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher M.

Volcanic Eruptions and Climate

Volcanic eruptions represent some of the most climatically important and societally disruptive short-term events in human history. Large eruptions inject ash, dust, sulfurous gases (e.g. SO2, H2S), halogens (e.g. Hcl and Hbr), and water vapor into the Earth's atmosphere. Sulfurous emissions principally interact with the climate by converting into sulfate aerosols that reduce incoming solar radiation, warming the stratosphere and altering ozone creation, reducing global mean surface temperature, and suppressing the hydrological cycle. In this issue, we focus on the history, processes, and consequences of these large eruptions that inject enough material into the stratosphere to significantly affect the climate system. In terms of the changes wrought on the energy balance of the Earth System, these transient events can temporarily have a radiative forcing magnitude larger than the range of solar, greenhouse gas, and land use variability over the last millennium. In simulations as well as modern and paleoclimate observations, volcanic eruptions cause large inter-annual to decadal-scale changes in climate. Active debates persist concerning their role in longer-term (multi-decadal to centennial) modification of the Earth System, however.

Volcanoes

Comparison of Forced ENSO-Like Hydrological Expressions in Simulations of the Preindustrial and Mid-Holocene

Using the water isotope- and vapor source distribution (VSD) tracer-enabled Goddard Institute for Space Studies ModelE-R, we examine changing El Nino-Southern Oscillation (ENSO)-like expressions in the hydrological cycle in a suite of model experiments. We apply strong surface temperature anomalies associated with composite observed El Nino and La Nina events as surface boundary conditions to preindustrial and mid-Holocene model experiments in order to investigate ENSO-like expressions in the hydrological cycle under varying boundary conditions. We find distinct simulated hydrological anomalies associated with El Nino-like ("ENSOWARM") and La Nina-like ("ENSOCOOL") conditions, and the region-specific VSD tracers show hydrological differences across the Pacific basin between El Nino-like and La Nina-like events. The application of ENSOCOOL forcings does not produce climatological anomalies that represent the equal but opposite impacts of the ENSOWARM experiment, as the isotopic anomalies associated with ENSOWARM conditions are generally stronger than with ENSOCOOL and the spatial patterns of change distinct. Also, when the same ENSO-like surface temperature anomalies are imposed on the mid-Holocene, the hydrological response is muted, relative to the preindustrial. Mid-Holocene changes in moisture sources to the analyzed regions across the Pacific reveal potentially complex relationships between ENSO-like conditions and boundary conditions. Given the complex impacts of ENSO-like conditions on various aspects of the hydrological cycle, we suggest that proxy record insights into paleo-ENSO variability are most likely to be robust when synthesized from a network of many spatially diverse archives, which can account for the potential nonstationarity of ENSO teleconnections under different boundary conditions.

boundary conditions

Configuration and Assessment of the GISS ModelE2 Contributions to the CMIP5 Archive

We present a description of the ModelE2 version of the Goddard Institute for Space Studies (GISS) General Circulation Model (GCM) and the configurations used in the simulations performed for the Coupled Model Intercomparison Project Phase 5 (CMIP5). We use six variations related to the treatment of the atmospheric composition, the calculation of aerosol indirect effects, and ocean model component. Specifically, we test the difference between atmospheric models that have noninteractive composition, where radiatively important aerosols and ozone are prescribed from precomputed decadal averages, and interactive versions where atmospheric chemistry and aerosols are calculated given decadally varying emissions. The impact of the first aerosol indirect effect on clouds is either specified using a simple tuning, or parameterized using a cloud microphysics scheme. We also use two dynamic ocean components: the Russell and HYbrid Coordinate Ocean Model (HYCOM) which differ significantly in their basic formulations and grid. Results are presented for the climatological means over the satellite era (1980-2004) taken from transient simulations starting from the preindustrial (1850) driven by estimates of appropriate forcings over the 20th Century. Differences in base climate and variability related to the choice of ocean model are large, indicating an important structural uncertainty. The impact of interactive atmospheric composition on the climatology is relatively small except in regions such as the lower stratosphere, where ozone plays an important role, and the tropics, where aerosol changes affect the hydrological cycle and cloud cover. While key improvements over previous versions of the model are evident, these are not uniform across all metrics.

ocean model

CMIP5 Historical Simulations (1850-2012) with GISS ModelE2

Observations of climate change during the CMIP5 extended historical period (1850-2012) are compared to trends simulated by six versions of the NASA Goddard Institute for Space Studies ModelE2 Earth System Model. The six models are constructed from three versions of the ModelE2 atmospheric general circulation model, distinguished by their treatment of atmospheric composition and the aerosol indirect effect, combined with two ocean general circulation models, HYCOM and Russell. Forcings that perturb the model climate during the historical period are described. Five-member ensemble averages from each of the six versions of ModelE2 simulate trends of surface air temperature, atmospheric temperature, sea ice and ocean heat content that are in general agreement with observed trends, although simulated warming is slightly excessive within the past decade. Only simulations that include increasing concentrations of long-lived greenhouse gases match the warming observed during the twentieth century. Differences in twentieth-century warming among the six model versions can be attributed to differences in climate sensitivity, aerosol and ozone forcing, and heat uptake by the deep ocean. Coupled models with HYCOM export less heat to the deep ocean, associated with reduced surface warming in regions of deepwater formation, but greater warming elsewhere at high latitudes along with reduced sea ice. All ensembles show twentieth-century annular trends toward reduced surface pressure at southern high latitudes and a poleward shift of the midlatitude westerlies, consistent with observations.

ocean heat content

Spring-Summer Temperatures Since AD 1780 Reconstructed from Stable Oxygen Isotope Ratios in White Spruce Tree-Rings from the Mackenzie Delta, Northwestern Canada

High-latitude delta(exp 18)O archives deriving from meteoric water (e.g., tree-rings and ice-cores) can provide valuable information on past temperature variability, but stationarity of temperature signals in these archives depends on the stability of moisture source/trajectory and precipitation seasonality, both of which can be affected by atmospheric circulation changes. A tree-ring delta(exp 18)O record (AD 1780-2003) from the Mackenzie Delta is evaluated as a temperature proxy based on linear regression diagnostics. The primary source of moisture for this region is the North Pacific and, thus, North Pacific atmospheric circulation variability could potentially affect the tree-ring delta(exp 18)O-temperature signal. Over the instrumental period (AD 1892-2003), tree-ring delta(exp 18)O explained 29% of interannual variability in April-July minimum temperatures, and the explained variability increases substantially at lower-frequencies. A split-period calibration/verification analysis found the delta(exp 18)O-temperature relation was time-stable, which supported a temperature reconstruction back to AD 1780. The stability of the delta(exp 18)O-temperature signal indirectly implies the study region is insensitive to North Pacific circulation effects, since North Pacific circulation was not constant over the calibration period. Simulations from the NASA-GISS ModelE isotope-enabled general circulation model confirm that meteoric delta(exp 18)O and precipitation seasonality in the study region are likely insensitive to North Pacific circulation effects, highlighting the paleoclimatic value of tree-ring and possibly other delta(exp 18)O records from this region. Our delta(exp 18)O-based temperature reconstruction is the first of its kind in northwestern North America, and one of few worldwide, and provides a long-term context for evaluating recent climate warming in the Mackenzie Delta region.

Canada

Linking the 8.2 ka Event and its Freshwater Forcing in the Labrador Sea

The 8.2 ka event was the last deglacial abrupt climate event. A reduction in the Atlantic meridional overturning circulation (AMOC) attributed to the drainage of glacial Lake Agassiz may have caused the event, but the freshwater signature of Lake Agassiz discharge has yet to be identified in (delta)18O of foraminiferal calcite records from the Labrador Sea, calling into question the connection between freshwater discharge to the North Atlantic and AMOC strength. Using Mg/Ca-paleothermometry, we demonstrate that approx. 3 C of near-surface ocean cooling masked an 1.0 % decrease in western Labrador Sea (delta)18O of seawater concurrent with Lake Agassiz drainage. Comparison with North Atlantic (delta)18O of seawater records shows that the freshwater discharge was transported to regions of deep-water formation where it could perturb AMOC and force the 8.2 ka event.

Hoffman, Jeremy S.

Sensitivity of Seawater Oxygen Isotopes to Climatic and Tectonic Boundary Conditions in an Early Paleogene Simulation with GISS ModelE-R

An isotope-enabled ocean-atmosphere general circulation model (GISS ModelE -R) is used to estimate the spatial gradients of the oxygen isotopic composition of seawater (delta O-18(sub sw), where delta is the deviation from a known reference material in per mil) during the early Paleogene (45.65 Ma). Understanding the response of delta O-18(sub sw) to changes in climatic and tectonic boundary conditions is important because records of carbonate delta O-18 document changes in hydrology, as well as changes in temperature and global ice -volume. We present results from an early Paleogene configuration of ModelE -R which indicate that spatial gradients of surface ocean delta O-18(sub sw) during this period could have been significantly different to those in the modern ocean. The differences inferred from ModelE -R are sufficient to change early Paleogene sea surface temperature estimates derived from primary carbonate delta O-18 signatures by more than +/-2 C in large areas of the ocean. In the North Atlantic, Indian, and Southern Oceans, the differences in d18Osw inferred from our simulation with ModelE -R are in direct contrast with those from another d18O ]tracing model study which used different, but equally plausible, early Paleogene boundary conditions. The large differences in delta O-18(sub sw) between preindustrial and early Paleogene simulations, and between models, emphasizes the sensitivity of d18Osw to climatic and tectonic boundary conditions. For this reason, absolute estimates of Eocene/ Paleocene temperature derived from carbonate delta O-18 alone are likely to have larger uncertainties than are usually assumed.

Roberts, Christopher D.

Model, Proxy and Isotopic Perspectives on the East African Humid Period

Both North and East Africa experienced more humid conditions during the early and mid-Holocene epoch (11,000-5000yr BP; 11-5 ka) relative to today. The North African Humid Period has been a major focus of paleoclimatic study, and represents a response of the hydrological cycle to the increase in boreal summer insolation and associated ocean, atmosphere and land surface feedbacks. Meanwhile, the mechanisms that caused the coeval East African Humid Period are poorly understood. Here, we use results from isotopeenabled coupled climate modeling experiments to investigate the cause of the East African Humid Period. The modeling results are interpreted alongside proxy records of both water balance and the isotopic composition of rainfall. Our simulations show that the orbitally-induced increase in dry season precipitation and the subsequent reduction in precipitation seasonality can explain the East African Humid Period, and this scenario agrees well with regional lake level and pollen paleoclimate data. Changes in zonal moisture flux from both the Atlantic and Indian Ocean account for the simulated increase in precipitation from June through November. Isotopic paleoclimate data and simulated changes in moisture source demonstrate that the western East African Rift Valley in particular experienced more humid conditions due to the influx of Atlantic moisture and enhanced convergence along the Congo Air Boundary. Our study demonstrates that zonal changes in moisture advection are an important determinant of climate variability in the East African region.

Tierney, Jessica E.

The Anomalous Winter of 1783-1784: Was the Laki Eruption or an Analog of the 2009-2010 Winter to Blame?

The multi ]stage eruption of the Icelandic volcano Laki beginning in June, 1783 is speculated to have caused unusual dry fog and heat in western Europe and cold in North America during the 1783 summer, and record cold and snow the subsequent winter across the circum-North Atlantic. Despite the many indisputable impacts of the Laki eruption, however, its effect on climate, particularly during the 1783.1784 winter, may be the most poorly constrained. Here we test an alternative explanation for the unusual conditions during this time: that they were caused primarily by a combined negative phase of the North Atlantic Oscillation (NAO) and an El Nino ]Southern Oscillation (ENSO) warm event. A similar combination of NAO ]ENSO phases was identified as the cause of record cold and snowy conditions during the 2009.2010 winter in Europe and eastern North America. 600-year tree-ring reconstructions of NAO and ENSO indices reveal values in the 1783.1784 winter second only to their combined severity in 2009.2010. Data sources and model simulations support our hypothesis that a combined, negative NAO ]ENSO warm phase was the dominant cause of the anomalous winter of 1783.1784, and that these events likely resulted from natural variability unconnected to Laki.

D'Arrigo, Rosanne