Mapping spectroscopic uncertainties into prospective methane retrieval errors from Sentinel-5 and its precursor

Sentinel-5 (S5) and its precursor (S5P) are future European satellite missions aiming at global monitoring of methane (CH4) column-average dry air mole fractions (XCH4). The spectrometers to be deployed onboard the satellites record spectra of sunlight backscattered from the Earth's surface and atmosphere. In particular, they exploit CH4 absorption in the shortwave infrared spectral range around 1.65 mu m (S5 only) and 2.35 mu m (both S5 and S5P) wavelength. Given an accuracy goal of better than 2% for XCH4 to be delivered on regional scales, assessment and reduction of potential sources of systematic error such as spectroscopic uncertainties is crucial. Here, we investigate how spectroscopic errors propagate into retrieval errors on the global scale. To this end, absorption spectra of a ground-based Fourier transform spectrometer (FTS) operating at very high spectral resolution serve as estimate for the quality of the spectroscopic parameters. Feeding the FTS fitting residuals as a perturbation into a global ensemble of simulated S5- and S5P-like spectra at relatively low spectral resolution, XCH4 retrieval errors exceed 0.6% in large parts of the world and show systematic correlations on regional scales, calling for improved spectroscopic parameters.

Retreat of top cliff of Kurungnakh Island, Lena Delta, Siberia, Russia, 2010-2014, with links to shapefiles

Thawing-induced cliff top retreat in permafrost landscapes is mainly due to thermo-erosion. Ground-ice-rich permafrost landscapes are specifically vulnerable to thermo-erosion and may show high degradation rates. Within the HGF Alliance Remote Sensing and the FP7 PAGE21 permafrost programs we investigated how SAR and optical remote sensing can contribute to the monitoring of erosion rates of ice-rich cliffs in Arctic Siberia (Lena Delta, Russia). We produced two different vector products: i) Intra-annual cliff top retreat based on TerraSAR-X (TSX) satellite data (2012-2014): High-temporal resolution time series of TSX satellite data allow the inter-annual and intra-annual monitoring of the upper cliff-line retreat also under bad weather conditions and continuous cloud coverage. This published SAR product contains the retreating upper cliff lines of a 1.5 km long part of eroding ice-rich coast of Kurungnakh Island in the central Lena Delta. The upper cliff line was mapped using a thresholding approach for images acquired in the years 2012, 2013 and 2014 for the months June (2013, 2014), July (2013, 2014), August (2012, 2013, 2014) and September (2013, 2014). The cliff top retreat vector product is called 'upper_cliff_TerraSAR-X'. While the 2014 cliff lines show a clear retreat of 2 to 3 m/month, the cliff top lines for 2012 and 2013 are not chronologically ordered. However, lines from the end of the season of a year are always close to the lines from the beginning of the next summer season, indicating low cliff retreat in winter. ii) 4-year cliff top retreat based on optical satellite data (2010-2014): Long-term cliff top retreat could be assessed with two high-spatial resolution optical satellite images (GeoEye-1, 2010-08-05 and Worldview-1, 2014-08-19). The cliff top retreat vector product is called 'upper_cliff_optical'. Results: The long-term cliff top retreat derived from optical satellite data are 35 m cliff retreat within 4 years. The higher-temporal resolution SAR data equivalently show long-term rates of 18 m within 2 years and nearly now degradation activities in winter but maximum erosion rates in summer months.The Intra-seasonal cliff top retreat lines from 2014 show a rate of 2 to 3 m per month.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Physical properties of sediment cores from the Labrador Sea

Sediment core logs from six sediment cores in the Labrador Sea show millennial-scale climate variability during the last glacial by recording all Heinrich events and several major Dansgaard-Oeschger cycles. The same millennial-scale climate change is documented for surface-water d18O records of Neogloboquadrina pachyderma (left coiled); hence the surface-water d18O record can be derived from sediment core logging by means of multiple linear regression, providing a paleoclimate proxy record at very high temporal resolution (70 yrs). For the Labrador Sea, sediment core logs contain important information about deep-water current velocities and also reflect the variable input of IRD from different sources as inferred from grain-size analysis, benthic d18O, the relation of density and p-wave velocity, and magnetic susceptibility. For the last glacial, faster deep-water currents which correspond to highs in sediment physical properties, occurred during iceberg discharge and lasted for a several centuries to a few millennia. Those enhanced currents might have contributed to increased production of intermediate waters during times of reduced production of North Atlantic Deep Water. Hudson Strait might have acted as a major supplier of detrital carbonate only during lowered sea level (greater ice extent). During coldest atmospheric temperatures over Greenland, deep-water currents increased during iceberg discharge in the Labrador Sea, then surface water freshened shortly after, while the abrupt atmospheric temperature rise happened after a larger time lag of >=1 kyr. The correlation implies a strong link and common forcing for atmosphere, sea surface, and deep water during the last glacial at millennial time scales but decoupling at orbital time scales.

Osmium isotope record of Late Pleistocene sediments

We report high temporal resolution osmium isotopes records of bulk sediment and sediment leachates from DSDP Site 480 (Gulf of California) over the last 30 ka; from ODP Site 849 (Eastern equatorial Pacific) from the last 200 ka and from ODP Site 1002C (Cariaco Basin) across the 9-17 ka time interval in order to critically evaluate claims of a global 10% shift in the 187Os/188Os of seawater from glacial to interglacial intervals. We use organic-rich continental margin sites and carbonate-rich pelagic sites to isolate the temporal variations of the osmium seawater isotopic composition. Our results reveal that variations in 187Os/188Os fail to correlate with global changes in temperature across glacials/interglacials cycles as previously claimed. Instead, these results indicate differences of a few percent in the measured 187Os/188Os between each oceanic basin. We argue that these differences strongly suggest that seawater is not well homogenized with respect to its Os isotope composition.

Meteorological measurements from the Fourcade and Polar Club Glacier, Warszawa Icefield, King George Island, West Antarctica

The South Shetland Islands are located at the northern tip of the AP which is among the fastest warming regions on Earth. The islands are especially vulnerable to climate change due to their exposure to transient low-pressure systems and their maritime climate. Surface air temperature increases (2.5K in 50 years) are concurrent with retreating glacier fronts, an increase in melt areas, ice surface lowering and rapid break-up and disintegration of ice shelves. We have compiled a unique meteorological data set for the King George Island (KGI)/Isla 25 de Mayo, the largest of the South Shetland Islands. It comprises high-temporal resolution and spatially distributed observations of surface air temperature, wind directions and wind velocities, as well as glacier ice temperatures in profile with a fully equipped automatic weather station on the Warszawa Icefield, from November 2010 and ongoing. In combination with two long-term synoptic datasets (40 and 10 years, respectively) and NCEP/NCAR reanalysis data, we have looked at changes in the climatological drivers of the glacial melt processes, and the sensitivity of the inland ice cap with regard to winter melting periods and pressure anomalies. The analysis has revealed, a positive trend of 5K over four decades in minimum surface air temperatures for winter months, clearly exceeding the published annual mean statistics, associated to a decrease in mean monthly winter sea level pressure. This concurs with a positive trend in the Southern Annular Mode (SAM) index, which gives a measure for the strength and extension of the Antarctic vortex. We connect this trend with a higher frequency of low-pressure systems hitting the South Shetland Islands during austral winter, bringing warm and moist air masses from lower latitudes. Due to its exposure to the impact of transient synoptic weather systems, the ice cap of KGI is especially vulnerable to changes during winter glacial mass accumulation period. A revision of seasonal changes in adiabatic air temperature lapse rates and their dependency on exposure and elevation has shown a clear decoupling of atmospheric surface layers between coastal areas and the higher-elevation ice cap, showing the higher sensitivity to free atmospheric flow and synoptic changes. Observed surface air temperature lapse rates show a high variability during winter months (standard deviations up to ±1.0K/100 m), and a distinct spatial variability reflecting the impact of synoptic weather patterns. The observed advective conditions bringing warm, moist air with high temperatures and rain, lead to melt conditions on the ice cap, fixating surface air temperatures to the melting point. This paper assesses the impact of large-scale atmospheric circulation variability and climatic changes on the atmospheric surface layer and glacier mass accumulation of the upper ice cap during winter season for the Warszawa Icefield on KGI.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

Radiocarbon dates and dinoflagellate cysts of sediment core DP30PC

To obtain insight into character and potential forcing of short-term climatic and oceanographic variability in the southern Italian region during the "Roman Classical Period" (60 BC-AD 200), climatic and environmental reconstructions based on a dinoflagelate cyst record from a well dated site in the Gulf of Taranto located at the distal end of the Po-river discharge plume have been established with high temporal resolution. Short-term fluctuations in accumulation rates of the Adriatic Surface Water species Lingulodinium machaerophorum, the freshwater algae Concentricystes and species resistant to aerobic degradation indicate that fluctuations in the trophic state of the upper waters are related to river discharge of northern and eastern Italian rivers which in turn are strongly related to precipitation in Italy. The dinoflagellate cyst association indicates that local sea surface temperatures which in this region are strongly linked to local air temperatures were slightly higher than today. We reconstruct that sea surface temperatures have been relatively high and stable between 60 BC-AD 90 and show a decreasing trend after AD 90. Fluctuations in temperature and river discharge rates have a strong cyclic character with main cyclicities of 7-8 and 11 years. We argue that these cycles are related to variations of the North Atlantic Oscillation climate mode. A strong correlation is observed with global variation in Delta14C anomalies suggesting that solar variability might be one of the major forcings of the regional climate. Apart from cyclic climate variability we observed a good correlation between non-cyclic temperature drops and global volcanic activity indicating that the latter forms an additional major forcing factor of the southern Italian climate during the Roman Classical Period.

Temperature reconstruction of sediment cores from the northeastern Arabian Sea

In order to reconstruct the monsoonal variability during the late Holocene we investigated a complete, annually laminated sediment record from the oxygen minimum zone (OMZ) off Pakistan for oxygen isotopes of planktic foraminifera and alkenone-derived sea surface temperatures (SST). Significant SST changes of up to 3°C which cannot be explained by changes in the alkenone-producing coccolithophorid species (inferred from the Gephyrocapsa oceanica / Emiliania huxleyi ratio) suggest that SST changes are driven by changes in the monsoon strength. Our high-(decadal)-resolution data indicate that the late Holocene in the northeastern Arabian Sea was not characterized by a stable uniform climate, as inferred from the Greenland ice cores, but by variations in the dominance of the SW monsoon conditions with significant effects on temperatures. Highest SST fluctuations of up to 3.0°C and 2.5°C were observed for the time interval from 4600 to 3300 years B.P. and during the past 500 years. The significant, short-term SST changes during the past 500 years might be related to climatic instabilities known from the northern latitudes ("Little Ice Age") and confirm global effects. Surface salinity values, reconstructed from delta18O records after correction for temperature-related oxygen isotope fractionation, suggest that in general, the past 5000 years were characterized by higher-than-recent evaporation and more intense SW monsoon conditions. However, between 4600 and 3700 years B.P., evaporation dropped, SW monsoon weakened, and NE monsoon conditions were comparatively enhanced. For the past 1500 years we infer strongly fluctuating monsoon conditions. Comparisons of reconstructed salinity records with ice accumulation data from published Tibetan ice core and Tibetan tree ring width data reveal that during the past 2000 years, enhanced evaporation in the northeastern Arabian Sea correlates with periods of increased ice accumulation in Tibet, and vice versa. This suggests a strong climatic relationship between both monsoon-controlled areas.

ERA-interim reanalysis debiased at FLUXNET sites

(preliminary) Exchanges of carbon, water and energy between the land surface and the atmosphere are monitored by eddy covariance technique at the ecosystem level. Currently, the FLUXNET database contains more than 500 sites registered and up to 250 of them sharing data (Free Fair Use dataset). Many modelling groups use the FLUXNET dataset for evaluating ecosystem model's performances but it requires uninterrupted time series for the meteorological variables used as input. Because original in-situ data often contain gaps, from very short (few hours) up to relatively long (some months), we develop a new and robust method for filling the gaps in meteorological data measured at site level. Our approach has the benefit of making use of continuous data available globally (ERA-interim) and high temporal resolution spanning from 1989 to today. These data are however not measured at site level and for this reason a method to downscale and correct the ERA-interim data is needed. We apply this method on the level 4 data (L4) from the LaThuile collection, freely available after registration under a Fair-Use policy. The performances of the developed method vary across sites and are also function of the meteorological variable. On average overall sites, the bias correction leads to cancel from 10% to 36% of the initial mismatch between in-situ and ERA-interim data, depending of the meteorological variable considered. In comparison to the internal variability of the in-situ data, the root mean square error (RMSE) between the in-situ data and the un-biased ERA-I data remains relatively large (on average overall sites, from 27% to 76% of the standard deviation of in-situ data, depending of the meteorological variable considered). The performance of the method remains low for the Wind Speed field, in particular regarding its capacity to conserve a standard deviation similar to the one measured at FLUXNET stations.

Radiocarbon ages and diatom abundance in cores GA306-BC4 and GA306-GC4 obtained during Galathea 3 cruise, West Greenland

Diatom assemblages from Holsteinsborg Dyb on the West Greenland shelf were analysed with high temporal resolution for the last 1200 years. A high degree of consistency between changes in frequency of selected diatom species and instrumental data from the same area during the last 70 years confirms the reliability of diatoms (particularly sea-ice species and warm-water species) for the study of palaeoceanographic changes in this area. A general cooling trend with some fluctuations is marked by an increase in sea-ice species throughout the last 1200 years. A relatively warm period with increased influence of Atlantic water masses of the Irminger Current (IC) is found at AD 750-1330, although with some oceanographic variability after AD 1000. A pronounced oceanographic shift occurred at AD 1330, corresponding in time to the transition from the so-called 'Medieval Warm Period' (MWP) to the 'Little Ice Age' (LIA). The LIA cold episode is characterized by three intervals with particularly cold sea-surface conditions at AD 1330-1350, AD 1400-1575 and AD 1660-1710 as a result of variable influence of Polar waters in the area. During the last 70 years, two relatively warm periods and one cold period (the early 1960s to mid-1990s) are indicated by changes in the diatom components. Our study demonstrates that sedimentary records on the West Greenland shelf provide valuable palaeoenvironment data that confirm a linkage between local and large-scale North Atlantic oceanographic and atmospheric oscillations.

Geochemistry of FeMn oxyhydroxides of IODP Site 303-U1302 and 303-1303

Regional weathering intensity must have changed dramatically at high latitudes during the Quaternary as a consequence of repeated continental glaciation. Investigation of these glacial/interglacial changes at high temporal resolution is possible with the recent development of Pb isotopes in FeMn oxyhydroxide phases as a proxy for region-specific weathering intensity, where increases in the radiogenic component are thought to correspond to increased continental weathering fluxes. Here we present a Pb isotope record sourced from the FeMn oxyhydroxide fraction in marine sediments from IODP Sites U1302/3 on Orphan Knoll (~3500 mbsl, NW Atlantic), spanning the last 37 ka. Located at the eastern edge of the Laurentide Ice Sheet (LIS), Site U1302/3 is well-placed to monitor changes in weathering intensity associated with LIS glacial history. Overall, the data show a close correspondence to local surface water d18O, with least radiogenic values during times of heavy d18O (glacial maximum) and most radiogenic values during times of light d18O (Holocene). This supports the prediction that weathering intensity in glaciated regions of the North Atlantic correlates with the exposure age of glacial debris. Superimposed on these background trends are extreme radiogenic excursions (e.g. variation in 206Pb/204Pb from ~19.2-21.0) contemporaneous with Heinrich events and the Younger Dryas. These data are substantially more radiogenic than existing records from the NW Atlantic, and most likely represent episodes of exceptionally high inputs of pre-formed FeMn oxyhydroxides during drainage of the LIS. Due to its extreme isotope composition, at least in the NW Atlantic region, Pb would appear to be a good proxy for the fluxes of weathered continental material and perhaps, by inference, nutrients to the surface ocean

Seagrass time-series analysis products, in Moreton Bay, Australia, with links to GeoTIFFs

The spatial and temporal dynamics of seagrasses have been well studied at the leaf to patch scales, however, the link to large spatial extent landscape and population dynamics is still unresolved in seagrass ecology. Traditional remote sensing approaches have lacked the temporal resolution and consistency to appropriately address this issue. This study uses two high temporal resolution time-series of thematic seagrass cover maps to examine the spatial and temporal dynamics of seagrass at both an inter- and intra-annual time scales, one of the first globally to do so at this scale. Previous work by the authors developed an object-based approach to map seagrass cover level distribution from a long term archive of Landsat TM and ETM+ images on the Eastern Banks (~200 km**2), Moreton Bay, Australia. In this work a range of trend and time-series analysis methods are demonstrated for a time-series of 23 annual maps from 1988 to 2010 and a time-series of 16 monthly maps during 2008-2010. Significant new insight was presented regarding the inter- and intra-annual dynamics of seagrass persistence over time, seagrass cover level variability, seagrass cover level trajectory, and change in area of seagrass and cover levels over time. Overall we found that there was no significant decline in total seagrass area on the Eastern Banks, but there was a significant decline in seagrass cover level condition. A case study of two smaller communities within the Eastern Banks that experienced a decline in both overall seagrass area and condition are examined in detail, highlighting possible differences in environmental and process drivers. We demonstrate how trend and time-series analysis enabled seagrass distribution to be appropriately assessed in context of its spatial and temporal history and provides the ability to not only quantify change, but also describe the type of change. We also demonstrate the potential use of time-series analysis products to investigate seagrass growth and decline as well as the processes that drive it. This study demonstrates clear benefits over traditional seagrass mapping and monitoring approaches, and provides a proof of concept for the use of trend and time-series analysis of remotely sensed seagrass products to benefit current endeavours in seagrass ecology.