Oxygen isotope and palaeotemperature records from six Greenland ice‐core stations: Camp Century,Dye‐3, GRIP,GISP2, Renland and NorthGRIP

Oxygen isotope variations spanning the last glacial cycle and the Holocene derived from ice‐core records for six sites in Greenland (Camp Century, Dye‐3, GRIP, GISP2, Renland and NorthGRIP) show strong similarities. This suggests that the dominant influence on oxygen isotope variations reflected in the ice‐sheet records was regional climatic change. Differences in detail between the records probably reflect the effects of basal deformation in the ice as well as geographical gradients in atmospheric isotope ratios. Palaeotemperature estimates have been obtained from the records using three approaches: (i) inferences based on the measured relationship between mean annual δ¹⁸O of snow and of mean annual surface temperature over Greenland; (ii) modelled inversion of the borehole temperature profile constrained either by the dated isotopic profile, or (iii) by using Monte Carlo simulation techniques. The third of these approaches was adopted to reconstruct Holocene temperature variations for the Dye 3 and GRIP temperature profiles, which yields remarkably compatible results. A new record of Holocene isotope variations obtained from the NorthGRIP ice‐core matches the GRIP short‐term isotope record, and also shows similar long‐term trends to the Dye‐3 and GRIP inverted temperature data. The NorthGRIP isotope record reflects: (i) a generally stronger isotopic signal than is found in the GRIP record; (ii) several short‐lived temperature fluctuations during the first 1500 yr of the Holocene; (iii) a marked cold event at ca. 8.2 ka (the ‘8.2 ka event’); (iv) optimum temperatures for the Holocene between ca. 8.6 and 4.3 ka, a signal that is 0.6‰ stronger than for the GRIP profile; (v) a clear signal for the Little Ice Age; and (vi) a clear signal of climate warming during the last century. These data suggest that the NorthGRIP stable isotope record responded in a sensitive manner to temperature fluctuations during the Holocene. Copyright © 2001 John Wiley & Sons, Ltd.     

3D DEM simulation of principal stress rotation in different planes of cross-anisotropic granular materials

Three-dimensional Discrete Element Method simulations have been performed to study the deformation of cross-anisotropic granular materials under principal stress rotation (PSR), for rotation planes oriented at different angles θ with respect to the bedding plane. The simulations have been conducted with a novel technique for applying specified stresses at three-dimensional boundaries. The results are qualitatively in agreement with experimental results from literature. Cumulative volume contraction is always observed under continuous PSR and increases with increasing θ. The dilatancy rate decreases with increasing number of PSR cycles, tending to zero. The noncoaxiality angle between the strain increment and the stress in the PSR plane increases with increasing number of cycles, reaching the same asymptotic value for samples of various densities and for various θ. Periodic oscillations of the dilatancy rate and noncoaxiality angle within each PSR cycle are observed with an increasing oscillation magnitude with increasing θ, due to the larger fabric anisotropy within the PSR plane. When θ = 30 or 60°, significant noncoaxial strain accumulation occurs in the plane perpendicular to the PSR plane due to the oblique angle between the PSR plane and the bedding plane, echoing the major principal fabric direction's being neither parallel nor perpendicular to the PSR plane. The macroscopic behavior of the samples is related to the microscopic parameters including coordination number and fabric anisotropy. With increasing number of cycles, the difference between normalized stress/strain/fabric increment tensors tends to become constant, with only a small lag between each pair, irrespective of θ.     

Reactions of Cl Atoms with Selected VOCs: Kinetics, Products and Mechanisms

The reactions of isoprene, MBO (2-methyl-3-buten-2-ol) and toluene with chlorine atoms have been studied at 298 ± 5 K and 740 ± 5 Torr with the use of FTIR spectroscopy. Major products of the isoprene-Cl reaction and of the MBO-Cl reaction have been identified and quantified, and reaction mechanisms have been tentatively proposed in order to explain the products formed. The reaction between isoprene and Cl atoms yields mainly HCl, formyl chloride, formic acid, methylglyoxal (pyruvic aldehyde), CO and CO2, while the MBO-Cl reaction forms acetone, HCl, formyl chloride, formic acid, CO, CO2. As products from the reaction between toluene and Cl we identified and quantified HCl and benzaldehyde. The rate constants for the reactions of isoprene and toluene with Cl atoms have also been determined using a relative rate method. The measured values are: kisoprene = (5.5 ± 1.0) × 10–10 cm3 molecule–1 s–1 and ktoluene = (5.6 ± 1.3) × 10–11 cm3 molecule–1 s–1. Atmospheric lifetimes have been estimated from these values.     

A Simple Model for the Vertical Transport of Reactive Species in the Convective Atmospheric Boundary Layer

We have developed a simple, steady-state, one-dimensional second-order closure model to obtain continuous profiles of turbulent fluxes and mean concentrations of non-conserved scalars in a convective boundary layer without shear. As a basic tool we first set up a model for conserved species with standard parameterizations. This leads to formulations for profiles of the turbulent diffusivity and the ratio of temperature-scalar covariance to the flux of the passive scalar. The model is then extended to solving, in terms of profiles of mean concentrations and fluxes, the NO _x –O₃ triad problem. The chemical reactions involve one first-order reaction, the destruction of NO₂ with decay time τ, and one second-order reaction, the destruction of NO and O₃ with the reaction constant k. Since the fluxes of the sum concentrations of NO _x = NO + NO₂ and O₃ + NO₂ turn out to be constant throughout the boundary layer, the problem reduces to solving two differential equations for the concentration and the flux of NO₂. The boundary conditions are the three surface fluxes and the fluxes at the top of the boundary layer, the last obtained from the entrainment velocity, and the concentration differences between the free troposphere and the top of the boundary layer. The equations are solved in a dimensionless form by using 1/(kτ) as the concentration unit, the depth h of the boundary layer as the length unit, the convective velocity scale w _* as the velocity unit, and the surface temperature flux divided by w _* as the temperature unit. Special care has been devoted to the inclusion of the scalar–scalar covariance between the concentrations of O₃ and NO. Sample calculations show that the fluxes of the reactive species deviate significantly from those of non-reactive species. Further, the diffusivities, defined by minus the flux divided by the concentration gradient may become negative for reactive species in contrast to those of non-reactive species, which in the present model are never negative.     

Morphodynamic effects of riparian vegetation growth after stream restoration

The prediction of the morphological evolution of renaturalized streams is important for the success of restoration projects. Riparian vegetation is a key component of the riverine landscape and is therefore essential for the natural rehabilitation of rivers. This complicates the design of morphological interventions, since riparian vegetation is influenced by and influences the river dynamics. Morphodynamic models, useful tools for project planning, should therefore include the interaction between vegetation, water flow and sediment processes. Most restoration projects are carried out in USA and Europe, where rivers are highly intervened and where the climate is temperate and vegetation shows a clear seasonal cycle. Taking into account seasonal variations might therefore be relevant for the prediction of the river morphological adaptation. This study investigates the morphodynamic effects of riparian vegetation on a re‐meandered lowland stream in the Netherlands, the Lunterse Beek. The work includes the analysis of field data covering 5 years and numerical modelling. The results allow assessment of the performance of a modelling tool in predicting the morphological evolution of the stream and the relevance of including the seasonal variations of vegetation in the computations. After the establishment of herbaceous plants on its banks, the Lunterse Beek did not show any further changes in channel alignment. This is here attributed to the stabilizing effects of plant roots together with the small size of the stream. It is expected that the morphological restoration of similarly small streams may result in important initial morphological adaptation followed by negligible changes after full vegetation establishment. Copyright © 2018 John Wiley & Sons, Ltd.     

Diffusive and total oxygen uptake of deep-sea sediments in the eastern South Atlantic Ocean:in situ and laboratory measurements

Total O₂ uptake rates were measured by the benthic flux chamber lander ELINOR, and O₂ microprofiles were measured by the profiling lander PROFILUR in the eastern South Atlantic. Diffusive O₂ fluxes through the diffusive boundary layer and the depth distribution of O₂ consumption rates within the sediment were calculated from the obtained microprofiles. The depth integrated O₂ consumption rate agreed closely with the diffusive O₂ uptake at all stations. Total O₂ uptake was 1.2–4.2 times the diffusive O₂ uptake, and the difference correlated with the abundance of macrofauna in the sediment. Diffusive O₂ uptake and O₂-penetration depths correlated with the organic content of the sediments and exhibited an inverse correlation with water depth. Total and diffusive rates of in situ O₂ uptake were higher than previously published data for shelf and abyssal sediments in the Atlantic, but were comparable to rates from upwelling areas in the eastern Pacific. Laboratory measurements on recovered sediment cores showed lower O₂ penetration depths and higher diffusive uptake rates than in situ measurements. The differences increased with increasing water depth. We primarily ascribe this compression of O₂ profiles to a transiently increased temperature during recovery and enhanced microbial activity in decompressed sediment cores. Total O₂ uptake rates measured in the laboratory on macrofauna-rich stations were, in contrast, lower than those measured in situ because of underrepresentation and disturbance of the macrofauna.     

An Afterthought: Ecosystem Metrics and Pressure Indicators

Editing the symposium contributions has been an interestingand instructive experience, because they provide a broad overviewof the well-nigh unlimited creativity of marine scientists totry and find building blocks for an ecosystem approach to fisherymanagement. A wealth of ideas appears to be emerging on howdifferent aspects of the state of exploited marine ecosystemscan be measured. Although the metrics produced are commonlyreferred to as "ecosystem indicators", however, it is generallyless clear what they actually indicate. The concept that management must take the wider ranging     

Characterization of the reactivity of riverine heterogeneous sediments using a facies-based approach; the Rhine–Meuse delta (The Netherlands)

Large groundwater resources are found in densely populated lowland areas, which consist often of young unconsolidated and reduced sediments. When anthropogenic activities lead to oxygenation of the aquifer, breakdown of the main reduced fractions, i.e. sedimentary organic matter (SOM) and pyrite, could lead to severe groundwater deterioration such as acidification, heavy metal mobilization, and increased hardness. The characterization of the reactive properties of these sediments is important in predicting groundwater deterioration, but is often complicated by the high degree of heterogeneity of these sediments. In this study, the potential reduction capacity (PRC, based on SOM and pyrite content), the potential buffer capacity (PBC, based on carbonate content), potential acidification capacity (PAC, based on the potential acid production by sulfide oxidation), and the measured reduction capacity (MRC) of five facies, which are typical of the riverine sediments in the Rhine–Meuse delta (The Netherlands) were determined. A universal facies-classification model was used to classify the deposits into more homogeneous sub-units based on lithologic and geogenic properties, with a further sub-division into oxic or anoxic redox environment based upon groundwater data and field observations. The bulk chemical data show strong variation across facies for the median values of PRC (186–9093 mmol O₂ kg⁻¹), PBC (17–132 mmol O₂ kg⁻¹), and PAC (36–1530 mmol H⁺ kg⁻¹). The MRC was measured as reactivity to molecular O₂ exposure and was 0.5–567.3 mmol O₂ kg⁻¹. Steady-state oxidation rates were in the wide range of 0.001–10.355 mmol O₂ kg⁻¹ day⁻¹ but were typically about 3–8 times faster in fine facies than in coarse facies. Both the PRC and MRC depend strongly on grain size, but also on the syn/post-depositional environment and redox conditions. The main part of the PRC consists of SOM, but pyrite reactivity is higher than SOM reactivity as shown by the relative depletion of pyrite in oxic subfacies and the preferential oxidation during the oxidation experiments. Some facies are very prone to acidification because the PAC is higher than the PBC, but the oxidation experiments also show that acidification could already start before the PRC is fully exhausted. This study, is one of the few that combines bulk chemical data, groundwater data, and reactivity measurements and shows that a facies-based approach is a practical tool in characterizing the reactivity of heterogeneous deposits.     

Palaeomagnetism of middle Proterozoic (c. 1.25 Ga) dykes from central North Greenland

From a nunatak in central North Greenland (81.5°N, 44.7°W) nine sites of Middle Proterozoic basic dykes, cutting Archaean basement, were palaeomagnetically investigated. After AF and thermal cleaning the nine dyke sites and three adjacently baked gneiss sites give a stable characteristic remanent mean direction of D = 265°, I = 21.5° ( N = 12, α ₉₅= 5.6°), the direction being confirmed by a detailed and positive baked contact test.
The polarity of the dykes in the nunatak area is opposite to that of the Zig-Zag Dal Basalts and the Midsommersø Dolerites in eastern North Greenland some 200–300 km away, the volcanics of which are assumed to be of similar age (about 1.25 Ga). The remanent directions of the two sets of data are antiparallel within the 95 per cent significance level of confidence.
When rotating Greenland 18° clockwise back to North America by the 'Bullard fit', the pole of the central North Greenland dolerites (NDL) falls at (14.3°N, 144.3°W). The reversed pole (14.3°S, 35.7°E) fits well on to the loop between 1.2 and 1.4 Ma on the apparent polar wander swath of Berger & York for cratonic North America.
The palaeomagnetic results from the Middle Proterozoic basic dykes from central North Greenland thus strengthen previous palaeomagnetic results from the Midsommersø Dolerites and Zig-Zag Dal Basalts from the Peary Land Region in eastern North Greenland, suggesting that Greenland was part of the North American craton at least for the period between c . 1.3 and 1 Ma (and probably up to the end of Cretaceous time). The major geographical meridian of Greenland was orientated approximately E–W, and the palaeo-latitude of Greenland was about 10°–15°.