Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.

Stable isotope stacks from GRIP and GISP ice cores

Recent efforts to link the isotopic composition of snow in Greenland with meteorological and climatic parameters have indicated that relatively local information such as observed annual temperatures from coastal Greenland sites, as well as more synoptic scale features such as the North Atlantic Oscillation (NAO) and the temperature seesaw between Jakobshaven, Greenland, and Oslo, Norway, are significantly correlated with d18O and dD values from the past few hundred years measured in ice cores. In this study we review those efforts and then use a new record of isotope values from the Greenland Ice Sheet Project 2 and Greenland Ice Core Project sites at Summit, Greenland, to compare with meteorological and climatic parameters. This new record consists of six individual annually resolved isotopic records which have been average to produce a Summit stacked isotope record. The stacked record is significantly correlated with local Greenland temperatures over the past century (r=0.471), as well as a number of other records including temperatures and pressures from specific locations as well as temperature and pressure patterns such as the temperature seesaw and the North Atlantic Oscillation. A multiple linear regression of the stacked isotope record with a number of meteorological and climatic parameters in the North Atlantic region reveals that five variables contribute significantly to the variance in the isotope record: winter NAO, solar irradiance (as recorded by sunspot numbers), average Greenland coastal temperature, sea surface temperature in the moisture source region for Summit (30°-20°N), and the annual temperature seesaw between Jakobshaven and Oslo. Combined, these variables yield a correlation coefficient of r=0.71, explaining half of the variance in the stacked isotope record.

Age determination and paleomagnetic investigations on a sediment profile from laguna Potrok Aike

Paleomagnetic inclination, declination and relative paleointensity were reconstructed from the sediments of Laguna Potrok Aike in the framework of the International Continental scientific Drilling Program (ICDP) Potrok Aike maar lake Sediment Archive Drilling prOject (PASADO). Here we present the u-channel-based full vector paleomagnetic field reconstruction since 51.2 ka cal BP. The relative paleointensity proxy (RPI) was built by normalising the natural remanent magnetisation with the anhysteretic remanent magnetisation using the average ratio at 4 demagnetisation steps part of the ChRM interval (NRM/ARM10e40 mT). A grain size influence on the RPI was removed using a correction based on the linear relationship between the RPI and the median destructive field of the natural remanent magnetisation (MDFNRM). The new record is compared with other lacustrine and marine records and stacks from the mid- to high-latitudes of the Southern Hemisphere, revealing consistent millennial-scale variability, the identification of the Laschamp and possibly the Mono Lake geomagnetic excursions, and a direction swing possibly associated to the Hilina Pali excursion at 20 ka cal BP. Nonetheless, a global-scale comparison with other high-resolution records located on the opposite side of the Earth and with various dipole field references hint at a different behaviour of the geomagnetic field around southern South America at 46 ka cal BP.

Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

Constructing a stacked benthic d18O record

Composite stacks were constructed by superimposing 6 to 13 benthic foraminiferal d18O records covering the period 0-850 ka. An initial timescale for each core was established using radioisotopic age control points and assuming constant sedimentation rates between these points. The average of these records is our 13-core "untuned" stack. Next, we matched the 41 kyr component of each record individually to variations in Earth's obliquity. Four of the 13 records produced timescales that were inconsistent with one or more of the known radioisotopic ages. The nine remaining cores were averaged to create a "minimally tuned" stack. Six of the minimally tuned cores were assembled into a "tropical" stack. For each stack we estimated the uncertainty envelope from the standard deviation of the constituents. Spectral analysis of the three stacks indicates that benthic d18O is dominated by a 100 kyr oscillation that has a narrow spectral peak. The contribution of precession to the total variance is small when compared to prior results from planktic stacks.

Pliocene-Pleistocene stack of globally distributed benthic stable oxygen isotope records

We present a 5.3-Myr stack (the ''LR04'' stack) of benthic d18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic delta18O stack composed of more than three records to extend beyond 850 ka, and we use its improved signal quality to identify 24 new marine isotope stages in the early Pliocene. We also present a new LR04 age model for the Pliocene-Pleistocene derived from tuning the delta18O stack to a simple ice model based on 21 June insolation at 65 N. Stacked sedimentation rates provide additional age model constraints to prevent overtuning. Despite a conservative tuning strategy, the LR04 benthic stack exhibits significant coherency with insolation in the obliquity band throughout the entire 5.3 Myr and in the precession band for more than half of the record. The LR04 stack contains significantly more variance in benthic delta18O than previously published stacks of the late Pleistocene as the result of higher resolution records, a better alignment technique, and a greater percentage of records from the Atlantic. Finally, the relative phases of the stack's 41- and 23-kyr components suggest that the precession component of delta18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of d18O precession response changed suddenly at 1.6 Ma.