Shoreline change on nine reef islands over a seven-decade period

Atoll islands are subject to a variety of processes that influence their geomorphological development. Analysis of historical shoreline changes using remotely sensed images has become an efficient approach to both quantify past changes and estimate future island response. However, the detection of long-term changes in beach width is challenging mainly for two reasons: first, data availability is limited for many remote Pacific islands. Second, beach environments are highly dynamic and strongly influenced by seasonal or episodic shoreline oscillations. Consequently, remote-sensing studies on beach morphodynamics of atoll islands deal with dynamic features covered by a low sampling frequency. Here we present a study of beach dynamics for nine islands on Takú Atoll, Papua New Guinea, over a seven-decade period. A considerable chronological gap between aerial photographs and satellite images was addressed by applying a new method that reweighted positions of the beach limit by identifying "outlier" shoreline positions. On top of natural beach variability observed along the reweighted beach sections, we found that one third of the analyzed islands show a statistically significant decrease in reweighted beach width since 1943. The total loss of beach area for all islands corresponds to 44% of the initial beach area. Variable shoreline trajectories suggest that changes in beach width on Takú Atoll are dependent on local control (that is, human activity and longshore sediment transport). Our results show that remote imagery with a low sampling frequency may be sufficient to characterize prominent morphological changes in planform beach configuration of reef islands.

Late Cenozoic carbonate accumulation rates of ODP Leg 115 samples

The principal paleoceanographic objective of Ocean Drilling Program Leg 115 was to collect a suite of materials that would allow reconstruction of the dynamic features of the late Cenozoic carbonate system in the equatorial Indian Ocean. This goal was achieved with the recovery of sediments from a closely spaced depth transect (1541-4428 m) of five sites (Sites 707 through 711) from on and around the Mascarene Plateau that record the last 50 m.y. of pelagic deposition. More than 2200 measurements of carbonate content are combined here with a highly resolved bio- and magnetostratigraphy to produce the first detailed compilation of bulk, carbonate, and noncarbonate mass accumulation rates (MARs) from the Indian Ocean. These results allow us to recognize three major depositional intervals, each characterized by a distinct depth-dependent pattern of carbonate accumulation: (1) the Paleogene, a time of moderate accumulation rates (0.4-0.7 g/cm**2/1000 yr) and reduced between-site accumulation differences; (2) the early and middle Miocene, a period characterized by greatly reduced carbonate MARs (typically <0.2 g/cm**2/1000 yr) at all sites and a shallow carbonate compensation depth; and (3) the late Miocene to Holocene, a time span marked by the highest bulk and carbonate accumulation rates of the last 50 Ma (1.6-1.8 g/cm**2/1000 yr), and the first appearance of substantial contrasts in carbonate accumulation as a function of the water depth of the drill site. The fundamentally different character of the carbonate system during each of these intervals must represent a regional response to the complex evolution of late Cenozoic oceans and climate.

Wetlands in central Siberia from ASAR WS 2003, link to shapefiles

An area in central Siberia (partial coverage of Turukhansky und Yeniseysky districts) was investigated using satellite data. It covers freshwater ecosystems of non-forested peatlands in boreal forests. The satellite data represent the growing seasons of 2003/2004. Microwave data were acquired by the Advanced Synthetic Aperture Radar (ASAR) instrument onboard ENVISAT. The multi-temporal capabilities and resolution (150mx150m in WS mode) of the ASAR wide swath mode enabled the detection of dynamic features >2ha over this vast area. Scatterometer (QuikScat) data could be employed to distinguish hydro-periods. Wetland types have been identified on the basis of seasonal changes in backscatter. Results for peatlands have been compared with Russian forest inventory data which contain information on wetland distribution.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Probabilistic models for the prediction of a ship performance in dynamic ice

We introduce two probabilistic, data-driven models that predict a ship's speed and the situations where a ship is probable to get stuck in ice based on the joint effect of ice features such as the thickness and concentration of level ice, ice ridges, rafted ice, moreover ice compression is considered. To develop the models to datasets were utilized. First, the data from the Automatic Identification System about the performance of a selected ship was used. Second, a numerical ice model HELMI, developed in the Finnish Meteorological Institute, provided information about the ice field. The relations between the ice conditions and ship movements were established using Bayesian learning algorithms. The case study presented in this paper considers a single and unassisted trip of an ice-strengthened bulk carrier between two Finnish ports in the presence of challenging ice conditions, which varied in time and space. The obtained results show good prediction power of the models. This means, on average 80% for predicting the ship's speed within specified bins, and above 90% for predicting cases where a ship may get stuck in ice. We expect this new approach to facilitate the safe and effective route selection problem for ice-covered waters where the ship performance is reflected in the objective function.

Nimrod (dynamic penetrometer) measurements of a sandbar at Raglan, New Zealand in 2009

Dynamic penetrometer data obtained with the Nimrod penetrometer (MARUM). Data is presented as (i) penetration depth (including for different layers if present), (ii) measured deceleration and (iv) estimated quasi-static bearing capacity including range of uncertainty due to the processing method. Lat/Long coordinates are given.