Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Pacific Ocean

Total organic carbon (TOC) was analyzed on four transects along 140°W in 1992 using a high temperature combustion/discrete injection (HTC/DI) analyzer. For two of the transects, the analyses were conducted on-board ship. Mixed-layer concentrations of organic carbon varied from about 80 µM C at either end of the transect (12°N and 12°S) to about 60 µM C at the equator. Total organic carbon concentrations decreased rapidly below the mixed-layer to about 38-40 µM C at 1000 m across the transect. Little variation was observed below this depth; deep water concentrations below 2000 m were virtually monotonic at about 36 µM C. Repeat measurements made on subsequent cruises consistently found the same concentrations at 1000 m or deeper, but substantial variations were observed in the mixed-layer and the upper water column above 400 m depth. Linear mixing models of total organic carbon versus sigmaT exhibited zones of organic carbon formation and consumption. TOC was found to be inversely correlated with apparent oxygen utilization (AOU) in the region between the mixed-layer and the oxygen minimum. In the mixed-layer, TOC concentrations varied seasonally. Part of the variations in TOC at the equator was driven by changes in the upwelling rate in response to variations in physical forcing related to an El Niño and to the passage of tropical instability waves. TOC export fluxes, calculated from simple box models, averaged 8±4 mmol C/m**2/day at the equator and also varied seasonally. These export fluxes account for 50-75% of the total carbon deficit and are consistent with other estimates and model predictions.

X-ray fluorescence scanning of discrete samples on the Schwalbenberg II loess-paleosol sequence, raw data, magnetic susceptibility, organic carbon and U-ratio as grainsize index

The Schwalbenberg II loess-paleosol sequence (LPS) denotes a key site for Marine Isotope Stage (MIS 3) in Western Europe owing to eight succeeding cambisols, which primarily constitute the Ahrgau Subformation. Therefore, this LPS qualifies as a test candidate for the potential of temporal high-resolution geochemical data obtained X-ray fluorescence (XRF) scanning of discrete samplesproviding a fast and non-destructive tool for determining the element composition. The geochemical data is first contextualized to existing proxy data such as magnetic susceptibility (MS) and organic carbon (Corg) and then aggregated to element log ratios characteristic for weathering intensity [LOG (Ca/Sr), LOG (Rb/Sr), LOG (Ba/Sr), LOG (Rb/K)] and dust provenance [LOG (Ti/Zr), LOG (Ti/Al), LOG (Si/Al)]. Generally, an interpretation of rock magnetic particles is challenged in western Europe, where not only magnetic enhancement but also depletion plays a role. Our data indicates leaching and top-soil erosion induced MS depletion at the Schwalbenberg II LPS. Besides weathering, LOG (Ca/Sr) is susceptible for secondary calcification. Thus, also LOG (Rb/Sr) and LOG (Ba/Sr) are shown to be influenced by calcification dynamics. Consequently, LOG (Rb/K) seems to be the most suitable weathering index identifying the Sinzig Soils S1 and S2 as the most pronounced paleosols for this site. Sinzig Soil S3 is enclosed by gelic gleysols and in contrast to S1 and S2 only initially weathered pointing to colder climate conditions. Also the Remagen Soils are characterized by subtle to moderate positive excursions in the weathering indices. Comparing the Schwalbenberg II LPS with the nearby Eifel Lake Sediment Archive (ELSA) and other more distant German, Austrian and Czech LPS while discussing time and climate as limiting factors for pedogenesis, we suggest that the lithologically determined paleosols are in-situ soil formations. The provenance indices document a Zr-enrichment at the transition from the Ahrgau to the Hesbaye Subformation. This is explained by a conceptual model incorporating multiple sediment recycling and sorting effects in eolian and fluvial domains.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.