Svetlana Botsyun*
Todd A. Ehlers- Department of Geosciences, University of Tübingen, Tübingen, Germany
Paleoelevation reconstructions derived from proxy data such as stable oxygen isotope records in terrestrial archives have been determined for Cenozoic mountain ranges around the world. Recent studies have highlighted that a variety of paleoclimate processes can contribute to the isotopic composition of a measured precipitation (δ18Op) signal used in elevation reconstructions. These processes can include: regional, global, and topographic variations in paleotemperature; environmental conditions of an air mass before orographic ascent; evapotranspiration; water vapor recycling; and changes in the vapor source. In some cases, these processes can overprint the elevation signal sought in proxy data and preclude robust elevation reconstructions. Recent advances in isotope tracking climate models allow us to estimate paleoclimate changes during orogen development and associated changes in paleo δ18Op due to both climate and topographic changes. These models account for adiabatic and non-adiabatic temperature changes, relative humidity variations, changing continental evapotranspiration, vapor recycling, vapor source changes, etc. Modeling strategies using high-resolution isotopes-enabled General Circulation Models (iGCMs) together with time-specific boundary conditions and variable topography provide a powerful tool for enhancing elevation reconstructions from δ18Op proxy data. In this review, we discuss the principles, benefits and caveats of using iGCMs for interpreting isotopic records from natural archives for paleoelevation reconstructions. We also highlight future challenges for the application of iGCMs to paleoaltimetry proxy data that open up new avenues for research on tectonic-climate interactions.
Introduction
Quantitative paleoelevation techniques based on proxy data include approaches such as stable isotopes (e.g., δ18O, δD) measured in terrestrial deposits, paleofloral and faunal findings and their physiognomic characteristics, and clumped isotope (Δ47) paleothermometry. These techniques (and others) have been extensively applied to various mountain ranges such as the Himalayas and Tibetan Plateau (Rowley and Currie, 2006; Gébelin et al., 2013; Ding et al., 2014), the North America Cordillera (Chamberlain et al., 2012; Gébelin et al., 2012; Cassel et al., 2014), the Andes and Andean Plateau (Mulch et al., 2010; Garzione et al., 2017) for elevation reconstructions throughout the Cenozoic, ranging from the Paleocene and Eocene (e.g., Ding et al., 2014) up to the Pleistocene (e.g., Hoke et al., 2014). Oxygen isotope paleoaltimetry is one of the most widely applied techniques and is based on a progressive decrease in the δ18O value of precipitation (δ18Op) with elevation during the topographic (i.e., orographic) ascent of an air mass. This technique assumes Rayleigh distillation of an air mass under equilibrium conditions and without local evaporation or advection of additional vapor from outside the system (Rozanski et al., 1993; Gat, 1996).
Oxygen isotope paleoaltimetry uses δ18O preserved in pedogenic or lacustrine carbonates (hereafter referred to as δ18Oc) as a proxy for paleoelevation. Interpretation of measured δ18Oc generally requires two steps. For the first step, the reconstruction of paleoprecipitation δ18O (δ18Opp) using the calcite-water fractionation coefficients depend on the carbonate formation temperature as described through the equilibrium between calcite and water δ18O (Kim and O’Neil, 1997). Paleoaltimetry estimates based on δ18O alone, without independent information on paleotemperatures from other proxies or a paleoclimate model, depend on assumptions made about the paleoelevation and temperature (global and regional) changes relative to the modern. Knowledge of these paleo-conditions is challenging to come by. Moreover, if a surface temperature is assumed, then the calculation of paleoelevation from the proxy involves circular reasoning. This is because the temperature used at the point of interest depends on the elevation through the temperature lapse rate and often the chosen coefficients from the equations in Kim and O’Neil (1997) already implies high elevations.
For the second step, the reconstructed δ18Opp must be attributed to a paleoelevation using a δ18Opp-elevation relationship (the “isotopic lapse rate”) in a reference to low-elevation δ18O. A common approach (largely due to it’s simplicity) is to assume that the modern isotopic lapse rate measured from precipitation in an elevation-transect of stations (e.g., Gonfiantini et al., 2001; Fiorella et al., 2015), or simulated using a one-dimensional thermodynamic model for δ18Op composition based on the Rayleigh distillation, can be used to predict changes in the δ18Op response to changes in temperature, and the relative humidity of air parcels during orographic ascent (Rowley et al., 2001; Rowley and Garzione, 2007). This approach assumes the isotopic lapse rate at the time of carbonate formation not to change over millions of years. Thus, both of the previous steps rely on several assumptions regarding the contribution (if any) of global and regional climate change to δ18O and the stationarity of the δ18Opp-elevation relationship throughout geologic time.
However, numerous recent studies have shown that a range of atmospheric processes are important for the interpretation of δ18Opp. For example, hemispheric-scale atmospheric circulation and associated teleconnections (Schneider and Noone, 2007; Takahashi and Battisti, 2007; Pausata et al., 2011), moisture transport changes (Ehlers and Poulsen, 2009; Poulsen et al., 2010), evapotranspiration and water vapor recycling within continental interiors (Risi et al., 2013; Chamberlain et al., 2014), and shifts in precipitation amount and ratios of convective-to-large-scale precipitation (Lee and Fung, 2008; Feng et al., 2013; Botsyun et al., 2019a) all can affect δ18Opp. The effect of these complicating processes on isotopic lapse rates range from insignificant to major, up to the extent of flattening or inverting the paleo isotopic lapse rate relative to the modern (e.g., Moran et al., 2007). Thus, despite the intriguing potential of proxy records to reconstruct the paleoelevation history of orogens, the interpretation of the data requires several assumptions about the congruence between modern and paleo atmospheric processes and climate which may, or may not, be correct (e.g., Starke et al., 2020). Therefore, a rigorous interpretation of paleoelevations requires evaluating these assumptions on a case-by-case basis, which in practice is not always simple. One promising methods to address these concerns is an application of climate models for reconstructing time-specific temperatures, δ18Op and isotopic lapse rates for adequate interpretation of δ18Oc signal.
Basic Principles of Isotope Tracking Climate Models… in a Nutshell
Climate (and paleoclimate) models are typically based on a set of governing primitive equations for conservation of mass, energy, and momentum. These equations are solved simultaneously by discretizing (gridding) the surface of the earth and atmosphere and applying a numerical model. If only an atmospheric General Circulation Model (GCM) is used (as opposed to a coupled ocean-atmosphere GCM) then solving these equations also requires application of boundary conditions such as sea-surface temperatures, solar radiation, and also material properties of the atmosphere and land surface such as greenhouse gas concentrations, topography, vegetation/ice cover, sea-ice extent, soil properties, etc. Given the physics-based approach of these models, they can be used for diverse purposes including weather prediction/forecasts, future climate predictions, or paleoclimate predictions so long as the appropriate boundary conditions, and land surface properties are known. Thus, paleoclimate modeling studies that are applied to paleoelevation proxy data interpretation are “simply” applying well used and tested climate models from the atmospheric sciences community, with the caveat that they require a means for predicting the isotopic composition of precipitation, discussed next.
Isotopic Modeling and (Non-)stationarity of δ18Op-Lapse Rate
Climate modeling has been extensively applied to study the impact of mountain surface uplift on atmospheric physics and dynamics (Ruddiman and Kutzbach, 1989; Broccoli and Manabe, 1992; Sepulchre et al., 2006). The implementation of water isotope tracking in simple climate models (Dansgaard, 1964) and further into GCMs (Joussaume et al., 1984) enabled calculation of the atmospheric, land surface, and topographic processes impacting δ18Op. Since the development of iGCMs, a diverse range of studies have been conducted documenting the sensitivity of δ18Op to climate change resulting from variations in pCO2 (Poulsen et al., 2007; Poulsen and Jeffery, 2011), sea surface temperatures (SSTs) (Sturm et al., 2007), sea level variations (Poulsen et al., 2007), and paleogeography (Sewall and Fricke, 2013; Roe et al., 2016; Botsyun et al., 2019a), and much more (Risi, 2009; Sturm et al., 2010 and references herein).
The previous iGCM based studies have shown, among other things, that present-day isotopic lapse rates are not always suitable for paleoelevation reconstructions. First, modeling studies have shown spatial non-stationarity of the isotopic lapse rates. For example, Galewsky (2009) uses an idealized, fully nonlinear atmospheric model to evaluate the effects of stratified atmospheric flows on orographic precipitation isotopic ratios, with the aim of improving the foundation for interpreting proxy records used for paleoelevation studies. He showed that existing Rayleigh condensation models are strictly applicable for paleoaltimetry studies only when the atmosphere is not stratified. Changes in climate or the horizontal terrain aspect ratio can change precipitation isotopic ratios at least as much as changes in surface elevation (Galewsky, 2009), suggesting that proxy records must be interpreted within a broad context of climate variability and landscape evolution.
Second, iGCMs have demonstrated the impact of global climate change on δ18Op and the non-stationarity of isotopic lapse rates through time. For example, in climate simulations with elevated CO2 concentrations, the middle troposphere undergoes preferential moistening and warming, which weakens vertical stratification and gives rise to shallower lapse rates (Poulsen and Jeffery, 2011). Li et al. (2016) identified spatial and temporal variations in simulated isotopic lapse rates for different high relief zones along the flanks of the Tibetan plateau during Late Quaternary climates. They found that the precipitation weighted annual mean δ18Op lapse rate at the Himalaya is about 0.4‰/km larger during the Middle Holocene and 0.2‰/km smaller during the Last Glacial Maximum than during pre-industrial times. These changes are large enough to impact the interpretation of proxy data (Li et al., 2016). Furthermore, results from Botsyun et al. (2020) predict a similar magnitude (4–5‰ relative to present day) of change in δ18Op resulting from the topographic development of the European Alps and paleoclimate change during the Pliocene and the Last Glacial Maximum.
Third, modeling studies show that regional (e.g. continental or subcontinental scale) climate changes result from mountain surface uplift and also have a potential impact on δ18Op. For example, previous work has shown that Andean Plateau surface uplift resulted in the onset of convective precipitation (Poulsen et al., 2010) following a switch of the prevalent moisture source and associated transport paths from the South Pacific to the Equatorial Atlantic and initiation/strengthening of the South American Low Level Jet (Ehlers and Poulsen, 2009; Insel et al., 2010). These studies suggest that high precipitation rates enhance the isotope amount effect, leading to more negative δ18Op at high elevations and, thus, increases in the isotopic lapse rate. In agreement with this result, experiments with reduced elevations of the Andean Plateau by Insel et al. 2012 show that use of present-day isotopic lapse rate can lead to an underestimation of surface elevation by 2,000 m.
Regional climate changes associated with surface uplift have also been shown to impact δ18Op across the North American Cordillera (Feng et al., 2013) and Tibetan Plateau (Botsyun et al., 2016; Shen and Poulsen, 2019). For both these areas it was suggested that neither the common assumption that isotopic fractionation occurs primarily through rainout following Rayleigh distillation, nor the application of modern empirical δ18Op lapse rates to past environments, are valid. Shifts in atmospheric processes, including shifts in local precipitation type between convective and large-scale rain and between rain and snow, intensification of low-level vapor recycling particularly on leeward slopes, development of air mass mixing and changes in wind direction and moisture source changes all contributed to the deviations of the isotopic signal.
Finally, the robustness of applying the present-day isotopic lapse rate to interpret Eocene Tibetan Plateau proxy data was recently tested in work of Botsyun et al. (2019a). In this study, an iGCM with time specific Eocene boundary conditions was applied to assess the influence of changing Eocene paleogeography on climate and δ18Op signals. The authors found that a combination of increased convective precipitation, a mixture of air masses of different origin, widespread aridity, and an intensified water recycling resulted in a reversed isotopic lapse rate across the southern flank of the Tibetan Plateau. These processes resulted in the most negative δ18Op occurring over northern India and increased δ18Op northward (Figure 1). Taken together, these results indicate that standard stable isotope paleoaltimetry methods are not applicable in Eocene Asia and an alternative method, which takes into account paleo isotopic lapse rates, should be developed.
FIGURE 1. Schematic representation of the processes in the atmosphere and over land that potentially contribute to a δ18Op signal resulting in: (A) an isotopic lapse rate according to a Rayleigh distillation model for an isolated air parcel; and (B) flattened or inversed isotopic lapse rate. Examples of proxy data interpretations based on iGCM results.
In recent years, iGCMs have been extensively applied for δ18O proxy interpretation by using model-derived: 1) δ18Op lapse rates for various elevation scenarios; and 2) δ18Oc and soil temperatures for time-specific simulations including various paleogeographies and paleotopographies. Fan et al. (2017) and Gao and Fan (2018) used an iGCM-derived δ18Op-elevation gradient from Feng et al. (2013) to constrain the paleoelevations of the Cordillera orogenic system (Uinta Mountains, United States). This approach allows them to consider vapor mixing in the mountain flank and the warm climate during the Paleogene. Sundell et al. (2019) interpreted observed δ18Op across the Peruvian central Andes using different lapse rates including the thermodynamically-derived non-linear model from Rowley et al. (2001), an empirical approach from Quade et al. (2007), and reconstructed isotopic lapse rate using a climate model. Experiments with 25, 50, 75, and 100% of Andean elevations from Insel et al. (2012) and experiments with 50 and 100% Andean elevation from Poulsen et al. (2010) together with climate correction are used for the data interpretation.
Botsyun et al. (2019a) designed experiments with variable Tibetan Plateau elevations and Eocene paleogeography to compare simulated δ18Oc values with previously published data. These sensitivity tests identified that the best model-data fit corresponds to a low Tibetan Plateau in the Eocene. The authors showed that quantitative estimates of paleoaltimetry from stable isotopes depend on the outcome of a climate model, which in turn used paleoaltimetry, as input. This means that paleoaltimetry and paleoclimate models need to be iteratively tuned until the input elevation and the output relations between isotopic ratio and elevation are mutually consistent (Figure 2).
FIGURE 2. Principal scheme showing application of a paleoclimate model for paleoelevation reconstructions. t.f.—transfer function.
Limits of climate modeling technique to paleoaltimetry studies
From a stable isotope paleoaltimetry point of view the application of iGCMs have several limitations. The key caveats of climate modeling approach are linked to: 1) the limited ability of the model to simulate δ18O patterns over regions with complex topography due to the complexity of natural processes and limited scope of model parameterizations; 2) uncertainties in precipitation and δ18O changes in monsoonal areas; 3) limited computational resources for high-resolution fully-coupled ocean-atmosphere simulations; 4) uncertainties in the choice of boundary conditions (e.g. Figure 2), 5) relatively low model spatial resolution, and 6) a “social” challenge whereby the most fruitful applications involve a collaboration between proxy-oriented and model-oriented scientists.
Examples supporting the first point are as follows. Convective processes are known to have a significant imprint on the stable oxygen isotope composition of precipitation, especially in the tropics (Risi et al., 2008). However, convective processes in a iGCM strongly depend on the parameterization applied (Hourdin et al., 2006). Large uncertainties also exist in iGCMs depending on the representation of water stable isotopes parameterizations used for the isotopic exchanges between vapor and rain droplets during descent and the degree of partial reevaporation (Lee and Fung, 2008; Risi et al., 2010). Risi et al. (2010) demonstrated a sensitivity of modeling results to evaporation, especially in dry regions where this process is prevalent. In addition, the isotopic equilibration time of raindrops is drop-size dependent, with smaller equilibration time for small raindrops (Lee and Fung, 2008). Given this, if raindrop size is not taken into account an underestimation of equilibration times for heavy rains with big raindrops could be expected. Additional potential limitations of using iGCMs include complexity associated with realistic representation of soils, evapotranspiration schemes, geochemistry, vegetation, and clouds (Risi et al., 2010; Haese et al., 2013; Rio and Hourdin, 2015). These challenges are active areas of research in the climate modeling community and improvements in model parameterizations are always occurring.
The response of the hydrologic cycle to the climate forcing in monsoon areas vary depending on the model applied and, even for future climates, monsoons amplification has not been shown as a strong feature (Wang et al., 2020). Furthermore, the interpretation of water isotope variations and their link to monsoon is not straightforward even for recent climates (e.g., Middle Holocene) when there is no problem of mountain elevation change (Lin et al., 2019). Further uncertainties are associated with the application of atmosphere-only climate model with no feedbacks between the atmosphere, vegetation and the ocean. Use of modeling techniques with SSTs and sea ice extension prescribed from coupled simulations to run a separate high resolution atmospheric model suffer from not taking into account ocean variability and its feedback to the atmosphere. Moreover, vegetation can also play an important role in isotopic changes through not only evapotranspiration changes, but also through changing surface albedo feedback. Given this coupled high-resolution ocean-atmosphere and dynamic vegetation modeling approaches are preferred, but extremely challenging to conduct. However, usage of high-resolution atmosphere-only iGCMs that use previous coarser resolution coupled ocean-atmosphere models for boundary conditions have been shown to be a good trade-off between the high-resolution necessary to capture observed spatial, seasonal, and daily variations of δ18Op and computation resources required (Botsyun et al., 2019b).
Furthermore, paleoclimate studies are sensitive to the paleogeographical reconstruction used (e.g., Baatsen et al., 2016). Producing paleogeographic reconstructions for climate models is time-consuming and elaborate. For example, paleoclimate models frequently use reconstructions where the latest state-of-the-art of plate tectonic, crustal and paleomagnetic reconstructions, paleotopography, paleobathymetry, or vegetation distribution have not yet been incorporated. The availability of realistic Cenozoic boundary conditions, such as sea-surface temperatures, is also limited, unless a time-consuming coupled ocean-atmosphere model is used. Finally, coarse model resolutions used in iGCMs to reduce simulation time often prevent accurate representation of smaller-scale topographic features such that explicit simulation of air masses ascending mountain topography are only grossly represented (Hourdin et al., 2006; Haese et al., 2013). Despite the above caveats, numerous studies have applied iGCMs to capture regional changes in δ18Op through time, and comparison of model predicted δ18Op to observations (e.g., Global Network of Isotopes in Precipitation (GNIP) stations data) to understand model limitations are commonly conducted when applying models to new locations.
When are state-of-the-art climate models useful and when can simpler approaches be applied?
Despite the previously mentioned limitations associated with iGCMs, applications of climate models to enhance proxy interpretations are gaining success (e.g., van Hinsbergen and Boschman, 2019). Modeling studies have shown that simpler (i.e., non-iGCM based) paleoaltimetry approaches are valid when the Rayleigh distillation process is dominant. However, iGCMs have provided several examples of when atmospheric processes lead to δ18Op variations that are different, and sometimes opposite, of those estimated using Rayleigh distillation models of moist adiabatic condensation. For example, Shen and Poulsen (2019) used the ECHAM5-wiso model to show that Rayleigh distillation is only prevailing in the monsoonal regions of the Himalayas when the mountains are high. In contrast, when orogen topography is lower, local surface recycling and convective processes become important. This is due to weakened forced ascent causing weaker Asian monsoons. In this example, applying a standard paleoaltimetry approach is not valid. Similarly, over the North America Cordillera changes in atmospheric processes have been shown to impact δ18Op and violate the common assumption that isotopic fractionation occurs through rainout following Rayleigh distillation (Feng et al., 2013). In this example, the atmospheric changes included shifts in local precipitation type between convective and large-scale rain and also between rain and snow, intensification of low-level vapor recycling on leeward slopes, and development of air mass mixing and changes in wind direction and moisture source. Given the uncertainty of when, or when not, a one-dimensional Rayleigh distillation model adequately represents moisture transport, the significance of lateral vapor mixing, or temporal changes in the isotopic lapse rate, the application of iGCMs is recommended (e.g., Li et al., 2016) as an initial step to characterize these processes rather than assuming they are insignificant.
There are, however, observational (non-iGCM based) techniques that have emerged to help circumvent the effects of changing atmospheric processes on paleoaltimetry data. One such technique is the δ-δ approach. This approach assumes that a baseline sets of samples collected from a low elevation stratigraphic succession provides a reference frame for interpreting higher elevation samples. The implicit assumption is that any regional or global climate change will affect stable isotopes records from both (low and high elevation) localities equally if looking at the same time interval. If true, then the difference between the two records can be interpreted as resulting from elevation change, thereby enabling a paleoelevation reconstruction (e.g., Campani et al., 2012; Gébelin et al., 2012; Mulch, 2016; Pingel et al., 2020). While the δ-δ approach holds promise and there are geographic locations where this may work well, the approach does require that a suitable baseline location has been chosen–which is difficult to do without prior knowledge. In recent work by Botsyun et al. (2020), an iGCM was used to explore the δ18Op signal of Cenozoic topographic and paleoclimatic change in the European Alps. The results of their study suggest that the δ-δ approach is applicable for the case of the European Alps under the condition that the low-elevation reference section measured hasn’t experienced climate change, and that it is located far enough away and upwind from the area of interest to not experience climate changes associated with orogenic uplift.
In summary, when working with proxy records from “non-simple” locations where Rayleigh distillation assumptions may not be valid, the modern paleoelevation community is well equipped to make a step forward to more advanced reconstructions involving isotope enabled climate modeling. However, doing this requires a collaboration between proxy and climate modeling communities to evaluate how simple, or complex, the δ18Op response is to changing atmospheric processes. In this review, we have highlighted how application of climate models in conjunction with geological/geochemical data provides a powerful tool to incorporate climatic change effects into the analysis of paleoaltimetry work. Moreover, paleobotanic data and stable oxygen isotope paleoaltimetry can be reconciled using climate models. Further model-data comparison studies that combine multiple proxy types, together with iGCM modeling efforts testing the impact of complex topography structures as well as atmospheric parameterizations on δ18O are paving the way for more accurate and robust paleoelevation estimates in the decades to come.
Author Contributions
SB and TAE conceived and wrote the manuscript.
Funding
This study was supported as part of the German priority research program Mountain Building Processes in 4D (MB4D) grant (EH329/18-1 and 23-1) to TAE.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Baatsen, M., van Hinsbergen, D. J. J., von der Heydt, A. S., Dijkstra, H. A., Sluijs, A., Abels, H. A., et al. (2016). Reconstructing geographical boundary conditions for palaeoclimate modelling during the Cenozoic. Clim. Past 12, 1635–1644. doi:10.5194/cp-12-1635-2016
Botsyun, S., Sepulchre, P., Risi, C., and Donnadieu, Y. (2016). Impacts of Tibetan Plateau uplift on atmospheric dynamics and associated precipitation δ18O. Clim. Past 12, 1401–1420. doi:10.5194/cp-12-1401-2016
Botsyun, S., Sepulchre, P., Donnadieu, Y., Risi, C., Licht, A., and Caves Rugenstein, J. K. (2019a). Revised paleoaltimetry data show low Tibetan Plateau elevation during the Eocene. Science 363 (6430), eaaq1436. doi:10.1126/science.aaq1436
Botsyun, S., Sepulchre, P., Donnadieu, Y., Risi, C., Licht, A., and Caves Rugenstein, J. K. (2019b). Response to comment on “Revised paleoaltimetry data show low Tibetan Plateau elevation during the Eocene”. Science, 365 (6459), eaax8990. doi:10.1126/science.aax8990
Botsyun, S., Ehlers, T. A., Mutz, S. G., Methner, K., Krsnik, E., and Mulch, A. (2020). Opportunities and challenges for paleoaltimetry in “small” orogens: Insights from the European Alps. Geophys. Res. Lett., e2019GL086046. doi:10.1029/2019GL086046
Broccoli, A. J., and Manabe, S. (1992). The effects of orography on midlatitude northern hemisphere dry climates. J. Clim. 5, 1181–1201. doi:10.1175/1520-0442(1992)005<1181:TEOOOM>2.0.CO;2
Campani, M., Mulch, A., Kempf, O., Schlunegger, F., and Mancktelow, N. (2012). Miocene paleotopography of the Central Alps. Earth Planet. Sci. Lett. 337–338, 174–185. doi:10.1016/j.epsl.2012.05.017
Cassel, E. J., Breecker, D. O., Henry, C. D., Larson, T. E., and Stockli, D. F. (2014). Profile of a paleo-orogen: High topography across the present-day basin and range from 40 to 23 Ma. Geology 42, 1007–1010. doi:10.1130/G35924.1
Chamberlain, C. P., Mix, H. T., Mulch, A., Hren, M. T., Kent-Corson, M. L., Davis, S. J., et al. (2012). The Cenozoic climatic and topographic evolution of the western North American Cordillera. Am. J. Sci. 312, 213–262. doi:10.2475/02.2012.05
Chamberlain, C. P., Winnick, M. J., Mix, H. T., Chamberlain, S. D., and Maher, K. (2014). The impact of Neogene grassland expansion and aridification on the isotopic composition of continental precipitation. Global Biogeochem. Cycles 28, 992–1004. doi:10.1002/2014GB004822
Dansgaard, W. (1964). Stable isotopes in precipitation. Tellus 16, 436–468. doi:10.1111/j.2153-3490.1964.tb00181.x
Ding, L., Xu, Q., Yue, Y., Wang, H., Cai, F., and Li, S. (2014). The Andean-type Gangdese Mountains: Paleoelevation record from the Paleocene-Eocene Linzhou Basin. Earth Planet Sci. Lett. 392, 250–264. doi:10.1016/j.epsl.2014.01.045
Ehlers, T. A., and Poulsen, C. J. (2009). Influence of Andean uplift on climate and paleoaltimetry estimates. Earth Planet Sci. Lett. 281, 238–248. doi:10.1016/j.epsl.2009.02.026
Fan, M., Constenius, K. N., and Dettman, D. L. (2017). Prolonged high relief in the northern Cordilleran orogenic front during middle and late Eocene extension based on stable isotope paleoaltimetry. Earth Planet Sci. Lett. 457, 376–384. doi:10.1016/j.epsl.2016.10.038
Feng, R., Poulsen, C. J., Werner, M., Chamberlain, C. P., Mix, H. T., and Mulch, A. (2013). Early Cenozoic evolution of topography, climate, and stable isotopes in precipitation in the North American Cordillera. Am. J. Sci. 313, 613–648. doi:10.2475/07.2013.01
Fiorella, R. P., Poulsen, C. J., Zolá, R. S. P., Barnes, J. B., Tabor, C. R., and Ehlers, T. A. (2015). Spatiotemporal variability of modern precipitation δ18O in the central Andes and implications for paleoclimate and paleoaltimetry estimates. J. Geophys. Res. Atmos., 1–27. doi:10.1002/2014JD022893
Galewsky, J. (2009). Orographic precipitation isotopic ratios in stratified atmospheric flows: Implications for paleoelevation studies. Geology 37, 791–794. doi:10.1130/G30008A.1
Gao, M., and Fan, M. (2018). Depositional environment, sediment provenance and oxygen isotope paleoaltimetry of the early Paleogene greater Green River Basin, southwestern Wyoming, U.S.A. Am. J. Sci. 318, 1018–1055. doi:10.2475/10.2018.02
Garzione, C. N., McQuarrie, N., Perez, N. D., Ehlers, T. A., Beck, S. L., Kar, N., et al. (2017). Tectonic evolution of the Central Andean plateau and implications for the growth of plateaus. Annu. Rev. Earth Planet Sci. 45, 529–559. doi:10.1146/annurev-earth-063016-020612
Gat, J. R. (1996). Oxygen and hydrogen isotopes in the hydrologic cycle. Annu. Rev. Earth Planet Sci. 24, 225–262. doi:10.1146/annurev.earth.24.1.225
Gébelin, A., Mulch, A., Teyssier, C., Page Chamberlain, C., and Heizler, M. (2012). Coupled basin-detachment systems as paleoaltimetry archives of the western North American Cordillera. Earth Planet Sci. Lett. 335–336, 36–47. doi:10.1016/j.epsl.2012.04.029
Gébelin, A., Mulch, A., Teyssier, C., Jessup, M. J., Law, R. D., and Brunel, M. (2013). The Miocene elevation of mount everest. Geology 41, 799–802. doi:10.1130/G34331.1
Gonfiantini, R., Roche, M. A., Olivry, J. C., Fontes, J. C., and Zuppi, G. M. (2001). The altitude effect on the isotopic composition of tropical rains. Chem. Geol. 181, 147–167. doi:10.1016/S0009-2541(01)00279-0
Haese, B., Werner, M., and Lohmann, G. (2013). Stable water isotopes in the coupled atmosphere-land surface model ECHAM5-JSBACH. Geosci. Model Dev. (GMD) 6, 1463–1480. doi:10.5194/gmd-6-1463-2013
Hoke, G. D., Liu-Zeng, J., Hren, M. T., Wissink, G. K., and Garzione, C. N. (2014). Stable isotopes reveal high southeast Tibetan Plateau margin since the Paleogene. Earth Planet Sci. Lett. 394, 270–278. doi:10.1016/j.epsl.2014.03.007
Hourdin, F., Musat, I., Bony, S., Braconnot, P., Codron, F., Dufresne, J. L., et al. (2006). The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection. Clim. Dynam. 27, 787–813. doi:10.1007/s00382-006-0158-0
Insel, N., Poulsen, C. J., and Ehlers, T. A. (2010). Influence of the Andes Mountains on South American moisture transport, convection, and precipitation. Clim. Dynam. 35, 1477–1492. doi:10.1007/s00382-009-0637-1
Insel, N., Poulsen, C. J., Ehlers, T. A., and Sturm, C. (2012). Response of meteoric δ18O to surface uplift ‐ implications for Cenozoic Andean Plateau growth. Earth Planet. Sci. Lett. 317–318, 262–272. doi:10.1016/j.epsl.2011.11.039
Joussaume, S., Sadourny, R., and Jouzel, J. (1984). A general circulation model of water isotope cycles in the atmosphere. Nature 311, 24–29. doi:10.1038/311024a0
Kim, S. T., and O’Neil, J. R. (1997). Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochem. Cosmochim. Acta 61, 3461–3475. doi:10.1016/S0016-7037(97)00169-5
Lee, J. E., and Fung, I. (2008). “Amount effect” of water isotopes and quantitative analysis of post-condensation processes. Hydrol. Process. 22, 1–8. doi:10.1002/hyp.6637
Li, J., Ehlers, T. A., Werner, M., Mutz, S. G., Steger, C., and Paeth, H. (2016). Late quaternary climate, precipitation δ18O, and Indian monsoon variations over the Tibetan Plateau. Earth Planet Sci. Lett. doi:10.1016/j.epsl.2016.09.031
Lin, Y., Ramstein, G., Wu, H., Rani, R., Braconnot, P., Kageyama, M., et al. (2019). Mid-Holocene climate change over China: model–data discrepancy. Clim. Past, 15 (4), 1223–1249. doi:10.5194/cp-15-1223-2019
Moran, T. A., Marshall, S. J., Evans, E. C., and Sinclair, K. E. (2007). Altitudinal gradients of stable isotopes in lee-slope precipitation in the Canadian rocky mountains. Arctic Antarct. Alpine Res. 39, 455–467. doi:10.1657/1523-0430(06-022)[MORAN]2.0.CO;2
Mulch, A. (2016). Stable isotope paleoaltimetry and the evolution of landscapes and life. Earth Planet. Sci. Lett. 433, 180–191. doi:10.1016/j.epsl.2015.10.034
Mulch, A., Uba, C. E., Strecker, M. R., Schoenberg, R., and Chamberlain, C. P. (2010). Late Miocene climate variability and surface elevation in the central Andes. Earth Planet Sci. Lett. 290, 173–182. doi:10.1016/j.epsl.2009.12.019
Pausata, F. S. R., Li, C., Wettstein, J. J., Kageyama, M., and Nisancioglu, K. H. (2011). The key role of topography in altering North Atlantic atmospheric circulation during the last glacial period. Clim. Past 7, 1089–1101. doi:10.5194/cp-7-1089-2011
Pingel, H., Strecker, M. R., Mulch, A., Alonso, R. N., Cottle, J., and Rohrmann, A., (2020). Late Cenozoic topographic evolution of the Eastern Cordillera and Puna Plateau margin in the southern Central Andes (NW Argentina). Earth Planet Sci. Lett. 535, 116112. doi:10.1016/j.epsl.2020.116112
Poulsen, C. J., and Jeffery, M. L. (2011). Climate change imprinting on stable isotopic compositions of high-elevation meteoric water cloaks past surface elevations of major orogens. Geology 39, 595–598. doi:10.1130/G32052.1
Poulsen, C. J., Pollard, D., and White, T. S. (2007). General circulation model simulation of the δ18O content of continental precipitation in the middle Cretaceous: a model-proxy comparison. Geology 35, 199–202. doi:10.1130/G23343A.1
Poulsen, C. J., Ehlers, T. A., and Insel, N. (2010). Onset of convective rainfall during gradual late Miocene rise of the central Andes. Science 328 (5977), 490–493. doi:10.1126/science.1185078
Quade, J., Garzione, C., and Eiler, J. (2007). Paleoelevation reconstruction using pedogenic carbonates. Rev. Mineral. Geochem. 66, 53–87. doi:10.2138/rmg.2007.66.3
Rio, C., and Hourdin, F. (2015). A thermal plume model for the convective boundary layer: Representation of cumulus clouds. J. Atmos. Sci. 72. doi:10.1175/2007JAS2256.1
Risi, C., Bony, S., and Vimeux, F. (2008). Influence of convective processes on the isotopic composition (δ18O and δD) of precipitation and water vapor in the tropics: 2. Physical interpretation of the amount effect. J. Geophys. Res. Atmos. 113, 1–12. doi:10.1029/2008JD009943
Risi, C., Bony, S., Vimeux, F., and Jouzel, J. (2010). Water-stable isotopes in the LMDZ4 general circulation model: Model evaluation for present-day and past climates and applications to climatic interpretations of tropical isotopic records. J. Geophys. Res. Atmos. 115, 1–27. doi:10.1029/2009JD013255
Risi, C., Noone, D., Frankenberg, C., and Worden, J. (2013). Role of continental recycling in intraseasonal variations of continental moisture as deduced from model simulations and water vapor isotopic measurements. Water Resour. Res. 49, 4136–4156. doi:10.1002/wrcr.20312
Risi, C. (2009). Les isotopes stables de l’eau: applications à l’étude du cycle de l’eau et des variations du climat. PhD thesis. Université Pierre et Marie Curie.
Roe, G. H., Ding, Q., Battisti, D. S., Molnar, P., Clark, M. K., and Garzione, C. N. (2016). A modeling study of the response of Asian summertime climate to the largest geologic forcings of the past 50 Ma. J. Geophys. Res. Atmos. 121, 5453–5470. doi:10.1002/2015JD024370
Rowley, D. B., and Currie, B. S. (2006). Palaeo-altimetry of the late Eocene to Miocene Lunpola basin, central Tibet. Nature 439, 677–681. doi:10.1038/nature04506
Rowley, D. B., and Garzione, C. N. (2007). Stable isotope-based paleoaltimetry. Annu. Rev. Earth Planet Sci. 35, 463–508. doi:10.1146/annurev.earth.35.031306.140155
Rowley, D. B., Pierrehumbert, R. T., and Currie, B. S. (2001). A new approach to stable isotope-based paleoaltimetry: implications for paleoaltimetry and paleohypsometry of the High Himalaya since the Late Miocene. Earth Planet Sci. Lett. 188, 253–268. doi:10.1016/S0012-821X(01)00324-7
Rozanski, K., Araguás-Araguás, L., and Gonfiantini, R. (1993). Isotopic patterns in modern global precipitation. Clim. Change Cont. Isot. Rec., 1–36. doi:10.1029/GM078p0001
Ruddiman, W. F., and Kutzbach, J. E. (1989). Forcing of late Cenozoic northern hemisphere climate by plateau uplift in southern Asia and the American West. J. Geophys. Res. Atmos. 94, 18409–18427.
Schneider, D. P., and Noone, D. C. (2007). Spatial covariance of water isotope records in a global network of ice cores spanning twentieth‐century climate change. J. Geophys. Res. Atmos. 112. doi:10.1029/2007JD008652
Sepulchre, P., Ramstein, G., Fluteau, F., Schuster, M., Tiercelin, J. J., and Brunet, M. (2006). Tectonic uplift and Eastern Africa aridification. Science 313, 1419–1423. doi:10.1126/science.1129158
Sewall, J. O., and Fricke, H. C. (2013). Andean-scale highlands in the Late Cretaceous Cordillera of the North American western margin. Earth Planet Sci. Lett. 362, 88–98. doi:10.1016/j.epsl.2012.12.002
Shen, H., and Poulsen, C. J. (2019). Precipitation δ18O on the Himalaya–Tibet orogeny and its relationship to surface elevation. Clim. Past 15, 169–187. doi:10.5194/cp-15-169-2019
Starke, J., Ehlers, T. A., and Schaller, M. (2020). Latitudinal effect of vegetation on erosion rates identified along western South America. Science 367, 1358–1361. doi:10.1126/science.aaz0840
Sturm, C., Hoffmann, G., and Langmann, B. (2007). Simulation of the stable water isotopes in precipitation over South America: Comparing regional to global circulation models. J. Clim. 20, 3730–3750. doi:10.1175/JCLI4194.1
Sturm, C., Zhang, Q., and Noone, D. (2010). An introduction to stable water isotopes in climate models: benefits of forward proxy modelling for paleoclimatology. Clim. Past 6, 115–129. doi:10.5194/cp-6-115-2010
Sundell, K. E., Saylor, J. E., Lapen, T. J., and Horton, B. K. (2019). Implications of variable late Cenozoic surface uplift across the Peruvian central Andes. Sci. Rep. 9, 1–12. doi:10.1038/s41598-019-41257-3
Takahashi, K., and Battisti, D. S. (2007). Processes controlling the mean tropical pacific precipitation pattern. Part I: The Andes and the eastern Pacific ITCZ. J. Clim. 20, 3434–3451. doi:10.1175/JCLI4198.1
van Hinsbergen, D. J. J., and Boschman, L. M. (2019). How high were these mountains? Science 363, 928–929. doi:10.1126/science.aaw7705
Keywords: paleoaltimetry, oxygen isotopes (δ18O), isotope-enabled GCM, paleoelevation reconstruction, stable oxygen isotope composition, tectonic-climatic interaction, proxy interpretation, atmosphere circulation modeling
Citation: Botsyun S and Ehlers TA (2021) How Can Climate Models Be Used in Paleoelevation Reconstructions?. Front. Earth Sci. 9:624542. doi: 10.3389/feart.2021.624542
Received: 31 October 2020; Accepted: 13 January 2021;
Published: 18 February 2021.
Edited by:
Heiko Pingel, University of Potsdam, GermanyReviewed by:
Gilles Ramstein, UMR8212 Laboratoire des Sciences du Climat et de l’Environnement, FranceMichael Hren, University of Connecticut, United States
Copyright © 2021 Botsyun and Ehlers. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Svetlana Botsyun, botsyun.svetlana@gmail.com