2016-2017

Abstract: The objective of this paper is to present the multi-orbit (MO) surface soil moisture (SM) and angle-binned brightness temperature (TB) products for the SMOS (Soil Moisture and Ocean Salinity) mission based on a new multi-orbit algorithm. The Level 3 algorithm at CATDS (Centre Aval de Traitement des Donnees SMOS) makes use of MO retrieval to enhance the robustness and quality of SM retrievals. The motivation of the approach is to make use of the longer temporal autocorrelation length of the vegetation optical depth (VOD) compared to the corresponding SM autocorrelation in order to enhance the retrievals when an acquisition occurs at the border of the swath. The retrieval algorithm is implemented in a unique operational processor delivering multiple parameters (e.g. SM and VOD) using multi-angular dual-polarisation TB from MO. A subsidiary angle-binned TB product is provided. In this study the Level 3 TB V310 product is showcased and compared to SMAP (Soil Moisture Active Passive) TB. The Level 3 SM V300 product is compared to the single-orbit (SO) retrievals from the Level 2 SM processor from ESA with aligned configuration. The advantages and drawbacks of the Level 3 SM product (L3SM) are discussed. The comparison is done on a global scale between the two datasets and on the local scale with respect to in situ data from AMMA-CATCH and USDA ARS Watershed networks. The results obtained from the global analysis show that the MO implementation enhances the number of retrievals: up to 9% over certain areas. The comparison with the in situ data shows that the increase in the number of retrievals does not come with a decrease in quality, but rather at the expense of an increased time lag in product availability from 6 h to 3.5 days, which can be a limiting factor for applications like flood forecast but reasonable for drought monitoring and climate change studies. The SMOS L3 soil moisture and L3 brightness temperature products are delivered using an open licence and free of charge using a web application(https://www.catds.fr/sipad/).TheRE04 products, versions 300 and 310, used in this paper are also available at ftp://ext-catds-cpdc:catds2010@ftp.ifremer.fr/Land_products/GRIDDED/L3SM/RE04/.

Abstract: An increasing number of studies pointed out the ubiquitous presence of medical residues in surface and ground water as well as in soil compartments. Not only antibiotics can be found in the environment but also their transformation products about which little information is generally available. The development of bacterial resistance to antibiotics is particularly worrying as it can lead to sanitary and health problems. Studies about the dissemination of antibiotics and associated resistances in the Bolivian Altiplano are scarce. We provide baseline information on the occurrence of Sulfamethoxazole (SMX) and Trimethoprim (TMP) antibiotics as well as on the most common human SMX transformation products (TP) and on the occurrence of sulfonamide resistance genes. The studied water and soil compartments presented high levels of antibiotic pollution. This situation was shown to be mainly linked with uncontrolled discharges of treated and untreated wastewaters, resulting on the presence of antibiotics in the Titicaca Lake. SMXTPs were detected in surface waters and on soil sampled next to the wastewater treatment plant (WWTP). SMX resistance genes sull and sulll were widely detected in the basin hydrological network, even in areas unpolluted with antibiotics. Mechanisms of co-selection of antibiotic-and metal-resistance may be involved in the prevalence of ARG's in pristine areas with no anthropogenic activity and free of antibiotic pollution. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Recent advances in both ocean modeling and melt parameterization in ice-sheet models point the way toward coupled ice sheet�ocean modeling, which is needed to quantify Antarctic mass loss and the resulting sea-level rise. The latest Antarctic ocean modeling shows that complex interactions between the atmosphere, sea ice, icebergs, bathymetric features, and ocean circulation on many scales determine which water masses reach ice-shelf cavities and how much heat is available to melt ice. Meanwhile, parameterizations of basal melting in standalone ice-sheet models have evolved from simplified, depth-dependent functions to more sophisticated models, accounting for ice-shelf basal topography, and the evolution of the sub-ice-shelf buoyant flow. The focus of recent work has been on better understanding processes or adding new model capabilities, but a broader community effort is needed in validating models against observations and producing melt-rate projections. Given time, community efforts in coupled ice sheet�ocean modeling, already underway, will tackle the considerable challenges involved in building, initializing, constraining, and performing projections with coupled models, leading to reduced uncertainties in Antarctica�s contribution to future sea-level rise.

Abstract: Urban subsoil provides space, groundwater, geothermal energy and geomaterials, and the importance underlying the management of this strategic space is well-known. In the context of vertical urban growth, groundwater management is a crucial issue. In particular, operational tools dedicated to urban planners are needed in urban areas prone to conflicts over use in order to aid the decision-making process. The aim of this paper is to propose guidelines dedicated to urban underground development in view to incorporating groundwater resource management. These guidelines are supported by analyses of the scientific literature dealing with: (1) the impacts caused by underground structures on urban groundwater, and (2) the quantitative overview of urban groundwater resources at the operational level in terms of quantity, quality and temperature. Finally, perspectives are discussed regarding the research work needed to ensure better management of urban groundwater resources confronted by underground development. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: The date of the definitive start to the Anthropocene is still under debate, and although a lot of progress has been made, currently available data is not precise enough to define the start of the human-dominated geological epoch. We know that during and after the industrial revolution, humans started having a much greater impact on the Earth's environment. Increases in population have led to increases in resource and fossil fuel use, leaving a marked impact on our planet. This impact, in the Northern hemisphere, is effectively recorded in the snow and ice of the Arctic. Human activity has changed various biogeochemical cycles to such an extent, that the climate has started to change. This has disturbed the biosphere, pushing it to adapt in response to these changes through evolutionary pressure. The Arctic is a particularly vulnerable environment and mankind is having a profound impact on its fragile equilibrium. Higher concentrations of heavy metals, organic compounds and radionuclides have been detected in ice cores as well as snow. Although climatic changes are evident on a global scale, in the Arctic these changes have been amplified. Advances in laboratory analysis methods have been applied to ice cores and surface snow samples to help us understand the mechanisms governing this fragile environment and to evaluate the impact and amplitude of human activity. Despite these advances, the fluxes and distributions over time of anthropogenic organic compounds is largely unknown. Hopefully advances in analytical methods will mean that this is not the case in the future. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Climate change projections still suffer from a limited representation of the permafrost-carbon feedback. Predicting the response of permafrost temperature to climate change requires accurate simulations of Arctic snow and soil properties. This study assesses the capacity of the coupled land surface and snow models ISBA-Crocus and ISBA-ES to simulate snow and soil properties at Bylot Island, a high Arctic site. Field measurements complemented with ERA-Interim reanalyses were used to drive the models and to evaluate simulation outputs. Snow height, density, temperature, thermal conductivity and thermal insulance are examined to determine the critical variables involved in the soil and snow thermal regime. Simulated soil properties are compared to measurements of thermal conductivity, temperature and water content. The simulated snow density profiles are unrealistic, which is most likely caused by the lack of representation in snow models of the upward water vapor fluxes generated by the strong temperature gradients within the snowpack. The resulting vertical profiles of thermal conductivity are inverted compared to observations, with high simulated values at the bottom of the snowpack. Still, ISBA-Crocus manages to successfully simulate the soil temperature in winter. Results are satisfactory in summer, but the temperature of the top soil could be better reproduced by adequately representing surface organic layers, i.e., mosses and litter, and in particular their water retention capacity. Transition periods (soil freezing and thawing) are the least well reproduced because the high basal snow thermal conductivity induces an excessively rapid heat transfer between the soil and the snow in simulations. Hence, global climate models should carefully consider Arctic snow thermal properties, and especially the thermal conductivity of the basal snow layer, to perform accurate predictions of the permafrost evolution under climate change.

Abstract: In the verification field, stratification is the process of dividing the sample of forecast-observation pairs into quasi-homogeneous subsets, in order to learn more on how forecasts behave under specific conditions. A general framework for stratification is presented for the case of ensemble forecasts of continuous scalar variables. Distinction is made between forecast-based, observation-based, and external-based stratification, depending on the criterion on which the sample is stratified. The formalism is applied to two widely used verification measures: the continuous ranked probability score (CRPS) and the rank histogram. For both, new graphical representations that synthesize the added information are proposed. Based on the definition of calibration, it is shown that the rank histogram should be used within a forecast-based stratification, while an observation-based stratification leads to significantly nonflat histograms for calibrated forecasts. Nevertheless, as previous studies have warned, statistical artifacts created by a forecast-based stratification may still occur, thus a graphical test to detect them is suggested. To illustrate potential insights about forecast behavior that can be gained from stratification, a numerical example with two different datasets of mean areal precipitation forecasts is presented.

Abstract: This paper investigates the combined use of extensive geotechnical, hydrogeological and geophysical techniques to assess a small earth dyke with a permanent hydraulic head, namely a canal embankment. The experimental site was chosen because of known issues regarding internal erosion and piping phenomena. Two leakages were visually located following the emptying of the canal prior to remediation works. The results showed a good agreement between the geophysical imaging techniques (Electrical Resistivity Tomography, P- and SH-waves Tomography) and the geotechnical data to detect the depth to the bedrock and its lateral variations. It appeared that surface waves might not be fully adapted for dyke investigation because of the particular geometry of the studied dyke, non-respectful of the 10 assumption, and which induced depth and velocity discrepancies retrieved from Rayleigh and Love waves inversion. The use of these classical prospecting techniques however did not allow to directly locate the two leakages within the studied earth dyke. The analysis of ambient vibration time series with a modified beam-forming algorithm allowed to localize the most energetic water flow prior to remediation works. It was not possible to detect the leakage after remediation works, suggesting that they efficiently contributed to significantly reduce the water flow. The second leakage was not detected probably because of a non-turbulent water flow, generating few energetic vibrations. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Atmospheric methane (CH4) records reconstructed from polar ice cores represent an integrated view on processes predominantly taking place in the terrestrial biogeosphere. Here, we present dual stable isotopic methane records [delta(CH4)-C-13 and delta D(CH4)] from four Antarctic ice cores, which provide improved constraints on past changes in natural methane sources. Our isotope data show that tropical wetlands and seasonally inundated floodplains are most likely the controlling sources of atmospheric methane variations for the current and two older interglacials and their preceding glacial maxima. The changes in these sources are steered by variations in temperature, precipitation, and the water table as modulated by insolation, (local) sea level, and monsoon intensity. Based on our delta D(CH4) constraint, it seems that geologic emissions of methane may play a steady but only minor role in atmospheric CH4 changes and that the glacial budget is not dominated by these sources. Superimposed on the glacial/interglacial variations is a marked difference in both isotope records, with systematically higher values during the last 25,000 y compared with older time periods. This shift cannot be explained by climatic changes. Rather, our isotopic methane budget points to a marked increase in fire activity, possibly caused by biome changes and accumulation of fuel related to the late Pleistocene megafauna extinction, which took place in the course of the last glacial.

Abstract: The widespread use of Aerodyne aerosol mass spectrometers (AMS) has greatly improved real-time organic aerosol (OA) monitoring, providing mass spectra that contain sufficient information for source apportionment. However, AMS field deployments remain expensive and demanding, limiting the acquisition of long-term datasets at many sampling sites. The offline application of aerosol mass spectrometry entailing the analysis of nebulized water extracted filter samples (offline-AMS) increases the spatial coverage accessible to AMS measurements, being filters routinely collected at many stations worldwide. PM1 (particulate matter with an aerodynamic diameter <1 μm) filter samples were collected during an entire year in Lithuania at three different locations representative of three typical environments of the southeast Baltic region: Vilnius (urban background), Rugsteliskis (rural terrestrial), and Preila (rural coastal). Aqueous filter extracts were nebulized in Ar, yielding the first AMS measurements of water-soluble atmospheric organic aerosol (WSOA) without interference from air fragments. This enables direct measurement of the CO+ fragment contribution, whose intensity is typically assumed to be equal to that of CO2+. Offline-AMS spectra reveal that the water-soluble CO2+:CO+ ratio not only shows values systematically >1 but is also dependent on season, with lower values in winter than in summer. AMS WSOA spectra were analyzed using positive matrix factorization (PMF), which yielded four factors. These factors included biomass burning OA (BBOA), local OA (LOA) contributing significantly only in Vilnius, and two oxygenated OA (OOA) factors, summer OOA (S-OOA) and background OOA (B-OOA), distinguished by their seasonal variability. The contribution of traffic exhaust OA (TEOA) was not resolved by PMF due to both low concentrations and low water solubility. Therefore, the TEOA concentration was estimated using a chemical mass balance approach, based on the concentrations of hopanes, specific markers of traffic emissions. AMS-PMF source apportionment results were consistent with those obtained from PMF applied to marker concentrations (i.e., major inorganic ions, OC/EC, and organic markers including polycyclic aromatic hydrocarbons and their derivatives, hopanes, long-chain alkanes, monosaccharides, anhydrous sugars, and lignin fragmentation products). OA was the largest fraction of PM1 and was dominated by BBOA during winter with an average concentration of 2 μg m(-3) (53% of OM), while S-OOA, probably related to biogenic emissions, was the prevalent OA component during summer with an average concentration of 1.2 μg m(-3) (45% of OM). PMF ascribed a large part of the CO+ explained variability ( 97 %) to the OOA and BBOA factors. Accordingly, we discuss a new CO+ parameterization as a function of CO2+ and C2H4O2+ fragments, which were selected to describe the variability of the OOA and BBOA factors.

Abstract: The transformation of snow into ice is a complex phenomenon that is difficult to model. Depending on surface temperature and accumulation rate, it may take several decades to millennia for air to be entrapped in ice. The air is thus always younger than the surrounding ice. The resulting gas-ice age difference is essential to documenting the phasing between CO2 and temperature changes, especially during deglaciations. The air trapping depth can be inferred in the past using a firn densification model, or using delta N-15 of air measured in ice cores. All firn densification models applied to deglaciations show a large disagreement with delta N-15 measurements at several sites in East Antarctica, predicting larger firn thickness during the Last Glacial Maximum, whereas delta N-15 suggests a reduced firn thickness compared to the Holocene. Here we present modifications of the LGGE firn densification model, which significantly reduce the model-data mismatch for the gas trapping depth evolution over the last deglaciation at the coldest sites in East Antarctica (Vostok, Dome C), while preserving the good agreement between measured and modelled modern firn density profiles. In particular, we introduce a dependency of the creep factor on temperature and impurities in the firn densification rate calculation. The temperature influence intends to reflect the dominance of different mechanisms for firn compaction at different temperatures. We show that both the new temperature parameterization and the influence of impurities contribute to the increased agreement between modelled and measured delta N-15 evolution during the last deglaciation at sites with low temperature and low accumulation rate, such as Dome C or Vostok. We find that a very low sensitivity of the densification rate to temperature has to be used in the coldest conditions. The inclusion of im-purity effects improves the agreement between modelled and measured delta N-15 at cold East Antarctic sites during the last deglaciation, but deteriorates the agreement between modelled and measured delta N-15 evolution at Greenland and Antarctic sites with high accumulation unless threshold effects are taken into account. We thus do not provide a definite solution to the firnification at very cold Antarctic sites but propose potential pathways for future studies.

Abstract: This paper quantifies the impacts of approximations and assumptions in the parameterization of bulk formulas on the exchange of momentum, heat, and freshwater computed between the ocean and atmosphere. An ensemble of sensitivity experiments is examined. Climatologies of wind stress, turbulent heat flux, and evaporation for the period 1982-2014 are computed using SST and surface meteorological state variables from ERA-Interim. Each experiment differs from the defined control experiment in only one aspect of the parameterization of the bulk formulas. The wind stress is most sensitive to the closure used to relate the neutral drag coefficient to the wind speed in the bulk algorithm, which mainly involves the value of the Charnock parameter. The disagreement between the state-of-the-art algorithms examined is typically on the order of 10%. The largest uncertainties in turbulent heat flux and evaporation are also related to the choice of the algorithm (typically 15%) but also emerge in experiments examining approximations related to the surface temperature and saturation humidity. Thus, approximations for the skin temperature and the saltrelated reduction of saturation humidity have a substantial impact on the heat flux and evaporation (typically 10%). Approximations such as the use of a fixed air density, sea level pressure, or simplified formula for the saturation humidity lead to errors no larger than 4% when tested individually. The impacts of these approximations combine linearly when implemented together, yielding errors up to 20% over mid-and subpolar latitudes.

Abstract: High Mountain Asia hosts the largest glacier concentration outside the polar regions. These glaciers are important contributors to streamflow in one of the most populated areas of the world. Past studies have used methods that can provide only regionally averaged glacier mass balances to assess the glacier contribution to rivers and sea level rise. Here we compute the mass balance for about 92% of the glacierized area of High Mountain Asia using time series of digital elevation models derived from satellite stereo-imagery. We calculate a total mass change of -16.3 +/- 3.5 Gt yr(-1) (-0.18 +/- 0.04 m w.e. yr(-1)) between 2000 and 2016, which is less negative than most previous estimates. Region-wide mass balances vary from 4.0 +/- 1.5 Gt yr(-1) (-0.62 +/- 0.23 m w.e. yr(-1)) in Nyainqentanglha to +1.4 +/- 0.8 Gt yr(-1) (+0.14 +/- 0.08 m w.e. yr(-1)) in Kunlun, with large intra-regional variability of individual glacier mass balances (standard deviation within a region similar to 0.20m w.e. yr(-1)). Specifically, our results shed light on the Nyainqentanglha and Pamir glacier mass changes, for which contradictory estimates exist in the literature. They provide crucial information for the calibration of the models used for projecting glacier response to climatic change, as these models do not capture the pattern, magnitude and intra-regional variability of glacier changes at present.

Abstract: The East Australian Current (EAC), like many other subtropical western boundary currents, is believed to be penetrating further poleward in recent decades. Previous observational and model studies have used steady state dynamics to relate changes in the westerly winds to changes in the separation behavior of the EAC. As yet, little work has been undertaken on the impact of forcing variability on the EAC and Tasman Sea circulation. Here using an eddy-permitting regional ocean model, we present a suite of simulations forced by the same time-mean fields, but with different atmospheric and remote ocean variability. These eddy-permitting results demonstrate the nonlinear response of the EAC to variable, nonstationary inhomogeneous forcing. These simulations show an EAC with high intrinsic variability and stochastic eddy shedding. We show that wind stress variability on time scales shorter than 56 days leads to increases in eddy shedding rates and southward eddy propagation, producing an increased transport and southward reach of the mean EAC extension. We adopt an energetics framework that shows the EAC extension changes to be coincident with an increase in offshore, upstream eddy variance (via increased barotropic instability) and increase in subsurface mean kinetic energy along the length of the EAC. The response of EAC separation to regional variable wind stress has important implications for both past and future climate change studies.

Abstract: The contact between two particles of ice is studied experimentally and modeled. Experiments on the gradual indentation of a single or on two monocrystalline cylinders of ice demonstrate the effect of the relative orientation of their basal planes with the loading axis. Using polarized light and Digital Image Correlation, the main deformation mechanisms are documented. Experimental observations show the strong anisotropy of deformation, with strain localization into shear bands, corresponding to basal gliding of dislocations. When the rotation of cylinders is hindered by a large initial contact size, a simple model is proposed. It is based on indentation theory and takes into account the preferential viscoplastic deformation on the basal plane. The model form is validated by experimental data and by finite element simulations that incorporate the transverse isotropic law of an ice crystal. (C) 2017 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Abstract: Snow fabric, defined as the distribution of the c-axis orientations of the ice crystals in snow, is poorly known. So far, only one study exits that measured snow fabric based on a statistically representative technique. This recent study has revealed the impact of temperature gradient metamorphism on the evolution of fabric in natural snow, based on cold laboratory experiments. On polar ice sheets, snow properties are currently investigated regarding their strong variability in time and space, notably because of their potential influence on firn processes and consequently on ice core analysis. Here, we present measurements of fabric and microstructure of snow from Point Barnola, East Antarctica (close to Dome C). We analyzed a snow profile from 0 to 3 m depth, where temperature gradients occur. The main contributions of the paper are (1) a detailed characterization of snow in the upper meters of the ice sheet, especially by providing data on snow fabric, and (2) the study of a fundamental snow process, never observed up to now in a natural snowpack, namely the role of temperature gradient metamorphism on the evolution of the snow fabric. Snow samples were scanned by micro-tomography to measure continuous profiles of microstructural properties (density, specific surface area and pore thickness). Fabric analysis was performed using an automatic ice texture analyzer on 77 representative thin sections cut out from the samples. Different types of snow fabric could be identified and persist at depth. Snow fabric is significantly correlated with snow microstructure, pointing to the simultaneous influence of temperature gradient metamorphism on both properties. We propose a mechanism based on preferential grain growth to explain the fabric evolution under temperature gradients. Our work opens the question of how such a layered profile of fabric and microstructure evolves at depth and further influences the physical and mechanical properties of snow and firn. More generally, it opens the way to further studies on the influence of the snow fabric in snow processes related to anisotropic properties of ice such as grain growth, mechanical response, electromagnetic behavior. (C) 2016 The Author(s). Published by Elsevier B.V.

Abstract: Nucleation mechanisms occurring during dynamic recrystallization play a crucial role in the evolution of microstructures and textures during high temperature deformation. In polycrystalline ice, the strong viscoplastic anisotropy induces high strain heterogeneities between grains which control the recrystallization mechanisms. Here, we study the nucleation mechanisms occurring during creep tests performed on polycrystalline columnar ice at high temperature and stress (T=-5 degrees C; sigma = 0.5 MPa) by post-mortem analyses of deformation microstructures using cryogenic electron backscatter diffraction. The columnar geometry of the samples enables discrimination of the nuclei from the initial grains. Various nucleation mechanisms are deduced from the analysis of the nuclei relations with the dislocation sub-structures within grains and at grain boundaries. Tilt sub-grain boundaries and kink bands are the main structures responsible for development of polygonization and mosaic sub-structures. Nucleation by bulging at serrated grain boundaries is also an efficient nucleation mechanism near the grain boundaries where strain incompatibilities are high. Observation of nuclei with orientations not related to the 'parent' ones suggests the possibility of 'spontaneous' nucleation driven by the relaxation of the dislocation-related internal stress field. The complexity of the nucleation mechanisms observed here emphasizes the impact of stress and strain heterogeneities on dynamic recrystallization mechanisms. This article is part of the themed issue 'Microdynamics of ice'.

Abstract: Prediction of ice mass flow and associated dynamics is pivotal at a time of climate change. Ice flow is dominantly accommodated by the motion of crystal defects – the dislocations. In the specific case of ice, their observation is not always accessible by means of the classical tools such as X-ray diffraction or transmission electron microscopy (TEM). Part of the dislocation population, the geometrically necessary dislocations (GNDs) can nevertheless be constrained using crystal orientation measurements via electron backscattering diffraction (EBSD) associated with appropriate analyses based on the Nye (1950) approach. The present study uses the Weighted Burgers Vectors, a reduced formulation of the Nye theory that enables the characterization of GNDs. Applied to ice, this method documents, for the first time, the presence of dislocations with non-basal [c] or < c + a > Burgers vectors. These [c] or (c + a) dislocations represent up to 35% of the GNDs observed in laboratory-deformed ice samples. Our findings offer a more complex and comprehensive picture of the key plasticity processes responsible for polycrystalline ice creep and provide better constraints on the constitutive mechanical laws implemented in ice sheet flow models used to predict the response of Earth ice masses to climate change. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: We present a new water percolation routine added to the one-dimensional snowpack model Crocus as an alternative to the empirical bucket routine. This routine solves the Richards equation, which describes flow of water through unsaturated porous snow governed by capillary suction, gravity and hydraulic conductivity of the snow layers. We tested the Richards routine on two data sets, one recorded from an automatic weather station over the winter of 2013-2014 at Filefjell, Norway, and the other an idealized synthetic data set. Model results using the Richards routine generally lead to higher water contents in the snow layers. Snow layers often reached a point at which the ice crystals' surface area is completely covered by a thin film of water (the transition between pendular and funicular regimes), at which feedback from the snow metamorphism and compaction routines are expected to be nonlinear. With the synthetic simulation 18% of snow layers obtained a saturation of > 10% and 0.57% of layers reached saturation of > 15 %. The Richards routine had a maximum liquid water content of 173.6 kg m(-3) whereas the bucket routine had a maximum of 42.1 kg m(-3). We found that wet-snow processes, such as wet-snow metamorphism and wet-snow compaction rates, are not accurately represented at higher water contents. These routines feed back on the Richards routines, which rely heavily on grain size and snow density. The parameter sets for the water retention curve and hydraulic conductivity of snow layers, which are used in the Richards routine, do not represent all the snow types that can be found in a natural snowpack. We show that the new routine has been implemented in the Crocus model, but due to feedback amplification and parameter uncertainties, meaningful applicability is limited. Updating or adapting other routines in Crocus, specifically the snow compaction routine and the grain metamorphism routine, is needed before Crocus can accurately simulate the snowpack using the Richards routine.

Abstract: While the impacts of climate change on individual species and communities have been well documented there is little evidence on climate-mediated changes for entire ecosystems. Pristine alpine environments can provide unique insights into natural, physical and ecological response to climate change yet broad scale and long-term studies on these potential 'ecosystem sentinels' are scarce. We addressed this issue by examining cover changes of 1689 high-elevation wetlands (temporarily or perennial water-saturated grounds) in the Bolivian Cordillera Real, a region that has experienced significant warming and glacier melting over the last 30 years. We combined high spatial resolution satellite images from PLEIADES with the long-term images archive from LANDSAT to 1) examine environmental factors (e.g., glacier cover, wetland and watershed size) that affected wetland cover changes, and 2) identify wetlands' features that affect their vulnerability (using habitat drying as a proxy) in the face of climate change. Over the (1984-2011) period, our data showed an increasing trend in the mean wetland total area and number, mainly related to the appearance of wet grassland patches during the wetter years. Wetland cover also showed high inter-annual variability and their area for a given year was positively correlated to precipitation intensities in the three months prior to the image date. Also, round wetlands located in highly glacierized catchments were less prone to drying, while relatively small wetlands with irregularly shaped contours suffered the highest rates of drying over the last three decades. High Andean wetlands can therefore be considered as ecosystem sentinels for climate change, as they seem sensitive to glacier melting. Beyond the specific focus of this study, our work illustrates how satellite-based monitoring of ecosystem sentinels can help filling the lack of information on the ecological consequences of current and changing climate conditions, a common and crucial issue especially in less-developed countries.

Abstract: Rainfall is an important driver of erosion processes. The mean rainfall rate is often used to account for the erosive impact of a particular climate. However, for some erosion processes, erosion rate is a nonlinear function of rainfall, e.g., due to a threshold for erosion. When this is the case, it is important to take into account the full distribution of rainfall, instead of just the mean. In light of this, we have characterized the variability of daily rainfall over the Himalayan orogen using high spatial and temporal resolution rainfall data sets. We find significant variations in rainfall variability over the Himalayan orogen, with increasing rainfall variability to the west and north of the orogen. By taking into account variability of rainfall in addition to mean rainfall rate, we find a pattern of rainfall that, from a geomorphological perspective, is significantly different from mean rainfall rate alone. Using these findings, we argue that short-term rainfall variability may help explain observed short and long-term erosion rates in the Himalayan orogen. Plain Language Summary An important topic in earth science is understanding how climate and tectonic forces interact to shape the surface of the Earth. One of the main influences that climate has on the Earth's surface is to cause erosion by delivering water to landscapes as rain and snow. Wetter climates should, in general, cause more erosion than drier climates, which is why the mean annual rainfall has traditionally been used as a measure of the erosive strength of climate. However, field evidence seems to suggest that the effect of climate is more sophisticated than can be captured with the mean annual rainfall. There is a considerable body of theory demonstrating that the intensity of rainfall, or storminess, is an equally important aspect of how climate causes erosion. In light of this, we have characterized the storminess in the Himalayan mountain range, a place with high erosion rates driven by rainfall. We show that there is large variation in storminess from place to place in the Himalayas, and that this can help explain observed rates of erosion, which do not match mean annual rainfall rates.

Abstract: Forecasting the occurrence of debris flows is fundamental for issuing hazard warnings, and often focuses on rainfall as a triggering agent and on the use of empirical rainfall thresholds based on rain gauge observations. A recognized component of the uncertainty associated with the use of rainfall thresholds is related to the sampling of strongly varying rainfall variability with sparse rain gauge networks. In this work we examine the spatial distribution of rainfall depth in areas up to 10 km from the debris flow initiation points as a function of return period, and we exploit this information to analyze the errors expected in the estimation of debris flow triggering rainfall when rain gauge data are used. In particular, we investigate the impact of rain gauge density and of the use of different interpolation methods. High-resolution, adjusted radar rainfall estimates, representing the best available spatially-distributed-rainfall estimates at the debris flows initiation point and in the surrounding area, are sampled by stochastically generated rain gauge networks characterized by varying densities. Debris flow triggering rainfall is estimated by means of three rainfall interpolation methods: nearest neighbor, inverse distance weighting and ordinary kriging. On average, triggering rainfall shows a local peak corresponding to the debris flow initiation point, with a decay of rainfall with distance which increases with the return period of the triggering rainfall. Interpolation of the stochastically generated rain gauge measurements leads to an underestimation of the triggering rainfall that, irrespective of the interpolation methods, increases with the return period and decreases with the rain gauge density. For small return period events and high rain gauge density, the differences among the methods are minor. With increasing the return period and decreasing the rain gauge density, the nearest neighbor method is less biased, because it makes use only of the closest rain gauge to the debris flow initiation point. On the contrary, inverse distance weighting and ordinary kriging, which are using rain gauges located farther from the debris flows in addition to the closest one, exhibit negative biases that increase with return period. The standard deviation of the interpolated values is larger when the nearest neighbor is used with respect to inverse distance weighting and ordinary kriging. For large return period and low rain gauge density, the differences among the methods are minor. (C) 2016 Elsevier B.V. All ri
gths reserved.

Abstract: In this study, the relationship between trains of African easterly waves (AEWs) and downstream tropical cyclogenesis is studied. Based on 19 summer seasons (July-September from 1990 to 2008) of ERA-Interim reanalysis fields and brightness temperature from the Cloud User Archive, the signature of AEW troughs and embedded convection are tracked from the West African coast to the central Atlantic. The tracked systems are separated into four groups: (i) systems originating from the north zone of the midtropospheric African easterly jet (AEJ), (ii) those coming from the south part of AEJ, (iii) systems that are associated with a downstream trough located around 2000 km westward (termed DUO systems), and (iv) those that are not associated with such a close downstream trough (termed SOLO systems). By monitoring the embedded 700-hPa-filtered relative vorticity and 850-hPa wind convergence anomaly associated with these families along their trajectories, it is shown that the DUO generally have stronger dynamical structure and statistically have a longer lifetime than the SOLO ones. It is suggested that the differences between them may be due to the presence of the previous intense downstream trough in DUO cases, enhancing the low-level convergence behind them. Moreover, a study of the relationship between system trajectories and tropical depressions occurring between the West African coast and 40 degrees W showed that 90% of tropical depressions are identifiable from the West African coast in tracked systems, mostly in the DUO cases originating from the south zone of the AEJ.

Abstract: The influence of discrete and randomly distributed fibers (DRDF) on the improvement of the resistance to cracking of clay-based landfill covers was studied by conducting a series of centrifuge model tests. A motor-based differential settlement simulator (MDSS) was used to induce differential settlements to the soil layers during centrifuge tests performed at 40g using the 4.5 m radius beam centrifuge at Indian Institute of Technology Bombay. A digital image cross-correlation (DIC) technique was used to capture deformations of the model barriers and to obtain strain field distributions. Thickness of the unreinforced barriers (URSB) and fiber-reinforced barriers (FRSB) were varied as 15 mm (similar to 0.6 m in the field) and 25 mm (1.0 m). Fiber lengths were varied as 30, 60, and 90 mm. The strain at water breakthrough of 15 mm (0.6 m) and 25 mm (1 m) FRSB were respectively 1.9 and 2.2 times higher compared to identical URSBs. A relatively crack-free surface was observed for FRSB with longer fibers even after subjected to differential settlement of 1 m. (C) 2016 American Society of Civil Engineers.

Abstract: Simulation of vegetation-climate feedbacks in high latitudes in the ORCHIDEE land surface model was improved by the addition of three new circumpolar plant functional types (PFTs), namely non-vascular plants representing bryophytes and lichens, Arctic shrubs and Arctic C-3 grasses. Non-vascular plants are assigned no stomatal conductance, very shallow roots, and can desiccate during dry episodes and become active again during wet periods, which gives them a larger phenological plasticity (i.e. adaptability and resilience to severe climatic constraints) compared to grasses and shrubs. Shrubs have a specific carbon allocation scheme, and differ from trees by their larger survival rates in winter, due to protection by snow. Arctic C3 grasses have the same equations as in the original ORCHIDEE version, but different parameter values, optimised from in situ observations of biomass and net primary productivity (NPP) in Siberia. In situ observations of living biomass and productivity from Siberia were used to calibrate the parameters of the new PFTs using a Bayesian optimisation procedure. With the new PFTs, we obtain a lower NPP by 31% (from 55 degrees N), as well as a lower roughness length (-41%), transpiration (-33%) and a higher winter albedo (by +3.6%) due to increased snow cover. A simulation of the water balance and runoff and drainage in the high northern latitudes using the new PFTs results in an increase of fresh water discharge in the Arctic ocean by 11% (+140 km(3) yr(-1)), owing to less evapotranspiration. Future developments should focus on the competition between these three PFTs and boreal tree PFTs, in order to simulate their area changes in response to climate change, and the effect of carbon-nitrogen interactions.

Abstract: Snow spectral albedo in the visible/near-infrared range has been continuously measured during a winter season at Col de Porte alpine site (French Alps; 45.30(circle) N, 5.77(circle) E; 1325ma.s.l.). The evolution of such alpine snowpack is complex due to intensive precipitation, rapid melt events and Saharan dust deposition outbreaks. This study highlights that the resulting intricate variations of spectral albedo can be successfully explained by variations of the following snow surface variables: specific surface area (SSA) of snow, effective light-absorbing impurities content, presence of liquid water and slope. The methodology developed in this study disentangles the effect of these variables on snow spectral albedo. The presence of liquid water at the snow surface results in a spectral shift of the albedo from which melt events can be identified with an occurrence of false detection rate lower than 3.5 %. Snow SSA mostly impacts spectral albedo in the near-infrared range. Impurity deposition mostly impacts the albedo in the visible range but this impact is very dependent on snow SSA and surface slope. Our work thus demonstrates that the SSA estimation from spectral albedo is affected by large uncertainties for a tilted snow surface and medium to high impurity contents and that the estimation of impurity content is also affected by large uncertainties, especially for low values below 50 ng g(-1) black carbon equivalent. The proposed methodology opens routes for retrieval of SSA, impurity content, melt events and surface slope from spectral albedo. However, an exhaustive accuracy assessment of the snow black properties retrieval would require more independent in situ measurements and is beyond the scope of the present study. This time series of snow spectral albedo nevertheless already provides a new insight into our understanding of the evolution of snow surface properties.

Abstract: Satellite sensors increasingly provide high-resolution (HR) observations of the ocean. They supply observations of sea surface height (SSH) and of tracers of the dynamics such as sea surface salinity (SSS) and sea surface temperature (SST). In particular, the Surface Water Ocean Topography (SWOT) mission will provide measurements of the surface ocean topography at very high-resolution (HR) delivering unprecedented information on the meso-scale and submeso-scale dynamics. This study investigates the feasibility to use these measurements to reconstruct meso-scale features simulated by numerical models, in particular on the vertical dimension. A methodology to reconstruct three-dimensional (3D) multivariate meso-scale scenes is developed by using a HR numerical model of the Solomon Sea region. An inverse problem is defined in the framework of a twin experiment where synthetic observations are used. A true state is chosen among the 3D multivariate states which is considered as a reference state. In order to correct a first guess of this true state, a two-step analysis is carried out. A probability distribution of the first guess is defined and updated at each step of the analysis: (i) the first step applies the analysis scheme of a reduced-order Kalman filter to update the first guess probability distribution using SSH observation; (ii) the second step minimizes a cost function using observations of HR image structure and a new probability distribution is estimated. The analysis is extended to the vertical dimension using 3D multivariate empirical orthogonal functions (EOFs) and the probabilistic approach allows the update of the probability distribution through the two-step analysis. Experiments show that the proposed technique succeeds in correcting a multivariate state using meso-scale and submeso-scale information contained in HR SSH and image structure observations. It also demonstrates how the surface information can be used to reconstruct the ocean state below the surface.

Abstract: Since the early 80s, a progressively increasing number of sinkholes appeared along the Dead Sea coastal line. It has been found that their appearance is strongly correlating with the lowering of the Dead Sea level taking place with the rate of approximately 1 m/yr. Location of areas affected by sinkhole development corresponds to location of the salt formation deposited during the latest Pleistocene, when the Lake Lisan receded to later become the Dead Sea. Water flowing to the Dead Sea from adjacent and underlying aquifers dissolves salt and creates caverns that cause ground subsidence and consequent formation of sinkholes. Before subsidence, these caverns are not visible on the surface but can be investigated with surface geophysical methods. For that, we applied Surface Nuclear Magnetic Resonance (SNMR), Transient Electromagnetic (TEM) Seismic refraction and reflection, Multichannel Analysis of Surface waves (MASW), microgravity and magnetic surveys and their combinations. Our geophysical results allowed us to locate the salt formation and to detect caverns in salt thus contributing to better understanding sinkhole development mechanisms. Comparison of sinkhole appearance along the western DS shore derived from the recent database (2017) shows that predictions made on the base of geophysical data (2005-2008) are now confirmed thus demonstrating efficiency of our study. In this paper, we briefly present a summary of up to date knowledge of the geology and hydrogeology of Dead Sea basin, of the physical properties of the salt rock and the most popular models explaining mechanisms of sinkhole development. We also share our experience gained during geophysical studies carried out in the framework of national and international research projects in this area for the last 20 years.

Abstract: We review recent literature on atmospheric, surface ocean and sea-ice observations and modeling results in the Antarctic sector and relate the observed climatic trends with the potential changes in the surface mass balance (SMB) of the ice sheet since 1900. Estimates of regional scale SMB distribution and trends remain subject to large uncertainties. Approaches combining and comparing multiple satellite and model-based assessments of ice sheet mass balance aim at reducing these knowledge gaps. During the last decades, significant changes in atmospheric circulation occurred around Antarctica, due to the exceptional positive trend in the Southern Annular Mode and to the climate variability observed in the tropical Pacific at the end of the twentieth century. Even though climate over the East Antarctic Ice-Sheet remained quite stable, a warming and precipitation increase was observed over the West Antarctic Ice-Sheet and over the West Antarctic Peninsula (AP) during the twentieth century. However, the high regional climate variability overwhelms climate changes associated to human drivers of global temperature changes, as reflected by a slight recent decadal cooling trend over the AP. Climate models still fail to accurately reproduce the multi-decadal SMB trends at a regional scale, and progress has to be achieved in reproducing atmospheric circulation changes related to complex ocean/ice/atmosphere interactions. Complex processes are also still insufficiently considered, such as (1) specific polar atmospheric processes (clouds, drifting snow, and stable boundary layer physics), (2) surface firn physics involved in the surface drag variations, or in firn air depletion and albedo feedbacks. Finally, progress in reducing the uncertainties relative to projections of the future SMB of Antarctica will largely depend on climate model capability to correctly consider teleconnections with low and mid-latitudes, and on the ability to correct them for biases, taking into account the coupling between ocean, ice, and atmosphere in high southern latitudes.

Abstract: The wide range of ecological goods and services provided by tropical coastal lagoons and wetlands are under considerable pressure due to the synergistic effects of local anthropogenic impact and global climate change. In transitional waters, salinity is a key driver of ecological processes mostly depending on the balance between marine and river inputs, a balance that can be significantly modified by climate change and by anthropogenic alteration of the watershed. Mesoamerica being considered as a climate change hot-spot and as an ecoregion strongly vulnerable to global change, our study aimed at analyzing the relationship between salinity, river runoff, and rainfall variability in a tropical coastal lagoon and to assess the respective influence of climate change and watershed management. The study focusing on the large and shallow coastal lagoon of Laguna de Terminos in south eastern Mexico established: (1) the variability in salinity distribution along the yearly cycle and the occurrence of a high salinity anomaly period during the wet season of 2009; (2) the relationship between lagoon waters salinity and river inputs further underlying the anomalous situation encountered in 2009; (3) a long term increase in river discharge during the past 60 years, indicating potential salinity decrease in the lagoon during that same period; (4) an absence of any change in rainfall linking the increase in runoff to watershed management rather than long term trend in climate change. Additionally, the specific context of the 2009-2010 Central-Equatorial Pacific El Nio is underlined and the potential relationship between river discharge and ENSO is discussed. Those results should be of significant practical value to decision-makers who are often keen to point the finger at global climate change when local environmental management is also and sometime most significantly responsible.

Abstract: Thanks to its huge water storage capacity, Norway has an excess of energy generation at annual scale, although significant regional disparity exists. On average, the Mid-Norway region has an energy deficit and needs to import more electricity than it exports. We show that this energy deficit can be reduced with an increase in wind generation and transmission line capacity, even in future climate scenarios where both mean annual temperature and precipitation are changed. For the considered scenarios, the deficit observed in winter disappears, i.e., when electricity consumption and prices are high. At the annual scale, the deficit behaviour depends more on future changes in precipitation. Another consequence of changes in wind production and transmission capacity is the modification of electricity exchanges with neighbouring regions which are also modified both in terms of average, variability and seasonality.

Abstract: The basal topography is largely unknown beneath most glaciers and ice caps, and many attempts have been made to estimate a thickness field from other more accessible information at the surface. Here, we present a two-step reconstruction approach for ice thickness that solves mass conservation over single or several connected drainage basins. The approach is applied to a variety of test geometries with abundant thickness measurements including marine-and landterminating glaciers as well as a 2400 km(2) ice cap on Svalbard. The input requirements are kept to a minimum for the first step. In this step, a geometrically controlled, non-local flux solution is converted into thickness values relying on the shallow ice approximation (SIA). In a second step, the thickness field is updated along fast-flowing glacier trunks on the basis of velocity observations. Both steps account for available thickness measurements. Each thickness field is presented together with an error-estimate map based on a formal propagation of input uncertainties. These error estimates point out that the thickness field is least constrained near ice divides or in other stagnant areas. Withholding a share of the thickness measurements, error estimates tend to overestimate mismatch values in a median sense. We also have to accept an aggregate uncertainty of at least 25% in the reconstructed thickness field for glaciers with very sparse or no observations. For Vestfonna ice cap (VIC), a previous ice volume estimate based on the same measurement record as used here has to be corrected upward by 22 %. We also find that a 13% area fraction of the ice cap is in fact grounded below sea level. The former 5% estimate from a direct measurement interpolation exceeds an aggregate maximum range of 6-23% as inferred from the error estimates here.

Abstract: Passive acoustic monitoring of the self-generated noise of particle impacts has been shown to be correlated to bedload flux and bedload size. However, few studies have concentrated on the role of acoustic wave propagation in a river. For the first time, the river environment is modeled as a Pekeris waveguide, where a wave number integration technique is used to predict the transformation of sounds through their propagation paths. Focusing on the distance of a hydrophone from the channel bed and cutting off the low frequencies produced by impacts between gravel particles, we demonstrate that acoustic propagation modifies the spectral content of bedload-generated sound. Acoustic signals analyzed with the proposed model are interpreted by comparison to Helley-Smith bedload data obtained during flood conditions on the large gravel-bedded Arc-en-Maurienne River, France. This study shows that careful attention to acoustic propagation effects is required when estimating bedload grain size distribution with hydrophones in rivers, especially for rivers with slopes higher than 1%. Bedload monitoring with a hydrophone is particularly appropriate for large gravel-bed rivers – especially so during large floods, when in situ sampling is difficult or impractical and the impact of acoustic propagation is weaker relative to the self-generated noise of bedload impacts. Copyright (c) 2017 John Wiley & Sons, Ltd.

Abstract: Supersaturation often occurs at the top of the troposphere where cirrus clouds form, but is comparatively unusual near the surface where the air is generally warmer and laden with liquid and/or ice condensation nuclei. One exception is the surface of the high Antarctic Plateau. One year of atmospheric moisture measurement at the surface of Dome C on the East Antarctic Plateau is presented. The measurements are obtained using commercial hygrometry sensors modified to allow air sampling without affecting the moisture content, even in the case of supersaturation. Supersaturation is found to be very frequent. Common unadapted hygrometry sensors generally fail to report supersaturation, and most reports of atmospheric moisture on the Antarctic Plateau are thus likely biased low. The measurements are compared with results from two models implementing cold microphysics parameterizations: the European Center for Medium-range Weather Forecasts through its operational analyses, and the Model Atmospherique Regional. As in the observations, supersaturation is frequent in the models but the statistical distribution differs both between models and observations and between the two models, leaving much room for model improvement. This is unlikely to strongly affect estimations of surface sublimation because supersaturation is more frequent as temperature is lower, and moisture quantities and thus water fluxes are small anyway. Ignoring supersaturation may be a more serious issue when considering water isotopes, a tracer of phase change and temperature, largely used to reconstruct past climates and environments from ice cores. Because observations are easier in the surface atmosphere, longer and more continuous in situ observation series of atmospheric supersaturation can be obtained than higher in the atmosphere to test parameterizations of cold microphysics, such as those used in the formation of high-altitude cirrus clouds in meteorological and climate models.

Abstract: Aerosol-cloud interaction contributes to the largest uncertainties in the estimation and interpretation of the Earth's changing energy budget. The present study explores experimentally the impacts of water condensation-evaporation events, mimicking processes occurring, in atmospheric clouds, on the molecular composition of secondary organic aerosol (SOA) from the photooxidation of methacrolein. A range of on and off-line mass spectrometry techniques were used to obtain a detailed chemical characterization of SOA formed in control experiments in dry conditions, in triphasic experiments simulating gas-particle-cloud droplet interactions (starting from dry conditions and from 60% relative humidity (RH)), and in bulk aqueous-phase experiments. We observed that cloud events trigger fast SOA formation accompanied by evaporative losses. These evaporative losses decreased SOA concentration in the simulation chamber by 25-32% upon RH increase, while aqueous SOA was found to be metastable and slowly evaporated after cloud dissipation. In the simulation chamber, SOA composition measured with a high-resolution time-of flight aerosol mass spectrometer, did not change during cloud events compared with high RH conditions (RH > 80%). In all experiments, off-line mass spectrometry techniques emphasize the critical role of 2-methylglyceric acid as a major product of isoprene chemistry, as an important contributor to the total SOA mass (15-20%) and as a key building block of oligomers found in the particulate phase. Interestingly, the comparison between the series of oligomers obtained from experiments performed under different conditions show a markedly different reactivity. In particular, long reaction times at high RH seem to create the conditions for aqueous-phase processing to occur in a more efficient manner than during two relatively short cloud events.

Abstract: A 22.4 m-long shallow firn core was extracted during the 2006/2007 field season from coastal Adelie Land. Annual layer counting based on subannual analyses of delta O-18 and major chemical components was combined with 5 reference years associated with nuclear tests and non-retreat of summer sea ice to build the initial ice-core chronology (19462006), stressing uncertain counting for 8 years. We focus here on the resulting delta O-18 and accumulation records. With an average value of 21.8 +/- 6.9 cmw. e. yr(-1), local accumulation shows multi-decadal variations peaking in the 1980s, but no long-term trend. Similar results are obtained for delta O-18, also characterised by a remarkably low and variable amplitude of the seasonal cycle. The ice-core records are compared with regional records of temperature, stake area accumulation measurements and variations in sea-ice extent, and outputs from two models nudged to ERA (European Re-analysis) atmospheric reanalyses: the high-resolution atmospheric general circulation model (AGCM), including stable water isotopes ECHAM5-wiso (European Centre Hamburg model), and the regional atmospheric model Modele Atmospherique Regional (AR). A significant linear correlation is identified between decadal variations in delta O-18 and regional temperature. No significant relationship appears with regional sea-ice extent. A weak and significant correlation appears with Dumont d'Urville wind speed, increasing after 1979. The model-data comparison highlights the inadequacy of ECHAM5-wiso simulations prior to 1979, possibly due to the lack of data assimilation to constrain atmospheric reanalyses. Systematic biases are identified in the ECHAM5-wiso simulation, such as an overestimation of the mean accumulation rate and its interannual variability, a strong cold bias and an underestimation of the mean delta O-18 value and its interannual variability. As a result, relationships between simulated delta O-18 and temperature are weaker than observed. Such systematic precipitation and temperature biases are not displayed by MAR, suggesting that the model resolution plays a key role along the Antarctic ice sheet coastal topography. Interannual variations in ECHAM5-wiso temperature and precipitation accurately capture signals from meteorological data and stake observations and are used to refine the initial ice-core chronology within 2 years. After this adjustment, remarkable positive (negative) delta O-18 anomalies are identified in the ice-core record and the ECHAM5-wiso simulation in 1986 and 20

Abstract: The first results of a campaign of intensive observation of precipitation in Dumont d'Urville, Antarctica, are presented. Several instruments collected data from November 2015 to February 2016 or longer, including a polarimetric radar (MXPol), a Micro Rain Radar (MRR), a weighing gauge (Pluvio(2)), and a Multi-Angle Snowflake Camera (MASC). These instruments collected the first ground-based measurements of precipitation in the region of Adelie Land (Terre Adelie), including precipitation microphysics. Microphysical observations during the austral summer 2015/2016 showed that, close to the ground level, aggregates are the dominant hydrometeor type, together with small ice particles (mostly originating from blowing snow), and that riming is a recurring process. Eleven percent of the measured particles were fully developed graupel, and aggregates had a mean riming degree of about 30 %. Spurious precipitation in the Pluvio2 measurements in windy conditions, leading to phantom accumulations, is observed and partly removed through synergistic use of MRR data. The yearly accumulated precipitation of snow (300 m above ground), obtained by means of a local conversion relation of MRR data, trained on the Pluvio2 measurement of the summer period, is estimated to be 815 mm of water equivalent, with a confidence interval ranging between 739.5 and 989 mm. Data obtained in previous research from satellite-borne radars, and the ERA-Interim reanalysis of the European Centre for Medium-Range Weather Forecasts (ECMWF) provide lower yearly totals: 655 mm for ERA-Interim and 679 mm for the climatological data over DDU. ERA-Interim overestimates the occurrence of low-intensity precipitation events especially in summer, but it compensates for them by underestimating the snowfall amounts carried by the most intense events. Overall, this paper provides insightful examples of the added values of precipitation monitoring in Antarctica with a synergistic use of in situ and remote sensing measurements.

Abstract: Land surface processes play an important role in the West African monsoon variability. In addition, the evolution of hydrological systems in this region, and particularly the increase of surface water and runoff coefficients observed since the 1950s, has had a strong impact on water resources and on the occurrence of floods events. This study addresses results from phase 2 of the African Monsoon Multidisciplinary Analysis (AMMA) Land Surface Model Intercomparison Project (ALMIP2), carried out to evaluate the capability of different state-of-the-art land surface models to reproduce surface processes at the mesoscale. Evaluation of runoff and water fluxes over the Mali site is carried out through comparison with runoff estimations over endorheic watersheds as well as evapotranspiration (ET) measurements. Three remote-sensing-based ET products [ALEXI, MODIS, and Global Land Evaporation Amsterdam Model (GLEAM)] are also analyzed. It is found that, over deep sandy soils, surface runoff is generally overestimated, but the ALMIP2 multimodel mean reproduces in situ measurements of ET and water stress events rather well. However, ALMIP2 models are generally unable to distinguish among the two contrasted hydrological systems typical of the study area. Employing as input a soil map that explicitly represents shallow soils improves the representation of water fluxes for the models that can account for their representation. Shallow soils are shown to be also quite challenging for remote-sensing-based ET products, even if their effect on evaporative loss was captured by the diagnostic thermal-based ALEXI. A better representation of these soils, in soil databases, model parameterizations, and remote sensing algorithms, is fundamental to improve the estimation of water fluxes in this part of the Sahel.

Abstract: Aquatic ecosystems of the Bolivian Altiplano (similar to 3800 m a.s.l.) are characterized by extreme hydro-climatic constrains (e.g., high UV-radiations and low oxygen) and are under the pressure of increasing anthropogenic activities, unregulated mining, agricultural and urban development. We report here a complete inventory of mercury (Hg) levels and speciation in the water column, atmosphere, sediment and key sentinel organisms (i.e., plankton, fish and birds) of two endorheic Lakes of the same watershed differing with respect to their size, eutrophication and contamination levels. Total Hg (THg) and mono-methylmercury (MMHg) concentrations in filtered water and sediment of Lake Titicaca are in the lowest range of reported levels in other large lakes worldwide. Downstream, Hg levels are 3-10 times higher in the shallow eutrophic Lake Uru-Uru than in Lake Titicaca due to high Hg inputs from the surrounding mining region. High percentages of MMHg were found in the filtered and unfiltered water rising up from <1 to similar to 50% THg from the oligoihetero-trophic Lake Titicaca to the eutrophic Lake Uru-Uru. Such high % MMHg is explained by a high in situ MMHg production in relation to the sulfate rich substrate, the low oxygen levels of the water column, and the stabilization of MMHg due to abundant ligands present in these alkaline waters. Differences in MMHg concentrations in water and sediments compartments between Lake Titicaca and Uru-Uru were found to mirror the offset in MMHg levels that also exist in their respective food webs. This suggests that in situ MMHg baseline production is likely the main factor controlling MMHg levels in fish species consumed by the local population. Finally, the increase of anthropogenic pressure in Lake Titicaca may probably enhance eutrophication processes which favor MMHg production and thus accumulation in water and biota. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Analog methods are based on a statistical relationship between synoptic meteorological variables (predictors) and local weather (predictand, to be predicted). This relationship is defined by several parameters, which are often calibrated by means of a semiautomatic sequential procedure. This calibration approach is fast, but has strong limitations. It proceeds through successive steps, and thus cannot handle all parameter dependencies. Furthermore, it cannot automatically optimize some parameters, such as the selection of pressure levels and temporal windows (hours of the day) at which the predictors are compared. To overcome these limitations, the global optimization technique of genetic algorithms is considered, which can jointly optimize all parameters of the method, and get closer to a global optimum, by taking into account the dependencies of the parameters. Moreover, it can objectively calibrate parameters that were previously assessed manually and can take into account new degrees of freedom. However, genetic algorithms must be tailored to the problem under consideration. Multiple combinations of algorithms were assessed, and new algorithms were developed (e.g., the chromosome of adaptive search radius, which is found to be very robust), in order to provide recommendations regarding the use of genetic algorithms for optimizing several variants of analog methods. A global optimization approach provides new perspectives for the improvement of analog methods, and for their application to new regions or new predictands.

Abstract: In order to systematically plan river restoration actions at a regional scale, this paper develops a multi-criteria analysis that classifies rivers, based on their priority for hydromorphological restoration. This priority is defined by severe human pressures within the erodible corridor of the river, drastic alteration of the stream channel, and low intensity of river pattern functioning. Based on relevant indicators for three groups of features (human pressures, channel changes, and river functionality), a Hydromorphological Restoration Priority Index (HRPI) was designed. The high values (>66%) of HRPI reflect an urgent need for hydromorphological restoration while low values (<33%) reveal a less immediate necessity for restoration. The proposed methodology was applied on braided sectors of rivers crossing the south-eastern (Curvature) Subcarpathians (Romania). The values of the total HRPI ranged between 21% (Zabrdut River) and almost 44% (Prahova River). According to our results, most of the analyzed sectors have a low need for hydromorphological restoration of the braided pattern, while some have a moderate necessity for restoration. Whereas the Prahova River has the highest HRPI, it should be given priority for restoration at a regional scale, which corresponds to the objectives of River Basin Management Plans for the interval beyond 2021. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Similar to many other regions in the world, glaciers in the southern sub-polar regions are currently retreating. In the Kerguelen Islands (49 degrees S, 69 degrees E), the mass balance of the Cook Ice Cap (CIC), the largest ice cap in this region, experienced dramatic shrinking between 1960 and 2013 with retreat rates among the highest in the world. This observation needs to be evaluated in a long-term context. However, data on the past glacier extents are sparse in the sub-Antarctic regions. To investigate the deglaciation pattern since the Last Glacial Maximum (LGM) period, we present the first 13 cosmogenic Cl-36 surface exposure ages from four sites in the Kerguelen Islands. The Cl-36 ages from erratic and moraine boulders span from 24.4 +/- 2.7 ka to 03 +/- 0.1 ka. We combined these ages with existing glacio-marine radiocarbon ages and bathymetric data to document the temporal and spatial changes of the island's glacial history. Ice began to retreat on the main island before 24.4 +/- 2.7 ka until around the time of the Antarctic Cold Reversal (ACR) period (-14.5-12.9 ka), during which the Bontemps moraine was formed by the advance of a CIC outlet glacier. Deglaciation continued during the Holocene probably until 3 ka with evidence of minor advances during the last millennium. This chronology is in pace with major changes in delta O-18 in a recent West Antarctica ice core record, showing that Kerguelen Islands glaciers are particularly sensitive and relevant to document climate change in the southern polar regions. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: A 1/128 ocean model configuration of the Amundsen Sea sector is developed to better understand the circulation induced by ice-shelf melt and the impacts on the surrounding ocean and sea ice. Eighteen sensitivity experiments to drag and heat exchange coefficients at the ice shelf/ocean interface are performed. The total melt rate simulated in each cavity is function of the thermal Stanton number, and for a given thermal Stanton number, melt is slightly higher for lower values of the drag coefficient. Sub-ice-shelf melt induces a thermohaline circulation that pumps warm circumpolar deep water into the cavity. The related volume flux into a cavity is 100-500 times stronger than the melt volume flux itself. Ice-shelf melt also induces a coastal barotropic current that contributes 45612% of the total simulated coastal transport. Due to the presence of warm circumpolar deep waters, the melt-induced inflow typically brings 4-20 times more heat into the cavities than the latent heat required for melt. For currently observed melt rates, approximately 6-31% of the heat that enters a cavity with melting potential is actually used to melt ice shelves. For increasing sub-ice-shelf melt rates, the transport in the cavity becomes stronger, and more heat is pumped from the deep layers to the upper part of the cavity then advected toward the ocean surface in front of the ice shelf. Therefore, more ice-shelf melt induces less sea-ice volume near the ice sheet margins. Plain Language Summary The ice-shelf cavities of the Amundsen Sea, Antarctica, act as very powerful pumps that create strong inflows of warm water under the ice-shelves, as well as significant circulation changes in the entire region. Such warm inflows bring more heat than required to melt ice, so that a large part of that heat exits ice-shelf cavities without being used. Due to mixing between warm deep waters and melt freshwater, melt-induced flows are warm and buoyant when they leave cavities. Therefore, they reach the ocean surface near ice-shelf fronts and can melt significant amounts of sea ice. It is thus suggested that climatic trends in sub ice-shelf melt could partly explain sea ice trends near the ice-sheet margins in the Amundsen Sea region.

Abstract: Polar ice cores are unique climate archives. Indeed, most of them have a continuous stratigraphy and present high temporal resolution of many climate variables in a single archive. While water isotopic records (delta D or delta O-18) in ice cores are often taken as references for past atmospheric temperature variations, their relationship to temperature is associated with a large uncertainty. Several reasons are invoked to explain the limitation of such an approach; in particular, post-deposition effects are important in East Antarctica because of the low accumulation rates. The strong influence of post-deposition processes highlights the need for surface polar research programs in addition to deep drilling programs. We present here new results on water isotopes from several recent surface programs, mostly over East Antarctica. Together with previously published data, the new data presented in this study have several implications for the climatic reconstructions based on ice core isotopic data: (1) The spatial relationship between surface mean temperature and mean snow isotopic composition over the first meters in depth can be explained quite straightforwardly using simple isotopic models tuned to d-excess vs. delta O-18 evolution in transects on the East Antarctic sector. The observed spatial slopes are significantly higher (similar to 0.7-0.8 parts per thousand. degrees C (1) for delta O-18 vs. temperature) than seasonal slopes inferred from precipitation data at Vostok and Dome C (0.35 to 0.46 parts per thousand. degrees C (1)). We explain these differences by changes in condensation versus surface temperature between summer and winter in the central East Antarctic plateau, where the inversion layer vanishes in summer. (2) Post-deposition effects linked to exchanges between the snow surface and the atmospheric water vapor lead to an evolution of delta O-18 in the surface snow, even in the absence of any precipitation event. This evolution preserves the positive correlation between the delta O-18 of snow and surface temperature, but is associated with a much slower delta O-18-vs-temperature slope than the slope observed in the seasonal precipitation. (3) Post-deposition effects clearly limit the archiving of high-resolution (seasonal) climatic variability in the polar snow, but we suggest that sites with an accumulation rate of the order of 40 kg.m (2).yr (1) may record a seasonal cycle at shallow depths. (C) 2017 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Abstract: L-band radiometer measurements collected over the Dome similar to C area from 2010 to 2015 indicated that the brightness temperature (T-B) was relatively stable at vertical (V) polarization (standard deviation lower than 1 K at annual scale), while it was slightly more variable at horizontal (H) polarization. During the 2014-2015 austral summer, an exceptional situation was recorded by both the DOMEX ground radiometer and the European Space Agency (ESA)'s Soil Moisture and Ocean Salinity (SMOS) satellite. From November 2014 to March 2015, T-B H showed a progressive and significant increase until 20 March 2015 when it sharply decreased by about 5 K (at 52.5 incidence angle) within a few days. In parallel to the increase in T-B H, glaciological and meteorological in situ measurements showed a wind speed that was lower than usual and a low-density snow layer being progressively set up on the surface. This was consistent with the exceptional hoar event observed, as well as with snow accumulation on the surface. On the other hand, the decrease in T-B H was related to the passing over Dome C of a storm that removed or compacted the layer of light snow on the surface. The WALOMIS (Wave Approach for LOw-frequency Microwave emission in Snow) snow-emission model was used with in situ measurements of the snowpack as inputs for evaluating the effect of changes observed on the snow surface in T-B H. The simulations indicated that the surface snow density variations were sufficient for predicting the increasing and decreasing trends of the T-B H. However, the thickness variations of the superficial layer were essential so as to obtain a better agreement with the SMOS observations. This result confirmed that the L-band T-B H was affected by the snow properties of the top centimeters of the snowpack, in spite of the large penetration depth (hundreds of meters). Both the surface snow density and the thickness of the superficial layer were relevant, due to coherent interference effects.

Abstract: Multiple year-round (2006-2015) records of the bulk and size-segregated composition of aerosol were obtained at the inland site of Concordia located in East Antarctica. The well-marked maximum of non-sea salt sulfate (nssSO(4)) in January (100 +/- 28 ng m(-3) versus 4.4 +/- 2.3 ng m(-3) in July) is consistent with observations made at the coast (280 +/- 78 ng M-3 in January versus 16 +/- 9 ng m(-3) in July at Dumont d'Urville, for instance). In contrast, the well-marked maximum of MSA at the coast in January (60 +/- 23 ng M-3 at Dumont d'Urville) is not observed at Concordia (5.2 +/- 2.0 ng M-3 in January). Instead, the MSA level at Concordia peaks in October (5.6 +/- 1.9 ng m(-3)) and March (14.9 +/- 5.7 ng m(-3)). As a result, a surprisingly low MSA-to-nssSO(4) ratio (RMSA) is observed at Concordia in mid-summer (0.0 +/- 5 0.02 in January versus 0.25 +/- 0.09 in March). We find that the low value of RMSA in mid-summer at Concordia is mainly driven by a drop of MSA levels that takes place in submicron aerosol (0.3 pm diameter). The drop of MSA coincides with periods of high photochemical activity as indicated by high ozone levels, strongly suggesting the occurrence of an efficient chemical destruction of MSA over the Antarctic plateau in mid-summer. The relationship between MSA and nssSO(4) levels is examined separately for each season and indicates that concentration of non-biogenic sulfate over the Antarctic plateau does not exceed 1 ng M-3 in fall and winter and remains close to 5 ng M-3 in spring. This weak non-biogenic sulfate level is discussed in the light of radionuclides (Pb-210, Be-10, and Be-7) also measured on bulk aerosol samples collected at Concordia. The findings highlight the complexity in using MSA in deep ice cores extracted from inland Antarctica as a proxy of past dimethyl sulfide emissions from the Southern Ocean.

Abstract: The Maritime Continent is a region of intense rainfall characterised by a strong diurnal cycle. This study investigates the sensitivity of rainfall characteristics to resolution in a coupled regional climate model configuration, in particular focusing on processes that modulate the diurnal cycle. Model biases are resolution dependent. Increasing resolution from 3/4A degrees to 1/4A degrees improves the mean state sea surface temperature and precipitation biases. However, at higher resolutions (1/12A degrees) rainfall becomes too strong in most areas. Daily maximum rainfall is simulated about 2-4 h earlier than in observations over both the land and the ocean, with only small improvements over high topography at higher resolution. We develop a technique to examine cross-coastal processes associated with the rainfall diurnal cycle along all coastlines. This is used to investigate the sensitivity of the diurnal cycle to resolution and to the direction of the prevailing wind. During offshore prevailing winds, most land rainfall is confined near the coastline and associated with a shallow land-sea breeze circulation at all resolution (though rainfall partly develops directly inland at 1/12A degrees). During onshore prevailing winds, rainfall propagates from the coastline to the island interior at 1/4A degrees and 1/12A degrees, whereas it appears directly over the island interior at 3/4A degrees, and this is associated with a deep convective cell across the coastline for all resolutions. Oceanic rainfall propagates far offshore during offshore prevailing winds at all resolutions, whereas it tends to remain confined near the coastline under onshore prevailing winds condition, particularly at higher resolution.

Abstract: Over Saint-Sorlin Glacier in the French Alps (45 degrees N, 6.1 degrees E; similar to 3 km(2)) in summer, we study the atmospheric surface-layer dynamics, turbulent fluxes, their uncertainties and their impact on surface energy balance (SEB) melt estimates. Results are classified with regard to largescale forcing. We use high-frequency eddy-covariance data and mean air-temperature and wind-speed vertical profiles, collected in 2006 and 2009 in the glacier's atmospheric surface layer. We evaluate the turbulent fluxes with the eddycovariance (sonic) and the profile method, and random errors and parametric uncertainties are evaluated by including different stability corrections and assuming different values for surface roughness lengths. For weak synoptic forcing, local thermal effects dominate the wind circulation. On the glacier, weak katabatic flows with a wind-speed maximum at low height (2-3 m) are detected 71% of the time and are generally associated with small turbulent kinetic energy (TKE) and small net turbulent fluxes. Radiative fluxes dominate the SEB. When the large-scale forcing is strong, the wind in the valley aligns with the glacier flow, intense downslope flows are observed, no wind-speed maximum is visible below 5 m, and TKE and net turbulent fluxes are often intense. The net turbulent fluxes contribute significantly to the SEB. The surface-layer turbulence production is probably not at equilibrium with dissipation because of interactions of largescale orographic disturbances with the flow when the forcing is strong or low-frequency oscillations of the katabatic flow when the forcing is weak. In weak forcing when TKE is low, all turbulent fluxes calculation methods provide similar fluxes. In strong forcing when TKE is large, the choice of roughness lengths impacts strongly the net turbulent fluxes from the profile method fluxes and their uncertainties. However, the uncertainty on the total SEB remains too high with regard to the net observed melt to be able to recommend one turbulent flux calculation method over another.

Abstract: In this study, we investigate changes in the glaciated surface and the formation of lakes in the headwater of the Querococha watershed in Cordillera Blanca (Peru) using 24 Landsat images from 1975 to 2014. Information of glacier retreat was integrated with available climate data, the first survey of recent depositional dynamics in proglacial Yanamarey Lake (4600 m a.s.l.), and a relatively short hydrological record (2002-2014) at the outlet of Yanamarey Lake. A statistically significant temperature warming (0.21 degrees C decade(-1) for mean annual temperature) has been detected in the region, and it caused a reduction of the glacierized area since 1975 from 3.5 to 1.4 km(-2). New small lakes formed in the deglaciated areas, increasing the flooded area from 1.8 ha in 1976 to 2.8 ha in 2014. A positive correlation between annual rates of glacier recession and runoff was found. Sediment cores revealed a high sedimentation rate (>1 cm yr(-1)) and two contrasted facies, suggesting a shift toward a reduction of meltwater inputs and higher hydrological variability likely due to an increasing role of precipitation on runoff during the last decades. Despite the age control uncertainties, the main transition likely occurred around 19982000, correlating with the end of the phase with maximum warming rates and glacier retreat during the 1980s and 1990s, and the slowing down of expansion of surface lake-covered surface. With this hydrological – paleolimnological approach we have documented the association between recent climate variability and glacier recession and the rapid transfer of hydroclimate signal to depositional and geochemical processes in high elevation Andean environments. This, study also alerts about water quality risks as proglacial lakes act as secondary reservoirs that trap trace and minor elements in high altitude basins. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: The nature and timing of environmental changes throughout the last glacial-interglacial transition in the South Caucasus, and more widely in eastern Europe, are still not fully understood. According to certain pollen records, forest expansion occurred in many areas several millennia after what is considered worldwide as the onset of the Holocene. The current problem we face is that the time lag in forest expansion varies from one sequence to another, sometimes with no delay at all. Moreover, the potential forcing/controlling factors behind this complex pattern, contrary to the almost synchronous global Holocene warming, are still a matter for debate. Accordingly, we revisit the issue of forest expansion through vegetation history obtained in the South Caucasus using a new pollen record, retrieved from the Nariani paleolake (South Georgia). These data attest to a steppic phase, initially dominated by Amaranthaceae-Chenopodiaceae (12,700-10,500 cal yr BP), then by Poaceae (10,500-9000 cal yr BP), culminating with a more forested phase (9000-5000 cal yr BP). Although some palaeoclimatic regional reconstructions show a wet early Holocene, we interpret the delay in forest expansion recorded in Nariani (2500 years) as the result of reduced spring precipitation, which would have limited forest development at that time.

Abstract: Climate change is considered as one of the main factors controlling sediment fluxes in mountain belts. However, the effect of El Nino, which represents the primary cause of inter-annual climate variability in the South Pacific, on river erosion and sediment transport in the Western Andes remains unclear. Using an unpublished dataset of Suspended Sediment Yield (SSY) in Peru (1968-2012), we show that the annual SSY increases by 3-60 times during Extreme El Nino Events (EENE) compared to normal years. During EENE, 82% to 97% of the annual SSY occurs from January to April. We explain this effect by a sharp increase in river water discharge due to high precipitation rates and transport capacity during EENE. Indeed, sediments accumulate in the mountain and piedmont areas during dry normal years, and are then rapidly mobilized during EENE years. The effect of EENE on SSY depends on the topography, as it is maximum for catchments located in the North of Peru (3-7 degrees S), exhibiting a concave up hypsometric curve, and minimum for catchments in the South (7-18 degrees S), with a concave down hypsometric curve. These findings highlight how the sediment transport of different topographies can respond in very different ways to large climate variability.

Abstract: Tephra isochrons offer considerable potential for correlating diverse palaeoarchives and highlighting regional climatic differences. They are especially useful when applied to polar ice records encompassing the last glacial, as these clearly portray the pronounced millennial-scale climate variability that characterised this period. Here we present the continuous record of primary fallout tephra layers in the East Antarctic Talos Dome ice core (72 degrees 49 ' S, 159 degrees 11 ' E), developed upon examination of the core sections spanning the glacial period 16.5 to 71 ka. A total of ca. 45 discrete tephra deposits precisely positioned stratigraphically relative to the temperature record for the core and dated using the AICC2012 timescale, were identified. Quantitative grain size, particle morphology, major and trace element composition using Coulter Counter, SEM, EPMA-WDS, and LA-ICP-MS analytical methods were studied as diagnostic features for tephra characterisation. The tephrostratigraphic framework provides a reference for future precise comparison between ice and sediment sequences across the Antarctic continent. Indeed, several potential markers characterised by distinct volcanic glass geochemistry and/or particular stratigraphic location (e.g., a 17.6-ka ash layer deposited during the well-known major acidity event) are now available for the direct linkage of palaeoclimatic archives. The Talos Dome tephra sequence, dominated by mid distal pyroclastic products from the nearby Northern Victoria Land volcanoes, also represents the most comprehensive and best time-constrained record of regional Antarctic volcanism yet developed. It documents nearly continuous sustained explosive activity during the considered time interval and, combined with previous ice-core tephra results for the last and the current interglacial periods, suggests progressive compositional shift through time. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: The NOx (NO and NO2) and HOx (OH and HO2) budgets of the atmosphere exert a major influence on atmospheric composition, controlling removal of primary pollutants and formation of a wide range of secondary products, including ozone, that can influence human health and climate. However, there remain large uncertainties in the changes to these budgets over recent decades. Due to their short atmospheric lifetimes, NOx and HOx are highly variable in space and time, and so the measurements of these species are of limited value for examining long-term, largescale changes to their budgets. Here, we take an alternative approach by examining long-term atmospheric trends of alkyl nitrates, the production efficiency of which is dependent on the atmospheric [NO] / [HO2] ratio. We derive long-term trends in the alkyl nitrates from measurements in firn air from the NEEM site, Greenland. Their mixing ratios increased by a factor of 3-5 between the 1970s and 1990s. This was followed by a steep decline to the sampling date of 2008. Moreover, we examine how the trends in the alkyl nitrates compare to similarly derived trends in their parent alkanes (i.e. the alkanes which, when oxidised in the presence of NOx, lead to the formation of the alkyl nitrates). The ratios of the alkyl nitrates to their parent alkanes increased from around 1970 to the late 1990s. This is consistent with large changes to the [NO] / [HO2] ratio in the Northern Hemisphere atmosphere during this period. Alternatively, they could represent changes to concentrations of the hydroxyl radical, OH, or to the transport time of the air masses from source regions to the Arctic.

Abstract: A comprehensive hydrometeorological dataset is presented spanning the period 1 January 201131 December 2014 to improve the understanding of the hydrological processes leading to flash floods and the relation between rainfall, runoff, erosion and sediment transport in a mesoscale catchment (Auzon, 116 km(2)) of the Mediterranean region. Badlands are present in the Auzon catchment and well connected to high-gradient channels of bedrock rivers which promotes the transfer of suspended solids downstream. The number of observed variables, the various sensors involved (both in situ and remote) and the space-time resolution (similar to km(2), similar to min) of this comprehensive dataset make it a unique contribution to research communities focused on hydrometeorology, surface hydrology and erosion. Given that rainfall is highly variable in space and time in this region, the observation system enables assessment of the hydrological response to rainfall fields. Indeed, (i) rainfall data are provided by rain gauges (both a research network of 21 rain gauges with a 5 min time step and an operational network of 10 rain gauges with a 5 min or 1 h time step), S-band Doppler dual-polarization radars (1 km(2), 5 min resolution), disdrometers (16 sensors working at 30 s or 1 min time step) and Micro Rain Radars (5 sensors, 100m height resolution). Additionally, during the special observation period (SOP-1) of the HyMeX (Hydrological Cycle in the Mediterranean Experiment) project, two X-band radars provided precipitation measurements at very fine spatial and temporal scales (1 ha, 5 min). (ii) Other meteorological data are taken from the operational surface weather observation stations of Meteo-France (including 2m air temperature, atmospheric pressure, 2 m relative humidity, 10m wind speed and direction, global radiation) at the hourly time resolution (six stations in the region of interest). (iii) The monitoring of surface hydrology and suspended sediment is multi-scale and based on nested catchments. Three hydrometric stations estimate water discharge at a 2-10 min time resolution. Two of these stations also measure additional physico-chemical variables (turbidity, temperature, conductivity) and water samples are collected automatically during floods, allowing further geochemical characterization of water and suspended solids. Two experimental plots monitor overland flow and erosion at 1 min time resolution on a hillslope with vineyard. A network of 11 sensors installed in the intermittent hydrographic network continuously me
asures water level and water temperature in headwater subcatchments (from 0.17 to 116 km(2)) at a time resolution of 2-5 min. A network of soil moisture sensors enables the continuous measurement of soil volumetric water content at 20 min time resolution at 9 sites. Additionally, concomitant observations (soil moisture measurements and stream gauging) were performed during floods between 2012 and 2014. Finally, this dataset is considered appropriate for understanding the rainfall variability in time and space at fine scales, improving areal rainfall estimations and progressing in distributed hydrological and erosion modelling.

Abstract: Great effort is made to address heat waves (HWs) in developed countries because of their devastating impacts on society, economy, and environment. However, HWs are still understudied over developing countries. This is particularly true in West Africa, and especially in the Sahel, where temperatures recurrently reach critical values, such as during the 2010 HW event in the western Sahel. This work aims at characterizing the Sahelian HWs during boreal spring seasons (April-May-June) and understanding the mechanisms associated with such extreme events. Over the last three decades, Sahelian HWs have been becoming more frequent, lasting longer, covering larger areas, and reaching higher intensities. The physical mechanisms associated with HWs are examined to assess the respective roles of atmospheric dynamics and radiative and turbulent fluxes by analyzing the surface energy budget. Results suggest that the greenhouse effect of water vapor is the main driver of HWs in the western Sahel, increasing minimum temperatures by enhanced downward longwave radiation. Atmospheric circulation plays an important role in sustaining these warm anomalies by advecting moisture from the Atlantic Ocean and the Guinean coasts into the Sahel. Maximum temperature anomalies are mostly explained by increased downward shortwave radiation due to a reduction in cloud cover. Interannual variability of HWs is affected by the delayed impact of El Nio-Southern Oscillation (ENSO), with anomalous temperature warming following warm ENSO events, resulting from an amplified water vapor feedback.

Abstract: On average, the models in the Fifth Climate Model Intercomparison Project archive predict an increase in Antarctic precipitation from 5.5 to 24.5 % between 1986-2005 and 2080-2099, depending on greenhouse gas emissions scenarios. This translates into a moderation of future sea level rise ranging from -19 to -71 mm between 2006 and 2099. However, comparison with CloudSat and ERA-Interim data show that almost all the models overestimate current Antarctic precipitation, some by more than 100 %. If only the models that agree with CloudSat data within 20 % of error are considered, larger precipitation changes (from 7.4 to 29.3 %) and impact on sea level (from -25 to -85 mm) are predicted. A common practice of averaging all models to evaluate climate projections thus leads to a significant underestimation of the contribution of Antarctic precipitation to future sea level. Models simulate, on average, a 7.4 %/A degrees C precipitation change with surface temperature warming. The models in better agreement with CloudSat observations for Antarctic snowfall predict, on average, larger temperature and Antarctic sea ice cover changes, which could explain the larger changes in Antarctic precipitation simulated by these models. The agreement between the models, CloudSat data and ERA-Interim is generally less in the interior of Antarctica than at the peripheries, but the interior is also where climate change will induce the smallest absolute change in precipitation. About three-quarters of the impact on sea level will result from precipitation change over the half most peripheral and lowest elevation part of the surface of Antarctica.

Abstract: Polar ice cores provide exceptional archives of past environmental conditions. The dating of ice cores and the estimation of the age-scale uncertainty are essential to interpret the climate and environmental records that they contain. It is, however, a complex problem which involves different methods. Here, we present IceChrono1, a new probabilistic model integrating various sources of chronological information to produce a common and optimized chronology for several ice cores, as well as its uncertainty. IceChrono1 is based on the inference of three quantities: the surface accumulation rate, the lock-in depth (LID) of air bubbles and the thinning function. The chronological information integrated into IceChrono1 are modeling scenarios of the sedimentation process (accumulation of snow, densification of snow into ice and air trapping, ice flow), ice- and air-dated horizons, ice and air depth intervals with known durations, Delta depth observations (depth shift between synchronous events recorded in the ice and in the air) and finally ice, air or mix stratigraphic links in between ice cores. the inference problem is formulated as a least squares optimisation, implying that all densities of probabilities are assumed to be Gaussian. It is numerically solved using the Levenberg-Marquardt algorithm (thus assuming that the model is almost linear in the vicinity of the solution) and a numerical evaluation of the model's Jacobian. IceChrono1 is freely available under the General Public License v3 open source license.

Abstract: Basal melt rate is the most important physical quantity to be evaluated when looking for an old-ice drilling site, and it depends to a great extent on the geothermal flux (GF), which is poorly known under the East Antarctic ice sheet. Given that wet bedrock has higher reflectivity than dry bedrock, the wetness of the ice-bed interface can be assessed using radar echoes from the bedrock. But, since basal conditions depend on heat transfer forced by climate but lagged by the thick ice, the basal ice may currently be frozen whereas in the past it was generally melting. For that reason, the risk of bias between present and past conditions has to be evaluated. The objective of this study is to assess which locations in the Dome C area could have been protected from basal melting at any time in the past, which requires evaluating GF. We used an inverse approach to retrieve GF from radar-inferred distribution of wet and dry beds. A 1-D heat model is run over the last 800 ka to constrain the value of GF by assessing a critical ice thickness, i.e. the minimum ice thickness that would allow the present local distribution of basal melting. A regional map of the GF was then inferred over a 80 km x 130 km area, with a N-S gradient and with values ranging from 48 to 60m Wm(-2). The forward model was then emulated by a polynomial function to compute a time-averaged value of the spatially variable basal melt rate over the region. Three main subregions appear to be free of basal melting, two because of a thin overlying ice and one, north of Dome C, because of a low GF.

Abstract: The convection or settling of matter in the deep Earth's interior is mostly constrained by density variations between the different reservoirs. Knowledge of the density contrast between solid and molten silicates is thus of prime importance to understand and model the dynamic behavior of the past and present Earth. SiO2 is the main constituent of Earth's mantle and is the reference model system for the behavior of silicate melts at high pressure. Here, we apply our recently developed x-ray absorption technique to the density of SiO2 glass up to 110 GPa, doubling the pressure range for such measurements. Our density data validate recent molecular dynamics simulations and are in good agreement with previous experimental studies conducted at lower pressure. Silica glass rapidly densifies up to 40 GPa, but the density trend then flattens to become asymptotic to the density of SiO2 minerals above 60 GPa. The density data present two discontinuities at similar to 17 and similar to 60 GPa that can be related to a silicon coordination increase from 4 to a mixed 5/6 coordination and from 5/6 to sixfold, respectively. SiO2 glass becomes denser than MgSiO3 glass at similar to 40 GPa, and its density becomes identical to that of MgSiO3 glass above 80 GPa. Our results on SiO2 glass may suggest that a variation of SiO2 content in a basaltic or pyrolitic melt with pressure has at most a minor effect on the final melt density, and iron partitioning between the melts and residual solids is the predominant factor that controls melt buoyancy in the lowermost mantle.

Abstract: The strong socio-economic growth of Quito has led to significant projects of inter-basin transfers, intensifying high-altitude resources mobilization in environmentally sensitive areas and with accelerated melting of glaciers. In order to study various future scenarios, we propose modelling of the climate/glacier/hydrology/water resources management continuum, applied to Andean basins. Using the tool Water Evaluation and Planning (WEAP21), we developed: (1) semi-distributed hydro-climate modelling with monthly data using regional homogenization (vector method); (2) glacier modelling of water production and inter-annual evolution; (3) modelling in hydrological units depending on land cover; and (4) water management modelling distinguishing water rights, allocation and use. In this paper, we mainly present results of monthly hydrological calibrations (1963-2006), studying, in particular, the equifinality of various parameter settings. We show flexibility, robustness and limits of the proposed modelling, contributing to identification of uncertainties for evaluation of future scenarios.

Abstract: Glaciers are one of the terrestrial essential climate variables (ECVs) as they respond very sensitively to climate change. A key driver of their response is the glacier surface mass balance that is typically derived from field measurements. It deserves to be quantified over long time scales to better understand the accumulation and ablation processes at the glacier surface and their relationships with inter-annual changes in meteorological conditions and long-term climate changes. Glaciers with in situ monitoring of surface mass balance are scarce at the global scale, and satellite remote sensing provides a powerful tool to increase the number of monitored glaciers. In this study, we present a review of three optical remote sensing methods developed to quantify seasonal and annual glacier surface mass balances. These methodologies rely on the multitemporal monitoring of the end-of-summer snow line for the equilibrium-line altitude (ELA) method, the annual cycle of glacier surface albedo for the albedo method and the mapping of the regional snow cover at the seasonal scale for the snow-map method. Together with a presentation of each method, an application is illustrated. The ELA method shows promising results to quantify annual surface mass balance and to reconstruct multi-decadal time series. The other two methods currently need a calibration on the basis of existing in situ data; however, a generalization of these methods (without calibration) could be achieved. The two latter methods show satisfying results at the annual and seasonal scales, particularly for the summer surface mass balance in the case of the albedo method and for the winter surface mass balance in the case of the snow-map method. The limits of each method (e.g., cloud coverage, debris-covered glaciers, monsoon-regime and cold glaciers), their complementarities and the future challenges (e.g., automating of the satellite images processing, generalization of the methods needing calibration) are also discussed.

Abstract: This study compares the multivariate predictions of daily temperature, temperature range, precipitation, surface wind and solar radiation of a single-model analogue approach with an original multi-model analogy over 12 regions in Europe and Maghreb. Both approaches are based on two-level analogue models where atmospheric predictors are either dynamic or thermodynamic. In the multi-model approach, independent analogue models with predictand-specific predictors are used. In the single-model one, a unique analogue model and its associated set of predictors is applied to all predictands. Testing numerous large-scale predictors, we first identify the best predictor sets for each modelling strategy. Those obtained for the single-model approach are significantly different from those of the predictand-specific models. This is especially the case for local temperature and wind speed. Both methods perform similarly for precipitation, temperature range and radiation. We next assess the ability of both approaches to simulate physically coherent multivariate weather scenarios. With the single-model method, weather scenarios are obtained for each prediction day from observations sampled simultaneously on one analogue day. The physical consistency between variables is thus automatically fulfilled each day. This allows the single-model method to reproduce well the observed inter-predictand correlations, especially the significant correlations between radiation and precipitation and between radiation and temperature range. These results suggest a hybrid analogue model using a single-model for radiation, temperature range and precipitation, combined with a univariate approach for wind. Two options are proposed for temperature for which either the predictand-specific method or a single-model approach with an additional correction are conceivable. This hybrid approach leads to a possible compromise between reasonable univariate prediction skills and realistic inter-predictands correlations, both classically required for many impact studies.

Abstract: The Soil Moisture Active Passive (SMAP) mission Level-4 Surface and Root-Zone Soil Moisture (L4SM) data product is generated by assimilating SMAP L-band brightness temperature observations into the NASA Catchment land surface model. The L4SM product is available from 31 March 2015 to present (within 3 days from real time) and provides 3-hourly, global, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and land surface conditions. This study presents an overview of the L4SM algorithm, validation approach, and product assessment versus in situ measurements. Core validation sites provide spatially averaged surface (root zone) soil moisture measurements for 43 (17) "reference pixels'' at 9- and 36-km gridcell scales located in 17 (7) distinct watersheds. Sparse networks provide point-scale measurements of surface (root zone) soil moisture at 406 (311) locations. Core validation site results indicate that the L4SM product meets its soil moisture accuracy requirement, specified as an unbiased RMSE (ubRMSE, or standard deviation of the error) of 0.04 m(3) m(-3) or better. The ubRMSE for L4SM surface (root zone) soil moisture is 0.038 m(3) m(-3) (0.030 m(3) m(-3)) at the 9-km scale and 0.035 m(3) m(-3) (0.026 m(3) m(-3)) at the 36-km scale. The L4SM estimates improve (significantly at the 5% level for surface soil moisture) over model-only estimates, which do not benefit from the assimilation of SMAP brightness temperature observations and have a 9-km surface (root zone) ubRMSE of 0.042 m(3) m(-3) (0.032 m(3) m(-3)). Time series correlations exhibit similar relative performance. The sparse network results corroborate these findings over a greater variety of climate and land cover conditions.

Abstract: In order to understand atmospheric methane (CH4) biogeochemistry now and in the future, we must apprehend its natural variability, without anthropogenic influence. Samples of ancient air trapped within ice cores provide the means to do this. Here we analyze the ultrahigh-resolution CH4 record of the West Antarctic Ice Sheet Divide ice core 67.2-9.8 ka and find novel, atmospheric CH4 variability at centennial time scales throughout the record. This signal is characterized by recurrence intervals within a broad 80-500 year range, but we find that age-scale uncertainties complicate the possible isolation of any periodic frequency. Lower signal amplitudes in the Last Glacial relative to the Holocene may be related to incongruent effects of firn-based signal smoothing processes. Within interstadial and stadial periods, the peak-to-peak signal amplitudes vary in proportion to the underlying millennial-scale oscillations in CH4 concentration-the relative amplitude change is constant. We propose that the centennial CH4 signal is related to tropical climate variability that influences predominantly low-latitude wetland CH4 emissions. Plain Language Summary Using a new method to measure methane concentrations of ancient air trapped in ice cores, we have detected variability in atmospheric methane concentration on centennial time scales in the Last Glacial Period for the first time. We know these signals represent past changes in atmospheric methane because they appear in several ice core records. We ropose that changes in methane emissions from tropical wetlands are responsible. How this new variability might be related to similar signals found in the late Holocene ice core records and the instrumental record of atmospheric methane is an open question.

Abstract: Global rainfall information is useful for many applications. However, real-time versions of satellite-based rainfall products are known to contain errors. Recent studies have demonstrated how the information about rainfall intrinsically contained in soil moisture data can be utilised for improving rainfall estimates. That is, soil moisture dynamics are impacted for several days by the accumulated amount of rainfall following within a particular event. In this context, soil moisture data from the Soil Moisture Ocean Salinity (SMOS) satellite is used in this study to correct rainfall accumulation estimates provided by satellite-based real-time precipitation products such as CMORPH, TRMM-3B42RT or PERSIANN. An algorithm based on the SMOS measurements data assimilation is tested in two land-surface models of different complexity: a simple hydrological model (Antecedent Precipitation Index (API)) and a more sophisticated state-of-the-art land-surface model (SURFEX (Surface Extemalisee)). We show how the assimilation technique, based on a particle filter method, generally leads to a significant improvement in rainfall estimates, with slightly better results for the simpler (and less computationally demanding) API model. This methodology has been evaluated for six years at ten sites around the world with different land use and climatological features. The results also show the limitations of the methodology in regions highly affected by mountainous terrain, forest or intense radio-frequency interference (RFI), which can notably affect the quality of the retrievals. The satisfactory results shown here invite the future operational application of the methodology in near-real time on a global scale.

Abstract: Hysteresis in the relationship between bed load transport and river stage is a well-documented phenomenon with multiple known causes. Consequently, numerous studies have interpreted hysteresis in the relationship between seismic ground motion near rivers and some measure of flow strength (i.e., discharge or stage) as the signature of bed load transport. Here we test this hypothesis in the Erlenbach stream (Swiss Prealps) using a metric to quantitatively compare hysteresis in seismic data with hysteresis recorded by geophones attached beneath steel plates within the streambed, a well-calibrated proxy for direct sediment transport measurements. We find that while both the geophones and seismometers demonstrate hysteresis, the magnitude and direction of hysteresis are not significantly correlated between these data, indicating that the seismic signal at this site is primarily reflecting hysteresis in processes other than sediment transport. Seismic hysteresis also does not correlate significantly with the magnitude of sediment transport recorded by the geophones, contrary to previous studies' assumptions. We suggest that hydrologic sources and changes in water turbulence, for instance due to evolving boundary conditions at the bed, rather than changes in sediment transport rates, may sometimes contribute to or even dominate the hysteresis observed in seismic amplitudes near steep mountain rivers. Plain Language Summary An increasing number of studies have recently observed changes in the amount of seismic shaking (hysteresis) recorded near a river at a given discharge during floods. Most studies have assumed that this hysteresis was caused by changes in the amount of sediment being transported in the river and have therefore used the hysteresis to assess sediment transport rates and patterns. We examine concurrent seismic and sediment transport data from a steep mountain stream in the Swiss Prealps and find that changes in seismic shaking are unrelated and even opposed (increasing versus decreasing) to changes in sediment transport rates for four out of five transport events. Water turbulence, rather than sediment transport, appears to be the strongest source of seismic shaking, and changes in seismic shaking are most likely caused by changes in turbulence or how turbulence transmits energy through the river bed. These effects may be due to rearrangement of sediment around large boulders on the bed or slight shifting of the boulders themselves. Our results have significant implications for the growing field of fluvial seismology and the evaluation of seismic data near rivers, as previous interpretations of seismic hysteresis as evidence for sediment transport may not always be accurate.

Abstract: The high-resolution sedimentological and geochemical analysis of a sediment sequence from Lake Savine (Western Mediterranean Alps, France) led to the identification of 220 event layers for the last 6000 years. 200 were triggered by flood events and 20 by underwater mass movements possibly related to earthquakes that occurred in 5 clusters of increase seismicity. Because human activity could influence the flood chronicle, the presence of pastures was reconstructed through ancient DNA, which suggested that the flood chronicle was mainly driven by hydroclimate variability. Weather reanalysis of historical floods allow to identify that mesoscale precipitation events called “East Return” events were the main triggers of floods recorded in Lake Savine. The first part of this palaeoflood record (6-4 kyr BP) was characterized by increases in flood frequency and intensity in phase with Northern Alpine palaeoflood records. By contrast, the second part of the record (i.e., since 4 kyr BP) was phased with Southern Alpine palaeoflood records. These results suggest a palaeohydrological transition at approximately 4 kyr BP, as has been previously described for the Mediterranean region. This may have resulted in a change of flood-prone hydro-meteorological processes, i.e., in the balance between occurrence and intensity of local convective climatic phenomena and their influence on Mediterranean mesoscale precipitation events in this part of the Alps. At a centennial timescale, increases in flood frequency and intensity corresponded to periods of solar minima, affecting climate through atmospheric changes in the Euro-Atlantic sector. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Weather type classifications (WTCs) are a simple tool to analyse variations in weather patterns. Long series of WTCs could be used to address decadal changes in weather as a basis for studying changes in variability or extremes or for addressing contributions of sea-surface temperature or external forcings using climate models. However, there is no long series of daily objective weather types (WTs). A new method (Shortest Mahalanobis Distance, SMD) using daily European weather data is developed to reconstruct WTCs back in time. Here the SMD method is applied on the Cluster Analysis of Principal Components (CAP9) classification used by MeteoSwiss. The CAP9 daily WT time series (computed with ERA-40) is used as reference over the 1958-1998 period. Daily data (temperature, mean sea level pressure and pressure tendency) from 13 European stations covering the period 1763-2009 are used for the reconstruction. The reference CAP9 is reduced from nine to seven types so the new daily WTC is called CAP7. As an assessment, CAP7 is compared to the original classification CAP9 and to the same WTs computed with the Twentieth Century Reanalysis (20CR and 20CRv2c). Over the reference period up to 90% of all the daily WTs can be correctly reproduced in the new WTC compared to the original series, with higher reliability in winter than in summer. In addition, the reliability of the classification is increasing from 1763 onward. The annual occurrence of each type reveals some trends, mostly a decrease in the number of cyclonic days and an increase of cyclonic days.

Abstract: A global 1/4 degrees ocean/sea ice 50-member ensemble simulation is used to disentangle the low-frequency imprints of the atmospherically forced oceanic variability and of the chaotic intrinsic oceanic variability (IOV) on the large-scale (10 degrees x10 degrees) ocean heat content (OHC) between 1980 and 2010. The IOV explains most of the interannual-to-decadal large-scale OHC variance over substantial fractions of the global ocean area that increase with depth: 9%, 22%, and 31% in the 0-700m, 700-2000m and 2000m bottom layers, respectively. Such areas concern principally eddy-active regions, mostly found in the Southern Ocean and in western boundary current extensions, and also concern the subtropical gyres at intermediate and deep levels. The oceanic chaos may also induce random multidecadal fluctuations so that large-scale regional OHC trends computed on the 1980-2010 period cannot be unambiguously attributed to the atmospheric forcing in several oceanic basins at various depths. These results are likely to raise detection and attribution issues from real observations.

Abstract: Three debris-free glaciers with strongly differing annual glaciological glacier-wide mass balances (MBs) are monitored in the Everest region (central Himalaya, Nepal). The mass budget of Mera Glacier (5.1 km(2) in 2012), located in the southern part of this region, was balanced during 2007-15, whereas Pokalde (0.1 km(2) in 2011) and West Changri Nup glaciers (0.9 km(2) in 2013), similar to 30 km further north, have been losing mass rapidly with annual glacier-wide MBs of -0.69 +/- 0.28 m w.e. a(-1) (2009-15) and -1.24 +/- 0.27 m w.e. a(-1) (2010-15), respectively. An analysis of high-elevation meteorological variables reveals that these glaciers are sensitive to precipitation, and to occasional severe cyclonic storms originating from the Bay of Bengal. We observe a negative horizontal gradient of annual precipitation in south-to-north direction across the range (<= -21 mm km(-1), i.e. -2% km(-1)). This contrasted mass-balance pattern over rather short distances is related (i) to the low maximum elevation of Pokalde and West Changri Nup glaciers, resulting in years where their accumulation area ratio is reduced to zero and (ii) to a steeper vertical gradient of MB for glaciers located in the inner arid part of the range.

Abstract: During the lifespan of a Municipal Solid Waste landfill, its leachate drainage system may get clogged. Then, as a consequence of rainfall, leachate generation and possibly leachate injection, the moisture content in the landfill increases to the point that a leachate mound could be created. Therefore, pumping the leachate becomes a necessary solution. This paper presents an original analysis of leachate pumping and injection in an instrumented well. The water table level around the well is monitored by nine piezometers which allow the leachate flow behaviour to be captured. A numerical model based on Richards equation and an exponential relationship between saturated hydraulic conductivity and depth is used to analyze the landfill response to pumping and injection. Decreasing permeability with depth appears to have a major influence on the behaviour of the leachate flow. It could have a drastic negative impact on the pumping efficiency with a maximum quasi-stationary pumping rate limited to approximately 1 m(3)/11 for the tested well and the radius of influence is less than 20 m. The numerical model provides a reasonable description of both pumping and injection tests. However, an anomalous behaviour observed at the transition between pumping and recovery phases is observed. This could be due to a limitation of the Richards model in that it neglects the gas phase behaviour and other double porosity heterogeneous effects. (C) 2016 Elsevier Ltd. All rights reserved.

Abstract: The atmospheric deposition of mercury (Hg) occurs via several mechanisms, including dry and wet scavenging by precipitation events. In an effort to understand the atmospheric cycling and seasonal depositional characteristics of Hg, wet deposition samples were collected for approximately 5 years at 17 selected GMOS monitoring sites located in the Northern and Southern hemispheres in the framework of the Global Mercury Observation System (GMOS) project. Total mercury (THg) exhibited annual and seasonal patterns in Hg wet deposition samples. Interannual differences in total wet deposition are mostly linked with precipitation volume, with the greatest deposition flux occurring in the wettest years. This data set provides a new insight into baseline concentrations of THg concentrations in precipitation worldwide, particularly in regions such as the Southern Hemisphere and tropical areas where wet deposition as well as atmospheric Hg species were not investigated before, opening the way for future and additional simultaneous measurements across the GMOS network as well as new findings in future modeling studies.

Abstract: Quantifying the effect of the seawater density changes on sea level variability is of crucial importance for climate change studies, as the sea level cumulative rise can be regarded as both an important climate change indicator and a possible danger for human activities in coastal areas. In this work, as part of the Ocean Reanalysis Intercomparison Project, the global and regional steric sea level changes are estimated and compared from an ensemble of 16 ocean reanalyses and 4 objective analyses. These estimates are initially compared with a satellite-derived (altimetry minus gravimetry) dataset for a short period (2003-2010). The ensemble mean exhibits a significant high correlation at both global and regional scale, and the ensemble of ocean reanalyses outperforms that of objective analyses, in particular in the Southern Ocean. The reanalysis ensemble mean thus represents a valuable tool for further analyses, although large uncertainties remain for the inter-annual trends. Within the extended intercomparison period that spans the altimetry era (1993-2010), we find that the ensemble of reanalyses and objective analyses are in good agreement, and both detect a trend of the global steric sea level of 1.0 and 1.1 ± 0.05 mm/year, respectively. However, the spread among the products of the halosteric component trend exceeds the mean trend itself, questioning the reliability of its estimate. This is related to the scarcity of salinity observations before the Argo era. Furthermore, the impact of deep ocean layers is non-negligible on the steric sea level variability (22 and 12 % for the layers below 700 and 1500 m of depth, respectively), although the small deep ocean trends are not significant with respect to the products spread.

Abstract: The Mekong Delta, is home to 17 million inhabitants and faces numerous challenges relating to climate change, environmental degradation and water issues. In this study, we assess trace metals concentrations (Al, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Mo, Cd, Hg, Pb) in the water, suspended particulate matter and surface sediments of the Tien River, the Northern branch of the Mekong Delta, during both dry and rainy seasons. Metal concentrations in the dissolved and suspended particle phases remain in the low concentration range of the main Asian Tropical River. During transportation in the riverine part, we evidenced that V, Cr, Co, As and Pb are dominant in the particulate phase while Mo, Ni and Cu dominate in the dissolved fraction. In the salinity gradient, dissolved U, V, Mo exhibit conservative behaviour while Ni, Cu, As, Co and Cd showed additive behaviour suggesting desorption processes. In the surface sediment, metal concentrations are controlled by the particle-size, POC contents and Fe, Al and Mn – oxy(hydr)oxides. Calculated Enrichment Factor and Geoaccumulation Index evidenced As enrichment while the calculated mean effect range median quotients evidenced a low to medium ecotoxicological potential effects range in the surface sediments. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Geosynthetics have been extensively used in landfills as a lining system to prevent leachate infiltration into groundwater. In piggy-back landfill expansion (PBLE), consisting of building a new landfill over an existing one, a lining system is implemented between the old and new waste. In this context, interface failure (stability) and deformation (integrity) of the lining system should be considered for the design. Such stability and integrity mainly depend on the PBLE geometry and the mechanical properties of the geosynthetics. Comprehensive numerical modelling simulations were performed to show how these factors influence the shear stresses, shear displacements, translational stability and the axial strains/ forces within the various geosynthetics. The numerical modelling was conducted using the finite difference code FLAC 2D, focusing on a typical PBLE and considering geosynthetic interface strain softening, the nonlinear stiffness of geosynthetics, and the differentiation between the compressive and tensile behaviours of geosynthetics. The simulations showed that the lateral length of the PBLE, the type of geomembrane (textured or smooth) and the level of the leachate table in new waste are the factors that most influenced the mechanical behaviour of the lining system and its stability. Finally, a parameter called the stability ratio was proposed as a complement to the traditional factor of safety, to analyse the progressive slippage along the geosynthetic interfaces. The numerical results indicated that interface failure concurrently begins at the rightmost part of the lower flat area of the PBLE and near the corner of the inner slope before spreading out to the left as backfilling progresses.

Abstract: The hydrological cycle is expected to intensify under global warming(1), with studies reporting more frequent extreme rain events in many regions of the world(2-4), and predicting increases in future flood frequency(5). Such early, predominantly mid-latitude observations are essential because of shortcomings within climate models in their depiction of convective rainfall(6,7). A globally important group of intense storms-mesoscale convective systems (MCSs) 8-poses a particular challenge, because they organize dynamically on spatial scales that cannot be resolved by conventional climate models. Here, we use 35 years of satellite observations from the West African Sahel to reveal a persistent increase in the frequency of the most intense MCSs. Sahelian storms are some of the most powerful on the planet(9), and rain gauges in this region have recorded a rise in 'extreme'(17) daily rainfall totals. We find that intense MCS frequency is only weakly related to the multidecadal recovery of Sahel annual rainfall, but is highly correlated with global land temperatures. Analysis of trends across Africa reveals that MCS intensification is limited to a narrow band south of the Sahara desert. During this period, wet-season Sahelian temperatures have not risen, ruling out the possibility that rainfall has intensified in response to locally warmer conditions. On the other hand, the meridional temperature gradient spanning the Sahel has increased in recent decades, consistent with anthropogenic forcing driving enhanced Saharan warming(10). We argue that Saharan warming intensifies convection within Sahelian MCSs through increased wind shear and changes to the Saharan air layer. The meridional gradient is projected to strengthen throughout the twenty-first century, suggesting that the Sahel will experience particularly marked increases in extreme rain. The remarkably rapid intensification of Sahelian MCSs since the 1980s sheds new light on the response of organized tropical convection to global warming, and challenges conventional projections made by general circulation models.

Abstract: Total individual concentrations (in both gaseous and particulate phases) of 80 polycyclic aromatic compounds (PACs) including 32 nitro-PAHs, 27 oxy-PAHs (polycyclic aromatic hydrocarbons) and 21 parent PAHs have been investigated over a year in the ambient air of Grenoble (France) together with an extended aerosol chemical characterization. The results indicated that their concentrations were strongly affected by primary emissions in cold period, especially from residential heating (i.e. biomass burning). Besides, secondary processes occurred in summer but also in cold period under specific conditions such as during long thermal inversion layer periods and severe PM pollution events. Different secondary processes were involved during both PM pollution events observed in March April and in December 2013. During the first one, long range transport of air masses, nitrate chemistry and secondary nitro-PAH formation seemed linked. During the second one, the accumulation of primary pollutants over several consecutive days enhanced secondary chemical processes notably highlighted by the dramatic increase of oxy-PAH concentrations. The study of the time trends of ratios of individual nitro- or oxy-PAHs to parent PAHs, in combination with key primary or secondary aerosol species and literature data, allowed the identification of potential molecular markers of PAH oxidation. Finally, 6H-dibenzo[b,d]pyran-6-one, biphenyl-2,2'-dicarboxaldehyde and 3-nitrophenanthrene have been selected to be the best candidates as markers of PAH oxidation processes in ambient air. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Current understanding of mercury (Hg) behavior in the atmosphere contains significant gaps. Some key characteristics of Hg processes, including anthropogenic and geogenic emissions, atmospheric chemistry, and air-surface exchange, are still poorly known. This study provides a complex analysis of processes governing Hg fate in the atmosphere involving both measured data from ground-based sites and simulation results from chemical transport models. A variety of long-term measurements of gaseous elemental Hg (GEM) and reactive Hg (RM) concentration as well as Hg wet deposition flux have been compiled from different global and regional monitoring networks. Four contemporary global-scale transport models for Hg were used, both in their state-of-the-art configurations and for a number of numerical experiments to evaluate particular processes. Results of the model simulations were evaluated against measurements. As follows from the analysis, the interhemispheric GEM gradient is largely formed by the prevailing spatial distribution of anthropogenic emissions in the Northern Hemisphere. The contributions of natural and secondary emissions enhance the south-to-north gradient, but their effect is less significant. Atmospheric chemistry has a limited effect on the spatial distribution and temporal variation of GEM concentration in surface air. In contrast, RM air concentration and wet deposition are largely defined by oxidation chemistry. The Br oxidation mechanism can reproduce successfully the observed seasonal variation of the RM/GEM ratio in the near-surface layer, but it predicts a wet deposition maximum in spring instead of in summer as observed at monitoring sites in North America and Europe. Model runs with OH chemistry correctly simulate both the periods of maximum and minimum values and the amplitude of observed seasonal variation but shift the maximum RM/GEM ratios from spring to summer. O-3 chemistry does not predict significant seasonal variation of Hg oxidation. Hence, the performance of the Hg oxidation mechanisms under study differs in the extent to which they can reproduce the various observed parameters. This variation implies possibility of more complex chemistry and multiple Hg oxidation pathways occurring concurrently in various parts of the atmosphere.

Abstract: Light-absorbing impurities (LAIs) decrease snow albedo, increasing the amount of solar energy absorbed by the snowpack. Its most intuitive and direct impact is to accelerate snowmelt. Enhanced energy absorption in snow also modifies snow metamorphism, which can indirectly drive further variations of snow albedo in the near-infrared part of the solar spectrum because of the evolution of the near-surface snow microstructure. New capabilities have been implemented in the detailed snowpack model SURFEX/ISBACrocus (referred to as Crocus) to account for impurities' deposition and evolution within the snowpack and their direct and indirect impacts. Once deposited, the model computes impurities' mass evolution until snow melts out, accounting for scavenging by meltwater. Taking advantage of the recent inclusion of the spectral radiative transfer model TARTES (Two-stream Analytical Radiative TransfEr in Snow model) in Crocus, the model explicitly represents the radiative impacts of light-absorbing impurities in snow. The model was evaluated at the Col de Porte experimental site (French Alps) during the 2013-2014 snow season against in situ standard snow measurements and spectral albedo measurements. In situ meteorological measurements were used to drive the snowpack model, except for aerosol deposition fluxes. Black carbon (BC) and dust deposition fluxes used to drive the model were extracted from simulations of the atmospheric model ALADIN-Climate. The model simulates snowpack evolution reasonably, providing similar performances to our reference Crocus version in terms of snow depth, snow water equivalent (SWE), near-surface specific surface area (SSA) and shortwave albedo. Since the reference empirical albedo scheme was calibrated at the Col de Porte, improvements were not expected to be significant in this study. We show that the deposition fluxes from the ALADIN-Climate model provide a reasonable estimate of the amount of light-absorbing impurities deposited on the snowpack except for extreme deposition events which are greatly underestimated. For this particular season, the simulated melt-out date advances by 6 to 9 days due to the presence of light-absorbing impurities. The model makes it possible to apportion the relative importance of direct and indirect impacts of light-absorbing impurities on energy absorption in snow. For the snow season considered, the direct impact in the visible part of the solar spectrum accounts for 85% of the total impact, while the indirect impact related to accelerated snow metamorphi

Abstract: Water availability is the primary constraint on the improvement of food security in rural areas in northwestern Cambodia. A 4-year study was carried out in the upper Stung Sreng watershed to assess water resources. Four sub-watersheds with different land cover types, ranging in size from 1.5 to 185km(2), were monitored using dedicated weather stations and rain- and streamgauges. Geophysics and observation boreholes were used to characterize aquifers. Rainwater is mostly split into evapotranspiration (annual mean of 54% rainfall) and streamflow components (49%), because groundwater recharge is low (1%). Thus, rainwater and streamflow are the main sources for irrigation development. Groundwater can be used only in specific locations for low water-demand crops. A total of 186 household ponds and three village-scale dams were built and 31 wells were installed. The household pond was determined to be the best solution for irrigation development because of its simple management.

Abstract: This study aims to find the typical growth rate of the temperature inversion during the onset of the stable boundary layer around sunset. The sunset transition is a very challenging period for numerical weather prediction, since neither accepted theories for the convective boundary layer nor those for the stable boundary layer appear to be applicable. To gain more insight in this period, a systematic investigation of the temperature inversion growth rate is conducted. A statistical procedure is used to analyze almost 16 years of observations from the Cabauw observational tower, supported by observations from two additional sites (Dome C and Karlsruhe). The results show that, on average, the growth rate of the temperature inversion (normalized by the maximum inversion during the night) weakly declines with increasing wind speed. The observed growth rate is quantitatively consistent among the sites, and it appears insensitive to various other parameters. The results were also insensitive to the afternoon decay rate of the net radiation except when this decay rate was very weak. These observations are compared to numerical solutions of three models with increasing complexity: a bulk model, an idealized single-column model (SCM), and an operational-level SCM. It appears only the latter could reproduce qualitative features of the observations using a first-order closure. Moreover, replacing this closure with a prognostic TKE scheme substantially improved the quantitative performance. This suggests that idealized models assuming instantaneous equilibrium flux-profile relations may not aid in understanding this period, since history effects may qualitatively affect the dynamics.

Abstract: An extensive meteorological observational dataset at Dome C, East Antarctic Plateau, enabled estimation of the sensitivity of surface momentum and sensible heat fluxes to aerodynamic roughness length and atmospheric stability in this region. Our study reveals that (1) because of the preferential orientation of snow micro-reliefs (sastrugi), the aerodynamic roughness length varies by more than two orders of magnitude depending on the wind direction; consequently, estimating the turbulent fluxes with a realistic but constant of 1 mm leads to a mean friction velocity bias of in near-neutral conditions; (2) the dependence of the ratio of the roughness length for heat to on the roughness Reynolds number is shown to be in reasonable agreement with previous models; (3) the wide range of atmospheric stability at Dome C makes the flux very sensitive to the choice of the stability functions; stability function models presumed to be suitable for stable conditions were evaluated and shown to generally underestimate the dimensionless vertical temperature gradient; as these models differ increasingly with increases in the stability parameter z / L, heat flux and friction velocity relative differences reached when ; (4) the shallowness of the stable boundary layer is responsible for significant sensitivity to the height of the observed temperature and wind data used to estimate the fluxes. Consistent flux results were obtained with atmospheric measurements at heights up to 2 m. Our sensitivity study revealed the need to include a dynamical parametrization of roughness length over Antarctica in climate models and to develop new parametrizations of the surface fluxes in very stable conditions, accounting, for instance, for the divergence in both radiative and turbulent fluxes in the first few metres of the boundary layer.

Abstract: Investigation of meteorological measurements along a 45m tower at Dome C on the high East Antarctic Plateau revealed two distinct stable boundary layer (SBL) regimes at this location. The first regime is characterized by strong winds and continuous turbulence. It results in full vertical coupling of temperature, wind magnitude and wind direction in the SBL. The second regime is characterized by weak winds, associated with weak turbulent activity and very strong temperature inversions reaching up to 25 K in the lowest 10 m. Vertical temperature profiles are generally exponentially shaped (convex) in the first regime and 'convex-concave-convex' in the second. The transition between the two regimes is particularly abrupt when looking at the near-surface temperature inversion and it can be identified by a 10 m wind-speed threshold. With winds under this threshold, the turbulent heat supply toward the surface becomes significantly lower than the net surface radiative cooling. The threshold value (including its range of uncertainty) appears to agree with recent theoretical predictions from the so-called 'minimum wind speed for sustainable turbulence' (MWST) theory. For the quasi-steady, clear-sky winter cases, the relation between the near-surface inversion amplitude and the wind speed takes a characteristic 'S' shape. Closer analysis suggests that this relation corresponds to a 'critical transition' between a steady turbulent and a steady 'radiative' regime, with a dynamically unstable branch in the transition zone. These fascinating characteristics of the Antarctic boundary layer challenge present and future numerical models to represent this region in a physically correct manner.

Abstract: Snow cover variability has significant effects on local and global climate evolution. By changing surface energy fluxes and hydrological conditions, changes in snow cover can alter atmospheric circulation and lead to remote climate effects. To document such multi-scale climate effects, atmospheric reanalysis and derived products offer the opportunity to analyze snow variability in great detail far back to the early 20th century. So far only little is know about their quality. Comparing snow depth in four long-term reanalysis datasets with Russian in situ snow depth data, we find a moderately high daily correlation (around 0.6-0.7), which is comparable to correlations for the recent era (1981-2010), and a good representation of sub-decadal variability. However, the representation of pre-1950 inter-decadal snow variability is questionable, since reanalysis products divert towards different base states. Limited availability of independent long-term snow data makes it difficult to assess the exact cause for this bifurcation in snow states, but initial investigations point towards representation of the atmosphere rather than differences in assimilated data or snow schemes. This study demonstrates the ability of long-term reanalysis to reproduce snow variability accordingly.

Abstract: Mediterranean Alpine populations are particularly exposed to natural hazards like floods and earthquakes because of both the close Mediterranean humidity source and the seismically active Alpine region. Knowledge of longterm variability in flood and earthquake occurrences is of high value since it can be useful to improve risk assessment and mitigation. In this context, we explore the potential of a lake-sediment sequence from Lago Inferiore de Laures in Valle d'Aosta (Northern Italy) as a long-term record of past floods and earthquakes. The high-resolution sedimentological study revealed 76 event layers over the last ca. 270 years; 8 are interpreted as most probably induced by earthquakes and 68 by flood events. Comparison to historical seismic data suggests that the recorded earthquakes are strong (epicentral Medvedev-Sponheuer-Karnik (MSK) intensity of VI-IX) and/or close to the lake (distance of 25-120 km). Compared to other lake-sediment sequences, Lago Inferiore de Laures sediments appear to be regionally the most sensitive to earthquake shaking, offering a great potential to reconstruct the past regional seismicity further back in time. Comparison to historical and palaeoflood records suggests that the flood signal reconstructed from Lago Inferiore de Laures sediments represents the regional and (multi-) decadal variability of summer-autumn floods well, in connection to Mediterranean mesoscale precipitation events. Overall, our results reveal the high potential of Lago Inferiore de Laures sediments to extend the regional earthquake and flood catalogues far back in time.

Abstract: In this work, the problem of intracellular currents in longilinear bacteria, such as Escherichia coli, suspended in a physiological medium and submitted to a harmonic voltage (AC), is analyzed using the Finite-Element-based software COMSOL Multiphysics. Bacterium was modeled as a cylindrical capsule, ended by semi-spheres and surrounded by a dielectric cell wall. An equivalent single-layer cell wall was defined, starting from the well-recognized three-shell modeling approach. The bacterium was considered immersed in a physiological medium, which was also taken into account in the modeling. A new complex transconductance was thus introduced, relating the complex ratio between current inside the bacterium and voltage applied between two parallel equipotential planes, separated by a realistic distance. When voltage was applied longitudinally relative to the bacterium main axis, numerical results in terms of frequency response in the 1-20 MHz range for E. coli cells revealed that transconductance magnitude exhibited a maximum at a frequency depending on the cell wall capacitance. This occurred in spite of the purely passive character of the model and could be explained by an equivalent electrical network giving very similar results and showing special conditions for lateral paths of the currents through the cell wall. It is shown that the main contribution to this behavior is due to the conductive part of the current. (C) 2016 Wiley Periodicals, Inc.

Abstract: In this paper, three incremental analysis update schemes (IAU 0, IAU 50 and IAU 100) are compared in the same assimilation experiments with a realistic eddy permitting primitive equation model of the North Atlantic Ocean using the Ensemble Kalman Filter. The difference between the three IAU schemes lies on the position of the increment update window. The relevance of each IAU scheme is evaluated through analyses on both thermohaline and dynamical variables. The validation of the assimilation results is performed according to both deterministic and probabilistic metrics against different sources of observations. For deterministic validation, the ensemble mean and the ensemble spread are compared to the observations. For probabilistic validation, the continuous ranked probability score (CRPS) is used to evaluate the ensemble forecast system according to reliability and resolution. The reliability is further decomposed into bias and dispersion by the reduced centred random variable (RCRV) score. The obtained results show that 1) the IAU 50 scheme has the same performance as the IAU 100 scheme 2) the IAU 50/100 schemes outperform the IAU 0 scheme in error covariance propagation for thermohaline variables in relatively stable region, while the IAU 0 scheme outperforms the IAU 50/100 schemes in dynamical variables estimation in dynamically active region 3) in case with sufficient number of observations and good error specification, the impact of IAU schemes is negligible. The differences between the IAU 0 scheme and the IAU 50/100 schemes are mainly due to different model integration time and different instability (density inversion, large vertical velocity, etc.) induced by the increment update. The longer model integration time with the IAU 50/100 schemes, especially the free model integration, on one hand, allows for better re-establishment of the equilibrium model state, on the other hand, smooths the strong gradients in dynamically active region. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: We present a comparison of the numerics and simulation results for two “full” Stokes ice sheet models, FELIX-S (Leng et al., 2012) and Elmer/Ice (Gagliardini et al., 2013). The models are applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MIS-MIP3d). For the diagnostic experiment (P75D) the two models give similar results (<2% difference with respect to along-flow velocities) when using identical geometries and computational meshes, which we interpret as an indication of inherent consistencies and similarities between the two models. For the standard (Stnd), P75S, and P75R prognostic experiments, we find that FELIX-S (Elmer/Ice) grounding lines are relatively more retreated (advanced), results that are consistent with minor differences observed in the diagnostic experiment results and that we show to be due to different choices in the implementation of basal boundary conditions in the two models. While we are not able to argue for the relative favorability of either implementation, we do show that these differences decrease with increasing horizontal (i.e., both along-and across-flow) grid resolution and that grounding-line positions for FELIX-S and Elmer/Ice converge to within the estimated truncation error for Elmer/Ice. Stokes model solutions are often treated as an accuracy metric in model intercomparison experiments, but computational cost may not always allow for the use of model resolution within the regime of asymptotic convergence. In this case, we propose that an alternative estimate for the uncertainty in the grounding-line position is the span of grounding-line positions predicted by multiple Stokes models.

Abstract: CO (2) fluxes were measured during 18 months in a forest and a savannah in northern Benin. Higher values of carbon fluxes were found during the wet season at each site. A strong dependency of carbon fluxes on water relations was found in two contrasting sites. The forest sequestered 640 +/- 50 and the savannah 190 +/- 40 g C m (-2) year (-1) . In West Africa, the main mechanisms or factors governing the dynamics of ecosystems, especially the dynamics of the carbon fluxes and productivity, still remain less known. This study reports the carbon fluxes over two contrasting ecosystems, notably a protected forest (lat 9.79A degrees N, long 1.72A degrees E, alt 414 m) and a cultivated savannah (lat 9.74A degrees N, long 1.60A degrees E, alt 449 m) in northern Benin. The two sites were among those equipped by the AMMA-CATCH observatory and Ou,m, 2025 project. Flux data were analyzed at the daily and seasonal scales in order to understand their controlling variables. We discussed the patterns of CO2 fluxes and the characteristics of the two ecosystems. The study also focused on the different water usage strategies developed by the two ecosystems since the alternation between dry and wet seasons highly influenced the seasonal dynamics. Finally, the annual carbon sequestration was estimated together with its uncertainty. The carbon fluxes were measured during 18 months (July 2008-December 2009) by an eddy-covariance system over two contrasting sites in northern Benin. Fluxes data were computed following the standard procedure. The responses of CO2 fluxes to the principal climatic and edaphic factors, and the canopy conductance were studied. A clear CO2 fluxes response to main environmental factors was observed, however with difference according to the seasons and vegetation types. The ecosystem respiration showed the highest values during the wet season and a progressive decrease from wet to dry periods. Also, the carbon uptake values were high during the wet period, but low during the dry period. However, the CO2 fluxes for the protected forest were always higher than that for the cultivated savannah within each defined period. This was due to the seasonal changes not only in phenology and physiology but also to the acclimation to environmental conditions, especially to the soil water availability. The water use efficiency was influenced by VPD during the day conditions for two ecosystems. However, the VPD response curve of water usage was relatively constant for the protected forest during the transitional and wet seasons. In contrary, for the cultivated savannah the VPD response decreased about 46 % from transitional dry-wet to wet seasons and remained relatively constant until transitional wet-dry season. The close relationships between the net CO2 assimilation and the canopy conductance were found for the two ecosystems. This suggests a regulation of the stomata by a partial stomatal closure besides the radiation control. Finally, the forest sequestered at the annual scale 640 +/- 50 and the savannah 190 +/- 40 g C m(-2) year(-1). We conclude with a strong relation between the carbon fluxes and water in the two investigated ecosystems. Apart from the radiation control, the stomata also play an important role in the regulation of the CO2 assimilation in the two ecosystems.

Abstract: In polar regions, sastrugi are a direct manifestation of drifting snow and form the main surface roughness elements. In turn, sastrugi alter the generation of atmospheric turbulence and thus modify the wind field and the aeolian snow mass fluxes. Little attention has been paid to these feedback processes, mainly because of experimental difficulties. As a result, most polar atmospheric models currently ignore sastrugi over snow-covered regions. This paper aims at quantifying the potential influence of sastrugi on the local wind field and on snow erosion over a sastrugi-covered snowfield in coastal Ad,lie Land, East Antarctica. We focus on two erosion events during which sastrugi responses to shifts in wind direction have been interpreted from temporal variations in drag and aeolian snow mass flux measurements during austral winter 2013. Using this data set, it is shown that (i) neutral stability, 10aEuro-m drag coefficient (C-DN10) values are in the range of 1.3-1.5 x 10(-3) when the wind is well aligned with the sastrugi, (ii) as the wind shifts by only 20-30A degrees away from the streamlined direction, C-DN10 increases (by 30-120aEuro-%) and the aeolian snow mass flux decreases (by 30-80aEuro-%), thereby reflecting the growing contribution of the sastrugi form drag to the total surface drag and its inhibiting effect on snow erosion, (iii) the timescale of sastrugi aerodynamic adjustment can be as short as 3aEuro-h for friction velocities greater than 1aEuro-maEuro-s(-1) and during strong drifting snow conditions and (iv) knowing C-DN10 is not sufficient to estimate the snow erosion flux that results from drag partitioning at the surface because C-DN10 includes the contribution of the sastrugi form drag.

Abstract: Under the framework of the Global Mercury Observation System (GMOS) project, a 3.5-year record of atmospheric gaseous elemental mercury (Hg(0)) has been gathered at Dumont d'Urville (DDU, 66 degrees 40'S, 140 degrees 01'E, 43m above sea level) on the East Antarctic coast. Additionally, surface snow samples were collected in February 2009 during a traverse between Concordia Station located on the East Antarctic plateau and DDU. The record of atmospheric Hg(0) at DDU reveals particularities that are not seen at other coastal sites: a gradual decrease of concentrations over the course of winter, and a daily maximum concentration around midday in summer. Additionally, total mercury concentrations in surface snow samples were particularly elevated near DDU (up to 194.4 ng L-1) as compared to measurements at other coastal Antarctic sites. These differences can be explained by the more frequent arrival of inland air masses at DDU than at other coastal sites. This confirms the influence of processes observed on the Antarctic plateau on the cycle of atmospheric mercury at a continental scale, especially in areas subject to recurrent katabatic winds. DDU is also influenced by oceanic air masses and our data suggest that the ocean plays a dual role on Hg(0) concentrations. The open ocean may represent a source of atmospheric Hg(0) in summer whereas the sea-ice surface may provide reactive halogens in spring that can oxidize Hg(0). This paper also discusses implications for coastal Antarctic ecosystems and for the cycle of atmospheric mercury in high southern latitudes.

Abstract: We review the scientific efforts over the last decades to reconstruct erosion from continuous alpine lake sediment records. We focused both on methodological issues, showing the growing importance of non-destructive high resolution approaches (XRF core-scanner) as well as progresses in the understanding of processes leading to the creation of an “erosion signal” in lakes. We distinguish “continuous records” from “event-records”. Both provide complementary information but need to be studied with different approaches. Continuous regionally-relevant records proved to be particularly pertinent to document regional erosion patterns throughout the Holocene, in particular applying the source to sink approach. Event-based approaches demonstrated and took advantage of the strong non-linearity of sediment transport in high altitude catchment areas. This led to flood frequency and intensity reconstructions, highlighting the influence of climate change upon flood dynamics in the mountain. The combination of different record types, both in terms of location (high vs. low elevation), sedimentology (high vs. low terrigenous contribution) and significance (local vs. regional) is one of the main outputs of this paper. It allows the establishment of comprehensive histories of NW French Alps erosion, but also and consequently, soil dynamics and hydrological patterns throughout the Holocene. We also discuss the influence of glacier dynamics, one of the major agents of erosion in the Alps. A major feature is the growing human influence upon erosion at a local scale since at least the middle of the Bronze Age (3500 cal. BP). However and according to the regional record from Lake Bourget, only few periods of rising erosion at local scales generated a regional record that can be discriminated from wetter climatic periods. Among them, the period between 200 BCE and 400 AD appeared to be marked by a generalised rise in human-triggered erosion at local scales in the northern French Alps. This review highlights the importance of modern high-resolution and interdisciplinary studies of lake sediments, in order to better understand the complex relationships between humans, climate and the Earth system in general. We strongly argue that regional integration of data is now required to move a step further. Such an integration is easier with cost- and time-effective methods as well as after a better definition of approaches and their limits. This should lead to a stronger collaboration between paleo-data producers and modellers in the near future. (C) 2016 Elsevier Ltd. All rights reserved.

Abstract: Coupled ice sheet-ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet-Ocean Model Intercomparison Project ( MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects ( MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet ( WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase ( MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase ( ISOMIP+) and coupled ice sheet-ocean experiments under the MISOMIP first phase ( MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations ( MISOMIP1) to be compared directly to the individual component simulations ( MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.

Abstract: Underground structures have been shown to have a great influence on subsoil resources in urban aquifers. A methodology to assess the actual and the potential state of the groundwater flow in an urban area is proposed. The study develops a three-dimensional modeling approach to understand the cumulative impacts of underground infrastructures on urban groundwater flow, using a case in the city of Lyon (France). All known underground structures were integrated in the numerical model. Several simulations were run: the actual state of groundwater flow, the potential state of groundwater flow (without underground structures), an intermediate state (without impervious structures), and a transient simulation of the actual state of groundwater flow. The results show that underground structures fragment groundwater flow systems leading to a modification of the aquifer regime. For the case studied, the flow systems are shown to be stable over time with a transient simulation. Structures with drainage systems are shown to have a major impact on flow systems. The barrier effect of impervious structures was negligible because of the small hydraulic gradient of the area. The study demonstrates that the definition of a potential urban groundwater flow and the depiction of urban flow systems, which involves understanding the impact of underground structures, are important issues with respect to urban underground planning.

Abstract: Property economics favours the vertical development of cities but flow of groundwater can be affected by the use of underground space in them. This review article presents the state of the art regarding the impact of disturbances caused by underground structures (tunnels, basements of buildings, deep foundations, etc.) on the groundwater flow in urban aquifers. The structures built in the underground levels of urban areas are presented and organised in terms of their impact on flow: obstacle to the flow or disturbance of the groundwater budget of the flow system. These two types of disturbance are described in relation to the structure area and the urban area. The work reviewed shows, on one hand, the individual impacts of different urban underground structures, and on the other, their cumulative impacts on flow, using real case studies. Lastly, the works are placed in perspective regarding the integration of underground structures with the aim of operational management of an urban aquifer. The literature presents deterministic numerical modelling as a tool capable of contributing to this aim, in that it helps to quantify the effect of an underground infrastructure project on groundwater flow, which is crucial for decision-making processes. It can also be an operational decision-aid tool for choosing construction techniques or for formulating strategies to manage the water resource.

Abstract: We investigate the problem of balancing model complexity and input data requirements in snow hydrology. For this purpose, we analyze the performance of two models of different complexity in estimating variables of interest in snow hydrology applications. These are snow depth, bulk snow density, snow water equivalent and snowmelt run-off. We quantify the differences between data and model prediction using 18years of measurements from an experimental site in the French Alps (Col de Porte, 1325m AMSL). The models involved in this comparison are a one-layer temperature-index model (HyS) and a multilayer model (Crocus). Results show that the expected loss in performance in the one-layer temperature-index model with respect to the multilayer model is low when considering snow depth, snow water equivalent and bulk snow density. As for run-off, the comparison returns less clear indications for identification of a balance. In particular, differences between the models' prediction and data with an hourly resolution are higher when considering the Crocus model than the HyS model. However, Crocus is better at reproducing sub-daily cycles in this variable. In terms of daily run-off, the multilayer physically based model seems to be a better choice, while results in terms of cumulative run-off are comparable. The better reproduction of daily and sub-daily variability of run-off suggests that use of the multilayer model may be preferable for this purpose. Variation in performance is discussed as a function of both the calibration solution chosen and the time of year. Copyright (c) 2016 John Wiley & Sons, Ltd.

Abstract: A temperature reconstruction in the glacierized Himalaya over the past centuries using glacial length fluctuation records is challenging due to the abundance of debris-covered glaciers and a scarcity of glacial length fluctuation data. Using idealized flowline model simulations, we show that supraglacial debris cover significantly alters the length fluctuations only when the debris cover is very thick. An expanded database of length fluctuation records for 43 glaciers in the Himalaya and Karakoram is compiled and a standard linear inversion procedure is applied to a subset of 34 glaciers in this database. The reconstructed temperature anomaly during 1860-2010 indicates a continued warming of the region with a total temperature change of similar to 1.6 K. A close resemblance of the regional temperature anomaly to global trends is seen.

Abstract: In this paper, we reanalyze the glacier mass balance on Glaciar Antisana 15 alpha over the 1995-2012 period. Annual glacier mass balances were quantified on the basis of monthly glaciological measurements using an adaptation of Lliboutry's statistical approach. The geodetic mass balance was computed between 1997 and 2009 giving a cumulative balance of -1.39 +/- 1.97 m w.e. and a slightly negative adjusted annual glaciological mass balance (-0.12 +/- 0.16 m w.e. a(-1)). Despite a careful analysis of uncertainties, we found a large discrepancy between the cumulative glaciological and the geodetic mass balances over the common period, of 4.66 m w.e. This discrepancy can mainly be explained by underestimated net accumulation in the glacier upper reaches, which could be due to the peculiar climate conditions of the equatorial zone with year round accumulation, thereby preventing clear identification of annual layers. An increase of similar to 70% in measured rates of net accumulation would be needed to balance the glaciological and geodetic mass balances; a hypothesis confirmed by estimated ice flux in the vicinity of the ELA. Consequently, the vertical gradient of precipitation may be higher than previously estimated and the accumulation processes (including the role of frost deposition) need to be carefully analyzed.

Abstract: This paper presents some improvements of a probabilistic quantitative precipitation forecasting method based on analogues, formerly developed on small basins located in South-Eastern France. The scope is extended to large scale basins mainly influenced by frontal systems, considering a case study area related to the Saone river, a large basin in eastern France. For a given target situation, this method consists in searching for the most similar situations observed in a historical meteorological archive. Precipitation amounts observed during analogous situations are then collected to derive an empirical predictive distribution function, i.e. the probabilistic estimation of the precipitation amount expected for the target day. The former version of this forecasting method (Bontron, 2004) has been improved by introducing two innovative variables: temperature, that allows taking seasonal effects into account and vertical velocity, which enables a better characterization of the vertical atmospheric motion. The new algorithm is first applied in a perfect prognosis context (target situations come from a meteorological reanalysis) and then in an operational forecasting context (target situations come from weather forecasts) for a three years period. Results show that this approach yields useful forecasts, with a lower false alarm rate and improved performances from the present day D to day D + 2. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Better constraining the continental climate response to anthropogenic forcing is essential to improve climate projections. In this study, pattern scaling is used to extract, from observations, the patterned response of sea surface temperature (SST) and sea ice concentration (SICE) to anthropogenically dominated long-term global warming. The SST response pattern includes a warming of the tropical Indian Ocean, the high northern latitudes, and the western boundary currents. The SICE pattern shows seasonal variations of the main locations of sea ice loss. These SST-SICE response patterns are used to drive an ensemble of an atmospheric general circulation model, the National Center for Atmospheric Research (NCAR) Community Atmosphere Model, version 5 (CAM5), over the period 1980-2010 along with a standard AMIP ensemble using observed SST-SICE. The simulations enable attribution of a variety of observed trends of continental climate to global warming. On the one hand, the warming trends observed in all seasons across the entire Northern Hemisphere extratropics result from global warming, as does the snow loss observed over the northern midlatitudes and northwestern Eurasia. On the other hand, 1980-2010 precipitation trends observed in winter over North America and in summer over Africa result from the recent decreasing phase of the Pacific decadal oscillation and the recent increasing phase of the Atlantic multidecadal oscillation, respectively, which are not part of the global warming signal. The method holds promise for near-term decadal climate prediction but as currently framed cannot distinguish regional signals associated with oceanic internal variability from aerosol forcing and other sources of short-term forcing.

Abstract: Snowpack is a multiphase (photo) chemical reactor that strongly influences the air composition in polar and snow-covered regions. Snowpack plays a special role in the nitrogen cycle, as it has been shown that nitrate undergoes numerous recycling stages (including photolysis) in the snow before being permanently buried in the ice. However, the current understanding of these physicochemical processes remains very poor. Several modelling studies have attempted to reproduce (photo) chemical reactions inside snow grains, but these have relied on strong assumptions to characterise snow reactive properties, which are not well defined. Air-snow exchange processes such as adsorption, solid-state diffusion, or co-condensation also affect snow chemical composition. Here, we present a physically based model of these processes for nitrate. Using as input a 1-year-long time series of atmospheric nitrate concentration measured at Dome C, Antarctica, our model reproduces with good agreement the nitrate measurements in the surface snow. By investigating the relative importance of the main exchange processes, this study shows that, on the one hand, the combination of bulk diffusion and co-condensation allows a good reproduction of the measurements (correlation coefficient r = 0.95), with a correct amplitude and timing of summer peak concentration of nitrate in snow. During winter, nitrate concentration in surface snow is mainly driven by thermodynamic equilibrium, whilst the peak observed in summer is explained by the kinetic process of co-condensation. On the other hand, the adsorption of nitric acid on the surface of the snow grains, constrained by an already existing parameterisation for the isotherm, fails to fit the observed variations. During winter and spring, the modelled concentration of adsorbed nitrate is respectively 2.5 and 8.3-fold higher than the measured one. A strong diurnal variation driven by the temperature cycle and a peak occurring in early spring are two other major features that do not match the measurements. This study clearly demonstrates that co-condensation is the most important process to explain nitrate incorporation in snow undergoing temperature gradient metamorphism. The parameterisation developed for this process can now be used as a foundation piece in snowpack models to predict the inter-relationship between snow physical evolution and snow nitrate chemistry.

Abstract: Atmospheric particulate matter (PM) affects the climate in various ways and has a negative impact on human health. In populated mountain valleys in Alpine regions, emissions from road traffic contribute to carbonaceous aerosols, but residential wood burning can be another source of PM during winter. We determine the contribution of fossil and non-fossil carbon sources by measuring radiocarbon in aerosols using the recently installed AixMICADAS facility. The accelerator mass spectrometer is coupled to an elemental analyzer (EA) by means of a gas interface system directly connected to the gas ion source. This system provides rapid and accurate radiocarbon measurements for small samples (10-100 μgC) with minimal preparation from the aerosol filters. We show how the contamination induced by the EA protocol can be quantified and corrected for. Several standards and synthetic samples are then used to demonstrate the precision and accuracy of aerosol measurements over the full range of expected C-14/C-12 ratios, ranging from modern carbon to fossil carbon depleted in C-14. Aerosols sampled in Chamonix and Passy (Arve River valley, French Alps) from November 2013 to August 2014 are analyzed for both radiocarbon (124 analyses in total) and levoglucosan, which is commonly used as a specific tracer for biomass burning. NOx concentration, which is expected to be associated with traffic emissions, is also monitored. Based on C-14 measurements, we can show that the relative fraction of non-fossil carbon is significantly higher in winter than in summer. In winter, non-fossil carbon represents about 85% of total carbon, while in summer this proportion is still 75% considering all samples. The largest total carbon and levoglucosan concentrations are observed for winter aerosols with values up to 50 and 8 μg m(-3), respectively. These levels are higher than those observed in many European cities, but are close to those for other polluted Alpine valleys. The non-fossil carbon concentrations are strongly correlated with the levoglucosan concentrations in winter samples, suggesting that almost all of the non-fossil carbon originates from wood combustion used for heating during winter. For summer samples, the joint use of C-14 and levoglucosan measurements leads to a new model to separately quantify the contributions of biomass burning and biogenic emissions in the non-fossil fraction. The comparison of the biogenic fraction with polyols (a proxy for primary soil biogenic emissions) and with the temperature suggests a major influence of the secondary biogenic aerosols. Significant correlations are found between the NOx concentration and the fossil carbon concentration for all seasons and sites, confirming the relation between road traffic emissions and fossil carbon. Overall, this dual approach combining radiocarbon and levoglucosan analyses strengthens the conclusion concerning the impact of biomass burning. Combining these geochemical data serves both to detect and quantify additional carbon sources. The Arve River valley provides the first illustration of aerosols of this model.

Abstract: In this study, we used an energy balance model and two simple methods based on readily available data to identify the processes driving the point-scale energy and mass balance of the snowpack. Data were provided from an experimental site located at 3200m. All models were evaluated by comparing observed and modelled snow water equivalents. Performances are variable from one season to the next and the energy balance model gives better results (mean of root mean square error, RMSE=25mm and r(2)=0.90) than the two simplified approaches (mean of RMSE=54mm and r(2)=0.70). There are significant amounts of snow sublimation but they are highly variable from season to season, depending on wind conditions (between 7 and 20% of the total). While the main source of energy for melting is net radiation, the amount of heat brought by sensible heat flux is significant for two of the most windy snow seasons.

Abstract: Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO2 remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO2 fertilization, land use, wildfire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO2 dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO2 and delta(CO2)-C-13 changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO2 dynamics from 8 ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO2 changes after 122 ka BP. This failure to simulate late-Eemian CO2 dynamics could be a result of the imposed forcings such as prescribed CaCO3 accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO2 dynamics – shallow water CaCO3 accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO2 dynamics. (C) 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Abstract: Estimates of flood event magnitudes with a certain return period are required for the design of hydraulic structures. While the return period is clearly defined in a univariate context, its definition is more challenging when the problem at hand requires considering the dependence between two or more variables in a multivariate framework. Several ways of defining a multivariate return period have been proposed in the literature, which all rely on different probability concepts. Definitions use the conditional probability, the joint probability, or can be based on the Kendall's distribution or survival function. In this study, we give a comprehensive overview on the tools that are available to define a return period in a multivariate context. We especially address engineers, practitioners, and people who are new to the topic and provide them with an accessible introduction to the topic. We outline the theoretical background that is needed when one is in a multivariate setting and present the reader with different definitions for a bivariate return period. Here, we focus on flood events and the different probability concepts are explained with a pedagogical, illustrative example of a flood event characterized by the two variables peak discharge and flood volume. The choice of the return period has an important effect on the magnitude of the design variable quantiles, which is illustrated with a case study in Switzerland. However, this choice is not arbitrary and depends on the problem at hand.

Abstract: We present a synthesis of the land-atmosphere carbon flux from land use and land cover change (LULCC) in Asia using multiple data sources and paying particular attention to deforestation and forest regrowth fluxes. The data sources are quasi-independent and include the U.N. Food andAgriculture Organization-Forest Resource Assessment (FAO-FRA2015; country-level inventory estimates), the Emission Database for Global Atmospheric Research (EDGARv4.3), the 'Houghton' bookkeeping model that incorporates FAO-FRA data, an ensemble of 8 state-of-the-art Dynamic Global Vegetation Models (DGVM), and 2 recently published independent studies using primarily remote sensing techniques. The estimates are aggregated spatially to Southeast, East, and South Asia and temporally for three decades, 1980-1989, 1990-1999 and 2000-2009. Since 1980, net carbon emissions from LULCC in Asia were responsible for 20%-40% of global LULCC emissions, with emissions from Southeast Asia alone accounting for 15%-25% of global LULCC emissions during the same period. In the 2000s and for all Asia, three estimates (FAO-FRA, DGVM, Houghton) were in agreement of a net source of carbon to the atmosphere, with mean estimates ranging between 0.24 to 0.41 Pg Cyr(-1), whereas EDGARv4.3 suggested a net carbon sink of -0.17 Pg Cyr(-1). Three of 4 estimates suggest that LULCC carbon emissions declined by at least 34% in the preceding decade (1990-2000). Spread in the estimates is due to the inclusion of different flux components and their treatments, showing the importance to include emissions from carbon rich peatlands and land management, such as shifting cultivation and wood harvesting, which appear to be consistently underreported.

Abstract: The merger of two identical surface temperature vortices is studied in the surface quasi-geostrophic model. The motivation for this study is the observation of the merger of submesoscale vortices in the ocean. Firstly, the interaction between two point vortices, in the absence or in the presence of an external deformation field, is investigated. The rotation rate of the vortices, their stationary positions and the stability of these positions are determined. Then, a numerical model provides the steady states of two finite-area, constant-temperature, vortices. Such states are less deformed than their counterparts in two-dimensional incompressible flows. Finally, numerical simulations of the nonlinear surface quasi-geostrophic equations are used to investigate the finite-time evolution of initially identical and symmetric, constant temperature vortices. The critical merger distance is obtained and the deformation of the vortices before or after merger is determined. The addition of external deformation is shown to favor or to oppose merger depending on the orientation of the vortex pair with respect to the strain axes. An explanation for this observation is proposed. Conclusions are drawn towards an application of this study to oceanic vortices.

Abstract: Since 1950, the Niger River basin has gone through three main climatic periods: a wet period (1950-1960), an extended drought (1970-1980) and since 1990 a recent partial recovery of annual rainfall. Hydrological changes co-occur with these rainfall fluctuations. In most of the basin, the rainfall deficit caused an enhanced discharge deficit, but in the Sahelian region the runoff increased despite the rainfall deficit. Since 2000 the Sahelian part of the Niger has been hit by an increase of flood hazards during the so-called red flood period. In Niamey city, the highest river levels and the longest flooded period ever recorded occurred in 2003, 2010, 2012 and 2013, with heavy casualties and property damage. The reasons for these changes, and the relative role of climate versus land use-land cover (LULC) changes are still debated and are investigated in this paper. The evolution of the Niger red flood in Niamey from 1950 to 2012 is analysed based on long-term records of rainfall (three data sets based on in situ and/or satellite data) and discharge, and a hydrological model. The model is first run with the present LULC conditions in order to analyse solely the effect of rainfall variability. The impact of LULC and drainage area modification is investigated in a second step. The simulations based on the current surface conditions are able to reproduce the observed trend in the red flood occurrence and intensity since the 1980s. This has been verified with three independent rainfall data sets and implies that rainfall variability is the main driver for the red flood intensification observed over the last 30 years. The simulation results since 1953 have revealed that LULC and drainage area changes need to be invoked to explain the changes over a 60-year period.

Abstract: Although particulate organic and elemental carbon (OC and EC) are important constituents of the suspended atmospheric particulate matter (PM), measurements of OC and EC are much less common and More uncertain than measurements of e.g. the ionic components of PM. In the framework of atmospheric research infrastructures supported by the European Union, actions have been undertaken to determine and mitigate sampling artefacts, and assess the comparability of OC and EC data obtained in a network of 10 atmospheric observatories across Europe. Positive sampling artefacts (from 0:4 to 2.8 μg C/m(3)) and analytical discrepancies (between -50% and +40% for the EC/TC ratio) have been taken into account to generate a robust data set, from which we established the phenomenology of carbonaceous aerosols at regional background sites in Europe. Across the network, TC and EC annual average concentrations range from 0.4 to 9 μg C/m(3), and from 0.1 to 2 μg C/m(3), respectively. TC/PM10 annual mean ratios range from 0.11 at a Mediterranean site to 0.34 at the most polluted continental site, and TC/PM2.5 ratios are slightly greater at all sites (0.15-0.42). EC/TC annual mean ratios range from 0.10 to 0.22, and do not depend much on PM concentration levels, especially in winter. Seasonal variations in PM and TC concentrations, and in TC/PM and EC/TC ratios, differ across the network, which can be explained by seasonal changes in PM source contributions at some sites. (C) 2016 The Authors. Published by Elsevier Ltd.

Abstract: The particle identification scheme developed by Dolan and Rutledge for X-band polarimetric radar is tested for the first time in Africa and compared with in situ measurements. The data were acquired during the Megha-Tropiques mission algorithm-validation campaign that occurred in Niger in 2010. The radar classification is compared with the in situ observations gathered by an instrumented aircraft for the 13 August 2010 squall-line case. An original approach has been developed for the radar-in situ comparison: it consists of simulating synthetic radar variables from the microphysical-probe information and comparing the two datasets in a common “radar space.” The consistency between the two types of observation is good considering the differences in sampling illustrated in the paper. The time evolution of the hydrometeor types and their relative proportion in the convective and stratiform regions are analyzed. The farther away from the convection one looks, the more aggregation dominates, riming diminishes, and hydrometeors are less dense. Particle identification based on the polarimetric radar will be applied to a 5-yr African dataset in the future.

Abstract: Grassland management type (grazed or mown) and intensity (intensive or extensive) play a crucial role in the greenhouse gas balance and surface energy budget of this biome, both at field scale and at large spatial scale. However, global gridded historical information on grassland management intensity is not available. Combining modelled grass-biomass productivity with statistics of the grass-biomass demand by livestock, we reconstruct gridded maps of grassland management intensity from 1901 to 2012. These maps include the minimum area of managed vs. maximum area of unmanaged grasslands and the fraction of mown vs. grazed area at a resolution of 0.5A degrees by 0.5A degrees. The grass-biomass demand is derived from a livestock dataset for 2000, extended to cover the period 1901-2012. The grass-biomass supply (i.e. forage grass from mown grassland and biomass grazed) is simulated by the process-based model ORCHIDEE-GM driven by historical climate change, risingaEuro-CO2 concentration, and changes in nitrogen fertilization. The global area of managed grassland obtained in this study increases from 6.1aEuro-aEuro parts per thousand x aEuro-10(6)aEuro-km(2) in 1901 to 12.3aEuro-aEuro parts per thousand x aEuro-10(6)aEuro-km(2) in 2000, although the expansion pathway varies between different regions. ORCHIDEE-GM also simulated augmentation in global mean productivity and herbage-use efficiency over managed grassland during the 20th century, indicating a general intensification of grassland management at global scale but with regional differences. The gridded grassland management intensity maps are model dependent because they depend on modelled productivity. Thus specific attention was given to the evaluation of modelled productivity against a series of observations from site-level net primary productivity (NPP) measurements to two global satellite products of gross primary productivity (GPP) (MODIS-GPP and SIF data). Generally, ORCHIDEE-GM captures the spatial pattern, seasonal cycle, and interannual variability of grassland productivity at global scale well and thus is appropriate for global applications presented here.

Abstract: This paper examines the ability of optical reflectance data assimilation to improve snow depth and snow water equivalent simulations from a chain of models with the SAFRAN meteorological model driving the detailed multilayer snowpack model Crocus now including a two-stream radiative transfer model for snow, TARTES. The direct use of reflectance data, allowed by TARTES, instead of higher level snow products, mitigates uncertainties due to commonly used retrieval algorithms.Data assimilation is performed with an ensemble-based method, the Sequential Importance Resampling Particle filter, to represent simulation uncertainties. In snowpack modeling, uncertainties of simulations are primarily assigned to meteorological forcings. Here, a method of stochastic perturbation based on an autoregressive model is implemented to explicitly simulate the consequences of these uncertainties on the snowpack estimates.Through twin experiments, the assimilation of synthetic spectral reflectances matching the MODerate resolution Imaging Spectroradiometer (MODIS) spectral bands is examined over five seasons at the Col du Lautaret, located in the French Alps. Overall, the assimilation of MODIS-like data reduces by 45aEuro-% the root mean square errors (RMSE) on snow depth and snow water equivalent. At this study site, the lack of MODIS data on cloudy days does not affect the assimilation performance significantly. The combined assimilation of MODIS-like reflectances and a few snow depth measurements throughout the 2010/2011 season further reduces RMSEs by roughly 70aEuro-%. This work suggests that the assimilation of optical reflectances has the potential to become an essential component of spatialized snowpack simulation and forecast systems. The assimilation of real MODIS data will be investigated in future works.

Abstract: Up till now, the diclofenac (DCF) transformation into its nitrogen-derivatives, N-nitroso-DCF (NO-DCF) and 5-nitro-DCF (NO2-DCF), has been mainly investigated in waste water treatment plant under nitrification or denitrification processes. This work reports, for the first time, an additional DCF microbial mediated nitrosation pathway of DCF in soil under strictly anoxic conditions probably involving codenitrification processes and fungal activities. This transformation pathway was investigated by using field observations data at a soil aquifer system (Katari watershed, Bolivia) and by carrying out soil slurry batch experiments. It was also observed for diphenylamine (DPA). Field measurements revealed the occurrence of NO-DCF, NO2-DCF and NO-DPA in groundwater samples at concentration levels in the 6-68 s/L range. These concentration levels are more significant than those previously reported in waste water treatment plant effluents taking into account dilution processes in soil. Interestingly, the p-benzoquinone imine of 5-OH-DCF was also found to be rather stable in surface water. In laboratory batch experiments under strictly anoxic conditions, the transformation of DCF and DPA into their corresponding N-nitroso derivatives was well correlated to denitrification processes. It was also observed that NO-DCF evolved into NO2-DCF while NO-DPA was stable. In vitro experiments showed that the Fisher-Hepp rearrangement could not account for NO2-DCF formation. One possible mechanism might be that NO-DCF underwent spontaneous NO loss to give the resulting intermediates diphenylaminyl radical or nitrenium cation which might evolve into NO2-DCF in presence of NO2 radical or nitrite ion, respectively. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Vortex merger is a phenomenon characterizing the whole class of geophysical vortices, from atmospheric storms and large oceanic eddies up to small scale turbulence. Here we focus on the merger of subsurface oceanic anticyclones in an idealized primitive equations model. This study has been motivated by past and recent observations of colliding lens-like anticyclones off of Gibraltar Strait. The critical conditions for merger (critical merger distance and time needed for merger) are determined. We will show that the predictions of classical two-dimensional merger are not verified for subsurface isolated vortices. For instance, critical merger distances will be reduced because of the vortex potential vorticity (PV) structure. The post-merger characteristics of the vortex (radius, extension and PV), are also determined. Merger-related effects, like production of peripheral filaments and small-scale eddies are also investigated and suggest the contribution of merger in both direct and inverse energy cascades.

Abstract: A new rheological model is developed that builds on an elasto-brittle (EB) framework used for sea ice and rock mechanics, with the intent of representing both the small elastic deformations associated with fracturing processes and the larger deformations occurring along the faults/leads once the material is highly damaged and fragmented. A viscous-like relaxation term is added to the linear-elastic constitutive law together with an effective viscosity that evolves according to the local level of damage of the material, like its elastic modulus. The coupling between the level of damage and both mechanical parameters is such that within an undamaged ice cover the viscosity is infinitely large and deformations are strictly elastic, while along highly damaged zones the elastic modulus vanishes and most of the stress is dissipated through permanent deformations. A healing mechanism is also introduced, counterbalancing the effects of damaging over large timescales. In this new model, named Maxwell-EB after the Maxwell rheology, the irreversible and reversible deformations are solved for simultaneously; hence drift velocities are defined naturally. First idealized simulations without advection show that the model reproduces the main characteristics of sea ice mechanics and deformation: strain localization, anisotropy, intermittency and associated scaling laws.

Abstract: In studies of ocean mixing, it is generally assumed that small-scale turbulent overturns lose 15%-20% of their energy in eroding the background stratification. Accumulating evidence that this energy fraction, or mixing efficiency R-f, significantly varies depending on flow properties challenges this assumption, however. Here, the authors examine the implications of a varying mixing efficiency for ocean energetics and deep-water mass transformation. Combining current parameterizations of internal wave-driven mixing with a recent model expressing R-f as a function of a turbulence intensity parameter Re-b = epsilon(nu)/nu N-2, the ratio of dissipation epsilon(nu) to stratification N-2 and molecular viscosity nu, it is shown that accounting for reduced mixing efficiencies in regions of weak stratification or energetic turbulence (high Re-b) strongly limits the ability of breaking internal waves to supply oceanic potential energy and drive abyssal upwelling. Moving from a fixed R-f = 1/6 to a variable efficiency R-f(Re-b) causes Antarctic Bottom Water upwelling induced by locally dissipating internal tides and lee waves to fall from 9 to 4 Sverdrups (Sv; 1 Sv equivalent to 10(6) m(3) s(-1)) and the corresponding potential energy source to plunge from 97 to 44 GW. When adding the contribution of remotely dissipating internal tides under idealized distributions of energy dissipation, the total rate of Antarctic Bottom Water upwelling is reduced by about a factor of 2, reaching 5-15 Sv, compared to 10-33 Sv for a fixed efficiency. The results suggest that distributed mixing, overflow-related boundary processes, and geothermal heating are more effective in consuming abyssal waters than topographically enhanced mixing by breaking internal waves. These calculations also point to the importance of accurately constraining R-f (Re-b) and including the effect in ocean models.