LMDZ5B: the atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection

Based on a decade of research on cloud processes, a new version of the LMDZ atmospheric general circulation model has been developed that corresponds to a complete recasting of the parameterization of turbulence, convection and clouds. This LMDZ5B version includes a mass-flux representation of the thermal plumes or rolls of the convective boundary layer, coupled to a bi-Gaussian statistical cloud scheme, as well as a parameterization of the cold pools generated below cumulonimbus by re-evaporation of convective precipitation. The triggering and closure of deep convection are now controlled by lifting processes in the sub-cloud layer. An available lifting energy and lifting power are provided both by the thermal plumes and by the spread of cold pools. The individual parameterizations were carefully validated against the results of explicit high resolution simulations. Here we present the work done to go from those new concepts and developments to a full 3D atmospheric model, used in particular for climate change projections with the IPSL-CM5B coupled model. Based on a series of sensitivity experiments, we document the differences with the previous LMDZ5A version distinguishing the role of parameterization changes from that of model tuning. Improvements found previously in single-column simulations of case studies are confirmed in the 3D model: (1) the convective boundary layer and cumulus clouds are better represented and (2) the diurnal cycle of convective rainfall over continents is delayed by several hours, solving a longstanding problem in climate modeling. The variability of tropical rainfall is also larger in LMDZ5B at intraseasonal time-scales. Significant biases of the LMDZ5A model however remain, or are even sometimes amplified. The paper emphasizes the importance of parameterization improvements and model tuning in the frame of climate change studies as well as the new paradigm that represents the improvement of 3D climate models under the control of single-column case studies simulations.     

New U-Pb zircon and Ar/Ar muscovite age constraints on the emplacement of the Lizio syn-tectonic granite (Armorican Massif, France)

LA-ICP-MS U-Pb analyses performed on zircon grains from the Lizio granite yielded an emplacement age of 316 ± 6 Ma. Typical S-C structures show that the Lizio granite was emplaced contemporaneously with dextral shearing along the northern branch of the South Armorican Shear Zone and that it was therefore active at that time. ⁴⁰Ar/³⁹Ar analyses performed on muscovite grains yielded plateau dates ranging between 311.5 and 308.2 Ma. Muscovite chemistry is typical of primary magmatic muscovite, which precludes a late fluids-induced resetting of the K-Ar isotopic system. ⁴⁰Ar/³⁹Ar dates thus likely correspond to the cooling ages below the argon closure temperature. Considering the uncertainties on the measured ages, we can propose that either the Lizio granite cooled down quickly in less than a million of years or that it remained in a hot environment for several millions of years after its emplacement. This latter scenario could have been sustained by shear heating during dextral shearing along the northern branch of the South Armorican Shear Zone.     

H and Cl isotope systematics of apatite in brecciated lunar meteorites Northwest Africa 4472, Northwest Africa 773, Sayh al Uhaymir 169, and Kalahari 009

We have investigated the H and Cl systematics in apatite from four brecciated lunar meteorites. In Northwest Africa (NWA) 4472, most of the apatites contain ～2000–6000 ppm H₂O with δD between ?200 and 0‰, except for one grain isolated in the matrix, which contains ～6000 ppm H₂O with δD of ～500–900‰. This low‐δD apatite contains ～2500–7500 ppm Cl associated with δ³⁷Cl of ～15–20‰, while the high‐δD grain contains ～2500 ppm Cl with δ³⁷Cl of ～7–15‰. In NWA 773, apatites in a first group contain ～700–2500 ppm H₂O with δD values averaging around ～0 ± 100‰, while apatites in a second group contain ～5500–16500 ppm H₂O with δD ～250 ± 50‰. In Sayh al Uhaymir (SaU) 169 and Kalahari (Kal) 009, apatites are similar in terms of their H₂O contents (～600–3000 ppm) and δD values (?100 to 200‰). In SaU 169, apatites contain ～6000–10,000 ppm Cl, characterized by δ³⁷Cl of ～5–12‰. Overall, most of the analyzed apatite grains have δD within the range reported for carbonaceous chondrites, similar to apatite analyzed in ancient (>3.9 Ga) lunar magmatic. One grain in NWA 4472 has H and Cl isotope compositions similar to apatite from mare basalts. With an age of 4.35 Ga, this grain could be a representative of the oldest known lunar volcanic activity. Finally, since numerous evolved clasts in NWA 773 formed through silicate liquid immiscibility, the apatite grains with extremely high H₂O contents, reaching pure hydroxylapatite composition, could provide insights into the effects of such process on the evolution of volatiles in lunar magmas.     

Ellipsoidal Harmonic expansions of the gravitational potential: Theory and application

Small bodies of the solar system are now the targets of space exploration. Many of these bodies have elongated, non-spherical shapes, and the usual spherical harmonic expansions of their gravity fields are not well suited for the modelling of spacecraft orbits around these bodies. An elegant remedy is to use ellipsoidal harmonic expansions instead of the usual spherical ones. In this paper, we present their mathematical theory as well as a real application: the simulation of a landing on the surface of a kilometer-sized comet. We show that with an ellipsoidal harmonic expansion up to degree 5, the error on the landing position is at the meter level, while the corresponding error for the spherical harmonic expansion can reach tens of meters.This revised version was published online in October 2005 with corrections to the Cover Date.     

The April 2010 North African heatwave: when the water vapor greenhouse effect drives nighttime temperatures

Climate Dynamics - North Africa experienced a severe heatwave in April 2010 with daily maximum temperatures ( $$T_{max}$$ ) frequently exceeding $$40\,^{\circ }\mathrm{C}$$ and daily minimum...     

Linear stability analysis of the explicit treatment of mobilities in non-Newtonian and non-Darcy porous media flow simulation

A von Neumann stability analysis of the discretized conservation equation for single-phase porous media flows is performed, where non-Newtonian and non-Darcy effects are accounted for using a velocity (or mass flux)-dependent mobility factor. Comprehensive results in three dimensions for two low-order finite-volume discretizations typically encountered in reservoir simulation are provided, based on edge-centered and upstream cell-centered mobility calculations. It is found that common semi-implicit schemes, where the pressure gradient driving the flow is taken implicitly while the velocity-dependent mobility is evaluated explicitly, are subject to restrictions on the logarithmic derivative of mobility with respect to velocity. A remarkable new result is nevertheless obtained: for any physically acceptable strength of non-Newtonian and non-Darcy effects, there exists a stable and explicit method to evaluate the mobility, rendering the need to implement costly fully implicit schemes more difficult to justify.     

Geodynamic evolution of the Tethyan lithosphere as recorded in the Spontang Ophiolite,South Ladakh ophiolites(NW Himalaya,India)

The Spontang Ophiolite complex represents the most complete ophiolite sequence amongst the South Ladakh ophiolites and comprises mantle rocks(depleted harzburgites,dunites and minor lherzolites)as well as crustal rocks(basalt,isotropic gabbros,layered gabbros etc.).In the present study,detailed geochemistry(whole rock as well as mineral chemistry)and Sr-Nd isotopic analyses of thirty-six ultramaficmafic samples have been attempted to constraint the evolution and petrogenetic history of the Tethyan oceanic crust.Major,trace-element and REE patterns of the peridotites and their minerals indicate that the lherzolites experienced lower degrees of partial melting resembling abyssal peridotites(at higher temperatures,TREE=$1216℃)than the harzburgites(6%–8%versus 15%–17%).Elevated eNd(t)and variable⁸⁷ Sr/⁸⁶ Sr(t)ratios along with REE patterns suggest that the Spontang mafic rocks display N-MORB affinity with negligible participation of oceanic sediments in their genesis are originated from a depleted upper mantle with little contribution from subduction-related fluids.MORB-type Neotethyan oceanic crust is associated with the earliest phase of subduction(of older Jurassic age)through which a younger intra-oceanic island arc(Spong arc)subsequently developed.Harzburgites REE display typical U-shaped patterns,suggesting that these rocks have been metasomatized by LREE-enriched fluids.On the other side,mafic rocks are characterized by heterogeneous(Nb/La)PMand(Hf/Sm)PMand relatively homogeneous eNd(t),indicating interaction of subduction-related melts with the upper mantle during the initiation of subduction,in Early Cretaceous times.     

On the Use of Hydrological Models and Satellite Data to Study the Water Budget of River Basins Affected by Human Activities: Examples from the Garonne Basin of France

Natural and anthropogenic forcing factors and their changes significantly impact water resources in many river basins around the world. Information on such changes can be derived from fine scale in situ and satellite observations, used in combination with hydrological models. The latter need to account for hydrological changes caused by human activities to correctly estimate the actual water resource. In this study, we consider the catchment area of the Garonne river (in France) to investigate the capabilities of space-based observations and up-to-date hydrological modeling in estimating water resources of a river basin modified by human activities and a changing climate. Using the ISBA–MODCOU and SWAT hydrological models, we find that the water resources of the Garonne basin display a negative climate trend since 1960. The snow component of the two models is validated using the moderate-resolution imaging spectroradiometer snow cover extent climatology. Crop sowing dates based on remote sensing studies are also considered in the validation procedure. Use of this dataset improves the simulated evapotranspiration and river discharge amounts when compared to conventional data. Finally, we investigate the benefit of using the MAELIA multi-agent model that accounts for a realistic agricultural and management scenario. Among other results, we find that changes in crop systems have significant impacts on water uptake for agriculture. This work constitutes a basis for the construction of a future modeling framework of the sociological and hydrological system of the Garonne river region.