The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Jet instability and the stabilizing effect of topography on jets in two-layer rotating systems

Abstract

Laboratory experiments concerning azimuthal jets in two-layer rotating systems in the absence and presence of bottom topography aligned along the jets have been conducted. The jets were forced by the selective withdrawal of fluid from the upper layer of a two-fluid system contained in a circular dishpan geometry. The principal parameters measured in the experiments were the jet Rossby number, Ro, and a stratification parameter F = r ₁/(λ₁λ₂)^1/2 where r ₁ is the radius of the circular disc used for the selective withdrawal (i.e., r ₁ is the approximate radius of curvature of the jet) and λ₁,λ₂ are the internal Rossby radii of deformation in the upper and lower fluids, respectively.

The no-topography experiments show that for a sufficiently small F, the particular value depending on Ro, the jet is stable for the duration of the experiment. For sufficiently large F, again as a function of Ro, the jet becomes unstable, exhibiting horizontal wave disturbances from modes three to seven. An Ro against F flow regime diagram is presented.

Experiments are then conducted in the presence of a bottom topography having constant cross-section and extending around a mid-radius of the dishpan. The axis of the topography is in the vicinity of the jet axis forced in the no-topography experiments and the crest of the topography is in the vicinity of the interface between the two fluids (i.e., the front associated with the jet). The experiments show that in all cases investigated the jet tends to be stabilized by the bottom topography. Experiments with the topography in place, but with the interface between the fluids being above the topography crest, are shown to be unstable but more irregular than their no-topography counterparts.

Various quantitative measurements of the jet are presented. It is shown, for example, that the jet Rossby number defined in terms of the fluid withdrawal rate from the tank. Q, can be well correlated with a dimensionless vorticity gradient, V_G , across the upper layer jet. This allows for an assessment of the stability characteristics of a jet based on a knowledge of V_G (which can be estimated given a jet profile) and F.     

Geochemistry, Geochronology and Isotopic Evolution of the Chewore-Rufunsa Terrane, Southern Irumide Belt: a Mesoproterozoic Continental Margin Arc

The southern Irumide Belt (SIB) is an ENE–WSW-trending,late Mesoproterozoic orogenic belt located between the Congo–Tanzania–Bangweulu(CTB) and Kalahari cratons in central southern Africa. It isseparated from the late Mesoproterozoic Irumide Belt (IB) tothe north by Permo-Triassic graben, raising the possibilitythat the younger rifts reactivated a suture between the twobelts that has been rendered cryptic as a result of youngerKaroo cover. Both belts are dominated by calc-alkaline gneisses,but in addition the SIB contains abundant metavolcanic and metasedimentaryrocks. In this study we present detailed geochemical, isotopicand geochronological data for volcanic and plutonic lithologiesfrom the southernmost part of the SIB, the Chewore–RufunsaTerrane. This terrane comprises a wide variety of supracrustalto mid-crustal rocks that have major- and trace-element compositionssimilar to magmas formed in present-day subduction zones. Chondrite-normalizedrare earth element (REE) profiles and whole-rock Sm–Ndisotope compositions indicate that the parental supra-subductionmelts interacted with, and were contaminated by sialic continentalcrust, implying a continental-margin-arc setting. Secondaryionization mass spectrometry dating of magmatic zircon has yieldedcrystallization ages between c. 1095 and 1040 Ma, similar toelsewhere in the SIB. U–Pb dating and in situ Lu–Hfisotopic analyses of abundant xenocrystic zircon extracted fromthe late Mesoproterozoic granitoids indicate that the contaminantcontinental basement was principally Palaeoproterozoic in ageand had a juvenile isotopic signature at the time of its formation.These data are in contrast to those for the IB, which is characterizedby younger, c. 1020 Ma, calc-alkaline gneisses that formed bythe direct recycling of Archaean crust without significant additionof any juvenile material. We suggest that the SIB developedby the subduction of oceanic crust under the margin of an unnamedcontinental mass until ocean closure at c. 1040 Ma. Subsequentcollision between the SIB and the CTB margin led to the cessationof magmatism in the SIB and the initiation of compression andcrustal melting in the IB. KEY WORDS: geochemistry; Mesoproterozoic; SHRIMP zircon U–Pb dating; Sm–Nd isotopes; Southern Irumide Belt     

Effects of Hydrogel Amendment to Different Soils on Plant Available Water and Survival of Trees under Drought Conditions

The effect of super absorbent polyacrylate (SAP) hydrogel amendment to different soil types on plant available water (PAW), evapotranspiration and survival of Eucalyptus grandis, Eucalyptus citriodora, Pinus caribaea, Araucaria cunninghamii, Melia volkensii, Grevillea robusta, Azadirachta indica, Maesopsis eminii and Terminalia superba was investigated. The seedlings were potted in 3 kg size polythene bags filled with sand, loam, silt loam, sandy loam and clay soils, amended at 0 (control), 0.2 and 0.4% w/w hydrogel. The tree seedlings were allowed to grow normally with routine uniform watering in a glass house set up for a period of eight weeks, after which they were subjected to drought conditions by not watering any further. The 0.4% hydrogel amendment significantly (p < 0.05) increased the PAW by a factor of about three in sand, two fold in silt loam and one fold in sandy loam, loam and clay soils compared to the control. Similarly, the addition of either 0.2 or 0.4% hydrogel to the five soil types resulted in prolonged tree survival compared to the controls. Araucaria cunninghammi survived longest at 153 days, while Maesopsis eminii survived least (95 days) in sand amended at 0.4% after subjection to desiccation. Evapotranspiration was reduced in eight of the nine tree species grown in sandy loam, loam, silt loam and clay soils amended at 0.4% hydrogel. It is probable that soil amendment with SAP decreased the hydraulic soil conductivity that might reduce plant transpiration and soil evaporation.     

A steepness index unit volume flood hydrograph approach for sub‐daily flow disaggregation

The estimation of sub‐daily flows from daily flood flows is important for many hydrological and hydraulic applications. Flows during flood events often vary significantly within sub‐daily time‐scales, and failure to capture the sub‐daily flood characteristic can result in an underestimation of the instantaneous flood peaks, with possible risk of design failure. It is more common to find a longer record of daily flow series (observed or modelled using daily rainfall series) than sub‐daily flow data. This paper describes a novel approach, known as the steepness index unit volume flood hydrograph approach, for disaggregating daily flood flows into sub‐daily flows that takes advantage of the strong relationship between the standardized instantaneous flood peak and the standardized daily flood hydrograph rising‐limb steepness index. The strength of this relationship, which is considerably stronger than the relationship between the standardized flood peak and the event flood volume, is shown using data from six rivers flowing into the Gippsland Lakes in southeast Australia. The results indicate that the steepness index unit volume flood hydrograph approach can be used to disaggregate modelled daily flood flows satisfactorily, but its reliability is dependent on a model's ability to simulate the standardized daily flood hydrograph rising‐limb steepness index and the event flood volume. Copyright © 2007 John Wiley & Sons, Ltd.     

The Wieliczka Salt Mine

Just a few kilometres south of Cracow, in southern Poland, the scenery changes abruptly from the fluvial plain of the Vistula river, mantled by tills and loess deposits of the Pleistocene glaciation, to a boldly undulating topography representing the outermost thrust sheet of the Carpathian mountains. These northern foothills of the Carpathians are known in Poland as the Beskid Mountains. The Carpathian thrust front can be traced from Upper Silesia through Cracow and Przemysl in southern Poland as far as the Ukraine and Romania. At several places along the thrust front, salt deposits have been found, which have been mined over many centuries. One of the best known mines is located at Wieliczka, 15 km south-east of Cracow. The mine is a popular tourist location, but it also illustrates some spectacular geological features.     

Rock physics analysis for time-lapse seismic at Schiehallion Field, North Sea

Rock physics analysis plays a vital role in time‐lapse seismic interpretation because it provides the link between changes in rock and fluid properties and the resulting seismic data response. In this case study of the Schiehallion Field, we discuss a number of issues that commonly arise in rock physics analyses for time‐lapse studies. We show that:

• 1 Logarithmic fits of dry bulk (K_dry) and shear (G_dry) moduli vs. effective pressure (P_eff) are superior to polynomial fits.
• 2 2D surface fits of K_dry and G_dry over porosity (φ) and effective pressure using all the core data simultaneously are more useful and accurate than separate 1D fits over φ and P_eff for each individual core.
• 3 One average set (facies) of K_dry(φ, P_eff) and G_dry(φ, P_eff) can be chosen to represent adequately the entire Schiehallion reservoir.
• 4 Saturated velocities and densities modelled by fluid substitution of K_dry(φ, P_eff), G_dry(φ, P_eff) and the dry bulk density ρ_dry(φ) compare favourably with well‐log velocities and densities.
• 5 P‐ and S‐wave impedance values resulting from fluid substitution of K_dry(φ, P_eff), G_dry(φ, P_eff) and ρ_dry(φ) show that the largest impedance changes occur for high porosities and low effective pressures.
• 6 Uncertainties in K_dry(φ, P_eff) and G_dry(φ, P_eff) derived for individual cores can be used to generate error surfaces for these moduli that represent bounds for quantifying uncertainties in seismic modelling or pressure–saturation inversion.

     

The Ebro Deep-Sea Fan system

The Ebro Fan System consists of en echelon channel-levee complexes, 50×20 km in area and 200-m thick. A few strong reflectors in a generally transparent seismic facies identify the sand-rich channel floors and levee crests. Numerous continuous acoustic reflectors characterize overbank turbidites and hemipelagites that blanket abandoned channel-levee complexes. The interlobe areas between channel complexes fill with homogeneous mud and sand from mass flow and overbank deposition; these exhibit a transparent seismic character. The steep continental rise and sediment “drainage” of Valencia Trough at the end of the channel-levee complexes prevent the development of distributary channels and midfan lobe deposits. Margin setting represents fan and/or source area