Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

Intrusion and its implication for subsidence: A case from the Baiyun Sag, on the northern margin of the South China Sea

We develop a series of simple numerical models to explain the anomalous subsidence and deposition phenomena on the northern continental margin of the South China Sea, in particular in the Baiyun Sag. The results suggest that a short-period high rate deposition of around 17 Ma is related to a rapid subsidence event, which may be due to episodic emplacement of a dense intrusion. Necking and gravity models indicate that in the basement of the Baiyun Sag, there is a dense zone that is 100–200 kg/m³ more dense than the surrounding country rock. Considering its high magnetic intensity and regional igneous activity, the dense zone is thought to be related to a phase of basalt intrusion that may have taken place around 17 Ma. Thermal and subsidence models indicate that a hot denser intrusion can cause significant subsidence immediately after the intrusion. The subsidence rate then slows down with cooling, thus becoming consistent with the observed subsidence curves at around 17 Ma. The results also indicate that the lithospheric strength under the Baiyun Sag is negligible, and that the high-velocity layer in the lowermost crust may be not an original part of the pre-rift crust. Instead, it is thought to be underplated intrusion emplaced at around 30 Ma when the continental margin broke up.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Comparative studies of the kinetic parameters of various algaenans and kerogens via open-system pyrolyses

Kinetic parameters were determined for the first time, via open-system pyrolyses, on algaenans (highly resistant biomacromolecules that are selectively preserved during kerogen formation) isolated from extant microalgae. Parallel studies were also carried out on 10 kerogens exhibiting, with one exception, a low level of maturity. These kerogens included samples chiefly derived from the selective preservation of the above algaenans and samples mainly, or almost exclusively, derived from the “natural vulcanization” pathway. Important differences in activation energy (E_a) distributions were observed between the four algaenans investigated and correlated with their chemical structures. The kerogens predominantly derived from algaenan-selective preservation (Pula alginite, NE 70 and BJ 248 Torbanites, Rundle Oil Shale) also exhibited pronounced differences in E_a distributions. These distributions provided: (i) information on the diversity of the source materials; and (ii) reflected the occurrence of important differences in chemical structures and thermal behaviour between three of the tested kerogens, even though they are all classified as low maturity type I. The Kimmeridge Clay samples and the Lorca Oil Shale showed broad E_a distributions shifted to low energies when compared with the above algaenans and kerogens. Such shifts reflect an important (or even almost exclusive for some of these kerogens) contribution of materials originating from sulphur incorporation into various lipids during early diagenesis. Finally, the kinetic data derived for the nine low maturity fossil samples were extrapolated to a very low, geological heating rate of 3°C Ma⁻¹ and the generation rate curves and cumulative yield curves thus obtained were compared.     

Organic matter quality in reclaimed boreal forest soils following oil sands mining

Following surface mining of the Athabasca Oil Sands deposits in northeastern Alberta, Canada, land reclamation entails the reconstruction of soil-like profiles using salvaged soil materials such as peat and mining by-products. The overall objective of this research was to assess soil organic matter (SOM) quality in different reclamation practices as compared to undisturbed soils found in the region. Soil samples (0–10 cm) were taken from 45 plots selected to represent undisturbed reference ecosites and reclamation treatments. Soil OM pools were isolated using a combination of acid hydrolysis and physical separation techniques. Chemical composition of the low density fractions was characterized using ramped cross polarization (RAMP-CP) ¹³C nuclear magnetic resonance (NMR). Differences between disturbed and undisturbed sites reflected the influence of different botanical inputs (peat vs. forest litter) to SOM composition. Reconstructed soils were characterized by significantly lower alkyl over O-alkyl carbon (ALK/O-ALK) ratios (0.3) than undisturbed sites (0.5). For these reconstructed soils, a significantly higher proportion of soil carbon was present in the sand associated (Heavy sand) pool (49.2 vs. 37.3) and in the acid-unhydrolyzable residue (AUR) fraction (61.3 vs. 54.7). These SOM parameters were significantly related to time since reconstruction, with the AUR (p value = 0.012) and Heavy sand (p value = 0.05) fractions decreasing with time since reconstruction, while the ALK/O-ALK ratio increased (p value = 0.006). These findings suggest that the ALK/O-ALK, AUR, and Heavy sand parameters are suitable indicators for monitoring SOM quality in these reconstructed soils following oil sands mining.     

Ontogenic variation and effect of collection procedure on leaf biomechanical properties of Mediterranean seagrass Posidonia oceanica (L.) Delile

Leaf mechanical traits are important to understand how aquatic plants fracture and deform when subjected to abiotic (currents or waves) or biotic (herbivory attack) mechanical forces. The likely occurrence of variation during leaf ontogeny in these traits may thus have implications for hydrodynamic performance and vulnerability to herbivory damage, and may be associated with changes in morphologic and chemical traits. Seagrasses, marine flowering plants, consist of shoot bundles holding several leaves with different developmental stages, in which outer older leaves protect inner younger leaves. In this study we examined the long‐lived seagrass Posidonia oceanica to determine ontogenic variation in mechanical traits across leaf position within a shoot, representing different developmental stages. Moreover, we investigated whether or not the collection procedure (classical uprooted shoot versus non‐destructive shoot method: cutting the shoot without a portion of rhizome) and time span after collection influence mechanical measurements. Neither collection procedure nor time elapsed within 48 h of collection affected measurements of leaf biomechanical traits when seagrass shoots were kept moist in dark cool conditions. Ontogenic variation in mechanical traits in P. oceanica leaves over intermediate and adult developmental stages was observed: leaves weakened and lost stiffness with aging, while mid‐aged leaves (the longest and thickest ones) were able to withstand higher breaking forces. In addition, younger leaves had higher nitrogen content and lower fiber content than older leaves. The observed patterns may explain fine‐scale within‐shoot ecological processes of leaves at different developmental stages, such as leaf shedding and herbivory consumption in P. oceanica.     

Applications of C and N stable isotopes to ecological and environmental studies in seagrass ecosystems

Stable isotopes of carbon and nitrogen are increasingly used in marine ecosystems, for ecological and environmental studies. Here, we examine some applications of stable isotopes as ecological integrators or tracers in seagrass ecosystem studies. We focus on both the use of natural isotope abundance as food web integrators or environmental tracers and on the use of stable isotopes as experimental tools. As ecosystem integrators, stable isotopes have helped to elucidate the general structure of trophic webs in temperate, Mediterranean and tropical seagrass ecosystems. As environmental tracers, stable isotopes have proven their utility in sewage impact measuring and mapping. However, to make such environmental studies more comprehensible, future works on understanding of basic reasons for variations of N and C stable isotopes in seagrasses should be encouraged. At least, as experimental tracers, stable isotopes allow the study of many aspects of N and C cycles at the scale of a plant or at the scale of the seagrass ecosystem.     

Gas-driven lava lake fluctuations at Erta ’Ale volcano (Ethiopia) revealed by MODIS measurements

The long-lived lava lake of Erta ’Ale volcano (Ethiopia) is remotely monitored by moderate resolution imaging spectroradiometers (MODIS) installed on satellites. The Normalised Thermal Index (NTI) (Wright et al. Remote Sens Environ 82:135–155 2002) is shown to be proportional to the volume of the lava lake based on visual observations. The lava lake’s variable level can be plausibly related to a stable foam, i.e. a mixture composed of densely packed non-coalescing bubbles in suspension within a liquid. This foam is trapped at the top of the magma reservoir, and its thickness changes in response to the gas flux feeding the foam being successively turned on and off. The temporal evolution of the foam thickness, and the resulting variation of the volume of the lava lake, is calculated numerically by assuming that the gas flux feeding the foam, initially constant and homogeneous since December 9, 2002, is suddenly stopped on December 13, 2002 and not restarted before May 2003. The best fit between the theoretical foam thickness and the level of the lava lake deduced from the NTI provides an estimate of both the reservoir radius, 155–170 m, and the gas flux feeding the foam, 5.5×10^?3–7.2×10^?3 m ³ s ^?1 when existing. This is in agreement with previous estimates from acoustic measurements (Bouche et al. Earth Planet Sci Lett 295:37–48 2010). The very good agreement between the theoretical foam thickness and that deduced from MODIS data shows for the first time the existence of a regime based on the behaviour of a stable foam, whose spreading towards the conduit (“wide” conduit condition), can explain the long-lived activity. Our predictive model, which links the gas flux at the vent to the foam spreading, could potentially be used on any volcano with a long-lived activity. The underlying gas flux and the horizontal surface area of the magma reservoir can then be deduced by combining modelling to continuous measurements of gas flux. The lava lake, when high, often shows regular rise and fall of its level. We have recognised a minimum of 26 very well marked cycles between January 2001 and December 13, 2002, corresponding to a typical return time of 10.8 ± 2.3 days and a gas volume of 8.3×10⁵ ± 2.0×10⁵ m ³. This corresponds to a gas volume fraction in the reservoir equal to 0.023–0.063 %. The yearly gas flux, estimated between December 13, 2002 and September 27, 2004, varies between 2.3×10^?6 and 5.9×10^?6 m ³ s ^?1 at the depth of the reservoir. The long-time series provided by infra-red sensors mounted on satellites could be used on any persistent volcano to detect potential periodic variations in the level of lava lakes or lava columns, providing that the vent has a funnel shape, as often, and is sufficiently large.