Quantifying internal variability in a regional climate model: a case study for Southern Africa

This study aims at presenting various methodologies to separate the reproducible and irreproducible components of seasonal and intraseasonal climate variability simulated by a regional climate model over Southern Africa (south of 15°S), during an austral summer rainy season representative of the climatology. To that end, a 30-member ensemble simulation is performed using WRF laterally forced by the ERA40 reanalyses. Retained metrics include the analysis of weather regimes, signal-to-noise ratio, inter-member standard deviation and coefficient of variation. At the seasonal timescale, simulated rainfall amounts generally show a strong reproducibility, except in the subtropics and over the southern part of the Mozambique Channel. There, the number of rainy days is roughly similar in all members, while their average intensity varies extensively. At the intraseasonal timescale, the chronology of weather regimes, derived from the 500?hPa geopotential height, is highly reproducible. Rainfall variability is much less reproducible, especially in the central parts of the domain and near its outflow boundaries. Analysis of a South African regional index nonetheless indicates that both wet and dry spells tend to be accurately simulated and occur in phase in most members, demonstrating that they are embedded in large-scale variability patterns. Internal variability is lastly related to the lateral forcings along the domain boundaries. An objective classification of inflow/outflow mass fluxes allows identification of the recurrent synoptic configurations that favor strong or weak regional reproducibility. The main uncertainties concern the basic morphological features of rain-bearing systems (i.e., their spatial extension, location and propagation speed). Consequences for tropical-temperate interactions are then discussed.     

Deep inside a neoproterozoic intra-oceanic arc: growth,differentiation and exhumation of the Amalaoulaou complex (Gourma,Mali)

We show here that the Amalaoulaou complex, in the Pan-African belt of West Africa (Gourma, Mali), corresponds to the lower and middle sections of a Neoproterozoic intra-oceanic arc. This complex records a 90–130-Ma-long evolution of magmatic inputs and differentiation above a subducting oceanic slab. Early c. 793 Ma-old metagabbros crystallised at lower crustal or uppermost mantle depths (25–30 km) and have geochemical characteristic of high-alumina basalts extracted from a depleted mantle source slightly enriched by slab-derived sedimentary components ((La/Sm)_N < 1; ε_Nd: +5.4–6.2; ⁸⁷Sr/⁸⁶Sr: 0.7027–0.7029). In response to crustal thickening, these mafic rocks were recrystallised into garnet-granulites (850–1,000°C; 10–12 kbar) and subject to local dehydration–melting reactions, forming trondhjemititic leucosomes with garnet–clinopyroxene–rutile residues. Slightly after the granulitic event, the arc root was subject to strong HT shearing during partial exhumation (detachment faults/rifting or thrusting), coeval with the emplacement of spinel- and garnet-pyroxenite dykes crystallised from a high-Mg andesitic parental magma. Quartz and hornblende-gabbros (700–660 Ma) with composition typical of hydrous volcanic rocks from mature arcs ((La/Sm)_N: 0.9–1.8; ε_Nd: +4.6 to +5.2; ⁸⁷Sr/⁸⁶Sr: 0.7028–0.7031) were subsequently emplaced at mid-arc crust levels (~15 km). Trace element and isotopic data indicate that magmas tapped a depleted mantle source significantly more enriched in oceanic sedimentary components (0.2%). Exhumation occurred either in two stages (700–660 and 623 Ma) or in one stage (623 Ma) with a final exhumation of the arc root along cold P-T path (550°C, 6–9 kbar; epidote–amphibolite and greenschist facies conditions) during the main Pan-African collision event (620–580 Ma). The composition of magmas forming the Cryogenian Amalaoulaou arc and the processes leading to intra-arc differentiation are strikingly comparable to those observed in the deep section of exposed Mezosoic oceanic arcs, namely the Kohistan and Talkeetna complex. This evolution of the Amalaoulaou oceanic arc and its accretion towards the West African craton belong to the life and closure of the Pharusian Ocean that eventually led to the formation of the Greater Gondwana supercontinent, a similar story having occurred on the other side of the Sahara with the Mozambique Ocean.     

U redox fronts and kaolinisation in basement-hosted unconformity-related U ores of the Athabasca Basin (Canada): late U remobilisation by meteoric fluids

Proterozoic basement-hosted unconformity-related uranium deposits of the Athabasca Basin (Saskatchewan, Canada) were affected by significant uranium redistribution along oxidation–reduction redox fronts related to cold and late meteoric fluid infiltration. These redox fronts exhibit the same mineralogical and geochemical features as the well-studied uranium roll-front deposits in siliclastic rocks. The primary hydrothermal uranium mineralisation (1.6–1.3 Ga) of basement-hosted deposits is strongly reworked to new disseminated ores comprising three distinctly coloured zones: a white-green zone corresponding to the previous clay-rich alteration halo contemporaneous with hydrothermal ores, a uranium front corresponding to the uranium deposition zone of the redox front (brownish zone, rich in goethite) and a hematite-rich red zone marking the front progression. The three zones directly reflect the mineralogical zonation related to uranium oxides (pitchblende), sulphides, iron minerals (hematite and goethite) and alumino-phosphate-sulphate (APS) minerals. The zoning can be explained by processes of dissolution–precipitation along a redox interface and was produced by the infiltration of cold (<50°C) meteoric fluids to the hydrothermally altered areas. U, Fe, Ca, Pb, S, REE, V, Y, W, Mo and Se were the main mobile elements in this process, and their distribution within the three zones was, for most of them, directly dependent on their redox potential. The elements concentrated in the redox fronts were sourced by the alteration of previously crystallised hydrothermal minerals, such as uranium oxides and light rare earth element (LREE)-rich APS. The uranium oxides from the redox front are characterised by LREE-enriched patterns, which differ from those of unconformity-related ores and clearly demonstrate their distinct conditions of formation. Uranium redox front formation is thought to be linked to fluid circulation episodes initiated during the 400–300 Ma period during uplift and erosion of the Athabasca Basin when it was near the Equator and to have been still active during the last million years. A major kaolinisation event was caused by changes in the fluid circulation regime, reworking the primary uranium redox fronts and causing the redistribution of elements originally concentrated in the uranium-enriched meteoric-related redox fronts.     

Ab initio calculations of the Fe(II) and Fe(III) isotopic effects in citrates,nicotianamine, and phytosiderophore,and new Fe isotopic measurements in higher plants

Iron is one of the most abundant transition metal in higher plants and variations in its isotopic compositions can be used to trace its utilization. In order to better understand the effect of plant-induced isotopic fractionation on the global Fe cycling, we have estimated by quantum chemical calculations the magnitude of the isotopic fractionation between different Fe species relevant to the transport and storage of Fe in higher plants: Fe(II)-citrate, Fe(III)-citrate, Fe(II)-nicotianamine, and Fe(III)-phytosiderophore. The ab initio calculations show firstly, that Fe(II)-nicotianamine is ～3‰ (⁵⁶Fe/⁵⁴Fe) isotopically lighter than Fe(III)-phytosiderophore; secondly, even in the absence of redox changes of Fe, change in the speciation alone can create up to ～1.5‰ isotopic fractionation. For example, Fe(III)-phytosiderophore is up to 1.5‰ heavier than Fe(III)-citrate₂ and Fe(II)-nicotianamine is up to 1‰ heavier than Fe(II)-citrate. In addition, in order to better understand the Fe isotopic fractionation between different plant components, we have analyzed the iron isotopic composition of different organs (roots, seeds, germinated seeds, leaves and stems) from six species of higher plants: the dicot lentil (Lens culinaris), and the graminaceous monocots Virginia wild rye (Elymus virginicus), Johnsongrass (Sorghum halepense), Kentucky bluegrass (Poa pratensis), river oat (Uniola latifolia), and Indian goosegrass (Eleusine indica). The calculations may explain that the roots of strategy-II plants (Fe(III)-phytosiderophore) are isotopically heavier (by about 1‰ for the δ⁵⁶Fe) than the upper parts of the plants (Fe transported as Fe(III)-citrate in the xylem or Fe(II)-nicotianamine in the phloem). In addition, we suggest that the isotopic variations observed between younger and older leaves could be explained by mixing of Fe received from the xylem and the phloem.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Uniwavelength lidar sensitivity to spherical aerosol microphysical properties for the interpretation of Lagrangian stratospheric observations

The determination of stratospheric particle microphysical properties from multiwavelength lidar, including Rayleigh and/or Raman detection, has been widely investigated. However, most lidar systems are uniwavelength operating at 532 nm. Although the information content of such lidar data is too limited to allow the retrieval of the full size distribution, the coupling of two or more uniwavelength lidar measurements probing the same moving air parcel may provide some meaningful size information. Within the ORACLE-O3 IPY project, the coordination of several ground-based lidars and the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) space-borne lidar is planned during measurement campaigns called MATCH-PSC (Polar Stratospheric Clouds). While probing the same moving air masses, the evolution of the measured backscatter coefficient (BC) should reflect the variation of particles microphysical properties. A sensitivity study of 532 nm lidar particle backscatter to variations of particles size distribution parameters is carried out. For simplicity, the particles are assumed to be spherical (liquid) particles and the size distribution is represented with a unimodal log-normal distribution. Each of the four microphysical parameters (i.e. log-normal size distribution parameters, refractive index) are analysed separately, while the three others are remained set to constant reference values. Overall, the BC behaviour is not affected by the initial values taken as references. The total concentration (N₀) is the parameter to which BC is least sensitive, whereas it is most sensitive to the refractive index (m). A 2% variation of m induces a 15% variation of the lidar BC, while the uncertainty on the BC retrieval can also reach 15%. This result underlines the importance of having both an accurate lidar inversion method and a good knowledge of the temperature for size distribution retrieval techniques. The standard deviation (σ) is the second parameter to which BC is most sensitive to. Yet, the impact of m and σ on BC variations is limited by the realistic range of their variations. The mean radius (r_m) of the size distribution is thus the key parameter for BC, as it can vary several-fold. BC is most sensitive to the presence of large particles. The sensitivity of BC to r_m and σ variations increases when the initial size distributions are characterized by low r_m and large σ. This makes lidar more suitable to detect particles growing on background aerosols than on volcanic aerosols.     

Petroleum accumulation and leakage in a deeply buried carbonate reservoir,Níspero field (Mexico)

The filling history of the Níspero deeply buried Lower Cretaceous carbonate reservoir (below 4000 m) from the south part of Gulf of Mexico was studied using a combination of data from petrography, stable isotopes and fluid inclusions and compared with a one-dimensional burial model to derive timing.     

Zooplankton community structure in a highly turbid environment (Charente estuary,France): Spatio-temporal patterns and environmental control

Zooplankton assemblages were studied from January 2007 to January 2008 along the salinity gradient of the Charente estuary (France). A Lagrangian survey was performed monthly at five sampling stations defined by salinity (freshwater, 0.5, 5, 15 and 25) in order to collect zooplankton and measure the main environmental parameters (concentrations of suspended particulate matter, particulate organic carbon, chlorophyll a and phaeopigments). A combination of multivariate cluster analysis, species indicator index and canonical correspondence analysis was used to relate the spatio-temporal patterns of the zooplankton assemblages with environmental drivers. The estuary was divided into three different zones by means of environmental parameters while four zooplankton assemblages were identified along the salinity gradient. The Charente estuary appeared as one of the most turbid systems in Europe, with suspended particulate matter (SPM) concentration reaching 3.5 g l⁻¹ in the Maximum Turbidity Zone (MTZ). Algal heterotrophy and microphytobenthos resuspension from the wide mudflats could be responsible for the relatively high chlorophyll a concentrations measured within this MTZ. Salinity and SPM affected significantly the spatial distribution of zooplankton species while temperature and river flow seemed to control their temporal variations. From a zooplanktonic viewpoint, the highly turbid Charente estuary seemed to match an “ecotone–ecocline” model: the succession of species assemblages along the salinity gradient matched the concept of ecocline while the MTZ, which is a stressful narrow area, could be considered as an ecotone. Although such ecoclinal characteristics seemed to be a general feature of estuarine biocenoses, the ecotone could be more system-specific and biological compartment-specific.     

Tides in a body librating about a spin–orbit resonance: generalisation of the Darwin–Kaula theory

The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin–Kaula theory for a nonlibrating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque, so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.