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.

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

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.

Mechanical properties and comparisons of dynamic and static Young's moduli at DSDP Holes 69-504B, 70-504B and 83-504B

From laboratory tests under simulated downhole conditions we tentatively conclude that the higher the triaxial-compressive strength, the lower the drilling rate of basalts from DSDP Hole 504B. Because strength is roughly proportional to Young's modulus of elasticity, which is related in turn to seismic-wave velocities, one may be able to estimate drilling rates from routine shipboard measurements. However, further research is needed to verify that P-wave velocity is a generally useful predictor of relative drilling rate.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

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.

Global Ocean Biomes: Mean and time-varying maps (NetCDF 7.8 MB)

Large-scale studies of ocean biogeochemistry and carbon cycling have often partitioned the ocean into regions along lines of latitude and longitude despite the fact that spatially more complex boundaries would be closer to the true biogeography of the ocean. Herein, we define 17 open-ocean biomes classified from four observational data sets: sea surface temperature (SST), spring/summer chlorophyll a concentrations (Chl a), ice fraction, and maximum mixed layer depth (maxMLD) on a 1° × 1° grid. By considering interannual variability for each input, we create dynamic ocean biome boundaries that shift annually between 1998 and 2010. Additionally we create a core biome map, which includes only the grid cells that do not change biome assignment across the 13 years of the time-varying biomes. These biomes can be used in future studies to distinguish large-scale ocean regions based on biogeochemical function.

Travel time and speed of gray whales and amphipod abundance along Chukotka Peninsula in 2006

The purpose of this study was to evaluate summer and fall residency and habitat selection by gray whales, Eschrichtius robustus, together with the biomass of benthic amphipod prey on the coastal feeding grounds along the Chukotka Peninsula. Thirteen gray whales were instrumented with satellite transmitters in September 2006 near the Chukotka Peninsula, Russia. Nine transmitters provided positions from whales for up to 81 days. The whales travelled within 5 km of the Chukotka coast for most of the period they were tracked with only occasional movements offshore. The average daily travel speeds were 23 km/day (range 9-53 km/day). Four of the whales had daily average travel speeds <1 km/day suggesting strong fidelity to the study area. The area containing 95% of the locations for individual whales during biweekly periods was on average 13,027 km**2 (range 7,097-15,896 km**2). More than 65% of all locations were in water <30 m, and between 45 and 70% of biweekly kernel home ranges were located in depths between 31 and 50 m. Benthic density of amphipods within the Bering Strait at depths <50 m was on average ~54 g wet wt/m**2 in 2006. It is likely that the abundant benthic biomass is more than sufficient forage to support the current gray whale population. The use of satellite telemetry in this study quantifies space use and movement patterns of gray whales along the Chukotka coast and identifies key feeding areas.

Time-series of seagrass and land cover maps from 1972 to 2010, in South East Queensland, Australia, with links to GeoTIFFs

Long term global archives of high-moderate spatial resolution, multi-spectral satellite imagery are now readily accessible, but are not being fully utilised by management agencies due to the lack of appropriate methods to consistently produce accurate and timely management ready information. This work developed an object-based remote sensing approach to map land cover and seagrass distribution in an Australian coastal environment for a 38 year Landsat image time-series archive (1972-2010). Landsat Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) imagery were used without in situ field data input (but still using field knowledge) to produce land and seagrass cover maps every year data were available, resulting in over 60 map products over the 38 year archive. Land cover was mapped annually using vegetation, bare ground, urban and agricultural classes. Seagrass distribution was also mapped annually, and in some years monthly, via horizontal projected foliage cover classes, sand and deep water. Land cover products were validated using aerial photography and seagrass maps were validated with field survey data, producing several measures of accuracy. An average overall accuracy of 65% and 80% was reported for seagrass and land cover products respectively, which is consistent with other studies in the area. This study is the first to show moderate spatial resolution, long term annual changes in land cover and seagrass in an Australian environment, created without the use of in situ data; and only one of a few similar studies globally. The land cover products identify several long term trends; such as significant increases in South East Queensland's urban density and extent, vegetation clearing in rural and rural-residential areas, and inter-annual variation in dry vegetation types in western South East Queensland. The seagrass cover products show that there has been a minimal overall change in seagrass extent, but that seagrass cover level distribution is extremely dynamic; evidenced by large scale migrations of higher seagrass cover levels and several sudden and significant changes in cover level. These mapping products will allow management agencies to build a baseline assessment of their resources, understand past changes and help inform implementation and planning of management policy to address potential future changes.

In-situ geotechnical measurements with the dynamic penetrometer Nimrod on the continental shelf off the Coromandel Peninsula, New Zealand

We tested the ability of a small dynamic penetrometer, Nimrod, to infer geotechnical properties of sediment mixtures in the inner shelf. The penetrometer is light and easy to operate, and its operation by scuba divers ensures a greater degree of precision than ship-based penetrometer deployments. We have studied selected positions along a sorted bedform (~ 100 m wide) on the continental shelf off the Coromandel Peninsula close to Tairua, North Island of New Zealand, and additionally took sediment samples at the exact positions of penetrometer impact, also by scuba divers. The derived dynamic penetrometer signatures (i) measured deceleration of the probe and estimated quasi-static bearing capacity as a measure of sediment strength, (ii) reflected changes in grain-size distribution ranging from very fine to very coarse sands, and (iii) revealed the uppermost seafloor stratification (top layer 2-6 cm) potentially being an indicator for sediment dynamics. In this manner, the device proved to be suitable for spatially fine-scaled surveys using divers' support and might deliver complementary information about sediment dynamics, in this case sorted-bedform maintenance.

Changes in dynamic viscosity of sea water with depth in the Pacific Ocean at Stations DM9-613W and DM9-618W

To check on the assumption that metabolic products of planktonic organisms can affect the coefficient of dynamic viscosity of seawater, viscosity was measured in water samples taken from depths of 0 to 1843 m, west of the Hawaiian Islands. Obtained results showed that plankton has no effect on viscosity of water in regions of low productivity and that viscosity can be determined with high degree of accuracy from the appropriate tables.

(Table 1) Comparisons of dynamic and static Young's moduli of basalt at DSDP Holes 69-504B, 70-504B and 83-504B

Compressional- and shear-wave velocity logs (Vp and Vs, respectively) that were run to a sub-basement depth of 1013 m (1287.5 m sub-bottom) in Hole 504B suggest the presence of Layer 2A and document the presence of layers 2B and 2C on the Costa Rica Rift. Layer 2A extends from the mudline to 225 m sub-basement and is characterized by compressional-wave velocities of 4.0 km/s or less. Layer 2B extends from 225 to 900 m and may be divided into two intervals: an upper level from 225 to 600 m in which Vp decreases slowly from 5.0 to 4.8 km/s and a lower level from 600 to about 900 m in which Vp increases slowly to 6.0 km/s. In Layer 2C, which was logged for about 100 m to a depth of 1 km, Vp and Vs appear to be constant at 6.0 and 3.2 km/s, respectively. This velocity structure is consistent with, but more detailed than the structure determined by the oblique seismic experiment in the same hole. Since laboratory measurements of the compressional- and shear-wave velocity of samples from Hole 504B at Pconfining = Pdifferential average 6.0 and 3.2 km/s respectively, and show only slight increases with depth, we conclude that the velocity structure of Layer 2 is controlled almost entirely by variations in porosity and that the crack porosity of Layer 2C approaches zero. A comparison between the compressional-wave velocities determined by logging and the formation porosities calculated from the results of the large-scale resistivity experiment using Archie's Law suggest that the velocity- porosity relation derived by Hyndman et al. (1984) for laboratory samples serves as an upper bound for Vp, and the noninteractive relation derived by Toksöz et al. (1976) for cracks with an aspect ratio a = 1/32 serves as a lower bound.

Mean dynamic topography of the North Indian Ocean

A satellite-only Mean Dynamic Topography (MDT) of the North Indian Ocean is estimated from the DIRR5 geoid and CNES_CLS11 mean sea surface (Schaffer et al. 2012). DIRR5 geoid is estimated from the latest release (Release 5) of GOCE gravity data according to previous studies (e.g., Johannessen et al. 2003; Raj, 2014). Note that this MDT estimated is referenced to a time period of 7 years (1993-1999). A correction data obtained from AVISO is later used to convert the MDT to a time reference of 20 years (1993-2012). More details are given in Raj (2016).