Comparisons of observed ice properties and Envisat Advanced Synthetic Aperture Radar (ASAR) data during Nathanial B. Palmer cruise NBP09-01 and Oden cruise OSO08/09

Envisat Advanced Synthetic Aperture Radar (ASAR) Wide Swath Mode (WSM) images are used to derive C-band HH-polarization normalized radar cross sections (NRCS). These are compared with ice-core analysis and visual ship-based observations of snow and ice properties observed according to the Antarctic Sea Ice Processes and Climate (ASPeCt) protocol during two International Polar Year summer cruises (Oden 2008 and Palmer 2009) in West Antarctica. Thick first-year (TFY) and multi-year (MY) ice were the dominant ice types. The NRCS value ranges between -16.3 ± 1.1 and -7.6 ± 1.0 dB for TFY ice, and is -12.6 ± 1.3 dB for MY ice; for TFY ice, NRCS values increase from ~-15 dB to -9 dB from December/January to mid-February. In situ and ASPeCt observations are not, however, detailed enough to interpret the observed NRCS change over time. Co-located Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) vertically polarized 37 GHz brightness temperatures (TB37V), 7 day and 1 day averages as well as the TB37V difference between ascending and descending AMSR-E overpasses suggest the low NRCS values (-15 dB) are associated with snowmelt being still in progress, while the change towards higher NRCS values (-9dB) is caused by commencement of melt-refreeze cycles after about mid-January.

(Table 1) Mean temperatures in january and july, and total annual snowfall at Great Bear Lake and Great Slave Lake, Canada

The sensitivity of brightness temperature (T(B)) at 6.9, 10.7, and 18.7 GHz from Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) observations is investigated over five winter seasons (2002-2007) on Great Bear Lake and Great Slave Lake, Northwest Territories, Canada. The T(B) measurements are compared to ice thicknesses obtained with a previously validated thermodynamic lake ice model. Lake ice thickness is found to explain much of the increase of T(B) at 10.7 and 18.7 GHz. T(B) acquired at 18.7 GHz (V-pol) and 10.7 GHz (H-pol) shows the strongest relation with simulated lake ice thickness over the period of study (R**2 > 0.90). A comparison of the seasonal evolution of T(B) for a cold winter (2003-2004) and a warm winter (2005-2006) reveals that the relationship between T(B) and ice growth is stronger in the cold winter (2003-2004). Overall, this letter shows the high sensitivity of T(B) to ice growth and, thus, the potential of AMSR-E mid-frequency channels to estimate ice thickness on large northern lakes.

Seasonal mean air temperatures and ice phenology characteristics of Great Bear Lake and Great Slave Lake (2002-2009)

Time series of brightness temperatures (T(B)) from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) are examined to determine ice phenology variables on the two largest lakes of northern Canada: Great Bear Lake (GBL) and Great Slave Lake (GSL). T(B) measurements from the 18.7, 23.8, 36.5, and 89.0 GHz channels (H- and V- polarization) are compared to assess their potential for detecting freeze-onset/melt-onset and ice-on/ice-off dates on both lakes. The 18.7 GHz (H-pol) channel is found to be the most suitable for estimating these ice dates as well as the duration of the ice cover and ice-free seasons. A new algorithm is proposed using this channel and applied to map all ice phenology variables on GBL and GSL over seven ice seasons (2002-2009). Analysis of the spatio-temporal patterns of each variable at the pixel level reveals that: (1) both freeze-onset and ice-on dates occur on average about one week earlier on GBL than on GSL (Day of Year (DY) 318 and 333 for GBL; DY 328 and 343 for GSL); (2) the freeze-up process or freeze duration (freeze-onset to ice-on) takes a slightly longer amount of time on GBL than on GSL (about 1 week on average); (3) melt-onset and ice-off dates occur on average one week and approximately four weeks later, respectively, on GBL (DY 143 and 183 for GBL; DY 135 and 157 for GSL); (4) the break-up process or melt duration (melt-onset to ice-off) lasts on average about three weeks longer on GBL; and (5) ice cover duration estimated from each individual pixel is on average about three weeks longer on GBL compared to its more southern counterpart, GSL. A comparison of dates for several ice phenology variables derived from other satellite remote sensing products (e.g. NOAA Interactive Multisensor Snow and Ice Mapping System (IMS), QuikSCAT, and Canadian Ice Service Database) show that, despite its relatively coarse spatial resolution, AMSR-E 18.7 GHz provides a viable means for monitoring of ice phenology on large northern lakes.

(Table 1) Iceberg characteristics and start and end date and location of iceberg buoys in the Weddell Sea

The drift of 52 icebergs tagged with GPS buoys in the Weddell Sea since 1999 has been investigated with respect to prevalent drift tracks, sea ice/iceberg interaction, and freshwater fluxes. Buoys were deployed on small- to medium-sized icebergs (edge lengths ? 5 km) in the southwestern and eastern Weddell Sea. The basin-scale iceberg drift of this size class was established. In the western Weddell Sea, icebergs followed a northward course with little deviation and mean daily drift rates up to 9.5 ± 7.3 km/d. To the west of 40°W the drift of iceberg and sea ice was coherent. In the highly consolidated perennial sea ice cover of 95% the sea ice exerted a steering influence on the icebergs and was thus responsible for the coherence of the drift tracks. The northward drift of buoys to the east of 40°W was interrupted by large deviations due to the passage of low-pressure systems. Mean daily drift rates in this area were 11.5 ± 7.2 km/d. A lower threshold of 86% sea ice concentration for coherent sea ice/iceberg movement was determined by examining the sea ice concentration derived from Special Sensor Microwave Imager (SSM/I) and Advanced Microwave Scanning Radiometer for EOS (AMSR-E) satellite data. The length scale of coherent movement was estimated to be at least 200 km, about half the value found for the Arctic Ocean but twice as large as previously suggested. The freshwater fluxes estimated from three iceberg export scenarios deduced from the iceberg drift pattern were highly variable. Assuming a transit time in the Weddell Sea of 1 year, the iceberg meltwater input of 31 Gt which is about a third of the basal meltwater input from the Filchner Ronne Ice Shelf but spreads across the entire Weddell Sea. Iceberg meltwater export of 14.2 × 103 m3 s?1, if all icebergs are exported, is in the lower range of freshwater export by sea ice.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Magnetic susceptibility dataset from the early Givetian to early Frasnian in the La Thure section

Recent advances in radiometric dating result in significant improvements in the geological timescale and provide better insight into the timing of various processes and evolutions within the Earth's system. However, no radiometric ages are contained within the Givetian. Consequently, the absolute ages of the Givetian Stage boundaries, as well as the stage's duration, remain poorly constrained. As an alternative, the analysis of sedimentary cycles allows for the estimation of the duration of this stage. We examined the high-resolution magnetic susceptibility signals of four Givetian outcrops in the Givet area for a possible astronomical imprint, to fully understand the rates of evolutionary and environmental change. All four sections are firmly correlated and wavelet analyses of the magnetic susceptibility signals reveal the imprint of astronomical eccentricity forcing. The highly stable 405 kyr cycles constrain the duration of the Givetian Stage at 4.35±0.45 Myr, which is in good agreement with the International Chronostratigraphic Chart (5.0 Myr). The studied sections also exhibit an imprint of obliquity, suggesting a climatic teleconnection between low and high latitudes. The corresponding microfacies curves demonstrate similar astronomical imprint, and thereby indicate that the observed 10**5 year-scale cyclicity is the result of climatic and environmental change.

Geochemistry of ODP Hole 135-841B

Two igneous rock units were recovered at Site 841. More than 200 m of island-arc rhyolites, rhyolitic tuffs, lapilli tuffs, and pumice breccias, divided into five units, compose the basement at the site. These rhyolitic volcanics are late middle Eocene or older and formed part of a subaerial rhyolitic volcano. These low-K rhyolites were produced by fractional crystallization of a more mafic arc-tholeiitic lava or by dehydration melting of lower crustal arc tholeiites. The Site 841 basement rocks are similar in composition to high-SiO2 lavas in the Eocene basement on 'Eua and crystallized from depleted island-arc-tholeiitic basalts like those exposed on 'Eua. No evidence is present in the rhyolites, or in the clasts enclosed within them, for boninite series magmas at Site 841. The Site 841 rhyolitic complex bears no resemblance to Cretaceous rhyolites from the Lord Howe Rise, which are enriched in K and incompatible elements. The volcanic rocks at Site 841 are part of a widely distributed Eocene volcanic episode that marked the earliest phases of subduction in the Tonga region; they are not part of an older crustal fragment. The second igneous sequence is a series of basaltic dikes and sills that intruded Miocene sediments. These basalts have trace element abundances and ratios identical to upper Miocene lavas from the Lau Ridge. The Site 841 basalts do not have any geochemical characteristics that suggest they were generated by unusual thermal conditions in the shallow sub-forearc mantle. They are most reasonably interpreted as intrusions fed by basement dikes propagated from the associated active arc. No evidence for local serpentinite exposures, like those that are common in the Mariana forearc, was found at Site 841. The results from Site 841 provide strong support for hypotheses of forearc evolution that have been advanced for the Izu-Bonin-Mariana system.

Global LAI uncertainty and relative uncertainty datasets derived from a Triple Collocation Error Model (TCEM)

Uncertainty information for global leaf area index (LAI) products is important for global modeling studies but usually difficult to systematically obtain at a global scale. Here, we present a new method that cross-validates existing global LAI products and produces consistent uncertainty information. The method is based on a triple collocation error model (TCEM) that assumes errors among LAI products are not correlated. Global monthly absolute and relative uncertainties, in 0.05° spatial resolutions, were generated for MODIS, CYCLOPES, and GLOBCARBON LAI products, with reasonable agreement in terms of spatial patterns and biome types. CYCLOPES shows the lowest absolute and relative uncertainties, followed by GLOBCARBON and MODIS. Grasses, crops, shrubs, and savannas usually have lower uncertainties than forests in association with the relatively larger forest LAI. With their densely vegetated canopies, tropical regions exhibit the highest absolute uncertainties but the lowest relative uncertainties, the latter of which tend to increase with higher latitudes. The estimated uncertainties of CYCLOPES generally meet the quality requirements (± 0.5) proposed by the Global Climate Observing System (GCOS), whereas for MODIS and GLOBCARBON only non-forest biome types have met the requirement. Nevertheless, none of the products seems to be within a relative uncertainty requirements of 20%. Further independent validation and comparative studies are expected to provide a fair assessment of uncertainties derived from TCEM. Overall, the proposed TCEM is straightforward and could be automated for the systematic processing of real time remote sensing observations to provide theoretical uncertainty information for a wider range of land products.

Probability sampling protocol of classification maps from spaceborne/airborne image

To deliver sample estimates provided with the necessary probability foundation to permit generalization from the sample data subset to the whole target population being sampled, probability sampling strategies are required to satisfy three necessary not sufficient conditions: (i) All inclusion probabilities be greater than zero in the target population to be sampled. If some sampling units have an inclusion probability of zero, then a map accuracy assessment does not represent the entire target region depicted in the map to be assessed. (ii) The inclusion probabilities must be: (a) knowable for nonsampled units and (b) known for those units selected in the sample: since the inclusion probability determines the weight attached to each sampling unit in the accuracy estimation formulas, if the inclusion probabilities are unknown, so are the estimation weights. This original work presents a novel (to the best of these authors' knowledge, the first) probability sampling protocol for quality assessment and comparison of thematic maps generated from spaceborne/airborne Very High Resolution (VHR) images, where: (I) an original Categorical Variable Pair Similarity Index (CVPSI, proposed in two different formulations) is estimated as a fuzzy degree of match between a reference and a test semantic vocabulary, which may not coincide, and (II) both symbolic pixel-based thematic quality indicators (TQIs) and sub-symbolic object-based spatial quality indicators (SQIs) are estimated with a degree of uncertainty in measurement in compliance with the well-known Quality Assurance Framework for Earth Observation (QA4EO) guidelines. Like a decision-tree, any protocol (guidelines for best practice) comprises a set of rules, equivalent to structural knowledge, and an order of presentation of the rule set, known as procedural knowledge. The combination of these two levels of knowledge makes an original protocol worth more than the sum of its parts. The several degrees of novelty of the proposed probability sampling protocol are highlighted in this paper, at the levels of understanding of both structural and procedural knowledge, in comparison with related multi-disciplinary works selected from the existing literature. In the experimental session the proposed protocol is tested for accuracy validation of preliminary classification maps automatically generated by the Satellite Image Automatic MapperT (SIAMT) software product from two WorldView-2 images and one QuickBird-2 image provided by DigitalGlobe for testing purposes. In these experiments, collected TQIs and SQIs are statistically valid, statistically significant, consistent across maps and in agreement with theoretical expectations, visual (qualitative) evidence and quantitative quality indexes of operativeness (OQIs) claimed for SIAMT by related papers. As a subsidiary conclusion, the statistically consistent and statistically significant accuracy validation of the SIAMT pre-classification maps proposed in this contribution, together with OQIs claimed for SIAMT by related works, make the operational (automatic, accurate, near real-time, robust, scalable) SIAMT software product eligible for opening up new inter-disciplinary research and market opportunities in accordance with the visionary goal of the Global Earth Observation System of Systems (GEOSS) initiative and the QA4EO international guidelines.

Water temperature and salinity data from southern elephant seals from Marion Island

Over the last decade, several hundred seals have been equipped with conductivity-temperature-depth sensors in the Southern Ocean for both biological and physical oceanographic studies. A calibrated collection of seal-derived hydrographic data is now available, consisting of more than 165,000 profiles. The value of these hydrographic data within the existing Southern Ocean observing system is demonstrated herein by conducting two state estimation experiments, differing only in the use or not of seal data to constrain the system. Including seal-derived data substantially modifies the estimated surface mixedlayer properties and circulation patterns within and south of the Antarctic Circumpolar Current. Agreement with independent satellite observations of sea ice concentration is improved, especially along the East Antarctic shelf. Instrumented animals efficiently reduce a critical observational gap, and their contribution to monitoring polar climate variability will continue to grow as data accuracy and spatial coverage increase.