Ocean temperature measured by elephant seal-CTDs in the Southern Ocean with links to datasets

The potential effects of ocean warming on marine predators are largely unknown, though the impact on the distribution of prey in vertical space may have far reaching impacts on diving predators such as southern elephant seals. We used data from satellite-tracked southern elephant seals from Marion Island to investigate the relationship between their dive characteristics (dive depths, dive durations and time-at-depth index values) and environmental variables (temperature at depth, depth of maximum temperature below 100 m, frontal zone and bathymetry) as well as other demographic and behavioural variables (migration stage, age-class, track day and vertical diel strategy). While other variables, such as bathymetry and vertical diel strategy also influenced dive depth, our results consistently indicated a significant influence of temperature at depth on dive depths. This relationship was positive for all groups of animals, indicating that seals dived to deeper depths when foraging in warmer waters. Female seals adjusted their dive depths proportionally more than males in warmer water. Dive durations were also influenced by temperature at depth, though to a lesser extent. Results from time-at-depth indices showed that both male and female seals spent less time at targeted dive depths in warmer water, and were presumably less successful foragers when diving in warmer water. Continued warming of the Southern Ocean may result in the distribution of prey for southern elephant seals shifting either poleward and/or to increasing depths. Marion Island elephant seals are expected to adapt their ranging and diving behaviour accordingly, though such changes may result in greater physiological costs associated with foraging.

Dive depth and dive duration data of subadult and adult, male and female southern elephant seals from Marion Island between 2004 and 2008 with links to datasets

Although numerous studies have addressed the migration and dive behaviour of southern elephant seals (Mirounga leonina), questions remain about their habitat use in the marine environment. We report on the vertical use of the water column in the species and the potential lifetime implications for southern elephant seals from Marion Island. Long-term mark-resight data were used to complement vertical habitat use for 35 known individuals tagged with satellite-relay data loggers, resulting in cumulative depth use extrapolated for each individual over its estimated lifespan. Seals spent on average 77.59% of their lives diving at sea, 7.06% at the sea surface, and 15.35% hauled out on land. Some segregation was observed in maximum dive depths and depth use between male and female animals-males evidently being physiologically more capable of exploiting increased depths. Females and males spent 86.98 and 80.89% of their lives at sea, respectively. While at sea, all animals spent more time between 300 and 400 m depth, than any other depth category. Males and females spent comparable percentages of their lifetimes below 100 m depth (males: 65.54%; females: 68.92%), though males spent 8.98% of their lives at depths in excess of 700 m, compared to females' 1.84% at such depths. Adult males often performed benthic dives in excess of 2,000 m, including the deepest known recorded dive of any air-breathing vertebrate (>2,133 m). Our results provide a close approximation of vertical habitat use by southern elephant seals, extrapolated over their lifespans, and we discuss some physiological and developmental implications of their variable depth use.

At surface behaviour data of juvenile southern elephant seals from Marion Island between 2001 and 2006 with links to datasets

Marine mammals forage in dynamic environments characterized by variables that are continuously changing in relation to large-scale oceanographic processes. In the present study, behavioural states of satellite-tagged juvenile southern elephant seals (n = 16) from Marion Island were assessed for each reliable location, using variation in turning angle and speed in a state-space modelling framework. A mixed modelling approach was used to analyse the behavioural response of juvenile southern elephant seals to sea-surface temperature and proximity to frontal and bathymetric features. The findings emphasised the importance of frontal features as potentially rewarding areas for foraging juvenile southern elephant seals and provided further evidence of the importance of the area west of Marion Island for higher trophic-level predators. The importance of bathymetric features during the transit phase of juvenile southern elephant seal migrations indicates the use of these features as possible navigational cues.

Hydrographic data from Atka Bay taken by Weddell seals during expedition NEU2008 with links to datasets

In order to gain insights into species-level behavioural responses to the physical environment, it is necessary to obtain information from various populations and at all times of year. We analysed the influences of physical environmental parameters on the mid-summer dive behaviour of Weddell seals (Leptonychotes weddellii) from a little-known population at Atka Bay, Antarctica. Dive depth distributions followed a typical bimodal pattern also exhibited by seals from other populations and seals targeted both shallow water layers of <50 m and depths near the seafloor. Increased stratification of temperature layers within the water column resulted in increased forage efforts by the seals through relatively high numbers of dives to the seafloor, as well as forage effort associated with shallow dives. We interpret these behavioural responses to be due to increased water temperature stratification resulting in the concentration of prey species in particular depth layers.

Geodetic and satellite derived ice surface velocities of Vestfonna ice cap

During 2007 we launched a geodetic campaign on the Svalbard ice cap Vestfonna in order to estimate the velocity field of the ice cap. This was done within the frame of the IPY project KINNVIKA. We present here the velocity measurements derived from our campaigns 2007-2010 and compare the geodetic measurements against InSAR velocity fields from satellite platforms from 1995/96 and 2008. We find the spatial distribution of ice speeds from the InSAR is in good agreement within the uncertainty limits with our geodetic measurements. We observe no clear indication of seasonal ice speed differences, but we find a speed-up of the outlet glacier Franklinbreen between the InSAR campaigns, and speculate the outlet is having a surge phase.

Dive depth and dive duration data of male southern elephant seals from Marion Island between 2004 and 2008 with links to datasets

We describe the habitat use of 22 male southern elephant seals (Mirounga leonina) satellite tagged at Marion Island between 2004 and 2008. While a few areas of increased utilization appeared to be associated with areas of shallower bathymetry (such as sea-floor ridges and fracture zones), seals in our study did not target other areas of shallow bathymetry within close proximity to Marion Island. Rather, most elephant seals foraged pelagically over very deep water where much variation was evident in diel vertical migration strategies. These strategies resulted in generally deeper and longer dives than what has been reported for male elephant seals from other colonies. No significant differences were recorded for dive durations or dive depths between adults and sub-adults. However, younger animals displayed a positive relationship between dive durations and age, as well as between dive depths and age, while these relationships became negative for older animals. Mixed model outputs suggested that seals increased their aerobic fitness as migrations progressed, enabling them to undertake longer dives. We conclude that Marion Island male elephant seals exhibit much variability in dive strategy and are seemingly capable of exploiting a range of different prey types occurring in various depth layers.

At surface behaviour, dive depth and dive duration and associated hydrographic data of adult male southern elephant seals from King George Island from expedition JUB2010 with links to datasets

Understanding the distribution and foraging ecology of major consumers within pelagic systems, specifically in relation to physical parameters, can be important for the management of bentho-pelagic systems undergoing rapid change associated with global climate change and other anthropogenic disturbances such as fishing (i.e., the Antarctic Peninsula and Scotia Sea). We tracked 11 adult male southern elephant seals (Mirounga leonina), during their five-month post-moult foraging migrations from King George Island (Isla 25 de Mayo), northern Antarctic Peninsula, using tags capable of recording and transmitting behavioural data and in situ temperature and salinity data. Seals foraged mostly within the Weddell-Scotia Confluence, while a few foraged along the western Antarctic Peninsula shelf of the Bellingshausen Sea. Mixed model outputs suggest that the at-sea behaviour of seals was associated with a number of environmental parameters, especially seafloor depth, sea-ice concentrations and the temperature structure of the water column. Seals increased dive bottom times and travelled at slower speeds in shallower areas and areas with increased sea-ice concentrations. Changes in dive depth and durations, as well as relative amount of time spent during the bottom phases of dives, were observed in relation to differences in overall temperature gradient, likely as a response to vertical changes in prey distribution associated with temperature stratification in the water column. Our results illustrate the likely complex influences of bathymetry, hydrography and sea ice on the behaviour of male southern elephant seals in a changing environment and highlight the need for region-specific approaches to studying environmental influences on behaviour.

Basic and other measurements of radiation from the Baseline Surface Radiation Network (BSRN) Station Ny-Ålesund in the years 1992 to 2013, reference list of 253 datasets

At Ny-Ålesund (78.9° N), Svalbard, surface radiation measurements of up- and downward short- and longwave radiation are operated since August 1992 in the frame of the Baseline Surface Radiation Network (BSRN), complemented with surface and upper air meteorology since August 1993. The long-term observations are the base for a climatological presentation of the surface radiation data. Over the 21-year observation period, ongoing changes in the Arctic climate system are reflected. Particularly, the observations indicate a strong seasonality of surface warming and related changes in different radiation parameters. The annual mean temperature at Ny-Ålesund has risen by +1.3 ± 0.7 K per decade, with a maximum seasonal increase during the winter months of +3.1 ± 2.6 K per decade. At the same time, winter is also the season with the largest long-term changes in radiation, featuring an increase of +15.6 ± 11.6 W/m**2 per decade in the downward longwave radiation. Furthermore, changes in the reflected solar radiation during the months of snow melt indicate an earlier onset of the warm season by about 1 week compared to the beginning of the observations. The online available dataset of Ny-Ålesund surface radiation measurements provides a valuable data source for the validation of satellite instruments and climate models.

Global Ocean Particulate Organic Carbon flux merged with satellite parameters

The efficiency of the biological pump of carbon to the deep ocean depends largely on the biologically mediated export of carbon from the surface ocean and its remineralization with depth. Global satellite studies have primarily focused on chlorophyll concentration and net primary production (NPP) to understand the role of phytoplankton in these processes. Recent satellite retrievals of phytoplankton composition now allow for the size of phytoplankton cells to be considered. Here, we improve understanding of phytoplankton size structure impacts on particle export, remineralization and transfer. Particulate organic carbon (POC) flux observations from sediment traps and 234Th are compiled across the global ocean. Annual climatologies of NPP, percent microplankton, and POC flux at four time series locations and within biogeochemical provinces are constructed, and sinking velocities are calculated to align surface variables with POC flux at depth. Parameters that characterize POC flux vs. depth (export flux ratio, labile fraction, remineralization length scale) are then fit to the aligned dataset. Times of the year dominated by different size compositions are identified and fit separately in regions of the ocean where phytoplankton cell size showed enough dynamic range over the annual cycle. Considering all data together, our findings support the paradigm of high export flux but low transfer efficiency in more productive regions and vice versa for oligotrophic regions. However, when parsing by dominant size class, we find periods dominated by small cells to have both greater export flux and lower transfer efficiency than periods when large cells comprise a greater proportion of the phytoplankton community.

Aircraft and satellite derived arctic BrO profiles of the tropospheric column

We derive tropospheric column BrO during the ARCTAS and ARCPAC field campaigns in spring 2008 using retrievals of total column BrO from the satellite UV nadir sensors OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a photochemical simulation. We conduct a comprehensive comparison of satellite-derived tropospheric BrO column to aircraft in-situ observations of BrO and related species. The aircraft profiles reveal that tropospheric BrO, when present during April 2008, was distributed over a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). Perturbations to the total column resulting from tropospheric BrO are the same magnitude as perturbations due to longitudinal variations in the stratospheric component, so proper accounting of the stratospheric signal is essential for accurate determination of satellite-derived tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radiances were obtained over bright surfaces (albedo >0.7), for solar zenith angle <80° and clear sky conditions. The rapid activation of BrO due to surface processes (the bromine explosion) is apparent in both the OMI and GOME-2 based tropospheric columns. The wide orbital swath of OMI allows examination of the evolution of tropospheric BrO on about hourly time intervals near the pole. Low surface pressure, strong wind, and high PBL height are associated with an observed BrO activation event, supporting the notion of bromine activation by high winds over snow.

Satellite-based daily lake areas in the Yangtze Basin, China, 2000-2011

The Yangtze River Basin downstream of China's Three Gorges Dam (TGD) (thereafter referred to as "downstream" basin) hosts the largest cluster of freshwater lakes in East Asia. These lakes are crucial water stocks to local biophysical environments and socioeconomic development. Existing studies document that individual lakes in this region have recently experienced dramatic changes under the context of enduring meteorological drought, continuous population growth, and extensive water regulation since TGD's initial impoundment (i.e., June, 2003). However, spatial and temporal patterns of lake dynamics across the complete downstream Yangtze basin remain poorly characterized. Using daily MODIS imagery and an advanced thematic mapping scheme, this study presents a comprehensive monitoring of area dynamics in the downstream lake system at a 10-day temporal resolution during 2000-2011. The studied lakes constitute ~76% (~11,400 km**2) of the total downstream lake area, including the entire +70 major lakes larger than 20 km**2. The results reveal a decadal net decline in lake inundation area across the downstream Yangtze Basin, with a cumulative decrease of 849 km**2 or 7.4% from 2000 to 2011. Despite an excessive precipitation anomaly in the year 2010, the decreasing trend was tested significant in all seasons. The most substantial decrease in the post-TGD period appears in fall (1.1%/yr), which intriguingly coincides with the TGD water storage season. Regional lake dynamics exhibit contrasting spatial patterns, manifested as evident decrease and increase of aggregated lake areas respectively within and beyond the Yangtze Plain. This contrast suggests a marked vulnerability of lakes in the Yangtze Plain, to not only local meteorological variability but also intensified human water regulations from both the upstream Yangtze main stem (e.g., the TGD) and tributaries (e.g., lakes/reservoirs beyond the Yangtze Plain). The produced lake mapping result and derived lake area dynamics across the downstream Yangtze Basin provides a crucial monitoring basis for continuous investigations of changing mechanisms in the Yangtze lake system.

Imaging Land Subsidence Induced by Groundwater Extraction in Beijing (China) Using Satellite Radar Interferometry

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.

Parallel geostatistics for sparse and dense datasets

Very large spatially-referenced datasets, for example, those derived from satellite-based sensors which sample across the globe or large monitoring networks of individual sensors, are becoming increasingly common and more widely available for use in environmental decision making. In large or dense sensor networks, huge quantities of data can be collected over small time periods. In many applications the generation of maps, or predictions at specific locations, from the data in (near) real-time is crucial. Geostatistical operations such as interpolation are vital in this map-generation process and in emergency situations, the resulting predictions need to be available almost instantly, so that decision makers can make informed decisions and define risk and evacuation zones. It is also helpful when analysing data in less time critical applications, for example when interacting directly with the data for exploratory analysis, that the algorithms are responsive within a reasonable time frame. Performing geostatistical analysis on such large spatial datasets can present a number of problems, particularly in the case where maximum likelihood. Although the storage requirements only scale linearly with the number of observations in the dataset, the computational complexity in terms of memory and speed, scale quadratically and cubically respectively. Most modern commodity hardware has at least 2 processor cores if not more. Other mechanisms for allowing parallel computation such as Grid based systems are also becoming increasingly commonly available. However, currently there seems to be little interest in exploiting this extra processing power within the context of geostatistics. In this paper we review the existing parallel approaches for geostatistics. By recognising that diffeerent natural parallelisms exist and can be exploited depending on whether the dataset is sparsely or densely sampled with respect to the range of variation, we introduce two contrasting novel implementations of parallel algorithms based on approximating the data likelihood extending the methods of Vecchia [1988] and Tresp [2000]. Using parallel maximum likelihood variogram estimation and parallel prediction algorithms we show that computational time can be significantly reduced. We demonstrate this with both sparsely sampled data and densely sampled data on a variety of architectures ranging from the common dual core processor, found in many modern desktop computers, to large multi-node super computers. To highlight the strengths and weaknesses of the diffeerent methods we employ synthetic data sets and go on to show how the methods allow maximum likelihood based inference on the exhaustive Walker Lake data set.

Satellite borne radar altimetry : empirical orbit refinement and ocean signal recovery techniques

Measurements of the sea surface obtained by satellite borne radar altimetry are irregularly spaced and contaminated with various modelling and correction errors. The largest source of uncertainty for low Earth orbiting satellites such as ERS-1 and Geosat may be attributed to orbital modelling errors. The empirical correction of such errors is investigated by examination of single and dual satellite crossovers, with a view to identifying the extent of any signal aliasing: either by removal of long wavelength ocean signals or introduction of additional error signals. From these studies, it was concluded that sinusoidal approximation of the dominant one cycle per revolution orbit error over arc lengths of 11,500 km did not remove a significant mesoscale ocean signal. The use of TOPEX/Poseidon dual crossovers with ERS-1 was shown to substantially improve the radial accuracy of ERS-1, except for some absorption of small TOPEX/Poseidon errors. The extraction of marine geoid information is of great interest to the oceanographic community and was the subject of the second half of this thesis. Firstly through determination of regional mean sea surfaces using Geosat data, it was demonstrated that a dataset with 70cm orbit error contamination could produce a marine geoid map which compares to better than 12cm with an accurate regional high resolution gravimetric geoid. This study was then developed into Optimal Fourier Transform Interpolation, a technique capable of analysing complete altimeter datasets for the determination of consistent global high resolution geoid maps. This method exploits the regular nature of ascending and descending data subsets thus making possible the application of fast Fourier transform algorithms. Quantitative assessment of this method was limited by the lack of global ground truth gravity data, but qualitative results indicate good signal recovery from a single 35-day cycle.

Combination of altimetry data from different satellite missions

Substantial altimetry datasets collected by different satellites have only become available during the past five years, but the future will bring a variety of new altimetry missions, both parallel and consecutive in time. The characteristics of each produced dataset vary with the different orbital heights and inclinations of the spacecraft, as well as with the technical properties of the radar instrument. An integral analysis of datasets with different properties offers advantages both in terms of data quantity and data quality. This thesis is concerned with the development of the means for such integral analysis, in particular for dynamic solutions in which precise orbits for the satellites are computed simultaneously. The first half of the thesis discusses the theory and numerical implementation of dynamic multi-satellite altimetry analysis. The most important aspect of this analysis is the application of dual satellite altimetry crossover points as a bi-directional tracking data type in simultaneous orbit solutions. The central problem is that the spatial and temporal distributions of the crossovers are in conflict with the time-organised nature of traditional solution methods. Their application to the adjustment of the orbits of both satellites involved in a dual crossover therefore requires several fundamental changes of the classical least-squares prediction/correction methods. The second part of the thesis applies the developed numerical techniques to the problems of precise orbit computation and gravity field adjustment, using the altimetry datasets of ERS-1 and TOPEX/Poseidon. Although the two datasets can be considered less compatible that those of planned future satellite missions, the obtained results adequately illustrate the merits of a simultaneous solution technique. In particular, the geographically correlated orbit error is partially observable from a dataset consisting of crossover differences between two sufficiently different altimetry datasets, while being unobservable from the analysis of altimetry data of both satellites individually. This error signal, which has a substantial gravity-induced component, can be employed advantageously in simultaneous solutions for the two satellites in which also the harmonic coefficients of the gravity field model are estimated.