Validation and calibration of soil δ²H and brGDGTs along (E-W) and strike (N-S) of the Himalayan climatic gradient

Reconstructing the timing of mountain range uplift and the evolution of high-altitude plateaus is important when attempting to understand potential feedbacks between tectonics and climate at geological timescales. This requires proxies that are able to accurately reconstruct elevation during different time periods in the past. Often, the sensitivity of climatic parameters to elevation gradients, recorded in geological archives such as soils, is used to estimate paleoelevations. However, most proxies reflect an indirect response to elevation change, adding uncertainties to reconstructions. In this study, we aim to identify those sources of uncertainty with respect to elevation reconstructions and test if the combined application of two such proxies, i.e., stable isotopes (δ²H) of plant waxes in modern soils and surface waters and bacterial membrane lipids (brGDGTs) in soils, which can potentially reduce uncertainties in the estimation of (paleo-) elevation. We performed this study in four Himalayan catchments (from west to east: Sutlej, Alaknanda, Khudi, and Arun), of which each individual catchment is subject to a unique precipitation regime, relative influences of moisture sources, and vegetation cover. In total, we analyzed 275 surface water samples, 9 precipitation samples, 131 xylem water samples, and 60 soil samples, which were collected between 2009 and 2014. The following key observations were made: Soil nC₃₁-alkane δ²H values (δ²H_wax) in the Sutlej, Alaknanda, Khudi, and Arun generally record surface water δ²H values, confirming that the first-order control on the plant wax isotopic signature is precipitation δ²H and, therefore, the elevation in orogenic settings. We identified aridity as the factor that introduces scatter to this relationship. BrGDGT-derived Mean Annual Temperature (MAT) correlates in a statistically significant manner with sample site elevation and a 14-year annual average of remotely sensed land-surface temperature, showing that the main process influencing the brGDGT distribution is the adiabatic cooling of air. In an effort to combine these proxies to improve uncertainties in elevation reconstruction, elevations were inferred from both the δ²H_wax and brGDGT distributions. Arid, high elevation sites appear to underestimate actual sample site elevations using δ²H values while sites subject to high (>23–25 °C) annual temperatures overestimate the actual sample site elevation using brGDGT distributions. Elevations inferred from both proxies under such paleoclimatic conditions should be interpreted with caution. Elevations derived from the brGDGT distribution appear to most accurately reconstruct elevation. However, we show that the difference in elevation between the two proxies, described by the proposed Δ_Elevation parameter, can provide information on the hydrological conditions of the soil's depositional environment. In conclusion, we emphasize that knowledge of the sample site's climatic conditions are essential to reconstruct elevation from paleoarchives. In particular, knowledge of moisture availability and annual air temperatures are important, as these have been found to cause the largest scatter in the observed data.