ALBERT

Description: Regional and local patterns in depth to water table, hydrochemistry and peat properties of bogs and their laggs in coastal British Columbia Hydrology and Earth System Sciences, 17, 3421-3435, 2013 Author(s): S. A. Howie and H. J. van Meerveld In restoration planning for damaged raised bogs, the lagg at the bog margin is often not given considerable weight and is sometimes disregarded entirely. However, the lagg is critical for the proper functioning of the bog, as it supports the water mound in the bog. In order to include the lagg in a restoration plan for a raised bog, it is necessary to understand the hydrological characteristics and functions of this rarely studied transition zone. We studied 13 coastal British Columbia (BC) bogs and identified two different gradients in depth to water table, hydrochemistry and peat properties: (1) a local bog expanse–bog margin gradient, and (2) a regional gradient related to climate and proximity to the ocean. Depth to water table generally increased across the transition from bog expanse to bog margin. In the bog expanse, pH was above 4.2 in the Pacific Oceanic wetland region (cooler and wetter climate) and below 4.3 in the Pacific Temperate wetland region (warmer and drier climate). Both pH and pH-corrected electrical conductivity increased significantly across the transition from bog expanse to bog margin, though not in all cases. Na + and Mg 2+ concentrations were generally highest in exposed, oceanic bogs and lower in inland bogs. Ash content in peat samples increased across the bog expanse–bog margin transition, and appears to be a useful abiotic indicator of the location of the bog margin. The observed variation in the hydrological and hydrochemical gradients across the bog expanse–bog margin transition highlights both local and regional diversity of bogs and their associated laggs.

Description: Bayesian networks modelling in support to cross-cutting analysis of water supply and sanitation in developing countries Hydrology and Earth System Sciences, 17, 3397-3419, 2013 Author(s): C. Dondeynaz, J. López Puga, and C. Carmona Moreno Despite the efforts made towards the Millennium Development Goals targets during the last decade, improved access to water supply or basic sanitation still remains unavailable for millions of people across the world. This paper proposes a set of models that use 25 key variables and country profiles from the WatSan4Dev data set involving water supply and sanitation (Dondeynaz et al., 2012). This paper suggests the use of Bayesian network modelling methods because they are more easily adapted to deal with non-normal distributions, and integrate a qualitative approach for data analysis. They also offer the advantage of integrating preliminary knowledge into the probabilistic models. The statistical performance of the proposed models ranges between 20 and 5% error rates, which are very satisfactory taking into account the strong heterogeneity of variables. Probabilistic scenarios run from the models allow an assessment of the relationships between human development, external support, governance aspects, economic activities and water supply and sanitation (WSS) access. According to models proposed in this paper, gaining a strong poverty reduction will require the WSS access to reach 75–76% through: (1) the management of ongoing urbanisation processes to avoid slums development; and (2) the improvement of health care, for instance for children. Improving governance, such as institutional efficiency, capacities to make and apply rules, or control of corruption is positively associated with WSS sustainable development. The first condition for an increment of the HDP (human development and poverty) remains of course an improvement of the economic conditions with higher household incomes. Moreover, a significant country commitment to the environment, associated with civil society freedom of expression constitutes a favourable setting for sustainable WSS services delivery. Intensive agriculture using irrigation practises also appears as a mean for sustainable WSS thanks to multi-uses and complementarities. With a WSS sector organised at national level, irrigation practices can support the structuring and efficiency of the agriculture sector. It may then induce rural development in areas where WSS access often is set back compared to urban areas 1 . External financial support, called Official Development Assistance (ODA CI), plays a role in WSS improvement but comes last in the sensitivity analyses of models. An overall 47% of the Official Development Assistance goes first to poor countries, and is associated to governance aspects: (1) political stability and (2) country commitment to the environment and civil society degree of freedom. These governance aspects constitute a good framework for aid implementation in recipient countries. Modelling is run with the five groups of countries as defined in Dondeynaz et al. (2012). Models for profile 4 (essential external support) and profile 5 (primary material consumption) are specifically detailed and analysed in this paper. For countries in profile 4, fighting against water scarcity and progressing desertification should be the priority. However, for countries in profile 5, efforts should first concentrate on consolidation of political stability while supporting diversification of the economic activities. Nevertheless, for both profiles, reduction of poverty should remain the first priority as previously indicated. 1 JMP statistics, 2004 http://www.wssinfo.org/data-estimates/table/ , last access: 22 July 2013.

Description: Technical Note: A comparison of model and empirical measures of catchment-scale effective energy and mass transfer Hydrology and Earth System Sciences, 17, 3389-3395, 2013 Author(s): C. Rasmussen and E. L. Gallo Recent work suggests that a coupled effective energy and mass transfer (EEMT) term, which includes the energy associated with effective precipitation and primary production, may serve as a robust prediction parameter of critical zone structure and function. However, the models used to estimate EEMT have been solely based on long-term climatological data with little validation using direct empirical measures of energy, water, and carbon balances. Here we compare catchment-scale EEMT estimates generated using two distinct approaches: (1) EEMT modeled using the established methodology based on estimates of monthly effective precipitation and net primary production derived from climatological data, and (2) empirical catchment-scale EEMT estimated using data from 86 catchments of the Model Parameter Estimation Experiment (MOPEX) and MOD17A3 annual net primary production (NPP) product derived from Moderate Resolution Imaging Spectroradiometer (MODIS). Results indicated positive and significant linear correspondence ( R 2 = 0.75; P 〈 0.001) between model and empirical measures with an average root mean square error (RMSE) of 4.86 MJ m −2 yr −1 . Modeled EEMT values were consistently greater than empirical measures of EEMT. Empirical catchment estimates of the energy associated with effective precipitation ( E PPT ) were calculated using a mass balance approach that accounts for water losses to quick surface runoff not accounted for in the climatologically modeled E PPT . Similarly, local controls on primary production such as solar radiation and nutrient limitation were not explicitly included in the climatologically based estimates of energy associated with primary production ( E BIO ), whereas these were captured in the remotely sensed MODIS NPP data. These differences likely explain the greater estimate of modeled EEMT relative to the empirical measures. There was significant positive correlation between catchment aridity and the fraction of EEMT partitioned into E BIO ( F BIO ), with an increase in F BIO as a fraction of the total as aridity increases and percentage of catchment woody plant cover decreases. In summary, the data indicated strong correspondence between model and empirical measures of EEMT with limited bias that agree well with other empirical measures of catchment energy and water partitioning and plant cover.

Description: Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model Geoscientific Model Development Discussions, 6, 4563-4601, 2013 Author(s): F. Cossu and K. Hocke This study examines how different microphysical parameterization schemes influence orographically-induced precipitation and the distributions of hydrometeors and water vapour for mid-latitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m −2 in the downwelling shortwave radiation and up to 33 W m −2 in the downwelling longwave radiation.

Description: Contribution of directly connected and isolated impervious areas to urban drainage network hydrographs Hydrology and Earth System Sciences, 17, 3473-3483, 2013 Author(s): Y. Seo, N.-J. Choi, and A. R. Schmidt This paper addresses the mass balance error observed in runoff hydrographs in urban watersheds by introducing assumptions regarding the contribution of infiltrated rainfall from pervious areas and isolated impervious area (IIA) to the runoff hydrograph. Rainfall infiltrating into pervious areas has been assumed not to contribute to the runoff hydrograph until Hortonian excess rainfall occurs. However, mass balance analysis in an urban watershed indicates that rainfall infiltrated to pervious areas can contribute directly to the runoff hydrograph, thereby offering an explanation for the long hydrograph tail commonly observed in runoff from urban storm sewers. In this study, a hydrologic analysis based on the width function is introduced, with two types of width functions obtained from both pervious and impervious areas, respectively. The width function can be regarded as the direct interpretation of the network response. These two width functions are derived to obtain distinct response functions for directly connected impervious areas (DCIA), IIA, and pervious areas. The results show significant improvement in the estimation of runoff hydrographs and suggest the need to consider the flow contribution from pervious areas to the runoff hydrograph. It also implies that additional contribution from flow paths through joints and cracks in sewer pipes needs to be taken into account to improve the estimation of runoff hydrographs in urban catchments.

Description: An original interpretation of the wet edge of the surface temperature–albedo space to estimate crop evapotranspiration (SEB-1S), and its validation over an irrigated area in northwestern Mexico Hydrology and Earth System Sciences, 17, 3623-3637, 2013 Author(s): O. Merlin The space defined by the pair surface temperature ( T ) and surface albedo (α), and the space defined by the pair T and fractional green vegetation cover ( f vg ) have been extensively used to estimate evaporative fraction (EF) from solar/thermal remote sensing data. In both space-based approaches, evapotranspiration (ET) is estimated as remotely sensed EF times the available energy. For a given data point in the T -α space or in the T - f vg space, EF is derived as the ratio of the distance separating the point from the line identified as the dry edge to the distance separating the dry edge and the line identified as the wet edge. The dry and wet edges are classically defined as the upper and lower limit of the spaces, respectively. When investigating side by side the T -α and the T - f vg spaces, one observes that the range covered by T values on the (classically determined) wet edge is different for both spaces. In addition, when extending the wet and dry lines of the T -α space, both lines cross at α ≈ 0.4 although the wet and dry edges of the T - f vg space never cross for 0 ≤ f vg 〈 1. In this paper, a new ET (EF) model (SEB-1S) is derived by revisiting the classical physical interpretation of the T -α space to make its wet edge consistent with that of the T - f vg space. SEB-1S is tested over a 16 km by 10 km irrigated area in northwestern Mexico during the 2007–2008 agricultural season. The classical T -α space-based model is implemented as benchmark to evaluate the performance of SEB-1S. Input data are composed of ASTER (Advanced Spaceborne Thermal Emission and Reflection radiometer) thermal infrared, Formosat-2 shortwave, and station-based meteorological data. The fluxes simulated by SEB-1S and the classical T -α space-based model are compared on seven ASTER overpass dates with the in situ measurements collected at six locations within the study domain. The ET simulated by SEB-1S is significantly more accurate and robust than that predicted by the classical T -α space-based model. The correlation coefficient and slope of the linear regression between simulated and observed ET is improved from 0.82 to 0.93, and from 0.63 to 0.90, respectively. Moreover, constraining the wet edge using air temperature data improves the slope of the linear regression between simulated and observed ET.

Description: On an improved sub-regional water resources management representation for integration into earth system models Hydrology and Earth System Sciences, 17, 3605-3622, 2013 Author(s): N. Voisin, H. Li, D. Ward, M. Huang, M. Wigmosta, and L. R. Leung Human influence on the hydrologic cycle includes regulation and storage, consumptive use and overall redistribution of water resources in space and time. Representing these processes is essential for applications of earth system models in hydrologic and climate predictions, as well as impact studies at regional to global scales. Emerging large-scale research reservoir models use generic operating rules that are flexible for coupling with earth system models. Those generic operating rules have been successful in reproducing the overall regulated flow at large basin scales. This study investigates the uncertainties of the reservoir models from different implementations of the generic operating rules using the complex multi-objective Columbia River Regulation System in northwestern United States as an example to understand their effects on not only regulated flow but also reservoir storage and fraction of the demand that is met. Numerical experiments are designed to test new generic operating rules that combine storage and releases targets for multi-purpose reservoirs and to compare the use of reservoir usage priorities and predictors (withdrawals vs. consumptive demands, as well as natural vs. regulated mean flow) for configuring operating rules. Overall the best performing implementation is with combined priorities rules (flood control storage targets and irrigation release targets) set up with mean annual natural flow and mean monthly withdrawals. The options of not accounting for groundwater withdrawals, or on the contrary, of assuming that all remaining demand is met through groundwater extractions, are discussed.

Description: On selection of the optimal data time interval for real-time hydrological forecasting Hydrology and Earth System Sciences, 17, 3639-3659, 2013 Author(s): J. Liu and D. Han With the advancement in modern telemetry and communication technologies, hydrological data can be collected with an increasingly higher sampling rate. An important issue deserving attention from the hydrological community is which suitable time interval of the model input data should be chosen in hydrological forecasting. Such a problem has long been recognised in the control engineering community but is a largely ignored topic in operational applications of hydrological forecasting. In this study, the intrinsic properties of rainfall–runoff data with different time intervals are first investigated from the perspectives of the sampling theorem and the information loss using the discrete wavelet transform tool. It is found that rainfall signals with very high sampling rates may not always improve the accuracy of rainfall–runoff modelling due to the catchment low-pass-filtering effect. To further investigate the impact of a data time interval in real-time forecasting, a real-time forecasting system is constructed by incorporating the probability distributed model (PDM) with a real-time updating scheme, the autoregressive moving-average (ARMA) model. Case studies are then carried out on four UK catchments with different concentration times for real-time flow forecasting using data with different time intervals of 15, 30, 45, 60, 90 and 120 min. A positive relation is found between the forecast lead time and the optimal choice of the data time interval, which is also highly dependent on the catchment concentration time. Finally, based on the conclusions from the case studies, a hypothetical pattern is proposed in three-dimensional coordinates to describe the general impact of the data time interval and to provide implications of the selection of the optimal time interval in real-time hydrological forecasting. Although nowadays most operational hydrological systems still have low data sampling rates (daily or hourly), the future is that higher sampling rates will become more widespread, and there is an urgent need for hydrologists both in academia and in the field to realise the significance of the data time interval issue. It is important that more case studies in different catchments with various hydrological forecasting systems are explored in the future to further verify and improve the proposed hypothetical pattern.

Description: Benchmark products for land evapotranspiration: LandFlux-EVAL multi-data set synthesis Hydrology and Earth System Sciences, 17, 3707-3720, 2013 Author(s): B. Mueller, M. Hirschi, C. Jimenez, P. Ciais, P. A. Dirmeyer, A. J. Dolman, J. B. Fisher, M. Jung, F. Ludwig, F. Maignan, D. G. Miralles, M. F. McCabe, M. Reichstein, J. Sheffield, K. Wang, E. F. Wood, Y. Zhang, and S. I. Seneviratne Land evapotranspiration (ET) estimates are available from several global data sets. Here, monthly global land ET synthesis products, merged from these individual data sets over the time periods 1989–1995 (7 yr) and 1989–2005 (17 yr), are presented. The merged synthesis products over the shorter period are based on a total of 40 distinct data sets while those over the longer period are based on a total of 14 data sets. In the individual data sets, ET is derived from satellite and/or in situ observations (diagnostic data sets) or calculated via land-surface models (LSMs) driven with observations-based forcing or output from atmospheric reanalyses. Statistics for four merged synthesis products are provided, one including all data sets and three including only data sets from one category each (diagnostic, LSMs, and reanalyses). The multi-annual variations of ET in the merged synthesis products display realistic responses. They are also consistent with previous findings of a global increase in ET between 1989 and 1997 (0.13 mm yr −2 in our merged product) followed by a significant decrease in this trend (−0.18 mm yr −2 ), although these trends are relatively small compared to the uncertainty of absolute ET values. The global mean ET from the merged synthesis products (based on all data sets) is 493 mm yr −1 (1.35 mm d −1 ) for both the 1989–1995 and 1989–2005 products, which is relatively low compared to previously published estimates. We estimate global runoff (precipitation minus ET) to 263 mm yr −1 (34 406 km 3 yr −1 ) for a total land area of 130 922 000 km 2 . Precipitation, being an important driving factor and input to most simulated ET data sets, presents uncertainties between single data sets as large as those in the ET estimates. In order to reduce uncertainties in current ET products, improving the accuracy of the input variables, especially precipitation, as well as the parameterizations of ET, are crucial.

Description: Spatial patterns in timing of the diurnal temperature cycle Hydrology and Earth System Sciences, 17, 3695-3706, 2013 Author(s): T. R. H. Holmes, W. T. Crow, and C. Hain This paper investigates the structural difference in timing of the diurnal temperature cycle (DTC) over land resulting from choice of measuring device or model framework. It is shown that the timing can be reliably estimated from temporally sparse observations acquired from a constellation of low Earth-orbiting satellites given record lengths of at least three months. Based on a year of data, the spatial patterns of mean DTC timing are compared between temperature estimates from microwave Ka-band, geostationary thermal infrared (TIR), and numerical weather prediction model output from the Global Modeling and Assimilation Office (GMAO). It is found that the spatial patterns can be explained by vegetation effects, sensing depth differences and more speculatively the orientation of orographic relief features. In absolute terms, the GMAO model puts the peak of the DTC on average at 12:50 local solar time, 23 min before TIR with a peak temperature at 13:13 (both averaged over Africa and Europe). Since TIR is the shallowest observation of the land surface, this small difference represents a structural error that possibly affects the model's ability to assimilate observations that are closely tied to the DTC. The equivalent average timing for Ka-band is 13:44, which is influenced by the effect of increased sensing depth in desert areas. For non-desert areas, the Ka-band observations lag the TIR observations by only 15 min, which is in agreement with their respective theoretical sensing depth. The results of this comparison provide insights into the structural differences between temperature measurements and models, and can be used as a first step to account for these differences in a coherent way.

Description: Integrated hydrological modeling of the North China Plain and implications for sustainable water management Hydrology and Earth System Sciences, 17, 3759-3778, 2013 Author(s): H. Qin, G. Cao, M. Kristensen, J. C. Refsgaard, M. O. Rasmussen, X. He, J. Liu, Y. Shu, and C. Zheng Groundwater overdraft has caused fast water level decline in the North China Plain (NCP) since the 1980s. Although many hydrological models have been developed for the NCP in the past few decades, most of them deal only with the groundwater component or only at local scales. In the present study, a coupled surface water–groundwater model using the MIKE SHE code has been developed for the entire alluvial plain of the NCP. All the major processes in the land phase of the hydrological cycle are considered in the integrated modeling approach. The most important parameters of the model are first identified by a sensitivity analysis process and then calibrated for the period 2000–2005. The calibrated model is validated for the period 2006–2008 against daily observations of groundwater heads. The simulation results compare well with the observations where acceptable values of root mean square error (RMSE) (most values lie below 4 m) and correlation coefficient ( R ) (0.36–0.97) are obtained. The simulated evapotranspiration (ET) is then compared with the remote sensing (RS)-based ET data to further validate the model simulation. The comparison result with a R 2 value of 0.93 between the monthly averaged values of simulated actual evapotranspiration (AET) and RS AET for the entire NCP shows a good performance of the model. The water balance results indicate that more than 70% of water leaving the flow system is attributed to the ET component, of which about 0.25% is taken from the saturated zone (SZ); about 29% comes from pumping, including irrigation pumping and non-irrigation pumping (net pumping). Sustainable water management analysis of the NCP is conducted using the simulation results obtained from the integrated model. An effective approach to improve water use efficiency in the NCP is by reducing the actual ET, e.g. by introducing water-saving technologies and changes in cropping.

Description: Simulation of a persistent medium-term precipitation event over the western Iberian Peninsula Hydrology and Earth System Sciences, 17, 3741-3758, 2013 Author(s): S. C. Pereira, A. C. Carvalho, J. Ferreira, J. P. Nunes, J. J. Keizer, and A. Rocha This study evaluated the performance of the WRF-ARW (Weather Research and Forecasting with Advanced Research) weather prediction model in simulating the spatial and temporal patterns of an extreme rainfall period over a complex orographic region in north-central Portugal. The analysis was performed during the rainy season and, more specifically, the month of December 2009. In this period, the region of interest was under the influence of a sequential passage of low-pressure systems associated with frontal surfaces. These synoptic weather patterns were responsible for long periods of rainfall, resulting in a high monthly precipitation. The WRF model results during the study period were furthermore evaluated with the specific objective to complement gaps in the precipitation recordings of a reference meteorological station (located in Pousadas), the data of which are fundamental for hydrological studies in nearby experimental catchments. Three distinct WRF model runs were forced with initial fields and boundary conditions obtained from a global domain model: (1) a reference experiment with no nudging (RunRef); (2) observational nudging for a specific location, i.e. the above-mentioned Pousadas reference station (RunObsN); and (3) nudging to the analysed field (RunGridN). Model performance was evaluated, using several statistical parameters, against a dataset of 27 rainfall stations that were grouped by elevation. The three model runs had similar performances, even though RunGridN resulted in a slight improvement. Regarding the other two experiments, this improvement justifies its use for complementing the surface measurements at the Pousadas reference station. Overall model accuracy, expressed in root mean square error (RMSE), of the three runs was comparable for the stations of the different elevations classes. Even so, it was slightly better for stations in the lowlands than the highlands. Furthermore, model predictions tended to be less accurate for stations located in rough terrain and deep valleys.

Description: Investigating the spatio-temporal variability in groundwater and surface water interactions: a multi-technique approach Hydrology and Earth System Sciences, 17, 3437-3453, 2013 Author(s): N. P. Unland, I. Cartwright, M. S. Andersen, G. C. Rau, J. Reed, B. S. Gilfedder, A. P. Atkinson, and H. Hofmann The interaction between groundwater and surface water along the Tambo and Nicholson rivers, southeast Australia, was investigated using 222 Rn, Cl, differential flow gauging, head gradients, electrical conductivity (EC) and temperature profiles. Head gradients, temperature profiles, Cl concentrations and 222 Rn activities all indicate higher groundwater fluxes to the Tambo River in areas of increased topographic variation where the potential to form large groundwater–surface water gradients is greater. Groundwater discharge to the Tambo River calculated by Cl mass balance was significantly lower (1.48 × 10 4 to 1.41 × 10 3 m 3 day −1 ) than discharge estimated by 222 Rn mass balance (5.35 × 10 5 to 9.56 × 10 3 m 3 day −1 ) and differential flow gauging (5.41 × 10 5 to 6.30 × 10 3 m 3 day −1 ) due to bank return waters. While groundwater sampling from the bank of the Tambo River was intended to account for changes in groundwater chemistry associated with bank infiltration, variations in bank infiltration between sample sites remain unaccounted for, limiting the use of Cl as an effective tracer. Groundwater discharge to both the Tambo and Nicholson rivers was the highest under high-flow conditions in the days to weeks following significant rainfall, indicating that the rivers are well connected to a groundwater system that is responsive to rainfall. Groundwater constituted the lowest proportion of river discharge during times of increased rainfall that followed dry periods, while groundwater constituted the highest proportion of river discharge under baseflow conditions (21.4% of the Tambo in April 2010 and 18.9% of the Nicholson in September 2010).

Description: Resolving structural errors in a spatially distributed hydrologic model using ensemble Kalman filter state updates Hydrology and Earth System Sciences, 17, 3455-3472, 2013 Author(s): J. H. Spaaks and W. Bouten In hydrological modeling, model structures are developed in an iterative cycle as more and different types of measurements become available and our understanding of the hillslope or watershed improves. However, with increasing complexity of the model, it becomes more and more difficult to detect which parts of the model are deficient, or which processes should also be incorporated into the model during the next development step. In this study, we first compare two methods (the Shuffled Complex Evolution Metropolis algorithm (SCEM-UA) and the Simultaneous parameter Optimization and Data Assimilation algorithm (SODA)) to calibrate a purposely deficient 3-D hillslope-scale model to error-free, artificially generated measurements. We use a multi-objective approach based on distributed pressure head at the soil–bedrock interface and hillslope-scale discharge and water balance. For these idealized circumstances, SODA's usefulness as a diagnostic methodology is demonstrated by its ability to identify the timing and location of processes that are missing in the model. We show that SODA's state updates provide information that could readily be incorporated into an improved model structure, and that this type of information cannot be gained from parameter estimation methods such as SCEM-UA. We then expand on the SODA result by performing yet another calibration, in which we investigate whether SODA's state updating patterns are still capable of providing insight into model structure deficiencies when there are fewer measurements, which are moreover subject to measurement noise. We conclude that SODA can help guide the discussion between experimentalists and modelers by providing accurate and detailed information on how to improve spatially distributed hydrologic models.

Description: AnaWEGE: a weather generator based on analogues of atmospheric circulation Geoscientific Model Development Discussions, 6, 4745-4774, 2013 Author(s): P. Yiou This paper presents a stochastic weather generator based on analogues of circulation (AnaWEGE). Analogues of circulation have been a promising paradigm to analyse climate variability and its extremes. The weather generator uses precomputed analogues of sea-level pressure over the North Atlantic. The stochastic rules of the generator constrain the continuity in time of the simulations. The generator then simulates spatially coherent time series of a climate variable, drawn from meteorological observations. The weather generator is tested for European temperatures, and for winter and summer seasons. The biases in temperature quantiles and autocorrelation are rather small compared to observed variability. The ability of simulating extremely hot summers and cold winters is also assessed.

Description: Methodological aspects of a pattern-scaling approach to produce global fields of monthly means of daily maximum and minimum temperature Geoscientific Model Development Discussions, 6, 4833-4882, 2013 Author(s): S. Kremser, G. E. Bodeker, and J. Lewis A Climate Pattern-Scaling Model (CPSM) that simulates global patterns of climate change, for a prescribed emissions scenario, is described. A CPSM works by quantitatively establishing the statistical relationship between a climate variable at a specific location (e.g. daily maximum surface temperature, T max ) and one or more predictor time series (e.g. global mean surface temperature, T global ) – referred to as the "training" of the CPSM. This training uses a regression model to derive fit-coefficients that describe the statistical relationship between the predictor time series and the target climate variable time series. Once that relationship has been determined, and given the predictor time series for any greenhouse gas (GHG) emissions scenario, the change in the climate variable of interest can be reconstructed – referred to as the "application" of the CPSM. The advantage of using a CPSM rather than a typical atmosphere-ocean global climate model (AOGCM) is that the predictor time series required by the CPSM can usually be generated quickly using a simple climate model (SCM) for any prescribed GHG emissions scenario and then applied to generate global fields of the climate variable of interest. The training can be performed either on historical measurements or on output from an AOGCM. Using model output from 21st century simulations has the advantage that the climate change signal is more pronounced than in historical data and therefore a more robust statistical relationship is obtained. The disadvantage of using AOGCM output is that the CPSM training might be compromised by any AOGCM inadequacies. For the purposes of exploring the various methodological aspects of the CPSM approach, AOGCM output was used in this study to train the CPSM. These investigations of the CPSM methodology focus on monthly mean fields of daily temperature extremes ( T max and T min ). Key conclusions are: (1) overall, the CPSM trained on simulations based on the Representative Concentration Pathway (RCP) 8.5 emissions scenario is able to reproduce AOGCM simulations of T max and T min based on predictor time series from an RCP 4.5 emissions scenario; (2) access to hemisphere average land and ocean temperatures as predictors improves the variance that can be explained, particularly over the oceans; (3) regression model fit-coefficients derived from individual simulations based on the RCP 2.6, 4.5 and 8.5 emissions scenarios agree well over most regions of the globe (the Arctic is the exception); (4) training the CPSM on concatenated time series from an ensemble of simulations does not result in fit-coefficients that explain significantly more of the variance than an approach that weights results based on single simulation fits; and (5) the inclusion of a linear time dependence in the regression model fit-coefficients improves the variance explained, primarily over the oceans.

Description: Modeling different freeze/thaw processes in heterogeneous landscapes of the Arctic polygonal tundra using an ecosystem model Geoscientific Model Development Discussions, 6, 4883-4932, 2013 Author(s): given_name prefix surname suffix, S. Yi, K. Wischnewski, M. Langer, S. Muster, and J. Boike Freeze/thaw (F/T) processes can be quite different under the various land surface types found in the heterogeneous polygonal tundra of the Arctic. Proper simulation of these different processes is essential for accurate prediction of the release of greenhouse gases under a warming climate scenario. In this study we have modified the dynamic organic soil version of the Terrestrial Ecosystem Model (DOS-TEM) to simulate F/T processes beneath the polygon rims, polygon centers (with and without water), and lakes that are common features in Arctic lowland regions. We first verified the F/T algorithm in the DOS-TEM against analytical solutions, and then compared the results with in situ measurements from Samoylov Island, Siberia. In the final stage, we examined the different responses of the F/T processes for different water levels at the various land surface types. The simulations revealed that (1) the DOS-TEM was very efficient and its results compared very well with analytical solutions for idealized cases, (2) the simulations compared reasonably well with in situ measurements although there were a number of model limitations and uncertainties, (3) the DOS-TEM was able to successfully simulate the differences in F/T dynamics under different land surface types, and (4) permafrost beneath water bodies was found to respond highly sensitive to changes in water depths between 1 and 2 m. Our results indicate that water is very important in the thermal processes simulated by the DOS-TEM; the heterogeneous nature of the landscape and different water depths therefore need to be taken into account when simulating methane emission responses to a warming climate.

Description: A standard test case suite for two-dimensional linear transport on the sphere: results from a collection of state-of-the-art schemes Geoscientific Model Development Discussions, 6, 4983-5076, 2013 Author(s): P. H. Lauritzen, P. A. Ullrich, C. Jablonowski, P. A. Bosler, D. Calhoun, A. J. Conley, T. Enomoto, L. Dong, S. Dubey, O. Guba, A. B. Hansen, E. Kaas, J. Kent, J.-F. Lamarque, M. J. Prather, D. Reinert, V. V. Shashkin, W. C. Skamarock, B. Sørensen, M. A. Taylor, and M. A. Tolstykh Recently, a standard test case suite for 2-D linear transport on the sphere was proposed to assess important aspects of accuracy in geophysical fluid dynamics with a "minimal" set of idealized model configurations/runs/diagnostics. Here we present results from 19 state-of-the-art transport scheme formulations based on finite-difference/finite-volume methods as well as emerging (in the context of atmospheric/oceanographic sciences) Galerkin methods. Discretization grids range from traditional regular latitude-longitude grids to more isotropic domain discretizations such as icosahedral and cubed-sphere tessellations of the sphere. The schemes are evaluated using a wide range of diagnostics in idealized flow environments. Accuracy is assessed in single- and two-tracer configurations using conventional error norms as well as novel diagnostics designed for climate and climate-chemistry applications. In addition, algorithmic considerations that may be important for computational efficiency are reported on. The latter is inevitably computing platform dependent, The ensemble of results from a wide variety of schemes presented here helps shed light on the ability of the test case suite diagnostics and flow settings to discriminate between algorithms and provide insights into accuracy in the context of global atmospheric/ocean modeling. A library of benchmark results is provided to facilitate scheme intercomparison and model development. Simple software and data-sets are made available to facilitate the process of model evaluation and scheme intercomparison.

Description: Spatial and seasonal variations in evapotranspiration over Canada's landmass Hydrology and Earth System Sciences, 17, 3561-3575, 2013 Author(s): S. Wang, Y. Yang, Y. Luo, and A. Rivera A 30 yr (1979–2008) dataset of actual evapotranspiration (ET) at 1 km resolution was generated over Canada's landmass by integrating remote sensing land surface data and gridded climate data using the EALCO model run at a 30 min time step. This long-term high-resolution dataset was used to characterize the spatiotemporal variations in ET across Canada. The results show that annual ET varied from 600 mm yr −1 over several regions in the south to less than 100 mm yr −1 in the northern Arctic. Nationally, ET in summer (i.e., June to August) comprised 65% of the annual total amount. ET in the cold season remained mostly below 10 mm month −1 over the country. Negative monthly ET was obtained over the Arctic region in winter, indicating EALCO simulated a larger amount of condensation than ET. Overall, the mean ET over the entire Canadian landmass for the 30 yr was 239 mm yr −1 , or 44% of its corresponding precipitation. Comparisons of available ET studies in Canada revealed large uncertainties in ET estimates associated with using different approaches. The scarcity of ET measurements for the diverse ecosystems in Canada remains a significant challenge for reducing the uncertainties; this gap needs to be addressed in future studies to improve capabilities in climate/weather modeling and water resource management.

Description: On the analytic approximation of bulk collision rates of non-spherical hydrometeors Geoscientific Model Development Discussions, 6, 5077-5116, 2013 Author(s): A. Seifert, U. Blahak, and R. Buhr Analytic approximations of the binary collision rates of hydrometeors are derived for use in bulk microphysical parameterizations. Special attention is given to non-spherical hydrometeors like raindrops and snowflakes. The terminal fall velocity of these particles cannot be sufficiently well approximated by power law relations which are used in most microphysical parameterizations and therefore an improved formulation is needed. The analytic approximations of the bulk collision rates given in this paper are an alternative to look-up tables and can replace the Wisner approximation which is used in many atmospheric models.

Description: Post-processing rainfall forecasts from numerical weather prediction models for short-term streamflow forecasting Hydrology and Earth System Sciences, 17, 3587-3603, 2013 Author(s): D. E. Robertson, D. L. Shrestha, and Q. J. Wang Sub-daily ensemble rainfall forecasts that are bias free and reliably quantify forecast uncertainty are critical for flood and short-term ensemble streamflow forecasting. Post-processing of rainfall predictions from numerical weather prediction models is typically required to provide rainfall forecasts with these properties. In this paper, a new approach to generate ensemble rainfall forecasts by post-processing raw numerical weather prediction (NWP) rainfall predictions is introduced. The approach uses a simplified version of the Bayesian joint probability modelling approach to produce forecast probability distributions for individual locations and forecast lead times. Ensemble forecasts with appropriate spatial and temporal correlations are then generated by linking samples from the forecast probability distributions using the Schaake shuffle. The new approach is evaluated by applying it to post-process predictions from the ACCESS-R numerical weather prediction model at rain gauge locations in the Ovens catchment in southern Australia. The joint distribution of NWP predicted and observed rainfall is shown to be well described by the assumed log-sinh transformed bivariate normal distribution. Ensemble forecasts produced using the approach are shown to be more skilful than the raw NWP predictions both for individual forecast lead times and for cumulative totals throughout all forecast lead times. Skill increases result from the correction of not only the mean bias, but also biases conditional on the magnitude of the NWP rainfall prediction. The post-processed forecast ensembles are demonstrated to successfully discriminate between events and non-events for both small and large rainfall occurrences, and reliably quantify the forecast uncertainty. Future work will assess the efficacy of the post-processing method for a wider range of climatic conditions and also investigate the benefits of using post-processed rainfall forecasts for flood and short-term streamflow forecasting.

Description: Comparison of climate change signals in CMIP3 and CMIP5 multi-model ensembles and implications for Central Asian glaciers Hydrology and Earth System Sciences, 17, 3661-3677, 2013 Author(s): A. F. Lutz, W. W. Immerzeel, A. Gobiet, F. Pellicciotti, and M. F. P. Bierkens Central Asian water resources largely depend on melt water generated in the Pamir and Tien Shan mountain ranges. To estimate future water availability in this region, it is necessary to use climate projections to estimate the future glacier extent and volume. In this study, we evaluate the impact of uncertainty in climate change projections on the future glacier extent in the Amu and Syr Darya river basins. To this end we use the latest climate change projections generated for the upcoming IPCC report (CMIP5) and, for comparison, projections used in the fourth IPCC assessment (CMIP3). With these projections we force a regionalized glacier mass balance model, and estimate changes in the basins' glacier extent as a function of the glacier size distribution in the basins and projected temperature and precipitation. This glacier mass balance model is specifically developed for implementation in large scale hydrological models, where the spatial resolution does not allow for simulating individual glaciers and data scarcity is an issue. Although the CMIP5 ensemble results in greater regional warming than the CMIP3 ensemble and the range in projections for temperature as well as precipitation is wider for the CMIP5 than for the CMIP3, the spread in projections of future glacier extent in Central Asia is similar for both ensembles. This is because differences in temperature rise are small during periods of maximum melt (July–September) while differences in precipitation change are small during the period of maximum accumulation (October–February). However, the model uncertainty due to parameter uncertainty is high, and has roughly the same importance as uncertainty in the climate projections. Uncertainty about the size of the decline in glacier extent remains large, making estimates of future Central Asian glacier evolution and downstream water availability uncertain.

Description: Prediction of dissolved reactive phosphorus losses from small agricultural catchments: calibration and validation of a parsimonious model Hydrology and Earth System Sciences, 17, 3679-3693, 2013 Author(s): C. Hahn, V. Prasuhn, C. Stamm, P. Lazzarotto, M. W. H. Evangelou, and R. Schulin Eutrophication of surface waters due to diffuse phosphorus (P) losses continues to be a severe water quality problem worldwide, causing the loss of ecosystem functions of the respective water bodies. Phosphorus in runoff often originates from a small fraction of a catchment only. Targeting mitigation measures to these critical source areas (CSAs) is expected to be most efficient and cost-effective, but requires suitable tools. Here we investigated the capability of the parsimonious Rainfall-Runoff-Phosphorus (RRP) model to identify CSAs in grassland-dominated catchments based on readily available soil and topographic data. After simultaneous calibration on runoff data from four small hilly catchments on the Swiss Plateau, the model was validated on a different catchment in the same region without further calibration. The RRP model adequately simulated the discharge and dissolved reactive P (DRP) export from the validation catchment. Sensitivity analysis showed that the model predictions were robust with respect to the classification of soils into "poorly drained" and "well drained", based on the available soil map. Comparing spatial hydrological model predictions with field data from the validation catchment provided further evidence that the assumptions underlying the model are valid and that the model adequately accounts for the dominant P export processes in the target region. Thus, the parsimonious RRP model is a valuable tool that can be used to determine CSAs. Despite the considerable predictive uncertainty regarding the spatial extent of CSAs, the RRP can provide guidance for the implementation of mitigation measures. The model helps to identify those parts of a catchment where high DRP losses are expected or can be excluded with high confidence. Legacy P was predicted to be the dominant source for DRP losses and thus, in combination with hydrologic active areas, a high risk for water quality.

Description: Ensemble initialization of the oceanic component of a coupled model through bred vectors at seasonal-to-interannual time scales Geoscientific Model Development Discussions, 6, 5189-5214, 2013 Author(s): J. Baehr and R. Piontek We evaluate the ensemble spread at seasonal-to-interannual time scales for two perturbation techniques implemented into the ocean component of a coupled model: (1) lagged initial conditions as commonly used for decadal predictions, (2) bred vectors as commonly used for weather and seasonal forecasting. We show that relative to an uninitialized reference simulation the implementation for bred vectors can improve the ensemble spread compared to lagged initialization at time scales from one months up to three years. As bred vectors have so far mostly been used at short time scales, we initially focus on the implementation of the bred vectors into the ocean component. We introduce a depth-dependent vertical rescaling norm, accounting for the vertical dependence of the variability, and extending the commonly used upper-ocean rescaling norm to the full water column. We further show that it is sufficient for the (sub-surface) ocean to breed temperature and salinity (i.e., scalar quantities), and rely on the governing physics to carry the temperature and salinity perturbations to the flow field. Using these bred vectors with a rescaling interval of 12 months, we initialize hindcast simulations and compare them to hindcast simulations initialized with lagged initial conditions. We quantify the ensemble spread by analyzing Talagrand diagrams and spread-error ratios. For both temperature and salinity, the lagged initialized ensemble is particularly under-dispersive for the first few months of predictable lead time. The ensemble initialized with bred vectors improves the spread for temperature and salinity for the 0–700 m and 1000–3500 m means, compared to the lagged ensemble at lead times of several months to one year. As the lead time increases to years, the differences between the two ensemble initialization techniques becomes more difficult to discern. While the results need to be confirmed in an initialized framework, the present analysis represents a first step towards an improved ensemble generation at the transition from seasonal-to-interannual time scales, in particular at lead times up to one year.

Description: Understanding the performance of the FLake model over the African Great Lakes Geoscientific Model Development Discussions, 6, 5141-5187, 2013 Author(s): W. Thiery, A. Martynov, F. Darchambeau, J.-P. Descy, P.-D. Plisnier, L. Sushama, and N. P. M. van Lipzig The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu, Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding Automatic Weather Stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration allowed to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and temperature structure. At Lake Kivu, mixolimnion temperatures predicted by FLake were found sensitive both to minimal variations in the external parameters (lake depth and water transparency) as to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may already lead to a regime switch from the correctly represented seasonal mixed layer deepening to either completely mixed (down to the model lake bottom) or permanently stratified (from ~10 m downwards) conditions. In contrast, model temperatures are found robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterize tropical lake water surface temperatures within atmospheric prediction models, but may be less appropriate, in its current form, to study complex limnological processes within tropical lakes. Furthermore, a study of different initial conditions showed that for tropical lakes lacking reliable initial data, a fully mixed, artificially warm initialisation is to be preferred, but only if the model is allowed to spin up until convergence is reached. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60 m during the main dry season is primarily due to enhanced lake evaporation and secondarily due to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Description: High-resolution Med-CORDEX regional climate model simulations for hydrological impact studies: a first evaluation of the ALADIN-Climate model in Morocco Hydrology and Earth System Sciences, 17, 3721-3739, 2013 Author(s): Y. Tramblay, D. Ruelland, S. Somot, R. Bouaicha, and E. Servat In the framework of the international CORDEX program, new regional climate model (RCM) simulations at high spatial resolutions are becoming available for the Mediterranean region (Med-CORDEX initiative). This study provides the first evaluation for hydrological impact studies of one of these high-resolution simulations in a 1800 km 2 catchment located in North Morocco. Different approaches are compared to analyze the climate change impacts on the hydrology of this catchment using a high-resolution RCM (ALADIN-Climate) from the Med-CORDEX initiative at two different spatial resolutions (50 and 12 km) and for two different Radiative Concentration Pathway scenarios (RCP4.5 and RCP8.5). The main issues addressed in the present study are: (i) what is the impact of increased RCM resolution on present-climate hydrological simulations and on future projections? (ii) Are the bias-correction of the RCM model and the parameters of the hydrological model stationary and transferable to different climatic conditions? (iii) What is the climate and hydrological change signal based on the new Radiative Concentration Pathways scenarios (RCP4.5 and RCP8.5)? Results indicate that high resolution simulations at 12 km better reproduce the seasonal patterns, the seasonal distributions and the extreme events of precipitation. The parameters of the hydrological model, calibrated to reproduce runoff at the monthly time step over the 1984–2010 period, do not show a strong variability between dry and wet calibration periods in a differential split-sample test. However the bias correction of precipitation by quantile-matching does not give satisfactory results in validation using the same differential split-sample testing method. Therefore a quantile-perturbation method that does not rely on any stationarity assumption and produces ensembles of perturbed series of precipitation was introduced. The climate change signal under scenarios 4.5 and 8.5 indicates a decrease of respectively −30 to −57% in surface runoff for the mid-term (2041–2062), when for the same period the projections for precipitation are ranging between −15 and −19% and for temperature between +1.3 and +1.9 °C.

Description: Errors in climate model daily precipitation and temperature output: time invariance and implications for bias correction Hydrology and Earth System Sciences, 17, 2147-2159, 2013 Author(s): E. P. Maurer, T. Das, and D. R. Cayan When correcting for biases in general circulation model (GCM) output, for example when statistically downscaling for regional and local impacts studies, a common assumption is that the GCM biases can be characterized by comparing model simulations and observations for a historical period. We demonstrate some complications in this assumption, with GCM biases varying between mean and extreme values and for different sets of historical years. Daily precipitation and maximum and minimum temperature from late 20th century simulations by four GCMs over the United States were compared to gridded observations. Using random years from the historical record we select a "base" set and a 10 yr independent "projected" set. We compare differences in biases between these sets at median and extreme percentiles. On average a base set with as few as 4 randomly-selected years is often adequate to characterize the biases in daily GCM precipitation and temperature, at both median and extreme values; 12 yr provided higher confidence that bias correction would be successful. This suggests that some of the GCM bias is time invariant. When characterizing bias with a set of consecutive years, the set must be long enough to accommodate regional low frequency variability, since the bias also exhibits this variability. Newer climate models included in the Intergovernmental Panel on Climate Change fifth assessment will allow extending this study for a longer observational period and to finer scales.

Description: Characterization of groundwater dynamics in landslides in varved clays Hydrology and Earth System Sciences, 17, 2171-2183, 2013 Author(s): J. E. van der Spek, T. A. Bogaard, and M. Bakker Groundwater dynamics may play a significant role in landslides. A detailed model is developed of the groundwater dynamics in landslides in varved clays in the Trièves area in the French Alps. The varved clays consist of a sequence of alternating silt and clay layers, covered by a colluvium layer and intersected by fissures. The hydraulic conductivity of the clay layers is negligible compared to the silt layers. It is conceptualized that fissures form a hydraulic connection between the colluvium and the varved clays. Groundwater recharge flows through the colluvium into the fissures, where water is exchanged horizontally between the fissure and the silt layers of the varved clays. Groundwater flow in the colluvium is simulated with the Boussinesq equation, while flow in the silt layers of the varved clays is simulated with the Richards equation. Longitudinal outflow from the fissure is simulated with a linear-reservoir model. Scattered data of relatively short monitoring periods is available for several landslides in the region. A good similarity between observed and simulated heads is obtained, especially when considering the lack of important physical parameters such as the fissure width and the distance between the monitoring point and the fissure. A simulation for the period 1959–2004 showed some correlation between peaks in the simulated heads and the recorded occurrence of landslides, while the bottom of the varved clays remained saturated during the entire simulation period.

Description: Assessing the hydrological effect of the check dams in the Loess Plateau, China, by model simulations Hydrology and Earth System Sciences, 17, 2185-2193, 2013 Author(s): Y. D. Xu, B. J. Fu, and C. S. He Check dams are commonly used for soil conservation. In the Loess Plateau of China, check dams have been widely constructed as the principal means to retain floodwater and intercept soil sediments since the 1970s. For instance, there are more than 6572 check dams in the Yanhe watershed with an area of 7725 km 2 in the Loess Plateau. However, little research has been done to quantify the hydrological effects of the check dams. In this research, the SWAT model (Soil and Water Assessment Tool) was applied to simulate the runoff and sediment in the Yanhe watershed. We treated the 1950s to 1960s as the reference period since there were very few check dams during the period. The model was firstly calibrated and validated in the reference period. The calibrated model was then used in the later periods to simulate the hydrological effects of the check dams. The results showed that the check dams had a regulation effect on runoff and a retention effect on sediment. From 1984 to 1987, the runoff in rainy season (from May to October) decreased by 1.54 m 3 s −1 (14.7%) to 3.13 m 3 s −1 (25.9%) due to the check dams; while in dry season (from November to the following April), runoff increased by 1.46 m 3 s −1 (60.5%) to 1.95 m 3 s −1 (101.2%); the sediment in rainy season decreased by 2.49 × 10 6 ton (34.6%) to 4.35 × 10 6 ton (48.0%). From 2006 to 2008, the runoff in rainy season decreased by 0.79 m 3 s −1 (15.5%) to 1.75 m 3 s −1 (28.9%), and the runoff in dry season increased by 0.51 m 3 s −1 (20.1%) to 0.97 m 3 s −1 (46.4%); the sediment in rainy season decreased by 2.03 × 10 6 ton (79.4%) to 3.12 × 10 6 ton (85.5%). Construction of the large number of check dams in the Loess Plateau has enhanced the region's capacity to control the runoff and sediment. In the Yanhe watershed, the annual runoff was reduced by less than 14.3% due to the check dams; and the sediment in rainy season was blocked by up to 85.5%. Thus, check dams are effective measures for soil erosion control in the Loess Plateau.

Description: Determining spatial variability of dry spells: a Markov-based method, applied to the Makanya catchment, Tanzania Hydrology and Earth System Sciences, 17, 2161-2170, 2013 Author(s): B. M. C. Fischer, M. L. Mul, and H. H. G. Savenije With a growing world population and a trend towards more resource-intensive diets, pressure on land and water resources for food production will continue to increase in the coming decades. Large parts of the world rely on rainfed agriculture for their food security. In Africa, 90% of the food production is from rainfed agriculture, generally with low yields and a high risk of crop failure. One of the main reasons for crop failure is the occurrence of dry spells during the growing season. Key indicators are the critical dry spell duration and the probability of dry spell occurrence. In this paper a new Markov-based framework is presented to spatially map the length of dry spells for fixed probabilities of non-exceedance. The framework makes use of spatially varying Markov coefficients that are correlated to readily available spatial information such as elevation and distance to the sea. The dry spell map thus obtained is compared to the spatially variable critical dry spell duration, based on soil properties and crop water requirements, to assess the probability of crop failure in different locations. The results show that in the Makanya catchment the length of dry spell occurrence is highly variable in space, even over relatively short distances. In certain areas the probability of crop failure reaches levels that make rainfed agricultural unsustainable, even close to areas where currently rainfed agriculture is successfully being practised. This method can be used to identify regions that are vulnerable to dry spells and, subsequently, to develop strategies for supplementary irrigation or rainwater harvesting.

Description: Estimating actual, potential, reference crop and pan evaporation using standard meteorological data: a pragmatic synthesis Hydrology and Earth System Sciences, 17, 1331-1363, 2013 Author(s): T. A. McMahon, M. C. Peel, L. Lowe, R. Srikanthan, and T. R. McVicar This guide to estimating daily and monthly actual, potential, reference crop and pan evaporation covers topics that are of interest to researchers, consulting hydrologists and practicing engineers. Topics include estimating actual evaporation from deep lakes and from farm dams and for catchment water balance studies, estimating potential evaporation as input to rainfall-runoff models, and reference crop evapotranspiration for small irrigation areas, and for irrigation within large irrigation districts. Inspiration for this guide arose in response to the authors' experiences in reviewing research papers and consulting reports where estimation of the actual evaporation component in catchment and water balance studies was often inadequately handled. Practical guides using consistent terminology that cover both theory and practice are not readily available. Here we provide such a guide, which is divided into three parts. The first part provides background theory and an outline of the conceptual models of potential evaporation of Penman, Penman–Monteith and Priestley–Taylor, as well as discussions of reference crop evapotranspiration and Class-A pan evaporation. The last two sub-sections in this first part include techniques to estimate actual evaporation from (i) open-surface water and (ii) landscapes and catchments (Morton and the advection-aridity models). The second part addresses topics confronting a practicing hydrologist, e.g. estimating actual evaporation for deep lakes, shallow lakes and farm dams, lakes covered with vegetation, catchments, irrigation areas and bare soil. The third part addresses six related issues: (i) automatic (hard wired) calculation of evaporation estimates in commercial weather stations, (ii) evaporation estimates without wind data, (iii) at-site meteorological data, (iv) dealing with evaporation in a climate change environment, (v) 24 h versus day-light hour estimation of meteorological variables, and (vi) uncertainty in evaporation estimates. This paper is supported by a Supplement that includes 21 sections enhancing the material in the text, worked examples of many procedures discussed in the paper, a program listing (Fortran 90) of Morton's WREVAP evaporation models along with tables of monthly Class-A pan coefficients for 68 locations across Australia and other information.

Description: Corrigendum to "Calibration of aerodynamic roughness over the Tibetan Plateau with Ensemble Kalman Filter analysed heat flux" published in Hydrol. Earth Syst. Sci., 16, 4291–4302, 2012 Hydrology and Earth System Sciences, 17, 1309-1310, 2013 Author(s): J. H. Lee, J. Timmermans, Z. Su, and M. Mancini No abstract available.

Description: Structural break or long memory: an empirical survey on daily rainfall data sets across Malaysia Hydrology and Earth System Sciences, 17, 1311-1318, 2013 Author(s): F. Yusof, I. L. Kane, and Z. Yusop A short memory process that encounters occasional structural breaks in mean can show a slower rate of decay in the autocorrelation function and other properties of fractional integrated I (d) processes. In this paper we employed a procedure for estimating the fractional differencing parameter in semiparametric contexts proposed by Geweke and Porter-Hudak (1983) to analyse nine daily rainfall data sets across Malaysia. The results indicate that all the data sets exhibit long memory. Furthermore, an empirical fluctuation process using the ordinary least square (OLS)-based cumulative sum (CUSUM) test for the break date was applied. Break dates were detected in all data sets. The data sets were partitioned according to their respective break date, and a further test for long memory was applied for all subseries. Results show that all subseries follows the same pattern as the original series. The estimate of the fractional parameters d 1 and d 2 on the subseries obtained by splitting the original series at the break date confirms that there is a long memory in the data generating process (DGP). Therefore this evidence shows a true long memory not due to structural break.

Description: Snow glacier melt estimation in tropical Andean glaciers using artificial neural networks Hydrology and Earth System Sciences, 17, 1265-1280, 2013 Author(s): V. Moya Quiroga, A. Mano, Y. Asaoka, S. Kure, K. Udo, and J. Mendoza Snow and glacier melt (SGM) estimation plays an important role in water resources management. Although melting process can be modelled by energy balance methods, such studies require detailed data, which is rarely available. Hence, new and simpler approaches are needed for SGM estimations. The present study aims at developing an artificial neural networks (ANN) based technique for estimating the energy available for melt (EAM) and SGM rates using available and easy to obtain data such as temperature, short-wave radiation and relative humidity. Several ANN and multiple linear regression models (MLR) were developed to represent the energy fluxes and estimate the EAM. The models were trained using measured data from the Zongo glacier located in the outer tropics and validated against measured data from the Antizana glacier located in the inner tropics. It was found that ANN models provide a better generalisation when applied to other data sets. The performance of the models was improved by including Antizana data into the training set, as it was proved to provide better results than other techniques like the use of a prior logarithmic transformation. The final model was validated against measured data from the Alpine glaciers Argentière and Saint-Sorlin. Then, the models were applied for the estimation of SGM at Condoriri glacier. The estimated SGM was compared with SGM estimated by an enhanced temperature method and proved to have the same behaviour considering temperature sensibility. Moreover, the ANN models have the advantage of direct application, while the temperature method requires calibration of empirical coefficients.

Description: Assessment of shallow subsurface characterisation with non-invasive geophysical methods at the intermediate hill-slope scale Hydrology and Earth System Sciences, 17, 1297-1307, 2013 Author(s): S. Popp, D. Altdorff, and P. Dietrich Hill-slopes of several hectares in size represent a difficult scale for subsurface characterisation, as these landscape units are well beyond the scope of traditional point-scale techniques. By means of electromagnetic induction (EMI) and gamma-ray spectroscopy, spatially distributed soil proxy data were collected from a heterogeneous hill-slope site. Results of EMI mapping using the EM38DD showed that soil electrical conductivity (ECa) is highly variable at both temporal and spatial scales. Calibration of the integral ECa signal to a specific target like soil moisture is hampered by the ambiguous response of EMI to the clay-rich hill-slope underground. Gamma-ray results were obtained during a single survey, along with EMI measurements and selected soil sampling. In contrast to ECa, a noticeable correlation between Total Count and K emission data and soil-water content seemed to be present. Relevant proxy variables from both methods were used for k means clustering in order to distinguish between hill-slope areas with different soil conditions. As a result, we obtained a suitable partition of hill-slope that was comparable with a previously obtained zonation model based on ecological factors.

Description: Multivariate return periods in hydrology: a critical and practical review focusing on synthetic design hydrograph estimation Hydrology and Earth System Sciences, 17, 1281-1296, 2013 Author(s): B. Gräler, M. J. van den Berg, S. Vandenberghe, A. Petroselli, S. Grimaldi, B. De Baets, and N. E. C. Verhoest Most of the hydrological and hydraulic studies refer to the notion of a return period to quantify design variables. When dealing with multiple design variables, the well-known univariate statistical analysis is no longer satisfactory, and several issues challenge the practitioner. How should one incorporate the dependence between variables? How should a multivariate return period be defined and applied in order to yield a proper design event? In this study an overview of the state of the art for estimating multivariate design events is given and the different approaches are compared. The construction of multivariate distribution functions is done through the use of copulas, given their practicality in multivariate frequency analyses and their ability to model numerous types of dependence structures in a flexible way. A synthetic case study is used to generate a large data set of simulated discharges that is used for illustrating the effect of different modelling choices on the design events. Based on different uni- and multivariate approaches, the design hydrograph characteristics of a 3-D phenomenon composed of annual maximum peak discharge, its volume, and duration are derived. These approaches are based on regression analysis, bivariate conditional distributions, bivariate joint distributions and Kendall distribution functions, highlighting theoretical and practical issues of multivariate frequency analysis. Also an ensemble-based approach is presented. For a given design return period, the approach chosen clearly affects the calculated design event, and much attention should be given to the choice of the approach used as this depends on the real-world problem at hand.

Description: The mid-Pliocene climate simulated by FGOALS-g2 Geoscientific Model Development Discussions, 6, 2403-2428, 2013 Author(s): W. Zheng, Z. Zhang, L. Chen, and Y. Yu Within the framework of Pliocene Model Intercomparison Project (PlioMIP), the mid-Pliocene (3.264–3.025 Ma) climate simulated by the Flexible Global Ocean-Atmosphere-Land System model grid-point version 2 (FGOALS-g2) are analyzed in this study. Results show that the model reproduces the large-scale features of the global warming over the land and ocean. The simulated mid-Pliocene global annual mean surface air temperature (TAS) and sea surface temperature (SST) are 4.17 and 2.62°C warmer than the pre-Industrial simulation, respectively. In particular, the feature of larger warming over mid-high latitudes is well captured. In the simulated warm mid-Pliocene climate, the Atlantic Meridional Overturning Circulation (AMOC) and El Niño-Southern Oscillation (ENSO) become weaker.

Description: The potential of an observational data set for calibration of a computationally expensive computer model Geoscientific Model Development Discussions, 6, 2369-2401, 2013 Author(s): D. J. McNeall, P. G. Challenor, J. R. Gattiker, and E. J. Stone We measure the potential of an observational data set to constrain a set of inputs to a complex and computationally expensive computer model. We use each member in turn of an ensemble of output from a computationally expensive model, corresponding to some observable part of a modelled system, as a proxy for an observational data set. We argue that our ability to constrain uncertain parameter inputs to a model using its own output as data, provides a maximum bound for our ability to constrain the model inputs using observations of the real system. The ensemble provides a set of known parameter input and model output pairs, which we use to build a computationally efficient statistical proxy for the full computer model, termed an emulator. We use the emulator to find and rule out ''implausible" values for the inputs of held-out ensemble members, given the computer model output. As we know the true values of the inputs for the ensemble, we can compare our constraint of the model inputs with the true value of the input for any ensemble member. Measures of the quality of constraint have the potential to inform strategy for data collection campaigns, before any real-world data is collected, as well as acting as an effective sensitivity analysis. We use an ensemble of the ice sheet model Glimmer to demonstrate our measures of quality of constraint. The ensemble has 250 model runs with 5 uncertain input parameters, and an output variable representing the pattern of the thickness of ice over Greenland. We have an observation of historical ice sheet thickness that directly matches the output variable, and offers an opportunity to constrain the model. We show that different ways of summarising our output variable (ice volume, ice surface area and maximum ice thickness) offer different potential constraints on individual input parameters. We show that combining the observational data gives increased power to constrain the model. We investigate the impact of uncertainty in observations or in model biases on our measures, showing that even a modest uncertainty can seriously degrade the potential of the observational data to constrain the model.

Description: Derivation of RCM-driven potential evapotranspiration for hydrological climate change impact analysis in Great Britain: a comparison of methods and associated uncertainty in future projections Hydrology and Earth System Sciences, 17, 1365-1377, 2013 Author(s): C. Prudhomme and J. Williamson Potential evapotranspiration (PET) is the water that would be lost by plants through evaporation and transpiration if water was not limited in the soil, and it is commonly used in conceptual hydrological modelling in the calculation of runoff production and hence river discharge. Future changes of PET are likely to be as important as changes in precipitation patterns in determining changes in river flows. However PET is not calculated routinely by climate models so it must be derived independently when the impact of climate change on river flow is to be assessed. This paper compares PET estimates from 12 equations of different complexity, driven by the Hadley Centre's HadRM3-Q0 model outputs representative of 1961–1990, with MORECS PET, a product used as reference PET in Great Britain. The results show that the FAO56 version of the Penman–Monteith equations reproduces best the spatial and seasonal variability of MORECS PET across GB when driven by HadRM3-Q0 estimates of relative humidity, total cloud, wind speed and linearly bias-corrected mean surface temperature. This suggests that potential biases in HadRM3-Q0 climate do not result in significant biases when the physically based FAO56 equations are used. Percentage changes in PET between the 1961–1990 and 2041–2070 time slices were also calculated for each of the 12 PET equations from HadRM3-Q0. Results show a large variation in the magnitude (and sometimes direction) of changes estimated from different PET equations, with Turc, Jensen–Haise and calibrated Blaney–Criddle methods systematically projecting the largest increases across GB for all months and Priestley–Taylor, Makkink, and Thornthwaite showing the smallest changes. We recommend the use of the FAO56 equation as, when driven by HadRM3-Q0 climate data, this best reproduces the reference MORECS PET across Great Britain for the reference period of 1961–1990. Further, the future changes of PET estimated by FAO56 are within the range of uncertainty defined by the ensemble of 12 PET equations. The changes show a clear northwest–southeast gradient of PET increase with largest (smallest) changes in the northwest in January (July and October) respectively. However, the range in magnitude of PET changes due to the choice of PET method shown in this study for Great Britain suggests that PET uncertainty is a challenge facing the assessment of climate change impact on hydrology mostly ignored up to now.

Description: On searching for optimized set of physical parameterization schemes in a multi-physics land surface process model Geoscientific Model Development Discussions, 6, 4511-4530, 2013 Author(s): S. Hong, X. Yu, S. K. Park, Y.-S. Choi, and B. Myoung Optimization of land surface models has been very challenging due to the increasing complexity of such models. Typical parameter calibration techniques often limit the solution of the spatiotemporal discrepancy in the modeling performance levels especially for regional applications. Thus, in this study, an attempt was made to perform scheme-based model optimization by designing a framework for coupling a micro-genetic algorithm (micro-GA) with the Noah land surface model that has multiple physics options (Noah-MP). Micro-GA controls the scheme selections in 10 different land surface parameterization fields in Noah-MP in order to extract the optimal scheme combination for a certain region. This coupling framework was successfully applied to the optimization of the surface water partitioning in the Korean Peninsula, promising not only the effectiveness of the scheme-based optimization but also model diagnosis capability by exploring the scheme sensitivity during the micro-GA evolution process. Then, the method was applied to four different regions in East Asia that have different climatic characteristics. The results indicate that (1) the optimal scheme combinations vary with the regions, (2) schemes related to the surface water partitioning are important for the modeling accuracy, and (3) specialized post-parameter optimization for each region may be required.

Description: Atmospheric inverse modeling with known physical bounds: an example from trace gas emissions Geoscientific Model Development Discussions, 6, 4531-4562, 2013 Author(s): S. M. Miller, A. M. Michalak, and P. J. Levi Many inverse problems in the atmospheric sciences involve parameters with known physical constraints. Examples include non-negativity (e.g., emissions of some urban air pollutants) or upward limits implied by reaction or solubility constants. However, probabilistic inverse modeling approaches based on Gaussian assumptions cannot incorporate such bounds and thus often produce unrealistic results. The atmospheric literature lacks consensus on the best means to overcome this problem, and existing atmospheric studies rely on a limited number of the possible methods with little examination of the relative merits of each. This paper investigates the applicability of several approaches to bounded inverse problems and is also the first application of Markov chain Monte Carlo (MCMC) to estimation of atmospheric trace gas fluxes. The approaches discussed here are broadly applicable. A common method of data transformations is found to unrealistically skew estimates for the examined example application. The method of Lagrange multipliers and two MCMC methods yield more realistic and accurate results. In general, the examined MCMC approaches produce the most realistic result but can require substantial computational time. Lagrange multipliers offer an appealing alternative for large, computationally intensive problems when exact uncertainty bounds are less central to the analysis. A synthetic data inversion of US anthropogenic methane emissions illustrates the strengths and weaknesses of each approach.

Description: A database and tool for boundary conditions for regional air quality modeling: description and evaluation Geoscientific Model Development Discussions, 6, 4665-4704, 2013 Author(s): B. H. Henderson, F. Akhtar, H. O. T. Pye, S. L. Napelenok, and W. T. Hutzell Transported air pollutants receive increasing attention as regulations tighten and global concentrations increase. The need to represent international transport in regional air quality assessments requires improved representation of boundary concentrations. Currently available observations are too sparse vertically to provide boundary information, particularly for ozone precursors, but global simulations can be used to generate spatially and temporally varying Lateral Boundary Conditions (LBC). This study presents a public database of global simulations designed and evaluated for use as LBC for air quality models (AQMs). The database covers the contiguous United States (CONUS) for the years 2000–2010 and contains hourly varying concentrations of ozone, aerosols, and their precursors. The database is complimented by a tool for configuring the global results as inputs to regional scale models (e.g., Community Multiscale Air Quality or Comprehensive Air quality Model with extensions). This study also presents an example application based on the CONUS domain, which is evaluated against satellite retrieved ozone vertical profiles. The results show performance is largely within uncertainty estimates for the Tropospheric Emission Spectrometer (TES) with some exceptions. The major difference shows a high bias in the upper troposphere along the southern boundary in January. This publication documents the global simulation database, the tool for conversion to LBC, and the fidelity of concentrations on the boundaries. This documentation is intended to support applications that require representation of long-range transport of air pollutants.

Description: Effects of vegetation structure on biomass accumulation in a Balanced Optimality Structure Vegetation Model (BOSVM v1.0) Geoscientific Model Development Discussions, 6, 4603-4663, 2013 Author(s): Z. Yin, S. C. Dekker, B. J. J. M. van den Hurk, and H. A. Dijkstra A myriad of interactions exist between vegetation and local climate for arid and semi-arid regions. Vegetation function, structure and individual behavior have large impacts on carbon-water-energy balances, which consequently influence local climate variability that, in turn, feeds back to the vegetation. In this study, a conceptual vegetation structure scheme is formulated and tested in a new carbon-water-energy coupled model to explore the importance of vegetation structure and vegetation adaptation to water stress on equilibrium biomass states. Surface energy, water and carbon fluxes are simulated for a range of vegetation structures across a precipitation gradient in West Africa and optimal vegetation structures that maximizes biomass for each precipitation regime are determined. Two different strategies of vegetation adaptation to water stress are included. Under dry conditions vegetation tries to maximize the Water Use Efficiency and Leaf Area Index as it tries to maximize carbon gain. However, an important negative feedback mechanism is found as the vegetation also tries to minimize its cover to optimize the surrounding bare ground area from which water can be extracted, thereby forming patches of vertical vegetation. Under larger precipitation, a positive feedback mechanism is found in which vegetation tries to maximize its cover as it then can reduce water loss from bare soil while having maximum carbon gain due to a large Leaf Area Index. The competition between vegetation and bare soil determines a transition between a "survival" state to a "growing" state.

Description: A distributed computing approach to improve the performance of the Parallel Ocean Program (v2.1) Geoscientific Model Development Discussions, 6, 4705-4744, 2013 Author(s): B. van Werkhoven, J. Maassen, M. Kliphuis, H. A. Dijkstra, S. E. Brunnabend, M. van Meersbergen, F. J. Seinstra, and H. E. Bal The Parallel Ocean Program (POP) is used in many strongly eddying ocean circulation simulations. Ideally one would like to do thousand-year long simulations, but the current performance of POP prohibits this type of simulations. In this work, using a new distributed computing approach, two innovations to improve the performance of POP are presented. The first is a new block partitioning scheme for the optimization of the load balancing of POP such that it can be run efficiently in a multi-platform setting. The second is an implementation of part of the POP model code on Graphics Processing Units. We show that the combination of both innovations leads to a substantial performance increase also when running POP simultaneously over multiple computational platforms.

Description: Impact of long-term drainage on summer groundwater flow patterns in the Mer Bleue peatland, Ontario, Canada Hydrology and Earth System Sciences, 17, 3485-3498, 2013 Author(s): B. J. Kopp, J. H. Fleckenstein, N. T. Roulet, E. Humphreys, J. Talbot, and C. Blodau Long-term impacts of drier conditions on the hydrology of northern peatlands are poorly understood. We used long-term drainage near a historic drainage ditch, separating an area from the main peatland, as an analogue for long-term drying in a northern temperate bog. The objective was to identify the impact of drier conditions on ecohydrological processes and groundwater flow patterns in an area now forested and an area that maintained a bog-like character. Groundwater flow patterns alternated between mostly downward flow and occasionally upward flow in the bog area and were mostly upward-orientated in the forested area, which suggested that there the flow pattern had shifted from bog- to fen-like conditions. Flow patterns were in agreement with changes in post-drainage hydraulic conductivities, storage capacity of the peat and water table levels. Compared to the bog, hydraulic conductivities in the forested area were one to three orders of magnitude lower in the uppermost 0.75 m of peat (paired t test, p 〈 0.05). Bulk density had increased and the water table level was lower and more strongly fluctuating in the forested area. Our findings suggest hydraulic gradients and flow patterns have changed due to increased evapotranspiration and interception with the emergence of a tree cover. The smaller size of the now-forested area relative to the remaining bog area appeared to be important for the hydrological change. With the main Mer Bleue bog as hinterland, enhanced runoff to the drainage channel had little effect on hydrologic and vegetation patterns. In the cut-off, smaller, now forested area pervasive changes in vegetation and hydrologic processes occurred. The difference in response to local drainage raises questions about tipping points with respect to the impact of drying on peatland ecosystems that need to be addressed in future research.

Description: Inherently mass-conservative version of the semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model dynamical core Geoscientific Model Development Discussions, 6, 4809-4832, 2013 Author(s): V. V. Shashkin and M. A. Tolstykh The semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model is the global semi-Lagrangian hydrostatic model used for operational medium-range and seasonal forecasts at Hydrometeorological centre of Russia. The distinct feature of SL-AV dynamical core is the semi-implicit semi-Lagrangian vorticity-divergence formulation on the unstaggered grid. Semi-implicit semi-Lagrangian approach allows for long time steps while violates the global and local mass-conservation. In particular, the total mass in simulations with semi-Lagrangian models can drift significantly if no aposteriori mass-fixing algorithms are applied. However, the global mass-fixing algorithms degrade the local mass conservation. The inherently mass-conservative version of SL-AV model dynamical core presented in the article ensures global and local mass conservation without mass-fixing algorithms. The mass conservation is achieved with the introduction of the finite-volume semi-Lagrangian discretization for continuity equation based on the 3-D extension of the conservative cascade semi-Lagrangian transport scheme (CCS). The numerical experiments show that the presented new version of SL-AV dynamical core combines the accuracy and stability of the standard SL-AV dynamical core with the mass-conservation properties. The results of the mountain induced Rossby wave test and baroclinic instability test for mass-conservative dynamical core are found to be in agreement with the results available in literature.

Description: Potential and limitations of multidecadal satellite soil moisture observations for selected climate model evaluation studies Hydrology and Earth System Sciences, 17, 3523-3542, 2013 Author(s): A. Loew, T. Stacke, W. Dorigo, R. de Jeu, and S. Hagemann Soil moisture is an essential climate variable (ECV) of major importance for land–atmosphere interactions and global hydrology. An appropriate representation of soil moisture dynamics in global climate models is therefore important. Recently, a first multidecadal, observation-based soil moisture dataset has become available that provides information on soil moisture dynamics from satellite observations (ECVSM, essential climate variable soil moisture). The present study investigates the potential and limitations of this new dataset for several applications in climate model evaluation. We compare soil moisture data from satellite observations, reanalysis and simulations from a state-of-the-art land surface model and analyze relationships between soil moisture and precipitation anomalies in the different dataset. Other potential applications like model parameter optimization or model initialization are not investigated in the present study. In a detailed regional study, we show that ECVSM is capable to capture well the interannual and intraannual soil moisture and precipitation dynamics in the Sahelian region. Current deficits of the new dataset are critically discussed and summarized at the end of the paper to provide guidance for an appropriate usage of the ECVSM dataset for climate studies.

Description: A fast input/output library for high resolution climate models Geoscientific Model Development Discussions, 6, 4775-4807, 2013 Author(s): X. Huang, W. Wang, H. Fu, G. Yang, B. Wang, and C. Zhang We describe the design and implementation of Climate Fast Input/Output (CFIO), a fast input/output (I/O) library for high resolution climate models. CFIO provides a simple method for modelers to overlap the I/O phase with the computing phase automatically, so as to shorten the running time of numerical simulations. To minimize the code modifications required for porting, CFIO provides similar interfaces and features to Parallel network Common Data Form (PnetCDF), which is one of the most widely used I/O libraries in climate models. We deployed CFIO in three high resolution climate models, including two ocean models (POP and LICOM) and one sea ice model (CICE). The experimental results show that CFIO improves the performance of climate models significantly versus the original serial I/O approach. When running with CFIO at 0.1° resolution with about 1000 CPU cores, we managed to reduce the running time by factors of 7.9, 4.6 and 2.0 for POP, CICE, and LICOM respectively. We also compared the performance of CFIO against PnetCDF in different scenarios. For scenarios with both data output operations and computations, CFIO decreases the I/O overhead by a factor of 5.1 compared to PnetCDF.

Description: Simultaneous estimation of model state variables and observation and forecast biases using a two-stage hybrid Kalman filter Hydrology and Earth System Sciences, 17, 3499-3521, 2013 Author(s): V. R. N. Pauwels, G. J. M. De Lannoy, H.-J. Hendricks Franssen, and H. Vereecken In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Description: Application and evaluation of McICA scheme with new radiation code in BCC_AGCM2.0.1 Geoscientific Model Development Discussions, 6, 4933-4982, 2013 Author(s): given_name prefix surname suffix, H. Zhang, X. Jing, and J. Li This research incorporates the Monte Carlo Independent Column Approximation (McICA) scheme with the correlated k-distribution BCC-RAD radiation model into the climate model BCC_AGCM2.0.1 and examines the impacts on modeled climate through several simulations with variations in cloud structures. Results from experiments with consistent sub-grid cloud structures show that both clear-sky radiation fluxes and cloud radiative forcings (CRFs) calculated by the new scheme are mostly improved relative to those calculated from the original one. The modeled atmospheric temperature and specific humidity are also improved due to changes in the radiative heating rates. The vertical overlap of fractional clouds and horizontal distribution of cloud condensation are important for computing CRFs. The maximum changes in seasonal CRF using the general overlap assumption (GenO) with different decorrelation depths ( L cf ) are larger than 10 and 20 Wm 2 for longwave (LW) CRF and shortwave (SW) CRF, respectively, mostly located in the Tropics and mid-latitude storm tracks. Larger (smaller) L cf in the Tropics (mid-latitude storm tracks) yield better cloud fraction and CRF compared with observations. The inclusion of an observation-based horizontal inhomogeneity of cloud condensation has a distinct impact on LW CRF and SW CRF, with global means of ∼1.2 Wm −2 and ∼3.7 Wm −2 at the top of atmosphere, respectively, making these much closer to observations. These results prove the reliability of the new model configuration to be used in BCC_AGCM2.0.1 for climate simulations, and also indicate that more detailed real-world information on cloud structures should be obtained to constrain cloud settings in McICA in the future.

Description: Land use change effects on runoff generation in a humid tropical montane cloud forest region Hydrology and Earth System Sciences, 17, 3543-3560, 2013 Author(s): L. E. Muñoz-Villers and J. J. McDonnell While tropical montane cloud forests (TMCF) provide critical hydrological services to downstream regions throughout much of the humid tropics, catchment hydrology and impacts associated with forest conversion in these ecosystems remain poorly understood. Here, we compare the annual, seasonal and event-scale streamflow patterns and runoff generation processes of three neighbouring headwater catchments in central Veracruz (eastern Mexico) with similar pedological and geological characteristics, but different land cover: old-growth TMCF, 20 yr-old naturally regenerating TMCF and a heavily grazed pasture. We used a 2 yr record of high resolution rainfall and stream flow data (2008–2010) in combination with stable isotope and chemical tracer data collected for a series of storms during a 6-week period of increasing antecedent wetness (wetting-up cycle). Our results showed that annual and seasonal streamflow patterns in the mature and secondary forest were similar. In contrast, the pasture showed a 10% higher mean annual streamflow, most likely because of a lower rainfall interception. During the wetting-up cycle, storm runoff ratios increased at all three catchments (from 11 to 54% for the mature forest, 7 to 52% for the secondary forest and 3 to 59% for the pasture). With the increasing antecedent wetness, hydrograph separation analysis showed progressive increases of pre-event water contributions to total stormflow (from 35 to 99% in the mature forest, 26 to 92% in the secondary forest and 64 to 97% in the pasture). At all three sites, rainfall-runoff responses were dominated by subsurface flow generation processes for the majority of storms. However, for the largest and most intense storm (typically occurring once every 2 yr), sampled under wet antecedent conditions, the event water contribution in the pasture (34% on average) was much higher than in the forests (5% on average), indicating that rainfall infiltration capacity of the pasture was exceeded. This result suggests that despite the high permeability of the volcanic soils and underlying substrate in this TMCF environment, the conversion of forest to pasture may lead to important changes in runoff generation processes during large and high intensity storms. On the other hand, our results also showed that 20 yr of natural regeneration may be enough to largely restore the original hydrological conditions of this TMCF.

Description: Extreme extension across Seram and Ambon, eastern Indonesia: evidence for Banda slab rollback Solid Earth, 4, 277-314, 2013 Author(s): J. M. Pownall, R. Hall, and I. M. Watkinson The island of Seram, which lies in the northern part of the 180°-curved Banda Arc, has previously been interpreted as a fold-and-thrust belt formed during arc-continent collision, which incorporates ophiolites intruded by granites thought to have been produced by anatexis within a metamorphic sole. However, new geological mapping and a re-examination of the field relations cause us to question this model. We instead propose that there is evidence for recent and rapid N–S extension that has caused the high-temperature exhumation of lherzolites beneath low-angle lithospheric detachment faults that induced high-temperature metamorphism and melting in overlying crustal rocks. These "Kobipoto Complex" migmatites include highly residual Al–Mg-rich garnet + cordierite + sillimanite + spinel + corundum granulites (exposed in the Kobipoto Mountains) which contain coexisting spinel + quartz, indicating that peak metamorphic temperatures likely approached 900 °C. Associated with these residual granulites are voluminous Mio-Pliocene granitic diatexites, or "cordierite granites", which crop out on Ambon, western Seram, and in the Kobipoto Mountains and incorporate abundant schlieren of spinel- and sillimanite-bearing residuum. Quaternary "ambonites" (cordierite + garnet dacites) emplaced on Ambon were also evidently sourced from the Kobipoto Complex migmatites as demonstrated by granulite-inherited xenoliths. Exhumation of the hot peridotites and granulite-facies Kobipoto Complex migmatites to shallower structural levels caused greenschist- to lower-amphibolite facies metapelites and amphibolites of the Tehoru Formation to be overprinted by sillimanite-grade metamorphism, migmatisation, and limited localised anatexis to form the Taunusa Complex. The extreme extension required to have driven Kobipoto Complex exhumation evidently occurred throughout Seram and along much of the northern Banda Arc. The lherzolites must have been juxtaposed against the crust at typical lithospheric mantle temperatures in order to account for such high-temperature metamorphism and therefore could not have been part of a cooled ophiolite. In central Seram, lenses of peridotites are incorporated with a major left-lateral strike-slip shear zone (the "Kawa Shear Zone"), demonstrating that strike-slip motions likely initiated shortly after the mantle had been partly exhumed by detachment faulting and that the main strike-slip faults may themselves be reactivated and steepened low-angle detachments. The geodynamic driver for mantle exhumation along the detachment faults and strike-slip faulting in central Seram is very likely the same; we interpret the extreme extension to be the result of eastward slab rollback into the Banda Embayment as outlined by the latest plate reconstructions for Banda Arc evolution.

Description: A regional climate modelling projection ensemble experiment – NARCliM Geoscientific Model Development Discussions, 6, 5117-5139, 2013 Author(s): J. P. Evans, F. Ji, C. Lee, P. Smith, D. Argüeso, and L. Fita Including the impacts of climate change in decision making and planning processes is a challenge facing many regional governments including the New South Wales (NSW) and Australian Capital Territory (ACT) governments in Australia. NARCliM (NSW/ACT Regional Climate Modelling project) is a regional climate modelling project that aims to provide a comprehensive and consistent set of climate projections that can be used by all relevant government departments when considering climate change. To maximise end user engagement and ensure outputs are relevant to the planning process, a series of stakeholder workshops were run to define key aspects of the model experiment including spatial resolution, time slices, and output variables. As with all such experiments, practical considerations limit the number of ensembles members that can be simulated such that choices must be made concerning which Global Climate Models (GCMs) to downscale from, and which Regional Climate Models (RCMs) to downscale with. Here a methodology for making these choices is proposed that aims to sample the uncertainty in both GCMs and RCMs, as well as spanning the range of future climate projections present in the full GCM ensemble. The created ensemble provides a more robust view of future regional climate changes.

Description: Using the nonlinear aquifer storage–discharge relationship to simulate the base flow of glacier- and snowmelt-dominated basins in northwest China Hydrology and Earth System Sciences, 17, 3577-3586, 2013 Author(s): R. Gan and Y. Luo Base flow is an important component in hydrological modeling. This process is usually modeled by using the linear aquifer storage–discharge relation approach, although the outflow from groundwater aquifers is nonlinear. To identify the accuracy of base flow estimates in rivers dominated by snowmelt and/or glacier melt in arid and cold northwestern China, a nonlinear storage–discharge relationship for use in SWAT (Soil Water Assessment Tool) modeling was developed and applied to the Manas River basin in the Tian Shan Mountains. Linear reservoir models and a digital filter program were used for comparisons. Meanwhile, numerical analysis of recession curves from 78 river gauge stations revealed variation in the parameters of the nonlinear relationship. It was found that the nonlinear reservoir model can improve the streamflow simulation, especially for low-flow period. The higher Nash–Sutcliffe efficiency, logarithmic efficiency, and volumetric efficiency, and lower percent bias were obtained when compared to the one-linear reservoir approach. The parameter b of the aquifer storage–discharge function varied mostly between 0.0 and 0.1, which is much smaller than the suggested value of 0.5. The coefficient a of the function is related to catchment properties, primarily the basin and glacier areas.

Description: Technical Note: Reducing the spin-up time of integrated surface water–groundwater models Hydrology and Earth System Sciences, 18, 5169-5179, 2014 Author(s): H. Ajami, J. P. Evans, M. F. McCabe, and S. Stisen One of the main challenges in the application of coupled or integrated hydrologic models is specifying a catchment's initial conditions in terms of soil moisture and depth-to-water table (DTWT) distributions. One approach to reducing uncertainty in model initialization is to run the model recursively using either a single year or multiple years of forcing data until the system equilibrates with respect to state and diagnostic variables. However, such "spin-up" approaches often require many years of simulations, making them computationally intensive. In this study, a new hybrid approach was developed to reduce the computational burden of the spin-up procedure by using a combination of model simulations and an empirical DTWT function. The methodology is examined across two distinct catchments located in a temperate region of Denmark and a semi-arid region of Australia. Our results illustrate that the hybrid approach reduced the spin-up period required for an integrated groundwater–surface water–land surface model (ParFlow.CLM) by up to 50%. To generalize results to different climate and catchment conditions, we outline a methodology that is applicable to other coupled or integrated modeling frameworks when initialization from an equilibrium state is required.

Description: Efficient performance of the Met Office Unified Model v8.2 on Intel Xeon partially used nodes Geoscientific Model Development Discussions, 7, 7395-7425, 2014 Author(s): I. Bermous The atmospheric Unified Model (UM) developed at the UK Met Office is used for weather and climate prediction by forecast teams at a number of international meteorological centres and research institutes on a wide variety of hardware and software environments. Over its 25 year history the UM sources have been optimised for a better application performance on a number of HPC systems including NEC SX vector architecture systems and recently the IBM Power6/Power7 platforms. Understanding the influence of the compiler flags, MPI libraries and run configurations is crucial to achieving the shortest elapsed times for a UM application on any particular HPC system. These aspects are very important for applications that must run within operational time frames. Driving the current study is the HPC industry trend since 1980 for processor arithmetic performance to increase at a faster rate than memory bandwidth. This gap has been growing especially fast for multicore processors in the past 10 years and it can have significant implication for the performance and performance scaling of memory bandwidth intensive applications, such as the UM. Analysis of partially used nodes on Intel Xeon clusters is provided in this paper for short and medium range weather forecasting systems using global and limited-area configurations. It is shown that on the Intel Xeon based clusters the fastest elapsed times and the most efficient system usage can be achieved using partially committed nodes.

Description: Hydrogeology of an Alpine rockfall aquifer system and its role in flood attenuation and maintaining baseflow Hydrology and Earth System Sciences, 18, 4437-4452, 2014 Author(s): U. Lauber, P. Kotyla, D. Morche, and N. Goldscheider The frequency and intensity of extreme hydrological events in Alpine regions is projected to increase with climate change. The goal of this study is to better understand the functioning of aquifers composed of complex alluvial and rockfall deposits in Alpine valleys and to quantify the role of these natural storage spaces in flood attenuation and baseflow maintenance. Geomorphological and hydrogeological mapping, tracer tests, and continuous flow measurements were conducted in the Reintal (German Alps), where runoff from a karst spring infiltrates a series of postglacial alluvial/rockfall aquifers. During high-flow conditions, groundwater velocities of 30 m h −1 were determined along 500 m; hydrograph analyses revealed short lag times (5 h) between discharge peaks upstream and downstream from the aquifer series; the maximum discharge ratio downstream (22) and the peak recession coefficient (0.196 d −1 ) are low compared with other Alpine catchments. During low-flow conditions, the underground flow path length increased to 2 km and groundwater velocities decreased to 13 m h −1 . Downstream hydrographs revealed a delayed discharge response after 101 h and peaks damped by a factor of 1.5. These results indicate that alluvial/rockfall aquifers might play an important role in the flow regime and attenuation of floods in Alpine regions.

Description: Observed groundwater temperature response to recent climate change Hydrology and Earth System Sciences, 18, 4453-4466, 2014 Author(s): K. Menberg, P. Blum, B. L. Kurylyk, and P. Bayer Climate change is known to have a considerable influence on many components of the hydrological cycle. Yet, the implications for groundwater temperature, as an important driver for groundwater quality, thermal use and storage, are not yet comprehensively understood. Furthermore, few studies have examined the implications of climate-change-induced groundwater temperature rise for groundwater-dependent ecosystems. Here, we examine the coupling of atmospheric and groundwater warming by employing stochastic and deterministic models. Firstly, several decades of temperature time series are statistically analyzed with regard to climate regime shifts (CRSs) in the long-term mean. The observed increases in shallow groundwater temperatures can be associated with preceding positive shifts in regional surface air temperatures, which are in turn linked to global air temperature changes. The temperature data are also analyzed with an analytical solution to the conduction–advection heat transfer equation to investigate how subsurface heat transfer processes control the propagation of the surface temperature signals into the subsurface. In three of the four monitoring wells, the predicted groundwater temperature increases driven by the regime shifts at the surface boundary condition generally concur with the observed groundwater temperature trends. Due to complex interactions at the ground surface and the heat capacity of the unsaturated zone, the thermal signals from distinct changes in air temperature are damped and delayed in the subsurface, causing a more gradual increase in groundwater temperatures. These signals can have a significant impact on large-scale groundwater temperatures in shallow and economically important aquifers. These findings demonstrate that shallow groundwater temperatures have responded rapidly to recent climate change and thus provide insight into the vulnerability of aquifers and groundwater-dependent ecosystems to future climate change.

Description: Evaluation of the satellite-based Global Flood Detection System for measuring river discharge: influence of local factors Hydrology and Earth System Sciences, 18, 4467-4484, 2014 Author(s): B. Revilla-Romero, J. Thielen, P. Salamon, T. De Groeve, and G. R. Brakenridge One of the main challenges for global hydrological modelling is the limited availability of observational data for calibration and model verification. This is particularly the case for real-time applications. This problem could potentially be overcome if discharge measurements based on satellite data were sufficiently accurate to substitute for ground-based measurements. The aim of this study is to test the potentials and constraints of the remote sensing signal of the Global Flood Detection System for converting the flood detection signal into river discharge values. The study uses data for 322 river measurement locations in Africa, Asia, Europe, North America and South America. Satellite discharge measurements were calibrated for these sites and a validation analysis with in situ discharge was performed. The locations with very good performance will be used in a future project where satellite discharge measurements are obtained on a daily basis to fill the gaps where real-time ground observations are not available. These include several international river locations in Africa: the Niger, Volta and Zambezi rivers. Analysis of the potential factors affecting the satellite signal was based on a classification decision tree (random forest) and showed that mean discharge, climatic region, land cover and upstream catchment area are the dominant variables which determine good or poor performance of the measure\-ment sites. In general terms, higher skill scores were obtained for locations with one or more of the following characteristics: a river width higher than 1km; a large floodplain area and in flooded forest, a potential flooded area greater than 40%; sparse vegetation, croplands or grasslands and closed to open and open forest; leaf area index 〉 2; tropical climatic area; and without hydraulic infrastructures. Also, locations where river ice cover is seasonally present obtained higher skill scores. This work provides guidance on the best locations and limitations for estimating discharge values from these daily satellite signals.

Description: Changes in soil organic carbon and nitrogen capacities of Salix cheilophila Schneid along a revegetation chronosequence in semi-arid degraded sandy land of the Gonghe Basin, Tibet Plateau Solid Earth, 5, 1045-1054, 2014 Author(s): Y. Yu and Z. Q. Jia The Gonghe Basin is a sandified and desertified region of China, but the distribution of soil organic carbon (SOC) and total nitrogen (TN) along the cultivation chronosequence across this ecologically fragile region is not well understood. This study was carried out to understand the effects of restoration with Salix cheilophila for different periods of time (6, 11, 16, 21 years) to test whether it enhanced C and N storage. Soil samples, in four replications from seven depth increments (0–10, 10–20, 20–30, 30–50, 50–100, 100–150 and 150–200 cm), were collected in each stand. Soil bulk density, SOC, TN, aboveground biomass and root biomass were measured. Results indicated that changes occurred in both the upper and deeper soil layers with an increase in revegetation time. The 0–200 cm soil showed that the 6-year stand gained 3.89 Mg C ha −1 and 1.00 Mg N ha −1 , which accounted for 40.82% of the original SOC and 11.06% of the TN of the 0-year stand. The 11-year stand gained 7.82 Mg C ha −1 and 1.98 Mg N ha −1 in the 0–200 cm soil layers, accounting for 58.06% of the SOC and 19.80% of the TN of the 0-year stand. The 16-year stand gained 11.32 Mg C ha −1 and 3.30 Mg N ha −1 in the 0–200 cm soil layers, accounting for 66.71% of the SOC and 21.98% of the TN of the 0-year stand. The 21-year stand gained 13.05 Mg C ha −1 and 5.45 Mg N ha −1 from the same soil depth, accounting for 69.79% of the SOC and 40.47% of the TN compared with the 0-year stand. The extent of these changes depended on soil depth and plantation age. The results demonstrated that, as stand age increased, the storage of SOC and TN increased. These results further indicated that restoration with S. cheilophila has positive impacts on the Gonghe Basin and has increased the capacity of SOC sequestration and N storage. The shrub's role as carbon sink is compatible with system management and persistence. The findings are significant for assessing C and N sequestration accurately in semi-arid degraded high, cold sandy regions in the future.

Description: Simulation of groundwater and surface water over the continental US using a hyperresolution, integrated hydrologic model Geoscientific Model Development Discussions, 7, 7317-7349, 2014 Author(s): R. M. Maxwell, L. E. Condon, and S. J. Kollet Interactions between surface and groundwater systems are well-established theoretically and observationally. While numerical models that solve both surface and subsurface flow equations in a single framework (matrix) are increasingly being applied, computational limitations have restricted their use to local and regional studies. Regional or watershed, scale simulations have been effective tools in understanding hydrologic processes, however there are still many questions, such as the adaptation of water resources to anthropogenic stressors and climate variability, that need to be answered across large spatial extents at high resolution. In response to this "grand challenge" in hydrology, we present the results of a parallel, integrated hydrologic model simulating surface and subsurface flow at high spatial resolution (1 km) over much of continental North America (~ 6 300 000 or 6.3 million km 2 ). These simulations provide predictions of hydrologic states and fluxes, namely water table depth and streamflow, at unprecedented scale and resolution. The physically-based modeling approach used here requires limited parameterizations and relies only on more fundamental inputs, such as topography, hydrogeologic properties and climate forcing. Results are compared to observations and provide mechanistic insight into hydrologic process interaction. This study demonstrates both the feasibility of continental scale integrated models and their utility for improving our understanding of large-scale hydrologic systems; the combination of high resolution and large spatial extent facilitates novel analysis of scaling relationships using model outputs.

Description: EDDA: integrated simulation of debris flow erosion, deposition and property changes Geoscientific Model Development Discussions, 7, 7267-7316, 2014 Author(s): H. X. Chen and L. M. Zhang Debris flow material properties change during the initiation, transportation and deposition processes, which influences the runout characteristics of the debris flow. A quasi-three-dimensional depth-integrated numerical model, EDDA, is presented in this paper to simulate debris flow erosion, deposition and induced material property changes. The model considers changes in debris flow density, yield stress and dynamic viscosity during the flow process. The yield stress of debris flow mixture is determined at limit equilibrium using the Mohr–Coulomb equation, which is applicable to clear water flow, hyper-concentrated flow and fully developed debris flow. To assure numerical stability and computational efficiency at the same time, a variable time stepping algorithm is developed to solve the governing differential equations. Four numerical tests are conducted to validate the model. The first two tests involve a one-dimensional dam-break water flow and a one-dimensional debris flow with constant properties. The last two tests involve erosion and deposition, and the movement of multi-directional debris flows. The changes in debris flow mass and properties due to either erosion or deposition are shown to affect the runout characteristics significantly. The model is also applied to simulate a large-scale debris flow in Xiaojiagou Ravine to test the performance of the model in catchment-scale simulations. The results suggest that the model estimates well the volume, inundated area, and runout distance of the debris flow. The model is intended for use as a module in a real-time debris flow warning system.

Description: Technical Note: On the Matt–Shuttleworth approach to estimate crop water requirements Hydrology and Earth System Sciences, 18, 4341-4348, 2014 Author(s): J. P. Lhomme, N. Boudhina, and M. M. Masmoudi The Matt–Shuttleworth method provides a way to make a one-step estimate of crop water requirements with the Penman–Monteith equation by translating the crop coefficients, commonly available in United Nations Food and Agriculture Organization (FAO) publications, into equivalent surface resistances. The methodology is based upon the theoretical relationship linking crop surface resistance to a crop coefficient and involves the simplifying assumption that the reference crop evapotranspiration (ET 0 ) is equal to the Priestley–Taylor estimate with a fixed coefficient of 1.26. This assumption, used to eliminate the dependence of surface resistance on certain weather variables, is questionable; numerical simulations show that it can lead to substantial differences between the true value of surface resistance and its estimate. Consequently, the basic relationship between surface resistance and crop coefficient, without any assumption, appears to be more appropriate for inferring crop surface resistance, despite the interference of weather variables.

Description: Variational assimilation of remotely sensed flood extents using a 2-D flood model Hydrology and Earth System Sciences, 18, 4325-4339, 2014 Author(s): X. Lai, Q. Liang, H. Yesou, and S. Daillet A variational data assimilation (4D-Var) method is proposed to directly assimilate flood extents into a 2-D dynamic flood model to explore a novel way of utilizing the rich source of remotely sensed data available from satellite imagery for better analyzing or predicting flood routing processes. For this purpose, a new cost function is specially defined to effectively fuse the hydraulic information that is implicitly indicated in flood extents. The potential of using remotely sensed flood extents for improving the analysis of flood routing processes is demonstrated by applying the present new data assimilation approach to both idealized and realistic numerical experiments.

Description: Translating aboveground cosmic-ray neutron intensity to high-frequency soil moisture profiles at sub-kilometer scale Hydrology and Earth System Sciences, 18, 4363-4379, 2014 Author(s): R. Rosolem, T. Hoar, A. Arellano, J. L. Anderson, W. J. Shuttleworth, X. Zeng, and T. E. Franz Above-ground cosmic-ray neutron measurements provide an opportunity to infer soil moisture at the sub-kilometer scale. Initial efforts to assimilate those measurements have shown promise. This study expands such analysis by investigating (1) how the information from aboveground cosmic-ray neutrons can constrain the soil moisture at distinct depths simulated by a land surface model, and (2) how changes in data availability (in terms of retrieval frequency) impact the dynamics of simulated soil moisture profiles. We employ ensemble data assimilation techniques in a "nearly-identical twin" experiment applied at semi-arid shrubland, rainfed agricultural field, and mixed forest biomes in the USA. The performance of the Noah land surface model is compared with and without assimilation of observations at hourly intervals, as well as every 2 days. Synthetic observations of aboveground cosmic-ray neutrons better constrain the soil moisture simulated by Noah in root-zone soil layers (0–100cm), despite the limited measurement depth of the sensor (estimated to be 12–20cm). The ability of Noah to reproduce a "true" soil moisture profile is remarkably good, regardless of the frequency of observations at the semi-arid site. However, soil moisture profiles are better constrained when assimilating synthetic cosmic-ray neutron observations hourly rather than every 2 days at the cropland and mixed forest sites. This indicates potential benefits for hydrometeorological modeling when soil moisture measurements are available at a relatively high frequency. Moreover, differences in summertime meteorological forcing between the semi-arid site and the other two sites may indicate a possible controlling factor to soil moisture dynamics in addition to differences in soil and vegetation properties.

Description: Large-scale 3-D modeling by integration of resistivity models and borehole data through inversion Hydrology and Earth System Sciences, 18, 4349-4362, 2014 Author(s): N. Foged, P. A. Marker, A. V. Christansen, P. Bauer-Gottwein, F. Jørgensen, A.-S. Høyer, and E. Auken We present an automatic method for parameterization of a 3-D model of the subsurface, integrating lithological information from boreholes with resistivity models through an inverse optimization, with the objective of further detailing of geological models, or as direct input into groundwater models. The parameter of interest is the clay fraction, expressed as the relative length of clay units in a depth interval. The clay fraction is obtained from lithological logs and the clay fraction from the resistivity is obtained by establishing a simple petrophysical relationship, a translator function, between resistivity and the clay fraction. Through inversion we use the lithological data and the resistivity data to determine the optimum spatially distributed translator function. Applying the translator function we get a 3-D clay fraction model, which holds information from the resistivity data set and the borehole data set in one variable. Finally, we use k -means clustering to generate a 3-D model of the subsurface structures. We apply the procedure to the Norsminde survey in Denmark, integrating approximately 700 boreholes and more than 100 000 resistivity models from an airborne survey in the parameterization of the 3-D model covering 156 km 2 . The final five-cluster 3-D model differentiates between clay materials and different high-resistivity materials from information held in the resistivity model and borehole observations, respectively.

Description: Regional parent flood frequency distributions in Europe – Part 2: Climate and scale controls Hydrology and Earth System Sciences, 18, 4391-4401, 2014 Author(s): J. L. Salinas, A. Castellarin, S. Kohnová, and T. R. Kjeldsen This study aims to better understand the effect of catchment scale and climate on the statistical properties of regional flood frequency distributions. A database of L-moment ratios of annual maximum series (AMS) of peak discharges from Austria, Italy and Slovakia, involving a total of 813 catchments with more than 25 yr of record length is presented, together with mean annual precipitation (MAP) and basin area as catchment descriptors surrogates of climate and scale controls. A purely data-based investigation performed on the database shows that the generalized extreme value (GEV) distribution provides a better representation of the averaged sample L-moment ratios compared to the other distributions considered, for catchments with medium to higher values of MAP independently of catchment area, while the three-parameter lognormal distribution is probably a more appropriate choice for drier (lower MAP) intermediate-sized catchments, which presented higher skewness values. Sample L-moment ratios do not follow systematically any of the theoretical two-parameter distributions. In particular, the averaged values of L-coefficient of skewness (L-Cs) are always larger than Gumbel 's fixed L-Cs. The results presented in this paper contribute to the progress in defining a set of "process-driven" pan-European flood frequency distributions and to assess possible effects of environmental change on its properties.

Description: Quantification of hydrologic impacts of climate change in a Mediterranean basin in Sardinia, Italy, through high-resolution simulations Hydrology and Earth System Sciences, 18, 5201-5217, 2014 Author(s): M. Piras, G. Mascaro, R. Deidda, and E. R. Vivoni Future climate projections robustly indicate that the Mediterranean region will experience a significant decrease of mean annual precipitation and an increase in temperature. These changes are expected to seriously affect the hydrologic regime, with a limitation of water availability and an intensification of hydrologic extremes, and to negatively impact local economies. In this study, we quantify the hydrologic impacts of climate change in the Rio Mannu basin (RMB), an agricultural watershed of 472.5 km 2 in Sardinia, Italy. To simulate the wide range of runoff generation mechanisms typical of Mediterranean basins, we adopted a physically based, distributed hydrologic model. The high-resolution forcings in reference and future conditions (30-year records for each period) were provided by four combinations of global and regional climate models, bias-corrected and downscaled in space and time (from ~25 km, 24 h to 5 km, 1 h) through statistical tools. The analysis of the hydrologic model outputs indicates that the RMB is expected to be severely impacted by future climate change. The range of simulations consistently predict (i) a significant diminution of mean annual runoff at the basin outlet, mainly due to a decreasing contribution of the runoff generation mechanisms depending on water available in the soil; (ii) modest variations in mean annual runoff and intensification of mean annual discharge maxima in flatter sub-basins with clay and loamy soils, likely due to a higher occurrence of infiltration excess runoff; (iii) reduction of soil water content and actual evapotranspiration in most areas of the basin; and (iv) a drop in the groundwater table. Results of this study are useful to support the adoption of adaptive strategies for management and planning of agricultural activities and water resources in the region.

Description: Northern Hemisphere storminess in the Norwegian Earth System Model (NorESM1-M) Geoscientific Model Development Discussions, 7, 8975-9015, 2014 Author(s): E. M. Knudsen and J. E. Walsh Metrics of storm activity in Northern Hemisphere high- and midlatitudes are evaluated from historical output and future projections by the Norwegian Earth System Model (NorESM1-M) coupled global climate model. The European Re-Analysis Interim (ERA-Interim) and the Community Climate System Model (CCSM4), a global climate model of the same vintage as NorESM1-M, provide benchmarks for comparison. The focus is on the autumn and early winter (September through December), the period when the ongoing and projected Arctic sea ice retreat is greatest. Storm tracks derived from a vorticity-based algorithm for storm identification are reproduced well by NorESM1-M, although the tracks are somewhat better resolved in the higher-resolution ERA-Interim and CCSM4. The tracks are projected to shift polewards in the future as climate changes under the Representative Concentration Pathway (RCP) forcing scenarios. Cyclones are projected to become generally more intense in the high-latitudes, especially over the Alaskan region, although in some other areas the intensity is projected to decrease. While projected changes in track density are less coherent, there is a general tendency towards less frequent storms in midlatitudes and more frequent storms in high-latitudes, especially the Baffin Bay/Davis Strait region. Autumn precipitation is projected to increase significantly across the entire high-latitudes. Together with the projected increases in storm intensity and sea level and the loss of sea ice, this increase in precipitation implies a greater vulnerability to coastal flooding and erosion, especially in the Alaskan region. The projected changes in storm intensity and precipitation (as well as sea ice and sea level pressure) scale generally linearly with the RCP value of the forcing and with time through the 21st century.

Description: Assessing winter cover crop nutrient uptake efficiency using a water quality simulation model Hydrology and Earth System Sciences, 18, 5239-5253, 2014 Author(s): I.-Y. Yeo, S. Lee, A. M. Sadeghi, P. C. Beeson, W. D. Hively, G. W. McCarty, and M. W. Lang Winter cover crops are an effective conservation management practice with potential to improve water quality. Throughout the Chesapeake Bay watershed (CBW), which is located in the mid-Atlantic US, winter cover crop use has been emphasized, and federal and state cost-share programs are available to farmers to subsidize the cost of cover crop establishment. The objective of this study was to assess the long-term effect of planting winter cover crops to improve water quality at the watershed scale (~ 50 km 2 ) and to identify critical source areas of high nitrate export. A physically based watershed simulation model, Soil and Water Assessment Tool (SWAT), was calibrated and validated using water quality monitoring data to simulate hydrological processes and agricultural nutrient cycling over the period of 1990–2000. To accurately simulate winter cover crop biomass in relation to growing conditions, a new approach was developed to further calibrate plant growth parameters that control the leaf area development curve using multitemporal satellite-based measurements of species-specific winter cover crop performance. Multiple SWAT scenarios were developed to obtain baseline information on nitrate loading without winter cover crops and to investigate how nitrate loading could change under different winter cover crop planting scenarios, including different species, planting dates, and implementation areas. The simulation results indicate that winter cover crops have a negligible impact on the water budget but significantly reduce nitrate leaching to groundwater and delivery to the waterways. Without winter cover crops, annual nitrate loading from agricultural lands was approximately 14 kg ha −1 , but decreased to 4.6–10.1 kg ha −1 with cover crops resulting in a reduction rate of 27–67% at the watershed scale. Rye was the most effective species, with a potential to reduce nitrate leaching by up to 93% with early planting at the field scale. Early planting of cover crops (~ 30 days of additional growing days) was crucial, as it lowered nitrate export by an additional ~ 2 kg ha −1 when compared to late planting scenarios. The effectiveness of cover cropping increased with increasing extent of cover crop implementation. Agricultural fields with well-drained soils and those that were more frequently used to grow corn had a higher potential for nitrate leaching and export to the waterways. This study supports the effective implementation of cover crop programs, in part by helping to target critical pollution source areas for cover crop implementation.

Description: Preface: Environmental benefits of biochar Solid Earth, 5, 1301-1303, 2014 Author(s): J. Paz-Ferreiro, A. Méndez, A. M. Tarquis, A. Cerdà, and G. Gascó

Description: Derivation and evaluation of landslide-triggering thresholds by a Monte Carlo approach Hydrology and Earth System Sciences, 18, 4913-4931, 2014 Author(s): D. J. Peres and A. Cancelliere Assessment of landslide-triggering rainfall thresholds is useful for early warning in prone areas. In this paper, it is shown how stochastic rainfall models and hydrological and slope stability physically based models can be advantageously combined in a Monte Carlo simulation framework to generate virtually unlimited-length synthetic rainfall and related slope stability factor of safety data, exploiting the information contained in observed rainfall records and field-measurements of soil hydraulic and geotechnical parameters. The synthetic data set, dichotomized in triggering and non-triggering rainfall events, is analyzed by receiver operating characteristics (ROC) analysis to derive stochastic-input physically based thresholds that optimize the trade-off between correct and wrong predictions. Moreover, the specific modeling framework implemented in this work, based on hourly analysis, enables one to analyze the uncertainty related to variability of rainfall intensity within events and to past rainfall (antecedent rainfall). A specific focus is dedicated to the widely used power-law rainfall intensity–duration ( I – D ) thresholds. Results indicate that variability of intensity during rainfall events influences significantly rainfall intensity and duration associated with landslide triggering. Remarkably, when a time-variable rainfall-rate event is considered, the simulated triggering points may be separated with a very good approximation from the non-triggering ones by a I – D power-law equation, while a representation of rainfall as constant–intensity hyetographs globally leads to non-conservative results. This indicates that the I – D power-law equation is adequate to represent the triggering part due to transient infiltration produced by rainfall events of variable intensity and thus gives a physically based justification for this widely used threshold form, which provides results that are valid when landslide occurrence is mostly due to that part. These conditions are more likely to occur in hillslopes of low specific upslope contributing area, relatively high hydraulic conductivity and high critical wetness ratio. Otherwise, rainfall time history occurring before single rainfall events influences landslide triggering, determining whether a threshold based only on rainfall intensity and duration may be sufficient or it needs to be improved by the introduction of antecedent rainfall variables. Further analyses show that predictability of landslides decreases with soil depth, critical wetness ratio and the increase of vertical basal drainage (leakage) that occurs in the presence of a fractured bedrock.

Description: Climate change and stream temperature projections in the Columbia River basin: habitat implications of spatial variation in hydrologic drivers Hydrology and Earth System Sciences, 18, 4897-4912, 2014 Author(s): D. L. Ficklin, B. L. Barnhart, J. H. Knouft, I. T. Stewart, E. P. Maurer, S. L. Letsinger, and G. W. Whittaker Water temperature is a primary physical factor regulating the persistence and distribution of aquatic taxa. Considering projected increases in air temperature and changes in precipitation in the coming century, accurate assessment of suitable thermal habitats in freshwater systems is critical for predicting aquatic species' responses to changes in climate and for guiding adaptation strategies. We use a hydrologic model coupled with a stream temperature model and downscaled general circulation model outputs to explore the spatially and temporally varying changes in stream temperature for the late 21st century at the subbasin and ecological province scale for the Columbia River basin (CRB). On average, stream temperatures are projected to increase 3.5 °C for the spring, 5.2 °C for the summer, 2.7 °C for the fall, and 1.6 °C for the winter. While results indicate changes in stream temperature are correlated with changes in air temperature, our results also capture the important, and often ignored, influence of hydrological processes on changes in stream temperature. Decreases in future snowcover will result in increased thermal sensitivity within regions that were previously buffered by the cooling effect of flow originating as snowmelt. Other hydrological components, such as precipitation, surface runoff, lateral soil water flow, and groundwater inflow, are negatively correlated to increases in stream temperature depending on the ecological province and season. At the ecological province scale, the largest increase in annual stream temperature was within the Mountain Snake ecological province, which is characterized by migratory coldwater fish species. Stream temperature changes varied seasonally with the largest projected stream temperature increases occurring during the spring and summer for all ecological provinces. Our results indicate that stream temperatures are driven by local processes and ultimately require a physically explicit modeling approach to accurately characterize the habitat regulating the distribution and diversity of aquatic taxa.

Description: Impact of changes in grain size and pore space on the hydraulic conductivity and spectral induced polarization response of sand Hydrology and Earth System Sciences, 15, 1785-1794, 2011 Author(s): K. Koch, A. Kemna, J. Irving, and K. Holliger Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. This underlines the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raises some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.

Description: WRF simulation of a precipitation event over the Tibetan Plateau, China – an assessment using remote sensing and ground observations Hydrology and Earth System Sciences, 15, 1795-1817, 2011 Author(s): F. Maussion, D. Scherer, R. Finkelnburg, J. Richters, W. Yang, and T. Yao Meteorological observations over the Tibetan Plateau (TiP) are scarce, and precipitation estimations over this remote region are difficult. The constantly improving capabilities of numerical weather prediction (NWP) models offer the opportunity to reduce this problem by providing precipitation fields and other meteorological variables of high spatial and temporal resolution. Longer time periods of years to decades can be simulated by NWP models by successive model runs of shorter periods, which can be described by the term "regional atmospheric reanalysis". In this paper, we assess the Weather Research and Forecasting (WRF) models capacity in retrieving rain- and snowfall on the TiP in such a configuration using a nested approach: the simulations are conducted with three nested domains at spatial resolutions of 30, 10, and 2 km. A validation study is carried out for a one-month period with a special focus on one-week (22–28 October 2008), during which strong rain- and snowfall was observed on the TiP. The output of the model in each resolution is compared to the Tropical Rainfall Measuring Mission (TRMM) data set for precipitation and to the Moderate Resolution Imaging Spectroradiometer (MODIS) data set for snow extent. TRMM and WRF data are then compared to weather-station measurements. Our results suggest an overall improvement from WRF over TRMM with respect to weather-station measurements. Various configurations of the model with different nesting and forcing strategies, as well as physical parameterisation schemes are compared to propose a suitable design for a regional atmospheric reanalysis over the TiP. The WRF model showed good accuracy in simulating snow- and rainfall on the TiP for a one-month simulation period. Our study reveals that there is nothing like an optimal model strategy applicable for the high-altitude TiP, its fringing high-mountain areas of extremely complex topography and the low-altitude land and sea regions from which much of the precipitation on the TiP is originating. The choice of the physical parameterisation scheme will thus be always a compromise depending on the specific purpose of a model simulation. Our study demonstrates the high importance of orographic precipitation, but the problem of the orographic bias remains unsolved since reliable observational data are still missing. The results are relevant for anyone interested in carrying out a regional atmospheric reanalysis. Many hydrological analyses and applications like rainfall-runoff modelling or the analysis of flood events require precipitation rates at daily or even hourly intervals. Thus, our study offers a process-oriented alternative for retrieving precipitation fields of high spatio-temporal resolution in regions like the TiP, where other data sources are limited.

Description: Spatio-temporal variations in soil moisture and physicochemical properties of a typical semiarid sand-meadow-desert landscape as influenced by land use Hydrology and Earth System Sciences, 15, 1865-1877, 2011 Author(s): L. Duan, T. Liu, X. Wang, G. Wang, L. Ma, and Y. Luo A good understanding of the interrelations between land cover alteration and changes in hydrologic conditions (e.g., soil moisture) as well as soil physicochemical properties (e.g., fine soil particles and nutrients) is crucial for maintaining the fragile hydrologic and environmental conditions of semiarid land, such as the Horqin Sandy Land in China, but is lacking in existing literature. The objectives of this study were to examine: (1) spatio-temporal variations of soil moisture and physicochemical properties in semiarid land; and (2) how those variations are influenced by land cover alteration. Using the data collected in a 9.71 km 2 well-instrumented area of the Horqin Sandy Land, this study examined by visual examination and statistical analyses the spatio-temporal variations of soil moisture and physicochemical properties. The results indicated that for the study area, the soil moisture and physicochemical properties were dependent on local topography, soil texture, vegetation density, and human activity. Long-term reclamation for agriculture was found to reduce soil moisture by over 23 % and significantly (p-value 〈 0.05) lower the contents of soil organic matter, fine soil particles, and nutrients.

Description: Exploring the potentials and limitations of the time-reversal imaging of finite seismic sources Solid Earth, 2, 95-105, 2011 Author(s): S. Kremers, A. Fichtner, G. B. Brietzke, H. Igel, C. Larmat, L. Huang, and M. Käser The characterisation of seismic sources with time-reversed wave fields is developing into a standard technique that has already been successful in numerous applications. While the time-reversal imaging of effective point sources is now well-understood, little work has been done to extend this technique to the study of finite rupture processes. This is despite the pronounced non-uniqueness in classic finite source inversions. The need to better constrain the details of finite rupture processes motivates the series of synthetic and real-data time reversal experiments described in this paper. We address questions concerning the quality of focussing in the source area, the localisation of the fault plane, the estimation of the slip distribution and the source complexity up to which time-reversal imaging can be applied successfully. The frequency band for the synthetic experiments is chosen such that it is comparable to the band usually employed for finite source inversion. Contrary to our expectations, we find that time-reversal imaging is useful only for effective point sources, where it yields good estimates of both the source location and the origin time. In the case of finite sources, however, the time-reversed field does not provide meaningful characterisations of the fault location and the rupture process. This result cannot be improved sufficiently with the help of different imaging fields, realistic modifications of the receiver geometry or weights applied to the time-reversed sources. The reasons for this failure are manifold. They include the choice of the frequency band, the incomplete recording of wave field information at the surface, the excitation of large-amplitude surface waves that deteriorate the depth resolution, the absence of a sink that should absorb energy radiated during the later stages of the rupture process, the invisibility of small slip and the neglect of prior information concerning the fault geometry and the inherent smoothness of seismologically inferred Earth models that prevents the beneficial occurrence of strong multiple-scattering. The condensed conclusion of our study is that the limitations of time-reversal imaging – at least in the frequency band considered here – start where the seismic source stops being effectively point-localised.

Description: Anatomy of extraordinary rainfall and flash flood in a Dutch lowland catchment Hydrology and Earth System Sciences, 15, 1991-2005, 2011 Author(s): C. C. Brauer, A. J. Teuling, A. Overeem, Y. van der Velde, P. Hazenberg, P. M. M. Warmerdam, and R. Uijlenhoet On 26 August 2010 the eastern part of The Netherlands and the bordering part of Germany were struck by a series of rainfall events lasting for more than a day. Over an area of 740 km 2 more than 120 mm of rainfall were observed in 24 h. This extreme event resulted in local flooding of city centres, highways and agricultural fields, and considerable financial loss. In this paper we report on the unprecedented flash flood triggered by this exceptionally heavy rainfall event in the 6.5 km 2 Hupsel Brook catchment, which has been the experimental watershed employed by Wageningen University since the 1960s. This study aims to improve our understanding of the dynamics of such lowland flash floods. We present a detailed hydrometeorological analysis of this extreme event, focusing on its synoptic meteorological characteristics, its space-time rainfall dynamics as observed with rain gauges, weather radar and a microwave link, as well as the measured soil moisture, groundwater and discharge response of the catchment. At the Hupsel Brook catchment 160 mm of rainfall was observed in 24 h, corresponding to an estimated return period of well over 1000 years. As a result, discharge at the catchment outlet increased from 4.4 × 10 −3 to nearly 5 m 3 s −1 . Within 7 h discharge rose from 5 × 10 −2 to 4.5 m 3 s −1 . The catchment response can be divided into four phases: (1) soil moisture reservoir filling, (2) groundwater response, (3) surface depression filling and surface runoff and (4) backwater feedback. The first 35 mm of rainfall were stored in the soil without a significant increase in discharge. Relatively dry initial conditions (in comparison to those for past discharge extremes) prevented an even faster and more extreme hydrological response.

Description: On the value of combined event runoff and tracer analysis to improve understanding of catchment functioning in a data-scarce semi-arid area Hydrology and Earth System Sciences, 15, 2007-2024, 2011 Author(s): M. Hrachowitz, R. Bohte, M. L. Mul, T. A. Bogaard, H. H. G. Savenije, and S. Uhlenbrook Hydrological processes in small catchments are not quite understood yet, which is true in particular for catchments in data scarce, semi-arid regions. This is in contrast with the need for a better understanding of water fluxes and the interactions between surface- and groundwater in order to facilitate sustainable water resources management in such environments, where both floods and droughts can result in severe crop loss. In this study, event runoff coefficient analysis and limited tracer data of four small, nested sub-catchments (0.4–25.3 km 2 ) in a data scarce, semi-arid region of Tanzania helped to characterize the distinct response of the study catchments and to gain insights into the dominant runoff processes. The estimated event runoff coefficients were very low and did not exceed 0.09. They were found to be significantly related to the 5-day antecedent precipitation totals as well as to base flow, indicating a close relation to changes in soil moisture and thus potential switches in runoff generation processes. The time scales of the "direct flow" reservoirs, used to compute the event runoff coefficients, were up to one order of magnitude reduced for extreme events, compared to "average" events, suggesting the activation of at least a third flow component, besides base- and direct flow, assumed to be infiltration overland flow. Analysis of multiple tracers highlighted the importance of pre-event water to total runoff, even during intense and high yield precipitation events. It further illustrated the distinct nature of the catchments, in particular with respect to the available water storage, which was suggested by different degrees of tracer damping in the individual streams. The use of multiple tracers subsequently allowed estimating uncertainties in hydrograph separations arising from the use of different tracers. The results highlight the presence of considerable uncertainties, emphasizing the need for multiple tracers in order to avoid misleading results. This study shows the value of hydrological data collection over one whole wet season using multi-tracers to improve the understanding of hydrological functioning and thus for water resources management in data scarce, semi-arid environments.

Description: Changes in land cover, rainfall and stream flow in Upper Gilgel Abbay catchment, Blue Nile basin – Ethiopia Hydrology and Earth System Sciences, 15, 1979-1989, 2011 Author(s): T. H. M. Rientjes, A. T. Haile, E. Kebede, C. M. M. Mannaerts, E. Habib, and T. S. Steenhuis In this study we evaluated changes in land cover and rainfall in the Upper Gilgel Abbay catchment in the Upper Blue Nile basin and how changes affected stream flow in terms of annual flow, high flows and low flows. Land cover change assessment was through classification analysis of remote sensing based land cover data while assessments on rainfall and stream flow data are by statistical analysis. Results of the supervised land cover classification analysis indicated that 50.9 % and 16.7 % of the catchment area was covered by forest in 1973 and 2001, respectively. This significant decrease in forest cover is mainly due to expansion of agricultural land. By use of a change detection procedure, three periods were identified for which changes in rainfall and stream flow were analyzed. Rainfall was analyzed at monthly base by use of the Mann-Kendall test statistic and results indicated a statistically significant, decreasing trend for most months of the year. However, for the wet season months of June, July and August rainfall has increased. In the period 1973–2005, the annual flow of the catchment decreased by 12.1 %. Low flow and high flow at daily base were analyzed by a low flow and a high flow index that is based on a 95 % and 5 % exceedance probability. Results of the low flow index indicated decreases of 18.1 % and 66.6 % for the periods 1982–2000 and 2001–2005 respectively. Results of high flows indicated an increase of 7.6 % and 46.6 % for the same periods. In this study it is concluded that over the period 1973–2005 stream flow has changed in the Gilgel Abbay catchment by changes in land cover and changes in rainfall.

Description: Native American lithic procurement along the international border in the boot heel region of southwestern New Mexico Solid Earth, 2, 75-93, 2011 Author(s): K. E. Zeigler, P. Hogan, C. Hughes, and A. Kurota Multidisciplinary field projects can be very useful to a more fundamental understanding of the world around us, though these projects are not as common as they should be. In particular, the combination of archeology and geology combines our understanding of human behavior and human use of the landscape with an intimate knowledge of geologic processes and the materials available for human use in order to gain a broader understanding of human-Earth interaction. Here we present data from a cross-disciplinary project that uses a common dataset, archeological artifacts, to explore the anthropological and geologic implications of useage patterns. Archeological excavations and surveys conducted by the Office of Contract Archeology in 2007 along the route of the proposed international border fence reveal patterns of use of geologic materials by Archaic, Formative and Protohistoric Native Americans in the Boot Heel of southwestern New Mexico. Thousands of artifacts were recorded in multiple sites from Guadalupe Pass in the southern Peloncillo Mountains to the Carrizalillo Hills west of Columbus. We identified the lithologies of artifacts, ranging from projectile points to groundstones, and then constructed material movement maps based on either known procurement sites ("quarries") or outcrops identified as the closest source to a given site for each lithology. Not unexpectedly, the majority of the rock types utilized by native peoples are local siliceous volcanic materials. However, several artifacts constructed from obsidian were transported into the region from northern Mexico and eastern Arizona, indicating long-distance travel and/or trade routes. We also examine useage pattern difference between Archaic, Formative and Protohistoric sites. Additionally, a dramatic change in distribution of sources for geologic materials occurs between one pre-Spanish site and one post-Spanish site that are adjacent to one another.

Description: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 1: The limited-area atmospheric chemistry model COSMO/MESSy Geoscientific Model Development Discussions, 4, 1305-1358, 2011 Author(s): A. Kerkweg and P. Jöckel The numerical weather prediction model of the Consortium for Small Scale Modelling (COSMO), maintained by the German weather service (DWD), is connected with the Modular Earth Submodel System (MESSy). This effort is undertaken in preparation of a~new, limited-area atmospheric chemistry model. This model is as consistent as possible, with respect to atmospheric chemistry and related processes, with a previously developed global atmospheric chemistry general circulation model: the ECHAM/MESSy Atmospheric Chemistry (EMAC) model. The combined system constitutes a new research tool, bridging the global to the meso-γ scale for atmospheric chemistry research. MESSy provides the infrastructure and includes, among others, the process and diagnostic submodels for atmospheric chemistry simulations. Furthermore, MESSy is highly flexible allowing model setups with tailor made complexity, depending on the scientific question. Here, the connection of the MESSy infrastructure to the COSMO model is documented. Previously published prototype submodels for simplified tracer studies are generalised to be plugged-in and used in the global and the limited-area model. They are used to evaluate the tracer transport characteristics of the new COSMO/MESSy model system, an important prerequisite for future atmospheric chemistry applications. A supplementary document with further details on the technical implementation of the MESSy interface into COSMO with a complete list of modifications to the COSMO code is provided.

Description: Hydropedological insights when considering catchment classification Hydrology and Earth System Sciences, 15, 1909-1919, 2011 Author(s): J. Bouma, P. Droogers, M. P. W. Sonneveld, C. J. Ritsema, J. E. Hunink, W. W. Immerzeel, and S. Kauffman Soil classification systems are analysed to explore the potential of developing classification systems for catchments. Soil classifications are useful to create systematic order in the overwhelming quantity of different soils in the world and to extrapolate data available for a given soil type to soils elsewhere with identical classifications. This principle also applies to catchments. However, to be useful, soil classifications have to be based on permanent characteristics as formed by the soil forming factors over often very long periods of time. When defining permanent catchment characteristics, discharge data would therefore appear to be less suitable. But permanent soil characteristics do not necessarily match with characteristics and parameters needed for functional soil characterization focusing, for example, on catchment hydrology. Hydropedology has made contributions towards the required functional characterization of soils as is illustrated for three recent hydrological catchment studies. However, much still needs to be learned about the physical behaviour of anisotropic, heterogeneous soils with varying soil structures during the year and about spatial and temporal variability. The suggestion is made therefore to first focus on improving simulation of catchment hydrology, possibly incorporating hydropedological expertise, before embarking on a catchment classification effort which involves major input of time and involves the risk of distraction. In doing so, we suggest to also define other characteristics for catchment performance than the traditionally measured discharge rates. Such characteristics may well be derived from societal issues being studied, as is illustrated for the Green Water Credits program.

Description: Quantifying spatial and temporal discharge dynamics of an event in a first order stream, using distributed temperature sensing Hydrology and Earth System Sciences, 15, 1945-1957, 2011 Author(s): M. C. Westhoff, T. A. Bogaard, and H. H. G. Savenije Understanding the spatial distribution of discharge can be important for water quality and quantity modeling. Non-steady flood waves can, particularly as a result of short high intensity summer rainstorms, influence small headwater streams significantly. The aim of this paper is to quantify the spatial and temporal dynamics of stream flow in a headwater stream during a summer rainstorm. These dynamics include gains and losses of stream water, the effect of bypasses that become active and hyporheic exchange fluxes that may vary over time as a function of discharge. We use an advection-dispersion model coupled with an energy balance model to simulate in-stream water temperature, which we compare with high resolution temperature observations obtained with Distributed Temperature Sensing. This model was used as a learning tool to stepwise unravel the complex puzzle of in-stream processes subject to varying discharge. Hypotheses were tested and rejected, which led to more insight in the spatial and temporal dynamics in discharge and hyporheic exchange processes. We showed that, for the studied stream infiltration losses increase during a small rain event, while gains of water remained constant over time. We conclude that, eventually, part of the stream water bypassed the main channel during peak discharge. It also seems that hyporheic exchange varies with varying discharge in the first 250 m of the stream; while further downstream it remains constant. Because we relied on solar radiation as the main energy input, we were only able to apply this method during a small summer storm and low flow conditions. However, when additional (artificial) energy is available, the presented method is also applicable in larger streams, during higher flow conditions or longer storms.

Description: Trend analysis of extreme precipitation in the Northwestern Highlands of Ethiopia with a case study of Debre Markos Hydrology and Earth System Sciences, 15, 1937-1944, 2011 Author(s): H. Shang, J. Yan, M. Gebremichael, and S. M. Ayalew Understanding extreme precipitation is very important for Ethiopia, which is heavily dependent on low-productivity rainfed agriculture but lacks structural and non-structural water regulating and storage mechanisms. There has been an increasing concern about whether there is an increasing trend in extreme precipitation as the climate changes. Existing analysis of this region has been descriptive, without taking advantage of the advances in extreme value modeling. After reviewing the statistical methodology on extremes, this paper presents an analysis based on the generalized extreme value modeling with daily time series of precipitation records at Debre Markos in the Northwestern Highlands of Ethiopia. We found no strong evidence to reject the null hypothesis that there is no increasing trend in extreme precipitation at this location.

Description: A framework for the quantitative assessment of climate change impacts on water-related activities at the basin scale Hydrology and Earth System Sciences, 15, 2025-2038, 2011 Author(s): D. Anghileri, F. Pianosi, and R. Soncini-Sessa While quantitative assessment of the climate change impact on hydrology at the basin scale is quite addressed in the literature, extension of quantitative analysis to impact on the ecological, economic and social sphere is still limited, although well recognized as a key issue to support water resource planning and promote public participation. In this paper we propose a framework for assessing climate change impact on water-related activities at the basin scale. The specific features of our approach are that: (i) the impact quantification is based on a set of performance indicators defined together with the stakeholders, thus explicitly taking into account the water-users preferences; (ii) the management policies are obtained by optimal control techniques, linking stakeholder expectations and decision-making; (iii) the multi-objective nature of the management problem is fully preserved by simulating a set of Pareto-optimal management policies, which allows for evaluating not only variations in the indicator values but also tradeoffs among conflicting objectives. The framework is demonstrated by application to a real world case study, Lake Como basin (Italy). We show that the most conflicting water-related activities within the basin (i.e. hydropower production and agriculture) are likely to be negatively impacted by climate change. We discuss the robustness of the estimated impacts to the climate natural variability and the approximations in modeling the physical system and the socio-economic system, and perform an uncertainty analysis of several sources of uncertainty. We demonstrate that the contribution of natural climate uncertainty is rather remarkable and that, among different modelling uncertainty sources, the one from climate modeling is very significant.

Description: Macroinvertebrate community responses to a dewatering disturbance gradient in a restored stream Hydrology and Earth System Sciences, 15, 1771-1783, 2011 Author(s): J. D. Muehlbauer, M. W. Doyle, and E. S. Bernhardt Dewatering disturbances are common in aquatic systems and represent a relatively untapped field of disturbance ecology, yet studying dewatering events along gradients in non-dichotomous (i.e. wet/dry) terms is often difficult. Because many stream restorations can essentially be perceived as planned hydrologic manipulations, such systems can make ideal test-cases for understanding processes of hydrological disturbance. In this study we used an experimental drawdown in a 440 ha stream/wetland restoration site to assess aquatic macroinvertebrate community responses to dewatering and subsequent rewetting. The geomorphic nature of the site and the design of the restoration allowed dewatering to occur predictably along a gradient and decoupled the hydrologic response from any geomorphic (i.e. habitat heterogeneity) effects. In the absence of such heterogeneous habitat refugia, reach-scale wetted perimeter and depth conditions exerted a strong control on community structure. The community exhibited an incremental response to dewatering severity over the course of this disturbance, which was made manifest not as a change in community means but as an increase in community variability, or dispersion, at each site. The dewatering also affected inter-species abundance and distributional patterns, as dewatering and rewetting promoted alternate species groups with divergent habitat tolerances. Finally, our results indicate that rapid rewetting – analogous to a hurricane breaking a summer drought – may represent a recovery process rather than an additional disturbance and that such processes, even in newly restored systems, may be rapid.

Description: A re-evaluation of the Italian historical geomagnetic catalogue: implications for paleomagnetic dating at active Italian volcanoes Solid Earth, 2, 65-74, 2011 Author(s): F. D'Ajello Caracciolo, A. Pignatelli, F. Speranza, and A. Meloni Paleomagnetism is proving to represent one of the most powerful dating tools of volcanics emplaced in Italy during the last few centuries/millennia. This method requires that valuable proxies of the local geomagnetic field (paleo)secular variation ((P)SV) are available. To this end, we re-evaluate the whole Italian geomagnetic directional dataset, consisting of 833 and 696 declination and inclination measurements, respectively, carried out since 1640 AD at several localities. All directions were relocated via the virtual geomagnetic pole method to Stromboli (38.8° N, 15.2° E), the rough centre of the active Italian volcanoes. For declination-only measurements, missing inclinations were derived (always by pole method) by French data (for period 1670–1789), and by nearby Italian sites/years (for periods 1640–1657 and 1790–1962). Using post-1825 declination values, we obtain a 0.46 ± 0.19° yr −1 westward drift of the geomagnetic field for Italy. The original observation years were modified, considering such drift value, to derive at a drift-corrected relocated dataset. Both datasets were found to be in substantial agreement with directions derived from the field models by Jackson et al. (2000) and Pavon-Carrasco et al. (2009). However, the drift-corrected dataset minimizes the differences between the Italian data and both field models, and eliminates a persistent 1.6° shift of 1933–1962 declination values from Castellaccio with respect to other nearly coeval Italian data. The relocated datasets were used to calculate two post-1640 Italian SV curves, with mean directions calculated every 30 and 10 years before and after 1790, respectively. The curve comparison suggests that both available field models yield the best available SV curve to perform paleomagnetic dating of 1600–1800 AD Italian volcanics, while the Italian drift-corrected curve is probably preferable for the 19th century. For the 20th century, the global model by Jackson et al. (2000) yields more accurate inclination values, while the declinations from our drift-corrected curve seem to better represent the local field evolution, at least for the first half of the century.

Description: Carbon monoxide as a tracer for tropical troposphere to stratosphere transport in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development Discussions, 4, 1185-1211, 2011 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, G. Günther, H.-C. Pumphrey, S. Viciani, F. D'Amato, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact on the composition of the tropical lower stratosphere of quasi-horizontal in-mixing into the tropical tropopause layer from the mid-latitude stratosphere. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F, and CO 2 ) in the lower tropical stratosphere. The boundary conditions at the ground are represented for the long-lived trace substances CH 4 , N 2 O, CCl 3 F, and CO 2 based on ground-based measurements. The boundary condition for CO in the free troposphere is deduced from MOPITT measurements. We find that the zonally averaged tropical CO anomaly patterns simulated by this model version of CLaMS are in good agreement with observations. The introduction of a new scheme in the ECMWF integrated forecast system (Tompkins et al., 2007) for the ice supersaturation after September 2006, results in a somewhat less good agreement between observed and simulated CO patterns in the tropical lower stratosphere after this date.

Description: The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations Geoscientific Model Development Discussions, 4, 1213-1271, 2011 Author(s): D. N. Walters, M. J. Best, A. C. Bushell, D. Copsey, J. M. Edwards, P. D. Falloon, C. M. Harris, A. P. Lock, J. C. Manners, C. J. Morcrette, M. J. Roberts, R. A. Stratton, S. Webster, J. M. Wilkinson, M. R. Willett, I. A. Boutle, P. D. Earnshaw, P. G. Hill, C. MacLachlan, G. M. Martin, W. Moufouma-Okia, M. D. Palmer, J. C. Petch, G. G. Rooney, A. A. Scaife, and K. D. Williams We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0. This paper provides a comprehensive technical and scientific description of the GA3.0 and GL3.0 (and related GA3.1 and GL3.1) configurations and presents the results of some initial evaluations of their performance in various applications. It is to be the first in a series of papers describing each subsequent Global Atmosphere release; this will provide a single source of reference for established users and developers as well as researchers requiring access to a current, but trusted, global MetUM setup.

Description: The geomorphic structure of the runoff peak Hydrology and Earth System Sciences, 15, 1853-1863, 2011 Author(s): R. Rigon, P. D'Odorico, and G. Bertoldi This paper develops a theoretical framework to investigate the core dependence of peak flows on the geomorphic properties of river basins. Based on the theory of transport by travel times, and simple hydrodynamic characterization of floods, this new framework invokes the linearity and invariance of the hydrologic response to provide analytical and semi-analytical expressions for peak flow, time to peak, and area contributing to the peak runoff. These results are obtained for the case of constant-intensity hyetograph using the Intensity-Duration-Frequency (IDF) curves to estimate extreme flow values as a function of the rainfall return period. Results show that, with constant-intensity hyetographs, the time-to-peak is greater than rainfall duration and usually shorter than the basin concentration time. Moreover, the critical storm duration is shown to be independent of rainfall return period as well as the area contributing to the flow peak. The same results are found when the effects of hydrodynamic dispersion are accounted for. Further, it is shown that, when the effects of hydrodynamic dispersion are negligible, the basin area contributing to the peak discharge does not depend on the channel velocity, but is a geomorphic propriety of the basin. As an example this framework is applied to three watersheds. In particular, the runoff peak, the critical rainfall durations and the time to peak are calculated for all links within a network to assess how they increase with basin area.

Description: FAMOUS, faster: using parallel computing techniques to accelerate the FAMOUS/HadCM3 climate model with a focus on the radiative transfer algorithm Geoscientific Model Development Discussions, 4, 1273-1303, 2011 Author(s): P. Hanappe, A. Beurivé, F. Laguzet, L. Steels, N. Bellouin, O. Boucher, Y. H. Yamazaki, T. Aina, and M. Allen We have optimised the atmospheric radiation algorithm of the FAMOUS climate model on several hardware platforms. The optimisation involved translating the Fortran code to C and restructuring the algorithm around the computation of a single air column. A task queue and a thread pool are used to distribute the computation to several processors. Finally, four air columns are packed together in a single data structure and computed simultaneously using Single Instruction Multiple Data operations. The modified algorithm runs more than 50 times faster on the CELL's Synergistic Processing Elements than on its main PowerPC processing element. On Intel-compatible processors, the new radiation code runs 4 times faster and on graphics processors, using OpenCL, more than 2.5 times faster, as compared to the original code. Because the radiation code takes more than 60 % of the total CPU time, FAMOUS executes more than twice as fast. Our version of the algorithm returns bit-wise identical results, which demonstrates the robustness of our approach.

Description: Downscaling of surface moisture flux and precipitation in the Ebro Valley (Spain) using analogues and analogues followed by random forests and multiple linear regression Hydrology and Earth System Sciences, 15, 1895-1907, 2011 Author(s): G. Ibarra-Berastegi, J. Saénz, A. Ezcurra, A. Elías, J. Diaz Argandoña, and I. Errasti In this paper, reanalysis fields from the ECMWF have been statistically downscaled to predict from large-scale atmospheric fields, surface moisture flux and daily precipitation at two observatories (Zaragoza and Tortosa, Ebro Valley, Spain) during the 1961–2001 period. Three types of downscaling models have been built: (i) analogues, (ii) analogues followed by random forests and (iii) analogues followed by multiple linear regression. The inputs consist of data (predictor fields) taken from the ERA-40 reanalysis. The predicted fields are precipitation and surface moisture flux as measured at the two observatories. With the aim to reduce the dimensionality of the problem, the ERA-40 fields have been decomposed using empirical orthogonal functions. Available daily data has been divided into two parts: a training period used to find a group of about 300 analogues to build the downscaling model (1961–1996) and a test period (1997–2001), where models' performance has been assessed using independent data. In the case of surface moisture flux, the models based on analogues followed by random forests do not clearly outperform those built on analogues plus multiple linear regression, while simple averages calculated from the nearest analogues found in the training period, yielded only slightly worse results. In the case of precipitation, the three types of model performed equally. These results suggest that most of the models' downscaling capabilities can be attributed to the analogues-calculation stage.

Description: The influence of constrained fossil fuel emissions scenarios on climate and water resource projections Hydrology and Earth System Sciences, 15, 1879-1893, 2011 Author(s): J. D. Ward, A. D. Werner, W. P. Nel, and S. Beecham Water resources planning requires long-term projections of the impact of climate change on freshwater resources. In addition to intrinsic uncertainty associated with the natural climate, projections of climate change are subject to the combined uncertainties associated with selection of emissions scenarios, GCM ensembles and downscaling techniques. In particular, unknown future greenhouse gas emissions contribute substantially to the overall uncertainty. We contend that a reduction in uncertainty is possible by refining emissions scenarios. We present a comprehensive review of the growing body of literature that challenges the assumptions underlying the high-growth emissions scenarios (widely used in climate change impact studies), and instead points to a peak and decline in fossil fuel production occurring in the 21st century. We find that the IPCC's new RCP 4.5 scenario (low-medium emissions), as well as the B1 and A1T (low emissions) marker scenarios from the IPCC's Special Report on Emissions Scenarios are broadly consistent with the majority of recent fossil fuel production forecasts, whereas the medium to high emissions scenarios generally depend upon unrealistic assumptions of future fossil fuel production. We use a simple case study of projected climate change in 2070 for the Scott Creek catchment in South Australia to demonstrate that even with the current suite of climate models, by limiting projections to the B1 scenario, both the median change and the spread of model results are reduced relative to equivalent projections under an unrealistic high emissions scenario (A1FI).

Description: Data-driven catchment classification: application to the pub problem Hydrology and Earth System Sciences, 15, 1921-1935, 2011 Author(s): M. Di Prinzio, A. Castellarin, and E. Toth A promising approach to catchment classification makes use of unsupervised neural networks (Self Organising Maps, SOM's), which organise input data through non-linear techniques depending on the intrinsic similarity of the data themselves. Our study considers ∼300 Italian catchments scattered nationwide, for which several descriptors of the streamflow regime and geomorphoclimatic characteristics are available. We compare a reference classification, identified by using indices of the streamflow regime as input to SOM, with four alternative classifications, which were identified on the basis of catchment descriptors that can be derived for ungauged basins. One alternative classification adopts the available catchment descriptors as input to SOM, the remaining classifications are identified by applying SOM to sets of derived variables obtained by applying Principal Component Analysis (PCA) and Canonical Correlation Analysis (CCA) to the available catchment descriptors. The comparison is performed relative to a PUB problem, that is for predicting several streamflow indices in ungauged basins. We perform an extensive cross-validation to quantify nationwide the accuracy of predictions of mean annual runoff, mean annual flood, and flood quantiles associated with given exceedance probabilities. Results of the study indicate that performing PCA and, in particular, CCA on the available set of catchment descriptors before applying SOM significantly improves the effectiveness of SOM classifications by reducing the uncertainty of hydrological predictions in ungauged sites.

Description: Analyzing numerics of bulk microphysics schemes in Community models: warm rain processes Geoscientific Model Development Discussions, 4, 1403-1434, 2011 Author(s): I. Sednev and S. Menon In the last decade there has been only one study that discussed time integration scheme (TIS) applied to advance governing differential equations in bulk microphysics (BLK) schemes. Recently, Morrison and Gettelman (2008) examine numerical aspects of double-moment BLK scheme with diagnostic treatment of precipitating hydrometeors implemented into Community Atmosphere Model, version 3 (CAM) to find an acceptable level of accuracy and numerical stability. However, stability condition for their explicit non-positive definite TIS was not defined. It is conventionally thought that the Weather Research and Forecasting (WRF) model can be applied for a broad range of spatial scales from large eddy up to global scale simulations if time steps used for model integration satisfy to a certain limit imposed mainly by dynamics. However, numerics used in WRF BLK schemes has never been analyzed in detail. To improve creditability of BLK schemes we derive a general analytical stability and positive definiteness criteria for explicit Eulerian time integration scheme used to advanced finite-difference equations that govern warm rain formation processes in microphysics packages in Community models (CAM and WRF) and define well-behaved, conditionally well-behaved, and non-well-behaved Explicit Eulerian Bulk Microphysics Code (EEBMPC) classes. We highlight that source codes of BLK schemes, originally developed for use in cloud-resolving models, implemented in Community models belong to conditionally well-behaved EEBMPC class and exhibit better performance for finer spatial resolutions when time steps do not exceed seconds or tenths of seconds. For coarser spatial resolutions used in regional and global scale simulations time steps are usually increased from hundredths up to thousands of seconds. This might lead to a degradation of conditionally well-behaved EEBMPCs ability to calculate the amount of precipitation as well as its spatial and temporal distribution since both stability and positive definiteness conditions are not met in the TIS. The correction through the so called “mass conservation” technique commonly used in many models with bulk microphysics is a main characteristic of non-well-behaved EEBMPC, whose utilization leads to erroneous conclusions regarding relative importance of different microphysical processes. Moreover, surface boundary conditions for ocean, land, lake, and sea ice models are dependent on the precipitation and its spatial and temporal distribution. Uncertainties in calculations of temporal and spatial patterns of accumulated precipitation influence the global water cycle. In fact, numerics in non-well-behaved EEBMPCs, which are used in Community Earth System Model, act as a hidden climate forcing agent, if relatively long time steps are used for the host model integration. By analyzing numerics of warm rain processes in EEBMPCs implemented in Community models we provide general guidelines regarding appropriate choice of integration time steps for use in these models.

Description: River flow time series using least squares support vector machines Hydrology and Earth System Sciences, 15, 1835-1852, 2011 Author(s): R. Samsudin, P. Saad, and A. Shabri This paper proposes a novel hybrid forecasting model known as GLSSVM, which combines the group method of data handling (GMDH) and the least squares support vector machine (LSSVM). The GMDH is used to determine the useful input variables which work as the time series forecasting for the LSSVM model. Monthly river flow data from two stations, the Selangor and Bernam rivers in Selangor state of Peninsular Malaysia were taken into consideration in the development of this hybrid model. The performance of this model was compared with the conventional artificial neural network (ANN) models, Autoregressive Integrated Moving Average (ARIMA), GMDH and LSSVM models using the long term observations of monthly river flow discharge. The root mean square error (RMSE) and coefficient of correlation ( R ) are used to evaluate the models' performances. In both cases, the new hybrid model has been found to provide more accurate flow forecasts compared to the other models. The results of the comparison indicate that the new hybrid model is a useful tool and a promising new method for river flow forecasting.

Description: Long term variability of the annual hydrological regime and sensitivity to temperature phase shifts in Saxony/Germany Hydrology and Earth System Sciences, 15, 1819-1833, 2011 Author(s): M. Renner and C. Bernhofer Recently, climatological studies report observational evidence of changes in the timing of the seasons, such as earlier timing of the annual cycle of surface temperature, earlier snow melt and earlier onset of the phenological spring season. Also hydrological studies report earlier timing and changes in monthly streamflows. From a water resources management perspective, there is a need to quantitatively describe the variability in the timing of hydrological regimes and to understand how climatic changes control the seasonal water budget of river basins. Here, the timing of hydrological regimes from 1930–2009 was investigated in a network of 27 river gauges in Saxony/Germany through a timing measure derived by harmonic function approximation of annual periods of runoff ratio series. The timing measure proofed to be robust and equally applicable to both mainly pluvial river basins and snow melt dominated regimes. We found that the timing of runoff ratio is highly variable, but markedly coherent across the basins analysed. Differences in average timing are largely explained by basin elevation. Also the magnitude of low frequent changes in the seasonal timing of streamflow and the sensitivity to the changes in the timing of temperature increase with basin elevation. This sensitivity is in turn related to snow storage and release, whereby snow cover dynamics in late winter explain a large part of the low- and high-frequency variability. A trend analysis based on cumulative anomalies revealed a common structural break around the year 1988. While the timing of temperature shifted earlier by 4 days, accompanied by a temperature increase of 1 K, the timing of runoff ratio within higher basins shifted towards occurring earlier about 1 to 3 weeks. This accelerated and distinct change indicates, that impacts of climate change on the water cycle may be strongest in higher, snow melt dominated basins.

Description: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 2: On-line coupling Geoscientific Model Development Discussions, 4, 1359-1402, 2011 Author(s): A. Kerkweg and P. Jöckel A new, highly flexible model system for the seamless dynamical down-scaling of meteorological and chemical processes from the global to the meso-γ scale is presented. A global model and a cascade of an arbitrary number of limited-area model instances run concurrently in the same parallel environment, in which the coarser grained instances provide the boundary data for the finer grained instances. Thus, disk-space intensive and time consuming intermediate and pre-processing steps are entirely avoided and the time interpolation errors of common off-line nesting approaches are minimised. More specifically, the regional model COSMO of the German Weather Service (DWD) is nested on-line into the atmospheric general circulation model ECHAM5 within the Modular Earth Submodel System (MESSy) framework. ECHAM5 and COSMO have previously been equipped with the MESSy infrastructure, implying that the same process formulations (MESSy submodels) are available for both models. This guarantees the highest degree of achievable consistency, between both, the meteorological and chemical conditions at the domain boundaries of the nested limited-area model, and between the process formulations on all scales. The on-line nesting of the different models is established by a client-server approach with the newly developed Multi-Model-Driver (MMD), an additional component of the MESSy infrastructure. With MMD an arbitrary number of model instances can be run concurrently within the same message passing interface (MPI) environment, the respective coarser model (either global or regional) is the server for the nested finer (regional) client model, i.e., it provides the data required to calculate the initial and boundary fields to the client model. On-line nesting means that the coupled (client-server) models exchange their data via the computer memory, in contrast to the data exchange via files on disk in common off-line nesting approaches. MMD consists of a library (Fortran95 and some parts in C) which is based on the MPI standard and two new MESSy submodels, MMDSERV and MMDCLNT (both Fortran95) for the server and client models, respectively. MMDCLNT contains a further sub-submodel, INT2COSMO, for the interpolation of the coarse grid data provided by the server models (either ECHAM5/MESSy or COSMO/MESSy) to the grid of the respective client model (COSMO/MESSy). INT2COSMO is based on the off-line pre-processing tool INT2LM provided by the DWD. The new achievements allow the setup of model cascades for zooming (down-scaling) from the global scale to the lower edge of the meso-γ scale (≈1 km) with a very high degree of consistency between the different models and between the chemical and meteorological boundary conditions.

Description: Modelling the statistical dependence of rainfall event variables through copula functions Hydrology and Earth System Sciences, 15, 1959-1977, 2011 Author(s): M. Balistrocchi and B. Bacchi In many hydrological models, such as those derived by analytical probabilistic methods, the precipitation stochastic process is represented by means of individual storm random variables which are supposed to be independent of each other. However, several proposals were advanced to develop joint probability distributions able to account for the observed statistical dependence. The traditional technique of the multivariate statistics is nevertheless affected by several drawbacks, whose most evident issue is the unavoidable subordination of the dependence structure assessment to the marginal distribution fitting. Conversely, the copula approach can overcome this limitation, by dividing the problem in two distinct parts. Furthermore, goodness-of-fit tests were recently made available and a significant improvement in the function selection reliability has been achieved. Herein the dependence structure of the rainfall event volume, the wet weather duration and the interevent time is assessed and verified by test statistics with respect to three long time series recorded in different Italian climates. Paired analyses revealed a non negligible dependence between volume and duration, while the interevent period proved to be substantially independent of the other variables. A unique copula model seems to be suitable for representing this dependence structure, despite the sensitivity demonstrated by its parameter towards the threshold utilized in the procedure for extracting the independent events. The joint probability function was finally developed by adopting a Weibull model for the marginal distributions.