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).     

Nonparametric Estimates of the Sea State Bias for the Jason-1 Radar Altimeter

The Jason-1 sea state bias (SSB) is analyzed in depth from the first year of GDR products. Compared to previous missions, this work benefits from two aspects of the empirical determination of the SSB from the altimetric data themselves. First, from a methodological point of view, a nonparametric technique (NP) has been developed and largely tested on TOPEX/Poseidon 1, GFO and Envisat data. The NP estimator has proven to be a useful tool in the SSB estimation, and it is now mature enough to be used for a refined analysis. On the other hand, the SSB can be extracted from three different data sets (crossovers, collinear data, and residuals) with different characteristics. It is then possible to cross calibrate various estimations of the SSB models and to determine the most accurate one. A systematic comparison is made between these different estimates for the Jason-1 altimeter. The collinear and crossover data sets yield very similar estimates despite their difference of spatial and temporal distributions. These SSB models assure consistency with the TOPEX mission when comparing Jason-1 and TOPEX residuals during the tandem phase. Thanks to the present work, the impact of the short wavelengths filtering on the SSB estimation is evidenced. More generally, our understanding of potential errors affecting the sea surface height and their impact onto the SSB estimation is also improved.     

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.     

Jets from Young Stellar Objects: Current Constraints and Challenges for the Future

I outline future directions of research that seem promising and necessary to identify the ejection process in YSO jets, including observational tests of ejection models, and challenging theoretical and observational issues.     

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.     

Intra-interstadial environmental changes in Last Glacial loess revealed by molluscan assemblages from the Upper Palaeolithic site of Amiens-Renancourt 1 (Somme,France)

The Amiens-Renancourt 1 site recently yielded one of the most important Upper Palaeolithic human occupations of northern France by the number of flint artefacts and especially by the presence of Venus figurines. All the material comes from a single archaeological layer located in a tundra gley bracketed by loess units. A multi-proxy study combining a detailed stratigraphy, luminescence and radiocarbon datings and high-resolution (5 cm per sample) grain size and molluscan analyses was therefore carried out to reconstruct and date the associated environmental changes and to determine the exact context of the human occupation. The chronological frame thus established supports the correlations of the archaeology-bearing tundra gley and of an underlying arctic brown soil with Greenland interstadials GI-4 and GI-3. Composition changes in the molluscan population enabled the identification of transitional and optimum phases and sub-phases within these two pedogenetic horizons. A conceptual correlation model linking molluscan phases with millennial-scale variations of Greenland ice-core and Sieben Hengste speleothem climate records is proposed. The Human occupation appears contemporaneous to the end of the stadial–interstadial transition of GI-3. Synchronous in Amiens-Renancourt 1 and Nussloch, subsequent micro-gleys may also result from a regional/global forcing. Such a level of detail is unprecedented in a loess sequence.     

Climate variability and trends in downscaled high-resolution simulations and projections over Metropolitan France

In order to fulfill the society demand for climate information at the spatial scale allowing impact studies, long-term high-resolution climate simulations are produced, over an area covering metropolitan France. One of the major goals of this article is to investigate whether such simulations appropriately simulate the spatial and temporal variability of the current climate, using two simulation chains. These start from the global IPSL-CM4 climate model, using two regional models (LMDz and MM5) at moderate resolution (15–20 km), followed with a statistical downscaling method in order to reach a target resolution of 8 km. The statistical downscaling technique includes a non-parametric method that corrects the distribution by using high-resolution analyses over France. First the uncorrected simulations are evaluated against a set of high-resolution analyses, with a focus on temperature and precipitation. Uncorrected downscaled temperatures suffer from a cold bias that is present in the global model as well. Precipitations biases have a season- and model-dependent behavior. Dynamical models overestimate rainfall but with different patterns and amplitude, but both have underestimations in the South-Eastern area (Cevennes mountains) in winter. A variance decomposition shows that uncorrected simulations fairly well capture observed variances from inter-annual to high-frequency intra-seasonal time scales. After correction, distributions match with analyses by construction, but it is shown that spatial coherence, persistence properties of warm, cold and dry episodes also match to a certain extent. Another aim of the article is to describe the changes for future climate obtained using these simulations under Scenario A1B. Results are presented on the changes between current and mid-term future (2021–2050) averages and variability over France. Interestingly, even though the same global climate model is used at the boundaries, regional climate change responses from the two models significantly differ.