Fission track analyses of ODP Hole 129-800A samples

Four models of fission track annealing in apatite are compared with measured fission track lengths in samples from Site 800 in the East Mariana Basin, Ocean Drilling Program Leg 129, given an independently determined temperature history. The temperature history of Site 800 was calculated using a one-dimensional, compactive, conductive heat flow model assuming two end-member thermal cases: one for cooling of Jurassic ocean crust that has experienced no subsequent heating, and one for cooling of Cretaceous ocean crust. Because the samples analyzed were only shallowly buried and because the tectonic history of the area since sample deposition is simple, resolution of the temperature history is high. The maximum temperature experienced by the sampled bed is between 16°-21°C and occurs at 96 Ma; temperatures since the Cretaceous have dropped in spite of continued pelagic sediment deposition because heat flow has continued to decay exponentially and bottom-water temperatures have dropped. Fission tracks observed within apatite grains from the sampled bed are 14.6 +/- 0.1 µm (1 sigma) long. Given the proposed temperature history of the samples, one unpublished and three published models of fission track annealing predict mean track lengths from 14.8 to 15.9 µm. These models require temperatures as much as 40°C higher than the calculated paleotemperature maximum of the sampled bed to produce the same degree of track annealing. Measured and predicted values are different because annealing models are based on extrapolation of high temperature laboratory data to geologic times. The model that makes the closest prediction is based on the greatest number of experiments performed at low temperature and on an apatite having composition closest to that of the core samples.

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.

ODP Site 162-980 age model and color reflectance record of sediments from ODP Leg 162 sites

Sediment spectral reflectance measurements were generated aboard the JOIDES Resolution during Ocean Drilling Program Leg 162 shipboard operations. The large size of the raw data set (over 1.3 gigabytes) and limited computer hard disk storage space precluded detailed analysis of the data at sea, although broad band averages were used as aids in developing splices and determining lithologic boundaries. This data report describes the methods used to collect these data and their shipboard and postcruise processing. These initial results provide the basis for further postcruise research.

Fission track age determinations of rocks from north Greenland

Apatite fission track (FT) ages and length characteristics of samples obtained from Cambrian to Paleocene-aged sandstones collected along the margin of Nares Strait in Ellesmere Island in the Canadian Arctic Archipelago are dominated by a thermal history related to Paleogene relative plate movements between Greenland and Ellesmere Island. A preliminary inverse FT thermal model for a Cambrian (Archer Fiord Formation) sandstone in the hanging wall of the Rawlings Bay thrust at Cape Lawrence is consistent with Paleocene exhumational cooling, likely as a result of erosion of the thrust. This suggests that thrusting at Cape Lawrence occurred prior to the onset of Eocene compression, likely due to transpression during earlier strikeslip along the strait. Models for samples from volcaniclastic sandstones of the Late Paleocene Pavy Formation (from Cape Back and near Pavy River), and a sandstone from the Late Paleocene Mount Lawson Formation (at Split Lake, near Makinson Inlet) are also consistent with minor burial heating following known periods of basaltic volcanism in Baffin Bay and Davis Strait (c. 61-59 Ma), or related tholeiitic volcanism and intrusive activity (c. 55-54 Ma). Thermal models for samples from sea level dykes from around Smith Sound suggest a period of Late Cretaceous - Paleocene heating prior to final cooling during Paleocene time. These model results imply that Paleocene tectonic movements along Nares Strait were significant, and provide limited support for the former existence of the Wegener Fault. Apatite FT data from central Ellesmere Island suggest however, that cooling there occurred during Early Eocene time (c. 50 Ma), which was likely a result of erosion of thrusts during Eurekan compression. This diachronous cooling suggests that Eurekan deformation was partitioned at discrete intervals across Ellesmere Island, and thus it is likely that displacements along the strait were much less than the 150 km that has been previously suggested for the Wegener Fault.

Age model for ODP Leg 175 sites

Upwelling along the western coast of Africa south of the equator may be partitioned into three major areas, each having its own dynamics and history: (1) the eastern equatorial region, comprising the Congo Fan and the area of Mid-Angola; (2) the Namibia upwelling system, extending from the Walvis Ridge to Lüderitz; and (3) the Cape Province region, where upwelling is subdued. The highest nutrient contents in thermocline waters are in the northern region, the lowest in the southern one. Wind effects are at a maximum near the southern end of the Namibia upwelling system, and maximum productivity occurs near Walvis Bay, where the product between upwelling rate and nutrient content of upwelled waters is at a maximum. In the Congo/Angola region, opal tends to follow organic carbon quite closely in the Quaternary record. However, organic carbon has a strong precessional component, while opal does not. Despite relatively low opal content, sediments off Angola show the same patterns as those off the Congo; thus, they are part of the same regime. The spectrum shows nonlinear interference patterns between high- and low-latitude forcing, presumably tied to thermocline fertility and wind. On Walvis Ridge, as in the Congo-Angola region, the organic matter record behaves normally; that is, supply is high during glacial periods. In contrast, interglacial periods are favorable for opal deposition. The pattern suggests reduction in silicate content of the thermocline during glacial periods. The reversed phase (opal abundant during interglacials) persists during the entire Pleistocene and can be demonstrated deep into the Pliocene, not just on Walvis Ridge but all the way to the Oranje River and off the Cape Province. From comparison with other regions, it appears that silicate is diminished in the global thermocline, on average, whenever winds become strong enough to substantially shorten the residence time of silicate in upper waters (Walvis Hypothesis, solving the Walvis Paradox of reversed phase in opal deposition). The central discovery during Leg 175 was the documentation of a late Pliocene opal maximum for the entire Namibia upwelling system (early Matuyama Diatom Maximum [MDM]). The maximum is centered on the period between the end of the Gauss Chron and the beginning of the Olduvai Chron. A rather sharp increase in both organic matter deposition and opal deposition occurs near 3 Ma in the middle of the Gauss Chron, in association with a series of major cooling steps. As concerns organic matter, high production persists at least to 1 Ma, when there are large changes in variability, heralding subsequent pulsed production periods. From 3 to 2 Ma, organic matter and opal deposition run more or less parallel, but after 2 Ma opal goes out of phase with organic matter. Apparently, this is the point when silicate becomes limiting to opal production. Thus, the MDM conundrum is solved by linking planetary cooling to increased mixing and upwelling (ramping up to the MDM) and a general removal of silicate from the upper ocean through excess precipitation over global supply (ramping down from the MDM). The hypothesis concerning the origin of the Namibia opal acme or MDM is fundamentally the same as the Walvis Hypothesis, stating that glacial conditions result in removal of silicate from the thermocline (and quite likely from the ocean as a whole, given enough time). The Namibia opal acme, and other opal maxima in the latest Neogene in other regions of the ocean, marks the interval when a cooling ocean selectively removes the abundant silicate inherited from a warm ocean. When the excess silicate is removed, the process ceases. According to the data gathered during Leg 175, major upwelling started in the late part of the late Miocene. Presumably, this process contributed to the drawing down of carbon dioxide from the atmosphere, helping to prepare the way for Northern Hemisphere glaciation.

Clay mineralogical assemblages, Sm-Nd concentrations, Sm and Nd isotopic ratios and oxygen isotope ratios, OIS number and the age model of ODP Site 105-646

We present 40 Sm-Nd isotope measurements of the clay-size (<2 µm) fractions of sediments from the Southern Greenland rise (ODP-646) that span the last 365 kyr. These data track changes in the relative supply of fine particles carried into the deep Labrador Sea by the Western Boundary Under Current (WBUC) back to the fourth glacial-interglacial cycles. Earlier studies revealed three general sources of particles to the core site: (i) Precambrian crustal material from Canada, Greenland, and/or Scandinavia (North American Shield - NAS), (ii) Palaeozoic or younger crustal material from East Greenland, NW Europe, and/or western Scandinavia (Young Crust - YC) and (iii) volcanic material from Iceland and the Mid-Atlantic Ridge (MAR). Clay-size fractions from glacial sediments have the lowest Nd isotopic ratios. Supplies of young crustal particles were similar during glacial oxygen isotope stages (OIS) 2, 6, and 10. In contrast the mean volcanic contributions decreased relative to old craton material from OIS 10 to OIS 6 and then from OIS 6 to OIS 2. The glacial OIS 8 interval displays a mean Sm/Nd ratio similar to those of interglacials OIS 1, 5, and 9. Compared with other interglacials, OIS 7 was marked by a higher YC contribution but a similar ~30% MAR supply. The overall NAS contribution dropped by a factor of 2 during each glacial/interglacial transition, with the MAR contribution broadly replacing it during interglacials. To decipher between higher supplies and/or dilution, particle fluxes from each end member were estimated. Glacial NAS fluxes were systematically higher than interglacial fluxes. During the time interval examined, fine particle supplies to the Labrador Sea were strongly controlled by proximal ice-margin erosion and thus echoed the glacial stage intensity. In contrast, the WBUC-carried MAR supplies from the eastern basins did not change significantly throughout the last 365 kyr, except for a marked increase in surface-sediments that suggests unique modern conditions. Distal WBUC-controlled inputs from the Northern and NE North Atlantic seem to have been less variable than proximal supplies linked with glacial erosion rate.

Meteorological measurements from 8 dropsondes released during POLAR 4 flight on 1993-03-04 along a track orthogonal to the pack ice edge north west of Svalbard

Different parameterizations of subgrid-scale fluxes are utilized in a nonhydrostatic and anelastic mesoscale model to study their influence on simulated Arctic cold air outbreaks. A local closure, a profile closure and two nonlocal closure schemes are applied, including an improved scheme, which is based on other nonlocal closures. It accounts for continuous subgrid-scale fluxes at the top of the surface layer and a continuous Prandtl number with respect to stratification. In the limit of neutral stratification the improved scheme gives eddy diffusivities similar to other parameterizations, whereas for strong unstable stratifications they become much larger and thus turbulent transports are more efficient. It is shown by comparison of model results with observations that the application of simple nonlocal closure schemes results in a more realistic simulation of a convective boundary layer than that of a local or a profile closure scheme. Improvements are due to the nonlocal formulation of the eddy diffusivities and to the inclusion of heat transport, which is independent of local gradients (countergradient transport).

Geochemistry of Fe-Mn crusts at seamounts in the Northwest Pacific

Hydrogenetic ferromanganese crusts were dredged from four seamounts in the western Pacific, OSM7, OSM2, Lomilik, and Lemkein, aligned in a NW-SE direction parallel to Pacific Plate movement. The crusts consist of four well-defined layers with distinct textural and geochemical properties. The topmost layer 1 is relatively enriched in Mn, Co, Ni, and Mo compared to the underlying layer 2, which is relatively enriched in Al, Ti, K, and Rb and Cu, Zn, and excess Ba. Textural and geochemical properties of layer 2 suggest growth conditions under high biogenic and detrital flux. Such conditions are met in the equatorial Pacific (i.e., between the Intertropical Convergence Zone (ITCZ) and equatorial high-productivity zone). Layer 2 likely formed when each seamount was beneath the equatorial Pacific along its back track path. On the other hand, layer 1 probably started to grow after seamounts moved northwest from the ITCZ. This interpretation is consistent with the thickness of layer 1 across the four crusts, which increases to the northwest. Ages of the layer 1-layer 2 boundary in each crust, a potential proxy for northern margin of the ITCZ, also increase to the northwest at 17, 11, 8, and 5 Ma for OSM7, OSM2, Lomilik, and Lemkein, respectively. Assuming Pacific Plate motion of 0.3°/Myr, the seamounts were located at 12°N, 11°N, 9°N, and 8°N at the time of boundary formation. This result suggests that the north edge of the ITCZ has shifted south since the middle Miocene in the western Pacific, which agrees with information from the eastern Pacific.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.

Underway measurements of dissolved gases along POLARSTERN cruise track ANT-XXVII/2

We present measurements of pCO2, O2 concentration, biological oxygen saturation (Delta O2/Ar) and N2 saturation (Delta N2) in Southern Ocean surface waters during austral summer, 2010-2011. Phytoplankton biomass varied strongly across distinct hydrographic zones, with high chlorophyll a (Chla) concentrations in regions of frontal mixing and sea-ice melt. pCO2 and Delta O2 /Ar exhibited large spatial gradients (range 90 to 450 µatm and -10 to 60%, respectively) and co-varied strongly with Chla. However, the ratio of biological O2 accumulation to dissolved inorganic carbon (DIC) drawdown was significantly lower than expected from photosynthetic stoichiometry, reflecting the differential time-scales of O2 and CO2 air-sea equilibration. We measured significant oceanic CO2 uptake, with a mean air-sea flux (~ -20 mmol m-2 d-1) that significantly exceeded regional climatological values. N2 was mostly supersaturated in surface waters (mean Delta N2 of +2.5 %), while physical processes resulted in both supersaturation and undersaturation of mixed layer O2 (mean Delta O2phys = 2.1 %). Box model calculations were able to reproduce much of the spatial variability of Delta N2 and Delta O2phys along the cruise track, demonstrating significant effects of air-sea exchange processes (e.g. atmospheric pressure changes and bubble injection) and mixed layer entrainment on surface gas disequilibria. Net community production (NCP) derived from entrainment-corrected surface Delta O2 /Ar data, ranged from ~ -40 to > 300 mmol O2 m-2 d-1 and showed good coherence with independent NCP estimates based on seasonal mixed layer DIC deficits. Elevated NCP was observed in hydrographic frontal zones and regions of sea-ice melt with shallow mixed layer depths, reflecting the importance of mixing in controlling surface water light and nutrient availability.