Abstract. In this study, a global land data assimilation system (LDAS-Monde) is applied over Europe and the Mediterranean basin to increase monitoring accuracy for land surface variables. LDAS-Monde is able to ingest information from satellite-derived surface soil moisture (SSM) and leaf area index (LAI) observations to constrain the interactions between soil–biosphere–atmosphere (ISBA, Interactions between Soil, Biosphere and Atmosphere) land surface model (LSM) coupled with the CNRM (Centre National de Recherches Météorologiques) version of the Total Runoff Integrating Pathways (ISBA-CTRIP) continental hydrological system. It makes use of the CO2-responsive version of ISBA which models leaf-scale physiological processes and plant growth. Transfer of water and heat in the soil rely on a multilayer diffusion scheme. SSM and LAI observations are assimilated using a simplified extended Kalman filter (SEKF), which uses finite differences from perturbed simulations to generate flow dependence between the observations and the model control variables. The latter include LAI and seven layers of soil (from 1 to 100 cm depth). A sensitivity test of the Jacobians over 2000–2012 exhibits effects related to both depth and season. It also suggests that observations of both LAI and SSM have an impact on the different control variables. From the assimilation of SSM, the LDAS is more effective in modifying soil moisture (SM) from the top layers of soil, as model sensitivity to SSM decreases with depth and has almost no impact from 60 cm downwards. From the assimilation of LAI, a strong impact on LAI itself is found. The LAI assimilation impact is more pronounced in SM layers that contain the highest fraction of roots (from 10 to 60 cm). The assimilation is more efficient in summer and autumn than in winter and spring. Results shows that the LDAS works well constraining the model to the observations and that stronger corrections are applied to LAI than to SM. A comprehensive evaluation of the assimilation impact is conducted using (i) agricultural statistics over France, (ii) river discharge observations, (iii) satellite-derived estimates of land evapotranspiration from the Global Land Evaporation Amsterdam Model (GLEAM) project and (iv) spatially gridded observation-based estimates of upscaled gross primary production and evapotranspiration from the FLUXNET network. Comparisons with those four datasets highlight neutral to highly positive improvement.

Abstract. Lakes and reservoirs are crucial elements of the hydrological and biochemical cycle and are a valuable resource for hydropower, domestic and industrial water use, and irrigation. Although their monitoring is crucial in times of increased pressure on water resources by both climate change and human interventions, publically available datasets of lake and reservoir levels and volumes are scarce. Within this study, a time series of variation in lake and reservoir volume between 1984 and 2015 were analysed for 137 lakes over all continents by combining the JRC Global Surface Water (GSW) dataset and the satellite altimetry database DAHITI. The GSW dataset is a highly accurate surface water dataset at 30 m resolution compromising the whole L1T Landsat 5, 7 and 8 archive, which allowed for detailed lake area calculations globally over a very long time period using Google Earth Engine. Therefore, the estimates in water volume fluctuations using the GSW dataset are expected to improve compared to current techniques as they are not constrained by complex and computationally intensive classification procedures. Lake areas and water levels were combined in a regression to derive the hypsometry relationship (dh ∕ dA) for all lakes. Nearly all lakes showed a linear regression, and 42 % of the lakes showed a strong linear relationship with a R2 > 0.8, an average R2 of 0.91 and a standard deviation of 0.05. For these lakes and for lakes with a nearly constant lake area (coefficient of variation < 0.008), volume variations were calculated. Lakes with a poor linear relationship were not considered. Reasons for low R2 values were found to be (1) a nearly constant lake area, (2) winter ice coverage and (3) a predominant lack of data within the GSW dataset for those lakes. Lake volume estimates were validated for 18 lakes in the US, Spain, Australia and Africa using in situ volume time series, and gave an excellent Pearson correlation coefficient of on average 0.97 with a standard deviation of 0.041, and a normalized RMSE of 7.42 %. These results show a high potential for measuring lake volume dynamics using a pre-classified GSW dataset, which easily allows the method to be scaled up to an extensive global volumetric dataset. This dataset will not only provide a historical lake and reservoir volume variation record, but will also help to improve our understanding of the behaviour of lakes and reservoirs and their representation in (large-scale) hydrological models.

The aim of this study is to compare the surface soil moisture (SSM) retrieved from ESA's Soil Moisture and Ocean Salinity mission (SMOS) with the output of the ORCHIDEE (ORganising Carbon and Hydrology In Dynamic EcosystEm) land surface model forced with two distinct atmospheric data sets for the period 2010 to 2012. The comparison methodology is first established over the REMEDHUS (Red de Estaciones de MEDición de la Humedad def Suelo) soil moisture measurement network, a 30 by 40. km catchment located in the central part of the Duero basin, then extended to the whole Iberian Peninsula (IP). The temporal correlation between the in-situ, remotely sensed and modelled SSM are satisfactory (r. >. 0.8). The correlation between remotely sensed and modelled SSM also holds when computed over the IP. Still, by using spectral analysis techniques, important disagreements in the effective inertia of the corresponding moisture reservoir are found. This is reflected in the spatial correlation over the IP between SMOS and ORCHIDEE SSM estimates, which is poor (¿. ~. 0.3). A single value decomposition (SVD) analysis of rainfall and SSM shows that the co-varying patterns of these variables are in reasonable agreement between both products. Moreover the first three SVD soil moisture patterns explain over 80% of the SSM variance simulated by the model while the explained fraction is only 52% of the remotely sensed values. These results suggest that the rainfall-driven soil moisture variability may not account for the poor spatial correlation between SMOS and ORCHIDEE products.Peer ReviewedPostprint (published version

Conventional wastewater treatment does not completely remove and/or inactive viruses; consequently, viruses excreted by the population can be detected in the environment. This study was undertaken to investigate the distribution and seasonality of human viruses and faecal indicator bacteria (FIB) in a river catchment located in a typical Mediterranean climate region and to discuss future trends in relation to climate change. Sample matrices included river water, untreated and treated wastewater from a wastewater treatment plant within the catchment area, and seawater from potentially impacted bathing water. Five viruses were analysed in the study. Human adenovirus (HAdV) and JC polyomavirus (JCPyV) were analysed as indicators of human faecal contamination of human pathogens; both were reported in urban wastewater (mean values of 10(6) and 10(5) GC/L, respectively), river water (10(3) and 10(2) GC/L) and seawater (10(2) and 10(1) GC/L). Human Merkel Cell polyomavirus (MCPyV), which is associated with Merkel Cell carcinoma, was detected in 75% of the raw wastewater samples (31/37) and quantified by a newly developed quantitative polymerase chain reaction (qPCR) assay with mean concentrations of 10(4) GC/L. This virus is related to skin cancer in susceptible individuals and was found in 29% and 18% of river water and seawater samples, respectively. Seasonality was only observed for norovirus genogroup II (NoV GGII), which was more abundant in cold months with levels up to 10(4) GC/L in river water. Human hepatitis E virus (HEV) was detected in 13.5% of the wastewater samples when analysed by nested PCR (nPCR). Secondary biological treatment (i.e., activated sludge) and tertiary sewage disinfection including chlorination, flocculation and UV radiation removed between 2.22 and 4.52 log10 of the viral concentrations. Climate projections for the Mediterranean climate areas and the selected river catchment estimate general warming and changes in precipitation distribution. Persistent decreases in precipitation during summer can lead to a higher presence of human viruses because river and sea water present the highest viral concentrations during warmer months. In a global context, wastewater management will be the key to preventing environmental dispersion of human faecal pathogens in future climate change scenarios.

Two data assimilation (DA) methods are compared for their ability to produce an accurate soil moisture analysis using the Météo-France land surface model: (i) SEKF, a simplified extended Kalman filter, which uses a climatological background-error covariance, and (ii) EnSRF, the ensemble square root filter, which uses an ensemble background-error covariance and approximates random rainfall errors stochastically. In situ soil moisture observations at 5 cm depth are assimilated into the surface layer and 30 cm deep observations are used to evaluate the root-zone analysis on 12 sites in south-western France (SMOSMANIA network). These sites differ in terms of climate and soil texture. The two methods perform similarly and improve on the open loop. Both methods suffer from incorrect linear assumptions which are particularly degrading to the analysis during water-stressed conditions: the EnSRF by a dry bias and the SEKF by an oversensitivity of the model Jacobian between the surface and the root-zone layers. These problems are less severe for the sites with wetter climates. A simple bias correction technique is tested on the EnSRF. Although this reduces the bias, it modifies the soil moisture fluxes and suppresses the ensemble spread, which degrades the analysis performance. However, the EnSRF flow-dependent background-error covariance evidently captures seasonal variability in the soil moisture errors and should exploit planned improvements in the model physics.Synthetic twin experiments demonstrate that when there is only a random component in the precipitation forcing errors, the correct stochastic representation of these errors enables the EnSRF to perform better than the SEKF. It might therefore be possible for the EnSRF to perform better than the SEKF with real data, if the rainfall uncertainty was ac-curately captured. However, the simple rainfall error model is not advantageous in our real experiments. More realistic rainfall error models are suggested.Published by Copernicus Publications on behalf of the European Geosciences Union.

No abstract

[1] We define a Markov-modulated autoregressive model with exogenous input (MARX) to generate runoff scenarios using climatic information. Runoff parameterization is assumed to be conditioned on a hidden climate state following a Markov chain, where state transition probabilities are functions of the climatic input. MARX allows stochastic modeling of nonstationary runoff, as runoff anomalies are described by a mixture of autoregressive models with exogenous input, each one corresponding to a climate state. We apply MARX to inflow time series of the Daule Peripa reservoir (Ecuador). El Niño-Southern Oscillation (ENSO) information is used to condition runoff parameterization. Among the investigated ENSO indexes, the NINO 1+2 sea surface temperature anomalies and the trans-Niño index perform best as predictors. In the perspective of reservoir optimization at various time scales, MARX produces realistic long-term scenarios and short-term forecasts, especially when intense El Niño events occur. Low predictive ability is found for negative runoff anomalies, as no climatic index correlating properly with negative inflow anomalies has yet been identified.

Abstract. Climate change is anticipated to alter the demand and supply of water at the earth's surface. Since many societal impacts from a lack of water happen under drought conditions, it is important to understand how droughts may develop with climate change. This study shows how hydrological droughts will change across Europe with increasing global warming levels (GWL of 1.5, 2 and 3 K above preindustrial temperature). We employ a low-flow index derived from river discharge simulations of a spatially-distributed physically-based hydrological and water use model, which was forced with a large ensemble of regional climate model projections under a high emissions (RCP8.5) and moderate mitigation (RCP4.5) pathway. Different traits of drought, including severity, duration and frequency, were investigated. The projected changes in these treats identify four main sub-regions in Europe that are characterized by somehow homogeneous and distinct behaviours with a clear southwest/northeast contrast. The Mediterranean and Boreal sub-regions of Europe show strong, but opposite, changes at all three GWLs, with the former area mostly interested by stronger droughts (with larger differences at 3 K) while the latter sees a reduction in droughts. In the Atlantic and Continental sub-regions the changes are less marked and characterized by a larger uncertainty, especially at the 1.5 and 2 K GWLs. Combining the projections in drought hazard with population and agricultural information shows that with 3 K global warming an additional 11 million people and 4.5 million ha of agricultural land will be exposed to droughts every year, on average. These are mostly located in the Mediterranean and Atlantic regions of Europe.