Surface chemistry and reactivity of biogenic silica

The surface chemistry of cultured diatoms was compared to that of biosiliceous material in Southern Ocean sediments, using potentiometric titrations and aluminum adsorption experiments. Aerosil 200, a well-studied synthetic amorphous silica, served as reference solid. Surface charge development and aluminum adsorption on cultured diatom shells were comparable to Aerosil. The surface chemical properties of biosiliceous material buried to depths of 15-25 cm in Southern Ocean sediments, however, deviated markedly from those of the cultured diatoms. In pH range 4-8.5, the surface charge density was systematically lower for biogenic silica from the sediments. In addition, the aluminum adsorption edge on the biosiliceous sediments was shifted to higher pHs by about 0.4 units. The results indicate that ionizable surface silanol groups on diagenetically altered diatom shells are less abundant and, possibly, less acidic than on freshly cultured diatoms. The observed differences in surface chemical structure are consistent with the progressive loss of reactivity, or aging, of biogenic silica which promotes its preservation in sediments.     

An Observational Study of the Mesoscale Mistral Dynamics

We investigate the mesoscale dynamics of the mistral through the wind profiler observations of the MAP (autumn 1999) and ESCOMPTE (summer 2001) field campaigns. We show that the mistral wind field can dramatically change on a time scale less than 3 hours. Transitions from a deep to a shallow mistral are often observed at any season when the lower layers are stable. The variability, mainly attributed in summer to the mistral/land–sea breeze interactions on a 10-km scale, is highlighted by observations from the wind profiler network set up during ESCOMPTE. The interpretations of the dynamical mistral structure are performed through comparisons with existing basic theories. The linear theory of R. B. Smith [Advances in Geophysics, Vol. 31, 1989, Academic Press, 1–41] and the shallow water theory [Schär, C. and Smith, R. B.: 1993a, J. Atmos. Sci. 50, 1373–1400] give some complementary explanations for the deep-to-shallow transition especially for the MAP mistral event. The wave breaking process induces a low-level jet (LLJ) downstream of the Alps that degenerates into a mountain wake, which in turn provokes the cessation of the mistral downstream of the Alps. Both theories indicate that the flow splits around the Alps and results in a persistent LLJ at the exit of the Rhône valley. The LLJ is strengthened by the channelling effect of the Rhône valley that is more efficient for north-easterly than northerly upstream winds despite the north–south valley axis. Summer moderate and weak mistral episodes are influenced by land–sea breezes and convection over land that induce a very complex interaction that cannot be accurately described by the previous theories.     

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

Aerosol polarized phase function and single-scattering albedo retrieved from ground-based measurements

Aerosol optical parameters, polarized phase function and single-scattering albdeo, have been retrieved from ground-based sun photometer measurements in Beijing 2003. The measured aerosol optical thickness varies from 0.12 to 0.77 with an average value of 0.39. The measured Ångström coefficient ranges from 0.75 to 1.47 with an average value of 1.21. The retrieved single-scattering albedo at 870 nm is within the 0.76–0.94 range and the average value is 0.85, suggests there are considerable aerosol absorptions in Beijing. The maximum value of retrieved polarized phase function at 870 nm ranges from 0.068 to 0.225 with an average value of 0.16, and it illustrates good correlations with the Ångström coefficient, i.e. the relative size of aerosol particles. Analyses of measurements and theoretical calculations show the polarized phase function is sensitive to aerosol size distribution and complex refractive index, especially the imaginary part of the refractive index which denotes aerosol light absorbing effects. These results suggest that the polarized phase function is an effective and unique aerosol optical parameter and is able to improve the retrieval of aerosol physical properties.     

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.     

Correction to: Future changes in Atlantic hurricanes with the rotated-stretched ARPEGE-Climat at very high resolution

Climate Dynamics - The original article can be found online.     

An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

Simulation of site-scale water fluxes in desert and natural oasis ecosystems of the arid region in Northwest China

The study of water fluxes is important to better understand hydrological cycles in arid regions. Data-driven machine learning models have been recently applied to water flux simulation. Previous studies have built site-scale simulation models of water fluxes for individual sites separately, requiring a large amount of data from each site and significant computation time. For arid areas, there is no consensus as to the optimal model and variable selection method to simulate water fluxes. Using data from seven flux observation sites in the arid region of Northwest China, this study compared the performance of random forest (RF), support vector machine (SVM), back propagation neural network (BPNN), and multiple linear regression (MLR) models in simulating water fluxes. Additionally, the study investigated inter-annual and seasonal variation in water fluxes and the dominant drivers of this variation at different sites. A universal simulation model for water flux was built using the RF approach and key variables as determined by MLR, incorporating data from all sites. Model performance of the SVM algorithm (R² = 0.25–0.90) was slightly worse than that of the RF algorithm (R² = 0.41–0.91); the BPNN algorithm performed poorly in most cases (R² = 0.15–0.88). Similarly, the MLR results were limited and unreliable (R² = 0.00–0.66). Using the universal RF model, annual water fluxes were found to be much higher than the precipitation received at each site, and natural oases showed higher fluxes than desert ecosystems. Water fluxes were highest during the growing season (May–September) and lowest during the non-growing season (October–April). Furthermore, the dominant drivers of water flux variation were various among different sites, but the normalized difference vegetation index (NDVI), soil moisture and soil temperature were important at most sites. This study provides useful insights for simulating water fluxes in desert and oasis ecosystems, understanding patterns of variation and the underlying mechanisms. Besides, these results can make a contribution as the decision-making basis to the water management in desert and oasis ecosystems.