An approximate energy cycle for inter-member variability in ensemble simulations of a regional climate model

The presence of internal variability (IV) in ensembles of nested regional climate model (RCM) simulations is now widely acknowledged in the community working on dynamical downscaling. IV is defined as the inter-member spread between members in an ensemble of simulations performed by a given RCM driven by identical lateral boundary conditions (LBC), where different members are being initialised at different times. The physical mechanisms responsible for the time variations and structure of such IV have only recently begun to receive attention. Recent studies have shown empirical evidence of a close parallel between the energy conversions associated with the time fluctuations of IV in ensemble simulations of RCM and the energy conversions taking place in weather systems. Inspired by the classical work on global energetics of weather systems, we sought a formulation of an energy cycle for IV that would be applicable for limited-area domain. We develop here a novel formalism based on local energetics that can be applied to further our understanding IV. Prognostic equations for ensemble-mean kinetic energy and available enthalpy are decomposed into contributions due to ensemble-mean variables (EM) and those due to deviations from the ensemble mean (IV). Together these equations constitute an energy cycle for IV in ensemble simulations of RCM. Although the energy cycle for IV was developed in a context entirely different from that of energetics of weather systems, the exchange terms between the various reservoirs have a rather similar mathematical form, which facilitates some interpretations of their physical meaning.     

Near‐infrared spectra of clinopyroxenes: Effects of calcium content and crystal structure

Abstract– Pyroxenes are among the most common minerals in the solar system and are ideally suited for remote geochemical analysis because of the sensitivity of their distinctive spectra to mineral composition. Fe²⁺ is responsible for the dominant pyroxene absorptions in the visible and near‐infrared, but substitutions of other cations such as Ca²⁺ change the crystal structure and site geometries and thus the crystal field splitting energies of the Fe cations. To define spectral systematics resulting from major pyroxene cations (Ca²⁺, Mg²⁺, and Fe²⁺), we focus on a suite of pyroxenes synthesized with only Ca²⁺, Mg²⁺, and Fe²⁺ in the two octahedral sites, specifically examining the effect of Ca²⁺ on pyroxene absorption bands. The modified Gaussian model is used to deconvolve pyroxene spectra into component bands that can then be linked directly to crystal field absorptions. In orthopyroxenes and low‐Ca clinopyroxenes, Ca²⁺‐content has a strong and predictable effect on the positions of the absorption bands. At a threshold of Wo₃₀, the crystal field environment stagnates and the M2 bands cease to change significantly as more Ca²⁺ is added. At Wo₅₀, when most of the M2 sites are filled by Ca²⁺, band positions do not change drastically, although the presence and strengths of the 1 and 2 μm bands are affected by even trace amounts of Fe²⁺ in the M2 site. It is thus apparent that next‐nearest neighbors and the distortions they impose on the pyroxene lattice affect the electronic states around the Fe²⁺ cations and control absorption band properties.     

Diagnosis of trace metal contamination in sediments: the example of Ensenada and El Sauzal, two harbors in Baja California, Mexico

Total metal concentrations in sediments from within Ensenada and El Sauzal Harbors are generally higher than at the mouths. Grain-size analyses suggested that this enrichment could be due to the presence of fine-grained sediments in the inner part of the harbors rather than to anthropogenic perturbations. The (Me/Al)(sample) ratios for Pb, Co, Ni and Fe were significantly higher for Ensenada Harbor relative to El Sauzal Harbor, whereas the ratios for Cd, Mn, Zn and Cu were statistically equivalent for both harbors. Calculated enrichment factors [EF(Me)=(Me/Al)(sample)/(Me/Al)(shale)] indicated that the metals showing slight enrichment were those associated with anthropogenic contamination (Pb, Zn), or probably related to primary productivity in the water column (Cd, Co). The levels of most of the metals were not greatly enriched, a consideration that is of the utmost importance when contamination issues are at stake.     

A seasonal study of the atmospheric dynamics over the Iberian Peninsula based on circulation types

A seasonal analysis of the atmospheric circulation over the Iberian Peninsula (IP) based on circulation types (CTs) obtained from sea level pressure and 500-hPa geopotential height is presented. The study covers the period of 1958–2008, when a high variability and important changes in winter and spring precipitation and temperature have been reported. Frequency, persistence, and the most probable transitions of the circulation types are analyzed. Among the clustering methods available in the literature, two of the most reliable classification methods have been tested, K-means and simulated annealing and diversified randomization. A comparison of both methods over the IP is presented for winter (DJF). The quality of the circulation types obtained through both methods as well as the better stability achieved by K-means suggest this method as more appropriated for our target area. Twelve CTs were obtained for each season and were analyzed. The patterns obtained were regrouped in five general situations: anticyclonic, cyclonic, zonal, summertime, and hybrid-mixed. The analysis of frequencies of these situations offers a similar characterization of the atmospheric circulation that others previously obtained by subjective methods. The analysis of the trends in frequency and persistence for each CT shows few significant trends, mainly in winter and spring with a general decrease of the cyclonic patterns and an increase of the anticyclonic situations. This can be related to the negative precipitation trends reported by other authors. Regarding the persistence, an interesting result is that there is a high interannual variability of the persistence in autumn and spring, when patterns can persist longer than in other seasons. An analysis of the most probable transitions between the CTs has been performed, revealing the existence of cyclic sequences in all seasons. These sequences are related to the high frequency of certain patterns such as the anticyclonic situations in winter. Finally, a clear seasonal dependence of the transitions between cyclonic situations associated with extratropical disturbances was found. This dependence suggests that the transitions of low-pressure systems towards the south of the IP are more likely in spring and autumn than in winter.     

Pole-to-pole moisture conditions for the IGY

Summary A study of the mean atmospheric humidity conditions on a planetary scale during the IGY covering the calendar year 1958 is presented. The fields of mean precipitable water content and of the zonal and meridional transports of water vapor are analyzed for the entire globe. Zonally averaged values of the various quantities at several levels are presented in tabular form and compared whenever possible with previous results or indirectly are analyzed on the basis of information obtained from different sources. The structure of these fields is studied and the corresponding implications for the general circulations of the atmosphere are discussed. Finally the zonal water balance for all the globe is discussed and its implications analyzed.     

Modelling of Ocean Shallow Convection under Fast Ice in the Arctic

Ocean convection in the Antarctic has been studied many times and has been revealed to be responsible for ice-cover reduction. In the Arctic, proof of that phenomenon has not been documented. It is believed that this phenomenon happens on a smaller scale in the Arctic when local circulation of deep warmer water melts and slows ice production. An example of this is the North Water (NOW) polynya in northern Baffin Bay. A polynya is an area of open water in an otherwise ice-covered area. As ice forms under the fast ice near the boundary of the polynya, ocean salts (brine) are ejected from the newly formed ice. This water, which has an increased concentration of salt, sinks and is replaced by warmer water from below, and this slows ice formation. In our study a coupled one-dimensional thermodynamic snow–fast ice model incorporating ocean heat flux input via a shallow convection model was used. Ice thickness was calculated using a thermodynamic model that included a current-induced entrainment model and a convection model to account for brine rejection during ice growth. Atmospheric observations from Grise Fiord and Thule and ocean profiles around the NOW polynya near these sites were used as input to the model. This purely thermodynamic study enables us to obtain ice thickness values that can be compared with qualitative observations. This modelling study compares two sites related to the NOW polynya. The results indicate that the shallow convection model simulates the reduction of fast ice near Thule but not near Grise Fiord.     

A paleoenvironmental reconstruction of Laguna Babícora, Chihuahua, Mexico based on ostracode paleoecology and trace element shell chemistry

Paleoecology of Laguna Babícora, Chihuahua, Mexico was reconstructed using ostracode faunal assemblages and shell chemistry. The paleolimnological record is used to show the magnitude of paleoclimatic changes in the area from 25,000 years to the present.Faunal assemblages consist of four species of the genus Limnocythere: L. sappaensis, L. ceriotuberosa, L. bradburyi and L. platyforma, all associated with Candona caudata, Candona patzcuaro and Cypridopsis vidua. A paleosalinity index developed from these assemblages indicates that the lake's salinity fluctuated frequently from oligo- to meso-haline conditions during the last 25,000 years. This pattern and low salinity range are in good agreement with modern TDS (here used as an indicator of salinity) values recorded from 26 wells and one spring from the area (258–975 mg l^–1). To estimate paleotemperature we examined the trace element content (Mg/Ca ratios) from individual valves of L. ceriotuberosa and L. platyforma, the two species most commonly recorded in Laguna Babícora.Shell Mg/Ca ratios of 204 specimens of these two species were used to estimate water temperature (Mg/Ca) by means of experimental standard coefficients. Our data show that paleowater temperature ranged from 5.6–21.3 °C (with 2 values ranging from 0.2–4.8 °C), which suggest a close correlation with atmospheric temperatures around the lake. These results are in good agreement with a modern mean winter temperature (3.5 °C) and mean summer temperature (20 °C) recorded in the area between 1970 and 1980.     

Angular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous Media

Wave-induced fluid flow generates a dominant attenuation mechanism in porous media. It consists of energy loss due to P-wave conversion to Biot (diffusive) modes at mesoscopic-scale inhomogeneities. Fractured poroelastic media show significant attenuation and velocity dispersion due to this mechanism. The theory has first been developed for the symmetry axis of the equivalent transversely isotropic (TI) medium corresponding to a poroelastic medium containing planar fractures. In this work, we consider the theory for all propagation angles by obtaining the five complex and frequency-dependent stiffnesses of the equivalent TI medium as a function of frequency. We assume that the flow direction is perpendicular to the layering plane and is independent of the loading direction. As a consequence, the behaviour of the medium can be described by a single relaxation function. We first consider the limiting case of an open (highly permeable) fracture of negligible thickness. We then compute the associated wave velocities and quality factors as a function of the propagation direction (phase and ray angles) and frequency. The location of the relaxation peak depends on the distance between fractures (the mesoscopic distance), viscosity, permeability and fractures compliances. The flow induced by wave propagation affects the quasi-shear (qS) wave with levels of attenuation similar to those of the quasi-compressional (qP) wave. On the other hand, a general fracture can be modeled as a sequence of poroelastic layers, where one of the layers is very thin. Modeling fractures of different thickness filled with CO₂ embedded in a background medium saturated with a stiffer fluid also shows considerable attenuation and velocity dispersion. If the fracture and background frames are the same, the equivalent medium is isotropic, but strong wave anisotropy occurs in the case of a frameless and highly permeable fracture material, for instance a suspension of solid particles in the fluid.