ALBERT

Description: The impact of dust on a six-day pulsation of the West African heat low (WAHL) in summertime (14–20 July 2006) is investigated, with convective rainfall and dust bursts being observed over the Sahel at the beginning and end of the episode. Three Meso-NH simulations were designed which differed in their dust representation. All the simulations capture the variation in the WAHL intensity well, including the simulation without any dust effects. This shows the primary role of large-scale forcing on the WAHL pulsation. In spite of additional daytime heating and night-time cooling effects over the Sahara, the simulation with dust climatology resembles the simulation without any dust effects. In contrast, the simulation using a prognostic dust scheme enhances alternating northward advection of warm and dry air and southward advection of cold and wet air associated with the propagation of an African easterly wave, leading to a strengthening of the WAHL variabilities. This study highlights two effects of dust on the WAHL over the Sahara: a so-called direct effect associated with dust radiative heating, which increases the WAHL thickness, and a so-called indirect effect that intensifies both the African easterly jet and a related African easterly wave. Copyright © 2011 Royal Meteorological Society

Description: The Geostationary Earth Radiation Budget Intercomparison of Longwave and Shortwave radiation (GERBILS) was an observational field experiment over North Africa during June 2007. The campaign involved 10 flights by the FAAM BAe-146 research aircraft over southwestern parts of the Sahara Desert and coastal stretches of the Atlantic Ocean. Objectives of the GERBILS campaign included characterisation of mineral dust geographic distribution and physical and optical properties, assessment of the impact upon radiation, validation of satellite remote sensing retrievals, and validation of numerical weather prediction model forecasts of aerosol optical depths (AODs) and size distributions. We provide the motivation behind GERBILS and the experimental design and report the progress made in each of the objectives. We show that mineral dust in the region is relatively non-absorbing (mean single scattering albedo at 550 nm of 0.97) owing to the relatively small fraction of iron oxides present (1–3%), and that detailed spectral radiances are most accurately modelled using irregularly shaped particles. Satellite retrievals over bright desert surfaces are challenging owing to the lack of spectral contrast between the dust and the underlying surface. However, new techniques have been developed which are shown to be in relatively good agreement with AERONET estimates of AOD and with each other. This encouraging result enables relatively robust validation of numerical models which treat the production, transport, and deposition of mineral dust. The dust models themselves are able to represent large-scale synoptically driven dust events to a reasonable degree, but some deficiencies remain both in the Sahara and over the Sahelian region, where cold pool outflow from convective cells associated with the intertropical convergence zone can lead to significant dust production. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Description: We studied the use of IASI data to improve the forecasts of extreme weather events in the Arctic region. For this purpose, the HARMONIE/Norway regional model was used. A set of 366 IASI channels was initially chosen from the ECMWF archived database. Active channels showing the best fit with the analysis system were selected by applying a multi-step monitoring technique. The IASI data were assimilated together with most of the available conventional and operational satellite observations using a three-dimensional variational data assimilation system. Four experiments with cyclic assimilations and subsequent 48-hour forecasts were performed during the IPY-THORPEX campaign period to evaluate the impact of the IASI data and the campaign observations on the hydrostatic HARMONIE/Norway analyses and forecasts. The assessment of the degrees of freedom for signals on the analysis showed that incorporating the IASI data in the assimilation system improved the contribution of the other observations. The utilization of an energy norm-based approach proved the sensitivity of the forecasts to the IASI channels in cases dominated by dynamic instabilities leading to quickly developing weather systems like, for example, polar lows. Comparison of the HARMONIE/Norway forecasts against independent observations and the ECMWF analyses showed a clear positive impact of the IASI data on geopotential fields in mid-troposphere and in the troposphere in general, respectively. We found small but significant positive impact on the temperature and humidity in the lower troposphere. A case-study showed positive impact of IASI radiances on the analysis and forecasts of a polar low. The impact on the forecasts lasted up to 24 hours when extra in situ campaign data were excluded from the analysis, and up to 36 hours when the campaign data were assimilated together with the IASI radiances. Copyright © 2011 Royal Meteorological Society

Description: Domestic wells have high probability of pumping septic tank leachate Hydrology and Earth System Sciences Discussions, 8, 5701-5732, 2011 Author(s): J. E. Horn and T. Harter Onsite wastewater treatment systems such as septic systems are common in rural and semi-rural areas around the world; in the US, about 25–30 % of households are served by a septic system and a private drinking water well. Site-specific conditions and local groundwater flow are often ignored when installing septic systems and wells. Particularly in areas with small lots, thus a high septic system density, these typically shallow wells are prone to contamination by septic system leachate. Typically, mass balance approaches are used to determine a maximum septic system density that would prevent contamination of the aquifer. In this study, we estimate the probability of a well pumping partially septic system leachate. A detailed groundwater and transport model is used to calculate the capture zone of a typical drinking water well. A spatial probability analysis is performed to assess the probability that a capture zone overlaps with a septic system drainfield depending on aquifer properties, lot and drainfield size. We show that a high septic system density poses a high probability of pumping septic system leachate. The hydraulic conductivity of the aquifer has a strong influence on the intersection probability. We conclude that mass balances calculations applied on a regional scale underestimate the contamination risk of individual drinking water wells by septic systems. This is particularly relevant for contaminants released at high concentrations, for substances which experience limited attenuation, and those being harmful even in low concentrations.

Description: Use of ENVISAT ASAR Global Monitoring Mode to complement optical data in the mapping of rapid broad-scale flooding in Pakistan Hydrology and Earth System Sciences Discussions, 8, 5769-5809, 2011 Author(s): D. O'Grady, M. Leblanc, and D. Gillieson Envisat ASAR Global Monitoring Mode (GM) data are used to produce maps of the extent of the flooding in Pakistan which are made available to the rapid response effort within 24 h of acquisition. The high temporal frequency and independence of the data from cloud-free skies makes GM data a viable tool for mapping flood waters during those periods where optical satellite data is unavailable, which may be crucial to rapid response disaster planning, where thousands of lives are affected. Image differencing techniques are used, with pre-flood baseline image backscatter values being deducted from target values to eliminate regions with a permanent flood-like radar response due to volume scattering and attenuation, and to highlight the low response caused by specular reflection by open flood water. The effect of local incidence angle on the received signal is mitigated by ensuring that the deducted image is acquired from the same orbit track as the target image. Poor separability of the water class with land in areas beyond the river channels is tackled using a region-growing algorithm which seeks threshold-conformance from seed pixels at the center of the river channels. The resultant mapped extents are tested against MODIS SWIR data where available, with encouraging results.

Description: Soil weathering rates in 21 catchments of the Canadian Shield Hydrology and Earth System Sciences Discussions, 8, 5743-5768, 2011 Author(s): D. Houle, P. Lamoureux, N. Bélanger, M. Bouchard, C. Gagnon, S. Couture, and A. Bouffard Soil mineral weathering represents an essential source of nutrient base cation (Ca, Mg and K) for forest growth in addition to provide a buffering power against precipitation acidity for soils and surface waters. Weathering rates of base cations were obtained for 21 catchments located within the temperate and the boreal forest of the Canadian Shield with the geochemical model PROFILE. Weathering rates ranged from 0.58 to 4.46 kmol c ha −1 yr −1 and their spatial variation within the studied area was mostly in agreement with spatial variations in soil mineralogy. Weathering rates of Ca and Mg were significantly correlated ( r = 0.80 and 0.64) with their respective lake concentrations. Weathering rates of K and Na did not correlate with lake concentrations of K and Na. The modeled weathering rates for each catchment were also compared with estimations of net catchment exportations. The result show that modeled weathering rates of Ca were not significantly different than the net catchment exportations while modeled weathering rates of Mg were higher by 51 %. Larger differences were observed for K and Na weathering rates that were significantly different than net catchment exportations being 6.9 and 2.2 times higher than net exportations, respectively. The results for K were expected given its high reactivity with biotic compartments and suggest that most of the K produced by weathering reactions was retained within soil catchments and/or above ground biomass. This explanation does not apply to Na, however, which is a conservative element in forest ecosystems because of the insignificant needs of Na for soil microorganisms and above ground vegetations. It raises concern about the liability of the PROFILE model to provide reliable values of Na weathering rates. Overall, we concluded that the PROFILE model is powerful enough to reproduce spatial geographical gradients in weathering rates for relatively large areas as well as adequately predict absolute weathering rates values for the sum of base cations, Ca and Mg.

Description: Extreme mesoscale weather in the Arctic region consists mainly of cases with shallow fronts that often form in the vicinity of the ice-edge and intense storms called polar lows. This article describes high-resolution numerical simulations of a severe weather event that occurred on 1 March 2008 over the Barents Sea. The event was recorded during the IPY–THORPEX field experiments carried out during February and March 2008. The numerical simulations indicated the formation of a low-pressure system over the Barents Sea on 29 February 2008 due to baroclinic instability. On 1 March, the surface low moved onto the sea-ice around Spitsbergen and decayed later on. The conditions that prevailed before the dissipation of the surface low were favourable for the formation of a polar low. Two experiments were performed to test the possibilities of triggering a polar low through certain modifications to the surface conditions. In the first experiment, the sea-ice around Spitsbergen was removed. No polar low developed in this case, since the static stability was too high. In the second experiment, an attempt to reduce the static stability was made by raising the sea-surface temperature by 5 K. The surface low persisted over the Barents Sea area due to the increased surface heating and led to a strong outbreak of Arctic air over the Norwegian Sea on 2 March. The Arctic-air outbreak formed a sharp baroclinic zone which was absent in the control simulation. A secondary mesoscale low was triggered near the baroclinic zone over the Norwegian Sea, which grew into an intense polar low with surface winds reaching hurricane force. Formation of the polar low was due to baroclinic instability, whereas convective instability was important during the growth of the low. Copyright © 2011 Royal Meteorological Society

Description: The role of catchment classification in rainfall-runoff modeling Hydrology and Earth System Sciences Discussions, 8, 6113-6153, 2011 Author(s): Y. He, A. Bárdossy, and E. Zehe A sound catchment classification scheme is a fundamental step towards improved catchment hydrology science and prediction in ungauged basins. Two categories of catchment classification methods are presented in the paper. The first one is based directly on physiographic properties and climatic conditions over a catchment and regarded as a Linnaean type or natural classification scheme. The second one is based on numerical clustering and regionalization methods and considered as a statistical or arbitrary classification scheme. This paper reviews each category including what has been done since recognition of the intrinsic value of catchment classification, what is being done in the current research, as well as what is to be done in the future.

Description: Infiltration-soil moisture redistribution under natural conditions: experimental evidence as a guideline for realizing simulation models Hydrology and Earth System Sciences Discussions, 8, 6199-6225, 2011 Author(s): R. Morbidelli, C. Corradini, C. Saltalippi, A. Flammini, and E. Rossi The evolution in time, t , of the experimental soil moisture vertical profile under natural conditions is investigated in order to address the corresponding simulation modelling. The measurements were conducted in a plot with a bare silty loam soil. The soil water content, θ, was continuously monitored at different depths, z , using a Time Domain Reflectometry (TDR) system. For each profile four buriable three-rod waveguides were inserted horizontally at different depths (5, 15, 25 and 35 cm). In addition, we used sensors of air temperature and relative humidity, wind speed, solar radiation, evaporation and rain as supports for the application of selected simulation models, as well as for the detection of elements leading to their improvement. The results indicate that, under natural conditions, very different trends of the θ( z , t ) function can be observed in the given fine-textured soil, where the formation of a sealing layer over the parent soil requires an adjustment of the simulation modelling commonly used for hydrological applications. In particular, because of the considerable variations in the shape of the moisture content vertical profile as a function of time, a generalization of the existing models should incorporate a representation of the variability in time of the saturated hydraulic conductivity of the uppermost soil. This conclusion is supported by the fact that the observed shape of θ( z ) can be appropriately reproduced by adopting this approach, however the observed rainfall rate and the occurrence of freeze-thaw cycles with high soil moisture contents have to be explicitly incorporated.

Description: Sand box experiments to evaluate the influence of subsurface temperature probe design on temperature based water flux calculation Hydrology and Earth System Sciences Discussions, 8, 6155-6197, 2011 Author(s): M. Munz, S. E. Oswald, and C. Schmidt Quantification of subsurface water fluxes based on the one dimensional solution to the heat transport equation depends on the accuracy of measured subsurface temperatures. The influence of temperature probe setup on the accuracy of vertical water flux calculation was systematically evaluated in this experimental study. Four temperature probe setups were installed into a sand box experiment to measure temporal highly resolved vertical temperature profiles under controlled water fluxes in the range of ±1.3 m d −1 . Pass band filtered time series provided amplitude and phase of the diurnal temperature signal varying with depth depending on water flux. Amplitude ratios of setups directly installed into the saturated sediment significantly varied with sand box hydraulic gradients. Amplitude ratios provided an accurate basis for the analytical calculation of water flow velocities, which matched measured flow velocities. Calculated flow velocities were sensitive to thermal properties of saturated sediment and to probe distance, but insensitive to thermal dispersivity equal to solute dispersivity. Amplitude ratios of temperature probe setups indirectly installed into piezometer pipes were influenced by thermal exchange processes within the pipes and significantly varied with water flux direction only. Temperature time lags of small probe distances of all setups were found to be insensitive to vertical water flux.

Description: Ground-based radar observations at three distinct geographical locations in West Africa along a common latitudinal band (Niamey, Niger (continental), Kawsara, Senegal (coastal), and Praia, Republic of Cape Verde (maritime)) are analyzed to determine convective system characteristics in each domain during a 29-day period in 2006. Ancillary datasets provided by the African Monsoon Multidisciplinary Analyses (AMMA) and NASA-AMMA (NAMMA) field campaigns are also used to place the radar observations in context. Results show that the total precipitation is dominated by propagating mesoscale convective systems. Convective characteristics vary according to environmental properties, such as vertical shear, CAPE, and the degree of synoptic forcing. Data are bifurcated based on the presence or absence of African easterly waves. In general, African easterly waves appear to enhance mesoscale convective system strength characteristics (e.g. total precipitation and vertical reflectivity profiles) at the inland and maritime sites. The wave regime also resulted in an increased population of the largest observed mesoscale convective systems observed near the coast, which led to an increase in stratiform precipitation. Despite this increase, differentiation of convective strength characteristics was less obvious between wave and no-wave regimes at the coast. Owing to the propagating nature of these advecting mesoscale convective systems, interaction with the regional thermodynamic and dynamic environment appears to result in more variability than enhancements due to the wave regime, independent of location. Copyright © 2011 Royal Meteorological Society

Description: A wavelet formulation on the sphere is being considered for modelling heterogeneous background-error correlations for the Météo-France global numerical weather prediction (NWP) model. This approach is compared with the operational spectral formulation, which is horizontally homogeneous to a large extent. Diagnostic studies have been conducted to examine geographical variations of three-dimensional correlations over the whole globe. Results indicate that the contrast between relatively broad horizontal correlations in the Tropics and sharp ones in midlatitudes is well represented by the wavelet formulation. Heterogeneities in vertical correlations are also better captured in the wavelet approach than in the spectral one, with visible changes as functions of e.g. latitude and land/sea contrasts. In addition, wavelet-based correlation estimates are shown to be partly sensitive to the choice of the calibration period. The impact of the wavelet formulation on the forecast quality has been investigated during a three-week calibration period, and also during the following three weeks. While the impact of the wavelet formulation is globally positive during the two periods, it tends to be more spectacular during the calibration time interval, as expected. These results indicate that an on-line calibration should be considered in the future, in order to exploit fully the ability of wavelets to extract correlation heterogeneities from ensemble data. Copyright © 2011 Royal Meteorological Society

Description: FY-3A, launched in May 2008, is the first in a series of seven polar-orbiting meteorological satellites due to be launched by China's Meteorological Administration in the period leading up to 2020. The FY-3A payload includes four instruments of particular interest for numerical weather prediction (NWP): microwave temperature and humidity sounders, a microwave imager, and an infrared sounder. The main features of these instruments are described. Data from the calibration–validation phase of the FY-3A mission were introduced into the ECMWF Integrated Forecasting System in order to assess the data quality and the influence of the data on analyses and forecasts. An analysis of first-guess departures has shown the data to be of good quality overall. Several issues with instrument performance and ground segment processing have been identified. The most serious of these are: uncertainties in the temperature sounder passbands on-orbit, orbital biases in the infrared instrument affecting the highest peaking channels, and scan biases in the microwave humidity sounder. Variational bias correction partially corrects for these errors, but more work remains to be done to correct the problems before the full benefit of the data is realised. In observing system experiments, the FY-3A instruments, both individually and as a package, show considerable skill when added to observation depleted control experiments. When added to a full observing system, the impacts are neutral to slightly positive, as expected. These initial results are encouraging and build confidence that the following series of FY-3 instruments will be widely used in NWP data assimilation systems. Copyright © 2011 Royal Meteorological Society

Description: The use of LIDAR as a data source for digital elevation models – a study of the relationship between the accuracy of digital elevation models and topographical attributes in northern peatlands Hydrology and Earth System Sciences Discussions, 8, 5497-5522, 2011 Author(s): A. Hasan, P. Pilesjö, and A. Persson It is important to study the factors affecting estimates of wetness since wetness is crucial in climate change studies. The availability of digital elevation models (DEMs) generated with high resolution data is increasing, and their use is expanding. LIDAR earth elevation data have been used to create several DEMs with different resolutions, using various interpolation parameters, in order to compare the models with collected surface data. The aim is to study the accuracy of DEMs in relation to topographical attributes such as slope and drainage area, which are normally used to estimate the wetness in terms of topographic wetness indices. Evaluation points were chosen from the high-resolution LIDAR dataset at a maximum distance of 10 mm from the cell center for each DEM resolution studied, 0.5, 1, 5, 10, 30 and 90 m. The interpolation method used was inverse distance weighting method with four search radii: 1, 2, 5 and 10 m. The DEM was evaluated using a quantile-quantile test and the normalized median absolute deviation. The accuracy of the estimated elevation for different slopes was tested using the DEM with 0.5 m resolution. Drainage areas were investigated at three resolutions, with coinciding evaluation points. The ability of the model to generate the drainage area at each resolution was obtained by pairwise comparison of three data subsets. The results show that the accuracy of the elevations obtained with the DEM model are the same for different resolutions, but vary with search radius. The accuracy of the values (NMAD of errors) varies from 29.7 mm to 88.9 mm, being higher for flatter areas. It was also found that the accuracy of the drainage area is highly dependent on DEM resolution. Coarse resolution yielded larger estimates of the drainage area but lower slope values. This may lead to overestimation of wetness values when using a coarse resolution DEM.

Description: Multivariate design via Copulas Hydrology and Earth System Sciences Discussions, 8, 5523-5558, 2011 Author(s): G. Salvadori, C. De Michele, and F. Durante Calculating return periods and design quantiles in a multivariate framework is a difficult problem: essentially, this is due to the lack of a natural total order in multi-dimensional Euclidean spaces. This paper tries to make the issue clear. First, we outline a possible way to introduce a coherent notion of multivariate total order, and discuss its consequences on the calculation of multivariate return period: in particular, the latter is based on Copulas and the Kendall's measure, which provides a consistent notion of multivariate quantile. Secondly, we introduce several approaches for the identification of critical design events: these latter quantities are of utmost importance in practical applications, but their calculation is yet limited, due to the lack of a suitable theoretical setting where to embed the problem. Throughout the paper, a case study involving the behavior of a dam is used to illustrate the new concepts outlined in this work.

Description: Influence of initial heterogeneities and recharge limitations on the evolution of aperture distributions in carbonate aquifers Hydrology and Earth System Sciences Discussions, 8, 5631-5666, 2011 Author(s): B. Hubinger and S. Birk Karst aquifers evolve where the dissolution of soluble rocks causes the enlargement of discrete pathways along fractures or bedding planes, thus creating highly conductive solution conduits. To identify general interrelations between hydrogeological conditions and the properties of the evolving conduit systems the aperture-size frequency distributions resulting from generic models of conduit evolution are analysed. For this purpose, a process-based numerical model coupling flow and rock dissolution is employed. Initial protoconduits are represented by tubes with log-normally distributed aperture sizes with a mean of 0.5 mm. Apertures are spatially uncorrelated and widen up to the metre range due to dissolution by chemically aggressive waters. Several examples of conduit development are examined focussing on influences of the initial heterogeneity and the available amount of recharge. If the available recharge is sufficiently high the evolving conduits compete for flow and those with large apertures and high hydraulic gradients attract more and more water. As a consequence, the positive feedback between increasing flow and dissolution causes the breakthrough of a conduit pathway connecting the recharge and discharge sides of the modelling domain. Under these competitive flow conditions dynamically stable bimodal aperture distributions are found to evolve, i.e. a certain percentage of tubes continues to be enlarged while the remaining tubes stay small-sized. The percentage of strongly widened tubes is found to be independent of the breakthrough time and decreases with increasing heterogeneity of the initial apertures and decreasing amount of available water. If the competition for flow is suppressed because the availability of water is strongly limited breakthrough of a conduit pathway is inhibited and the conduit pathways widen very slowly. The resulting aperture distributions are found to be unimodal covering some orders of magnitudes in size. Under these suppressed flow conditions the entire range of apertures continues to be enlarged. Hence, the number of tubes reaching aperture sizes in the order of centimetres or decimetres continues to increase with time and in the long term may exceed the number of large-sized tubes evolving under competitive flow conditions. This suggests that conduit development under suppressed flow conditions may significantly enhance the permeability of the formation e.g. in deep-seated carbonate settings.

Description: Skewness as measure of the invariance of instantaneous renormalized drop diameter distributions – Part 1: Convective vs. stratiform precipitation Hydrology and Earth System Sciences Discussions, 8, 5605-5629, 2011 Author(s): M. Ignaccolo and C. De Michele We investigate the variability of the instantaneous distribution shape of the renormalized drop diameter making use of the third order central moment: the skewness . Disdrometer data, collected at Darwin Australia, are considered either as whole or as divided in convective and stratiform precipitation intervals. We show that in all cases the distribution of the skewness is strongly peaked around 0.64. This allows to identify a most common distribution of renormalized drop diameters and two main variations, one with larger and one with smaller skewness. The distributions' shapes are independent from the stratiform vs. convective classification.

Description: The sensitivity of land emissivity estimates from AMSR-E at C and X bands to surface properties Hydrology and Earth System Sciences Discussions, 8, 5667-5699, 2011 Author(s): H. Norouzi, M. Temimi, W. B. Rossow, C. Pearl, M. Azarderakhsh, and R. Khanbilvardi Microwave observations at low frequencies exhibit more sensitivity to surface and subsurface properties with little interference from the atmosphere. The objective of this study is to develop a global land emissivity product using passive microwave observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and to investigate its sensitivity to land surface properties. The developed product complements existing land emissivity products from SSM/I and AMSU by adding land emissivity estimates at two lower frequencies, 6.9 and 10.65 GHz (C- and X-band, respectively). Observations at these low frequencies penetrate deeper into the soil layer. Ancillary data used in the analysis, such as surface skin temperature and cloud mask, are obtained from International Satellite Cloud Climatology Project (ISCCP). Atmospheric properties are obtained from the TIROS Operational Vertical Sounder (TOVS) observations to determine the small upwelling and downwelling atmospheric emissions as well as the atmospheric transmission. A sensitivity test confirms the small effect of the atmosphere but shows that skin temperature accuracy can significantly affect emissivity estimates. Retrieved emissivities at C- and X-bands and their polarization differences exhibit similar patterns of variation with changes in land cover type, soil moisture, and vegetation density as seen at SSM/I-like frequencies (Ka and Ku bands). The emissivity maps from AMSR-E at these higher frequencies agree reasonably well with the existing SSM/I-based product. The inherent but small discrepancy introduced by the difference between SSM/I and AMSR-E frequencies and incidence angles has been examined and found to be small. Large differences between emissivity estimates from ascending and descending overpasses were found at the lower frequencies due to the inconsistency between the thermal IR skin temperatures and passive microwave brightness temperatures which can come from below the surface. This issue must be addressed in future studies to improve the accuracy of the emissivity estimates at lower frequencies.

Description: Influence of soil parameters on the skewness coefficient of the annual maximum flood peaks Hydrology and Earth System Sciences Discussions, 8, 5559-5604, 2011 Author(s): A. Gioia, V. Iacobellis, S. Manfreda, and M. Fiorentino Understanding the spatial variability of key parameters of flood probability distributions represents a strategy to provide insights on hydrologic similarity and building probabilistic models able to reduce the uncertainty in flood prediction in ungauged basins. In this work, we exploited the theoretically derived distribution of floods TCIF (Gioia et al., 2008), based on two different threshold mechanisms associated respectively to ordinary and extraordinary events. The model is based on the hypotheses that ordinary floods are generally due to rainfall events exceeding a threshold infiltration rate in a small source area, while the so-called outlier events, responsible of the high skewness of flood distributions, are triggered when severe rainfalls exceed a storage threshold over a large portion of the basin. Within this scheme, a sensitivity analysis was performed in order to analyze the effects of climatic and geomorphologic parameters on the skewness coefficient. In particular, the analysis was conducted investigating the influence on flood distribution of physical factors such as rainfall intensity, soil infiltration capacity, and basin area, in order to provide insights in catchment classification and process conceptualization.

Description: Snow particle size distributions (particle size 〉400 µm) in the western Arctic measured with in situ aircraft instrumentation during the Surface Heat Budget of the Arctic/First ISCCP Regional Experiment - Arctic Clouds Experiment and Mixed-Phase Arctic Cloud Experiment are analysed. Three cases of shallow, precipitating mixed-phase boundary-layer clouds and two cases of deep, precipitating frontal clouds are examined. Overall, the shallow cases had much lower values of particle concentration and ice water content than the deep cases, indicating large differences in ice initiation and growth between these regimes. Within a given case for both the shallow and deep frontal systems, and for the dataset as a whole, crystal concentration had little correlation with temperature (height), despite an active aggregation process that was indicated by large aggregates (〉5 mm) observed in four out of the five cases. Exponential size distributions are fitted to the observations, allowing a direct comparison with the snow particle size distributions that are represented with exponential functions in many bulk microphysics schemes used in weather and climate models. Values of the fitted intercept parameter N 0 are generally 2–10 times smaller for the shallow compared to the deep frontal cases as a result of differences in crystal concentration between these regimes. Values of N 0 ∼ 10 7 m −4 specified for snow in many bulk microphysics schemes are broadly consistent with fitted N 0 for the deep cases but larger than values for the shallow cases. The deep frontal cases also exhibit a relationship between N 0 and temperature consistent with previous observations of midlatitude frontal systems. However, there are no consistent differences in N 0 between the shallow and deep cases when partitioned by ice water content. Fitted values of slope parameter λ for the shallow and deep cases are generally consistent with previous studies of lower-latitude cloud systems. Copyright © 2011 Royal Meteorological Society

Description: A bulk flux algorithm predicts the turbulent surface fluxes of momentum and sensible and latent heat from mean measured or modelled meteorological variables. The bulk flux algorithm resulting from data collected over winter sea ice during SHEBA, the experiment to study the Surface Heat Budget of the Arctic Ocean, failed, however, in its first trial to predict the turbulent momentum flux over sea ice in the Antarctic. This result suggests that the main parameter for predicting the momentum flux, the aerodynamics roughness length z 0 , does not respond just to the friction velocity, as in the SHEBA algorithm, but is closely related to the physical roughness of snow-covered sea ice and may need to be site-specific. I investigate this idea with simultaneous measurements of z 0 and the physical roughness of the surface, ξ , at Ice Station Weddell. The metric ξ derives from surveys of surface elevation and is related to but always less than the standard deviation in surface elevation. On combining the z 0 – ξ pairs from Ice Station Weddell with similar data obtained over Arctic sea ice, I show that the Arctic and Antarctic z 0 – ξ data lie along a continuum such that measuring ξ could provide a means for estimating a site-specific z 0 for any global sea ice surface. Backscatter data from satellite-borne synthetic aperture radar might provide a remotely sensed estimate of ξ . Copyright © 2011 Royal Meteorological Society

Description: Technical note: Towards a continuous classification of climate using bivariate colour mapping Hydrology and Earth System Sciences Discussions, 8, 5733-5742, 2011 Author(s): A. J. Teuling Climate is often defined in terms of discrete classes. Here I use bivariate colour mapping to show that the global distribution of Köppen-Geiger climate classes can largely be reproduced by combining the simple means of two key states of the climate system (i.e., air temperature and relative humidity). This allows for a classification that is not only continuous in space, but can be applied at and transferred between timescales ranging from minutes to decades.

Description: The importance of parameter resampling for soil moisture data assimilation into hydrologic models using the particle filter Hydrology and Earth System Sciences Discussions, 8, 5849-5890, 2011 Author(s): D. A. Plaza, R. De Keyser, G. J. M. De Lannoy, L. Giustarini, P. Matgen, and V. R. N. Pauwels The Ensemble Kalman filter (EnKF) and the Sequential Importance Resampling (SIR) particle filter are evaluated for their performance in soil moisture assimilation and the consequent effect on discharge. With respect to the resulting soil moisture time series, both filters perform similarly. However, both filters have a negative effect on the discharge due to inconsistency between the parameter values and the states after the assimilation. In order to overcome this inconsistency, parameter resampling is applied along with the SIR filter, to obtain consistent parameter values with the analyzed soil moisture state. Extreme parameter replication, which could lead to a particle collapse, is avoided by the perturbation of the parameters with white noise. Both the modelled soil moisture and discharge are improved if the complementary parameter resampling is applied. The SIR filter with parameter resampling offers an efficient way to deal with biased observations. The robustness of the methodology is evaluated for 3 model parameter sets and 3 assimilation frequencies. Overall, the results in this paper indicate that the particle filter is a promising tool for hydrologic modelling purposes, but that an additional parameter resampling may be necessary to consistently update all state variables and fluxes within the model.

Description: Spatial moments of catchment rainfall: rainfall spatial organisation, basin morphology, and flood response Hydrology and Earth System Sciences Discussions, 8, 5811-5847, 2011 Author(s): D. Zoccatelli, M. Borga, A. Viglione, G. B. Chirico, and G. Blöschl This paper provides a general analytical framework for assessing the dependence existing between spatial rainfall organisation, basin morphology and runoff response. The analytical framework builds upon a set of spatial rainfall statistics (termed " spatial moments of catchment rainfall ") which describe the spatial rainfall organisation in terms of concentration and dispersion statistics as a function of the distance measured along the flow routing coordinate. The introduction of these statistics permits derivation of a simple relationship for the quantification of storm velocity at the catchment scale. The paper illustrates the development of the analytical framework and explains the conceptual meaning of the statistics by means of application to five extreme flash floods occurred in various European regions in the period 2002–2007. High resolution radar rainfall fields and a distributed hydrologic model are employed to examine how effective are these statistics in describing the degree of spatial rainfall organisation which is important for runoff modelling. This is obtained by quantifying the effects of neglecting the spatial rainfall variability on flood modelling, with a focus on runoff timing. The size of the study catchments ranges between 36 to 982 km 2 . The analysis reported here shows that the spatial moments of catchment rainfall can be effectively employed to isolate and describe the features of rainfall spatial organization which have significant impact on runoff simulation. These statistics provide essential information on what space-time scales rainfall has to be monitored, given certain catchment and flood characteristics, and what are the effects of space-time aggregation on flood response modeling.

Description: Evaluation of the transferability of hydrological model parameters for simulations under changed climatic conditions Hydrology and Earth System Sciences Discussions, 8, 5891-5915, 2011 Author(s): S. Bastola, C. Murphy, and J. Sweeney Conceptual hydrological models are widely used for climate change impact assessment. The implicit assumption in most such work is that the parameters estimated from observations remain valid for future climatic conditions. This paper evaluates a simple threshold based approach for testing this assumption, where a set of behavioural simulators are identified for different climatic conditions for the future simulation i.e. wet, average and dry conditions. These simulators were derived using three different data sets that are generated by sampling a block of one year of data without replacement from the observations such that they define the different climatic conditions. The simulators estimated from the wet climatic data set showed the tendency to underestimate flow when applied to dry data set and vice versa. However, the performances of the three sets of basin simulators on chronologically coherent data are identical to the simulators identified from a sufficiently long data series that contains both wet and dry climatic conditions. The results presented suggest that the issue of time invariance in the value of parameters has a minimal effect on the simulation if the change in precipitation is less than 10 % of the data used for calibration.

Description: Mass transfer effects in 2-D dual-permeability modeling of field preferential bromide leaching with drain effluent Hydrology and Earth System Sciences Discussions, 8, 5917-5967, 2011 Author(s): H. H. Gerke, J. Dusek, and T. Vogel Subsurface drained experimental fields are frequently used for studying preferential flow (PF) in structured soils. Considering two-dimensional (2-D) transport towards the drain, however, the relevance of mass transfer coefficients, apparently reflecting small-scale soil structural properties, for the water and solute balances of the entire drained field is largely unknown. This paper reviews and analyzes effects of mass transfer reductions on Br − leaching for a subsurface drained experimental field using a numerical 2-D dual-permeability model (2D-DPERM). The sensitivity of the "diffusive" mass transfer component on bromide (Br − ) leaching patterns is discussed. Flow and transport is simulated in a 2-D vertical cross-section using parameters, boundary conditions (BC), and data of a Br − tracer irrigation experiment on a subsurface drained field (5000 m 2 area) at Bokhorst (Germany), where soils have developed from glacial till sediments. The 2D-DPERM simulation scenarios assume realistic irrigation and rainfall rates, and Br-application in the soil matrix (SM) domain. The mass transfer reduction controls preferential tracer movement and can be related to physical and chemical properties at the interface between flow path and soil matrix in structured soil. A reduced solute mass transfer rate coefficient allows a better match of the Br − mass flow observed in the tile drain discharge. The results suggest that coefficients of water and solute transfer between PF and SM domains have a clear impact on Br − effluent from the drain. Amount and composition of the drain effluent is analyzed as a highly complex interrelation between temporally and spatially variable mass transfer in the 2-D vertical flow domain that depends on varying "advective" and "diffusive" transfer components, the spatial distribution of residual tracer concentrations, and the lateral flow fields in both domains from plots of the whole subsurface drained field. The local-scale soil structural effects (e.g., such as macropore wall coatings), here conceptualized as changes in mass transfer coefficients, can have a clear effect on leaching at the plot and field-scales.

Description: Evaluation and bias correction of satellite rainfall data for drought monitoring in Indonesia Hydrology and Earth System Sciences Discussions, 8, 5969-5997, 2011 Author(s): R. R. E. Vernimmen, A. Hooijer, Mamenun, and E. Aldrian The accuracy of satellite rainfall data from different sources, TRMM 3B42RT, CMORPH and PERSIANN, was investigated through comparison with reliable ground station rainfall data in Indonesia, with a focus on their ability to detect patterns of low rainfall that may lead to drought conditions. It was found that all sources underestimated rainfall in dry season months. The CMORPH and PERSIANN data differed most from ground station data and are also very different from the TRMM data. However, it proved possible to improve TRMM data to yield sufficiently accurate estimates, both for dry periods ( R 2 0.65–0.92) and annually ( R 2 0.84–0.96), applying a single parameterized bias correction equation that is constant in space and time. It is proposed that these bias corrected TRMM data be used in real-time drought monitoring, in Indonesia and probably in other countries where similar conditions exist. This will yield major advantages, in terms of accuracy, spatial coverage, timely availability and cost efficiency, over drought monitoring with only ground stations.

Description: Effects of antecedent soil moisture on runoff modeling in small semiarid watersheds of southeastern Arizona Hydrology and Earth System Sciences Discussions, 8, 6227-6256, 2011 Author(s): Y. Zhang, H. Wei, and M. A. Nearing Antecedent soil moisture prior to a rain event influences the rainfall-runoff relationship. Very few studies have looked at the effects of antecedent soil moisture on runoff modeling sensitivities in arid and semi-arid areas. This study examines the influence of initial soil moisture on model runoff prediction capability in small semiarid watersheds using model sensitivity and by comparing the use of antecedent vs. average long term soil water content for defining the model initial conditions for the modified Green-Ampt Mein-Larson model within the Rangeland Hydrology and Erosion Model (RHEM). Measured rainfall, runoff, and soil moisture data from four semiarid rangeland watersheds ranging in size from 0.34 to 4.53 ha on the Walnut Gulch Experimental Watershed in southeastern Arizona, USA, were used. Results showed that: (a) there were no significant correlations between measured runoff ratio and antecedent soil moisture in any of the four watersheds; (b) average sensitivities of simulated runoff amounts and peaks to antecedent soil moisture were 0.05 mm and 0.18 mm h −1 , respectively, with each 1 % change in antecedent soil moisture; (c) runoff amounts and peaks simulated with long term average soil moisture were statistically equivalent to those simulated with measured antecedent soil moisture. The relative lack of sensitivity of modeled runoff to antecedent soil moisture in this case is contrary to results reported in other studies, and is largely due to the fact that the surface soil is nearly always very dry in this environment.

Description: During 3–4 March 2008, the Norwegian IPY-THORPEX field campaign successfully carried out three flight missions that observed the full life cycle of a polar low over the Norwegian Sea. Here the three-dimensional structure of the polar low has been investigated using dropsonde data from the three flights. The polar low developed in a cold air outbreak, with temperature differences between the sea surface and 500 hPa of about 45–50°C. Cross-sections show that the horizontal gradients of potential temperature weakened as the polar low matured, suggesting that baroclinic energy conversion took place. Dropsonde data of potential temperature and relative humidity show evidence of a tropopause fold, which is possibly a manifestation of upper-level forcing. This is corroborated by potential vorticity inversion, which shows a dominant role of upper-level forcing throughout the polar low's lifetime. During the cyclogenesis stage the polar low circulation was confined below 700 hPa, with a northerly low-level jet of 26 m s −1 . In the mature stage, its circulation reached up to the tropopause (∼450 hPa), with maximum wind speed between 700 and 900 hPa of about 26–28 m s −1 . At this stage the polar low warm core was about 3 K warmer than surrounding air masses. The deep moist towers at the eye-like structure of the polar low extended up to the tropopause with relative humidity values above 70%, indicating a possibly important role for condensational heating in the development. Estimates of surface fluxes of sensible and latent heat using temperature and moisture from the dropsonde data show latent heat fluxes west of the polar low increasing from 175 to 300 W m −2 as the low matured, while the sensible heat fluxes rose from 200 to 280 W m −2 , suggesting a gradually increasing contribution of surface fluxes to the energetics of the polar low with time. Copyright © 2011 Royal Meteorological Society

Description: We present simulations of London's meteorology using the Met Office Unified Model with a new, sophisticated surface energy-balance scheme to represent the urban surfaces, called MORUSES. Simulations are performed with the urban surfaces represented and with the urban surfaces replaced with grass in order to calculate the urban increment on the local meteorology. The local urban effects were moderated to some extent by the passage of an onshore flow that propagated up the Thames estuary and across the city, cooling London slightly in the afternoon. Validations of screen-level temperature show encouraging agreement to within 1–2 K, when the urban increment is up to 5 K. The model results are then used to examine factors shaping the spatial and temporal structure of London's atmospheric boundary layer. The simulations reconcile the differences in the temporal evolution of the urban heat island (UHI) shown in various studies and demonstrate that the variation of UHI with time depends strongly on the urban fetch. The UHI at a location downwind of the city centre shows a decrease in UHI during the night, while the UHI at the city centre stays constant. Finally, the UHI at a location upwind of the city centre increases continuously. The magnitude of the UHI by the time of the evening transition increases with urban fetch. The urban increments are largest at night, when the boundary layer is shallow. The boundary layer experiences continued warming after sunset, as the heat from the urban fabric is released, and a weakly convective boundary layer develops across the city. The urban land-use fraction is the dominant control on the spatial structure in the sensible heat flux and the resulting urban increment, although even the weak advection present in this case study is sufficient to advect the peak temperature increments downwind of the most built-up areas. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Description: Warm conveyor belts (WCBs) are key flow structures associated with extratropical cyclones. They transport moist air from the cyclone's warm sector poleward and upward close to the tropopause level, leading to the formation of elongated cloud bands, intense latent heating and surface precipitation. In this study a comprehensive dataset of airborne lidar observations of moisture and wind from different campaigns has been investigated with a trajectory-based approach to identify ‘lucky encounters’ with WCBs. On 19 July 2007, an upstream flight over the Iberian Peninsula during the European THORPEX Regional Campaign (ETReC 2007) in Central Europe intersected two WCBs: one in the upper tropospheric outflow region about 3 days after starting the ascent, and the other one in the boundary layer inflow region over Spain just prior to the strong ascent. Comparison of the lidar humidity measurements with analysis fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) reveals significant positive deviations, equivalent to an overestimation of the modelled humidity, in this low-tropospheric WCB inflow region (of about 1 g kg −1 (14%) on average and with peak deviations up to 7 g kg −1 ). It is noteworthy that this substantial bias occurs in a potentially dynamically highly relevant air mass that will be subsequently lifted within a WCB to the upper troposphere. A Lagrangian moisture source diagnostic reveals that these large moisture deviations occur within air masses that, according to the ECMWF analyses, are coherently transported from the western Mediterranean towards Spain and experience intense moisture uptake over the Ebro valley. It is suggested that inaccuracies in surface evapotranspiration, horizontal moisture advection, and turbulent vertical transport of moisture in the atmospheric boundary layer potentially contribute to the erroneous low-tropospheric humidity in the inflow region of this particular summertime WCB over Spain in the ECMWF analyses. Copyright © 2011 Royal Meteorological Society

Description: The role of the representation of deep convection on key elements of the West African summer monsoon climate is addressed using the Regional Climate Model RegCM3. Two simulations in which a scheme of deep convection is activated and then turned off are performed and intercompared. Results show that the presence of deep convective heating along the intertropical convergence zone sustains increased lower-level baroclinicity favoring intensification of the jet core and leading to a more realistic African easterly jet. In addition, although the isentropic potential vorticity (IPV) is lower when the convection scheme is switched off, African easterly waves (AEWs) are still generated and propagate westwards but they dissipate around the west coast. Substantial differences between the two simulations occur mainly at the 6- to 9-day time-scale over land, when much weaker activity is simulated in the absence of convection. This indicates that orographic friction and low-level large-scale moisture convergence, generating high values of latent heat and IPV, may play the dominant role in the genesis and growth of AEWs and that deep convection acts to strengthen the overall wave activity and to favor their west coast development. Copyright © 2011 Royal Meteorological Society

Description: A unique airborne differential absorption lidar (DIAL) for water vapour observations was developed at the Deutsches Zentrum für Luft- und Raumfahrt (DLR). Installed on board the DLR Falcon 20 aircraft, the system measured a dataset of about 3900 water vapour profiles during the T-PARC field campaign. These high-resolution humidity observations were assimilated into the European Centre for Medium-Range Weather Forecasts (ECMWF) global model using a version of the operational four-dimensional variational data assimilation system. The assimilation system is able to extract the information for DIAL observations, and verification with independent dropsonde observations shows a reduction in the analysis error when DIAL water vapour observations are assimilated. The forecast influence of the humidity observations is found to be small in most cases, but the observations are able to affect the forecast considerably under certain conditions. Systematic errors are investigated by comparison between humidity model fields, DIAL and dropsonde observations. Overall, DIAL observations are roughly 7–10% drier than model fields throughout the troposphere. Comparison with dropsonde observations suggests that the DIAL observations are too dry in the lower troposphere but not above it. Copyright © 2011 Royal Meteorological Society

Description: In this paper we provide an overview of various satellite products over the Sahara Desert that were available during the GERBILS field campaign. Our results indicate that all mid-visible satellite aerosol optical depth (AOD) products match well with AERONET retrievals. For low AOD (AOD 〈 1), the satellite AODs compare well with aircraft AOD values but they tend to underestimate at high AOD values. We then assessed the satellite products in 0.5 × 0.5 degree grids for the entire study region (10–30°N and 20°W–10°E). If we use a multi-angle imaging spectroradiometer (MISR) as a benchmark for AOD retrievals over bright targets, the estimated AOD derived from the ozone-monitoring instrument aerosol index–MISR relationship performs best when compared with MISR for the entire study region. Although differences exist among satellite products, the advancement in satellite retrieval techniques now provide AOD retrievals over bright targets such as deserts, which are useful for numerical modeling simulation comparisons and other studies. Furthermore, the in situ information from aircraft and the ground continue to provide valuable information for validating satellite products and for assessing their strengths and weaknesses. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, theMet Office

Description: Integrated versus isolated scenario for prediction dissolved oxygen at progression of water quality monitoring stations Hydrology and Earth System Sciences Discussions, 8, 6069-6112, 2011 Author(s): A. A. Najah, A. El-Shafie, O. A. Karim, and O. Jaafar This study examined the potential of Multi-layer Perceptron Neural Network (MLP-NN) in predicting dissolved oxygen (DO) at Johor River Basin. The river water quality parameters were monitored regularly each month at four different stations by the Department of Environment (DOE) over a period of ten years, i.e. from 1998 to 2007. The following five water quality parameters were selected for the proposed MLP-NN modelling, namely; temperature (Temp), water pH, electrical conductivity (COND), nitrate (NO 3 ) and ammonical nitrogen (NH 3 –NL). In this study, two scenarios were introduced; the first scenario (Scenario 1) was to establish the prediction model for DO at each station based on five input parameters, while the second scenario (Scenario 2) was to establish the prediction model for DO based on the five input parameters and DO predicted at previous station (upstream). The model needs to verify when output results and the observed values are close enough to satisfy the verification criteria. Therefore, in order to investigate the efficiency of the proposed model, the verification of MLP-NN based on collection of field data within duration 2009–2010 is presented. To evaluate the effect of input parameters on the model, the sensitivity analysis was adopted. It was found that the most effective inputs were oxygen-containing (NO 3 ) and oxygen demand (NH 3 –NL). On the other hand, Temp and pH were found to be the least effective parameters, whereas COND contributed the lowest to the proposed model. In addition, 17 neurons were selected as the best number of neurons in the hidden layer for the MLP-NN architecture. To evaluate the performance of the proposed model, three statistical indexes were used, namely; Coefficient of Efficiency (CE), Mean Square Error (MSE) and Coefficient of Correlation (CC). A relatively low correlation between the observed and predicted values in the testing data set was obtained in Scenario 1. In contrast, high coefficients of correlation were obtained between the observed and predicted values for the test sets of 0.98, 0.96 and 0.97 for all stations after adopting Scenario 2. It appeared that the results for Scenario 2 were more adequate than Scenario 1, with a significant improvement for all stations ranging from 4 % to 8 %.

Description: The causes of flow regime shifts in the semi-arid Hailiutu River, Northwest China Hydrology and Earth System Sciences Discussions, 8, 5999-6030, 2011 Author(s): Z. Yang, Y. X. Zhou, J. Wenninger, and S. Uhlenbrook Identifying the causes (climate vs. human activities) for hydrological variability is a major challenge in hydrology. This paper examines the flow regime shifts, changes in the climatic variables such as precipitation, evaporation, temperature, and crop area in the semi-arid Hailiutu catchment in the middle section of the Yellow River by performing several statistical analyses. The Pettitt test, cumulative sum charts (CUSUM), regime shift index (RSI) method, and harmonic analysis were carried out on annual, monthly, and daily discharges. Four major shifts in the flow regime have been detected in 1968, 1986, 1992 and 2001. Characteristics of the flow regime were analyzed in the five periods: 1957–1967, 1968–1985, 1986–1991, 1992–2000, and 2001–2007. From 1957 to 1967, the flow regime reflects quasi natural conditions with high variability and larger amplitude of 6 months periodic fluctuations. The river flow had been affected by the construction of two reservoirs in the period 1968–1985. In the period of 1986–1991, the river discharge decreased due to the combined influence of river diversions and increase of groundwater extractions for irrigation. In the fourth period of 1992–2000, the river discharge reached lowest flow values and variations corresponding to a large increase in crop area. The flow regime recovered, but not yet to natural status in the fifth period of 2001–2007. Climatic factors are not responsible for all these changes in the flow regime, but the changes are corresponding well to human activities.

Description: Integrating coarse-scale uncertain soil moisture data into a fine-scale hydrological modelling scenario Hydrology and Earth System Sciences Discussions, 8, 6031-6067, 2011 Author(s): H. Vernieuwe, B. De Baets, J. Minet, V. R. N. Pauwels, S. Lambot, M. Vanclooster, and N. E. C. Verhoest In a hydrological modelling scenario, often the modeller is confronted with external data, such as remotely-sensed soil moisture observations, that become available to update the model output. However, the scale triplet (spacing, extent and support) of these data is often inconsistent with that of the model. Furthermore, the external data can be cursed with epistemic uncertainty. Hence, a method is needed that not only integrates the external data into the model, but that also takes into account the difference in scale and the uncertainty of the observations. In this paper, a synthetic hydrological modelling scenario is set up in which a high-resolution distributed hydrological model is run over an agricultural field. At regular time steps, coarse-scale field-averaged soil moisture data, described by means of possibility distributions (epistemic uncertainty), are retrieved by synthetic aperture radar and assimilated into the model. A method is presented that allows to integrate the coarse-scale possibility distribution of soil moisture content data with the fine-scale model-based soil moisture data. To this end, a scaling relationship between field-averaged soil moisture content data and its corresponding standard deviation is employed.

Description: The impact of the high-frequency (HF, 〈90 days) variability on low-frequency (LF) interannual sea surface temperature (SST) variations associated with the El Niño-Southern Oscillation (ENSO) is investigated by conducting a series of oceanic general circulation model experiments. Two nonlinear rectification mechanisms are examined. The first is the internal oceanic nonlinear dynamics and the second is the nonlinear rectification of the LF surface wind stress by the HF wind. Numerical simulations show that the latter is dominant in modulating the LF SST variability. The HF wind increases both the amplitude and skewness of the LF wind stress anomaly. As a result, it increases both the amplitude and skewness of the SST anomaly (SSTA) in the eastern equatorial Pacific. For strong El Niño events in 1982/83 and 1997/98, such a nonlinear rectification effect may result in a SSTA increase of 1°C. A mixed-layer heat budget analysis reveals that whereas meridional and vertical advections primarily contribute to the strengthening of the warm and cold episodes, the nonlinear zonal advection is responsible for the increase of the SSTA skewness. Including the nonlinear rectification of the HF wind on both the surface wind stress and heat flux anomalies leads to a positively (negatively) skewed SSTA in the eastern (central) Pacific. Thus the combined dynamic and thermodynamic effect reshapes the ENSO zonal structure in such a way that it makes the maximum SSTA confined further to the eastern equatorial Pacific. Copyright © 2011 Royal Meteorological Society

Description: Scaling and trends of hourly precipitation extremes in two different climate zones – Hong Kong and the Netherlands Hydrology and Earth System Sciences Discussions, 8, 4701-4719, 2011 Author(s): G. Lenderink, H. Y. Mok, T. C. Lee, and G. J. van Oldenborgh Hourly precipitation extremes in very long time series from the Hong Kong Observatory and the Netherlands are investigated. Using the 2 m dew point temperature from 4 h before the rainfall event as a measure of near surface absolute humidity, hourly precipitation extremes closely follow a 14 % per degree dependency – a scaling twice as large as following from the Clausius-Clapeyron relation. However, for dew point temperatures above 23 °C no significant dependency on humidity was found. Strikingly, in spite of the large difference in climate, results are almost identical in Hong Kong and the Netherlands for the dew point temperature range where both observational sets have sufficient data. Trends in hourly precipitation extremes show substantial increases over the last century for both De Bilt (the Netherlands) and Hong Kong. For De Bilt, not only the long term trend, but also variations in hourly precipitation extremes on a inter-decadal timescale of 30 yr and longer, can be linked very well to the above scaling; there is a very close resemblance between variations in dew point temperature and precipitation intensity with an inferred dependency of hourly precipitation extremes of 10 to 14 % per degree. For Hong Kong there is no connection between variations in humidity and those in precipitation intensity in the wet season, May to September, consistent with the found zero-dependency of preciptation intensity on humidity for dew points above 23 °C. Yet, outside the wet season humidity changes do appear to explain the positive trend in hourly precipitation extremes, again following a dependency close to twice the Clausius-Clapeyron relation.

Description: In situ , satellite and model analyses data in April–July 2006 are used to investigate the links between sea surface temperature (SST) and the marine atmospheric boundary layer (MABL) in the Gulf of Guinea. The study region between 10°W and 6°E is divided into three areas with different characteristics: the North Area (4.5°N to 1°N) with the wettest atmosphere and the warmest SST, the Upwelling Area (1°N to 4°S) with the strongest SST decrease, and the South Area (4°S to 8°S) with a drier atmosphere and a more slowly decreasing SST than in the Upwelling Area. The key zone of the air-sea interactions in this region seems to be the SST front between North and Upwelling Areas. On the one hand, the study of the MABL on either side of the front shows a well-mixed layer between the surface and about 500 m high, sensitive to surface variations, which gets shallower (deeper) when the SST decreases (increases). The MABL height (about 1500 m) follows the same variations but is not exactly collocated with the SST variations. On the other hand, the observation of the MABL across the SST front shows a strengthening of southeasterlies in the South Area coinciding with a strong SST decrease in the Upwelling Area. In the latter, the wind weakens above the colder SST. Besides, in the North Area, the wind strengthens above the warmer SST. However, the wind acceleration spans from the equator to 2°N in April and as far as 4°N in June. A convergence zone is observed in the vicinity of 2°N in April, suggesting a convection activity there, favoured by the SST front. Copyright © 2011 Royal Meteorological Society

Description: This article presents and assesses an algorithm that constructs 3D distributions of cloud from passive satellite imagery and collocated 2D nadir profiles of cloud properties inferred synergistically from lidar, cloud radar and imager data. It effectively widens the active–passive retrieved cross-section (RXS) of cloud properties, thereby enabling computation of radiative fluxes and radiances that can be compared with measured values in an attempt to perform radiative closure experiments that aim to assess the RXS. For this introductory study, A-train data were used to verify the scene-construction algorithm and only 1D radiative transfer calculations were performed. The construction algorithm fills off-RXS recipient pixels by computing sums of squared differences (a cost function F ) between their spectral radiances and those of potential donor pixels/columns on the RXS. Of the RXS pixels with F lower than a certain value, the one with the smallest Euclidean distance to the recipient pixel is designated as the donor, and its retrieved cloud properties and other attributes such as 1D radiative heating rates are consigned to the recipient. It is shown that both the RXS itself and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery can be reconstructed extremely well using just visible and thermal infrared channels. Suitable donors usually lie within 10 km of the recipient. RXSs and their associated radiative heating profiles are reconstructed best for extensive planar clouds and less reliably for broken convective clouds. Domain-average 1D broadband radiative fluxes at the top of the atmosphere (TOA) for (21 km) 2 domains constructed from MODIS, CloudSat and Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) data agree well with coincidental values derived from Clouds and the Earth's Radiant Energy System (CERES) radiances: differences between modelled and measured reflected shortwave fluxes are within ±10 W m −2 for ∼35% of the several hundred domains constructed for eight orbits. Correspondingly, for outgoing longwave radiation ∼65% are within ±10 W m −2 . Copyright © 2011 Royal Meteorological Society and Crown in the right of Canada

Description: Forest cover influence on flood assessment in Italian catchments Hydrology and Earth System Sciences Discussions, 8, 4891-4926, 2011 Author(s): F. Preti, G. Forzieri, and G. B. Chirico The paper aims at evaluating to what extent forest cover could affect the flood peak frequency and magnitude in Italian catchments. The analysis is restricted to evaluating the component of the runoff coefficient which cannot be captured by the catchment lithology alone. A preliminary data mining is performed on data of 75 catchments distributed from South to Central Italy. Cluster and correlation structure analyses are conducted for distinguishing forest cover effects within sample sets of catchments characterized by hydro-morphological similarities. We propose a method for correcting the bias of the runoff coefficient estimated from the catchment lithology only, by accounting for the effect of forest cover. The bias correction becomes significant for small mountainous catchments, characterised by larger forest cover fraction and lower critical rainfall depth. Consistently with what suggested in previous studies, the bias correction decreases as the rainfall depth and return period increase.

Description: Technical Note: Development of an automated lysimeter for the calculation of peat soil actual evapotranspiration Hydrology and Earth System Sciences Discussions, 8, 5009-5033, 2011 Author(s): S. Proulx-McInnis, A. St-Hilaire, A. N. Rousseau, S. Jutras, G. Carrer, and G. Levrel A limited number of publications in the literature deal with the measurement of actual evapotranspiration (AET) from a peat soil. AET is an important parameter in the description of water pathways of an ecosystem. In peatlands, where the water table is near the surface and the vegetation is composed of nonvascular plants without stomatal resistance, the AET measurement represents a challenge. This paper discusses the development of an automated lysimeter installed between 12 and 27 July 2010, at a 11-ha bog site, Pont-Rouge (42 km west of Quebec City, Canada). This system was made of an isolated block of peat, maintained at the same water level as the surrounding water table by a system of submersible pressure transmitters and pumps. The change in water level in millimetres in the isolated block of peat was used to calculate the water lost through evapotranspiration (ET) while accounting the precipitation. The rates of AET were calculated for each day of the study period. Temperature fluctuated between 17.2 and 23.3 °C and total rainfall was 43.76 mm. AET rates from 0.6 to 6.9 mm day −1 were recorded, with a ΣAET/ΣP ratio of 1.38. The estimated potential ET (PET) resulting from Thornthwaite's semi-empirical formula suggested values between 2.8 and 3.9 mm day −1 . The average AET/PET ratio was 1.13. According to the literature, the results obtained are plausible. This system, relatively inexpensive and simple to install, may eventually be used to calculate AET on peaty soils in the years to come.

Description: In geophysical data assimilation, the control space is by definition the set of parameters which are estimated through the assimilation of observations. It has recently been proposed to design the discretizations of control space in order to assimilate observations optimally. The present paper describes the embedding of that formalism in a consistent Bayesian framework. General background errors are now accounted for. Scale-dependent errors, such as aggregation errors (that lead to representativeness errors) are consistently introduced. The optimal adaptive discretizations of control space minimize a criterion on a dictionary of grids. New criteria are proposed: degrees of freedom for the signal (DFS) built on the averaging kernel operator, and an observation-dependent criterion. These concepts and results are applied to atmospheric transport of pollutants. The algorithms are tested on the European tracer experiment (ETEX), and on a prototype of CO 2 flux inversion over Europe using a simplified CarboEurope-IP network. New types of adaptive discretization of control space are tested such as quaternary trees or factorised trees. Quaternary trees are proven to be both economical, in terms of storage and CPU time, and efficient on the test cases. This sets the path for the application of this methodology to high-dimensional and noisy geophysical systems. Part II of this article will develop asymptotic solutions for the design of control space representations that are obtained analytically and are contenders to exact numerical optimizations. Copyright © 2011 Royal Meteorological Society

Description: Technical Note: Demonstrating a 24/7 solution for monitoring water quality loads in rivers Hydrology and Earth System Sciences Discussions, 8, 5035-5050, 2011 Author(s): P. Jordan and R. Cassidy Quantifying nutrient and sediment loads in catchments is difficult owing to diffuse controls related to storm hydrology. Coarse sampling and interpolation methods are prone to very high uncertainties due to under-representation of high discharge, short duration events. Additionally, important low-flow processes such as diurnal signals linked to point source impacts are missed. Here we demonstrate a solution based on a time-integrated approach to sampling with a standard 24 bottle autosampler configured to take a sample every 7 h over a week. This is evaluated with a number of other sampling strategies using a two-year dataset of sub-hourly discharge and phosphorus concentration data. The 24/7 solution is shown to be the least uncertain in estimating load (inter-quartile range is 96 % to 110 % of actual load in year 1 and 97 % to 104 % in year 2) due to the increased frequency raising the probability of sampling storm events and point source signals. The 24/7 solution would appear to be most parsimonious in terms of technology requirements, the ability to be widely deployed and to represent important nutrient transfer processes in complex catchments.

Description: Spatial and temporal characteristics of a nocturnal low-level jet (LLJ) on the east side of the Western Ghat mountain range over India's west coast and processes leading to the formation of the jet are discussed. The boundary-layer jet has a regional scale extent, as revealed by high-resolution Advanced Research Weather Research and Forecasting (ARW) model simulations, and contributes to the formation of ‘atmospheric streams’ of water vapor over the selected land regions. Simulations indicate that the formation of LLJ is mainly attributed to the baroclinicity of the valley atmosphere due to the gently rolling terrain, which is assisted by the persistence of an unstable residual layer above the developing stable boundary layer in the valley and cooling over the slopes. Prior to the formation of LLJ, the boundary layer is dominated by deep roll circulations. The LLJ followed a gust front zone associated with a mountain wave. The low-level flow below the jet is decoupled from the upper-level flow as a result of strong vorticity below the jet and suppression of turbulence at the jet core. A conceptual model for the boundary layer interactions, dynamics of the mountain wave, LLJ, etc. are proposed for Western Ghat region. Copyright © 2011 Royal Meteorological Society

Description: This article examines the first-guess (FG) departures of microwave imager radiances assimilated in all-sky conditions (i.e. clear, cloudy and precipitating). Agreement between FG and observations is good in clear skies, with error standard deviations around 2 K, but in heavy cloud or precipitation errors increase to 20 K. The forecast model is not good at predicting cloud and precipitation with exactly the right intensity or location. This leads to apparently non-Gaussian behaviour, both heteroscedasticity, i.e. an increase in error with cloud amount, and boundedness, i.e. the size of errors is close to the geophysical range of the observations, which runs from clear to fully cloudy. However, the dependence of FG departure standard deviations on the mean cloud amount is predictable. Using this dependence to normalise the FG departures gives an error distribution that is close to Gaussian. Thus if errors are treated correctly, all-sky observations can be assimilated successfully under the assumption of Gaussianity on which assimilation systems are based. This ‘symmetric’ error model can be used to provide a robust threshold quality-control check and to determine the size of observation errors for all-sky assimilation. In practice, however, this ‘observation’ error is being used to account for the model's difficulty in forecasting cloud, which really comes from errors in the background and in the forecast model. Hence in future it will be necessary to improve the representation of background and model error. Separately, symmetric cloud amount is recommended as a predictor for bias correction schemes, avoiding the sampling problems associated with ‘asymmetric’ predictors like the FG cloud amount. Copyright © 2011 Royal Meteorological Society

Description: Bayesian inverse modelling of in situ soil water dynamics: using prior information about the soil hydraulic properties Hydrology and Earth System Sciences Discussions, 8, 2019-2063, 2011 Author(s): B. Scharnagl, J. A. Vrugt, H. Vereecken, and M. Herbst In situ observations of soil water state variables under natural boundary conditions are often used to estimate field-scale soil hydraulic properties. However, many contributions to the soil hydrological literature have demonstrated that the information content of such data is insufficient to reliably estimate all the soil hydraulic parameters. In this case study, we tested whether prior information about the soil hydraulic properties could help improve the identifiability of the van Genuchten-Mualem (VGM) parameters. Three different prior distributions with increasing complexity were formulated using the ROSETTA pedotransfer function (PTF) with input data that constitutes basic soil information and is readily available in most vadose zone studies. The inverse problem was posed in a formal Bayesian framework and solved using Markov chain Monte Carlo (MCMC) simulation with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm. Synthetic and real-world soil water content data were used to illustrate our approach. The results of this study corroborate and explicate findings previously reported in the literature. Indeed, soil water content data alone contained insufficient information to reasonably constrain all VGM parameters. The identifiability of these soil hydraulic parameters was substantially improved when an informative prior distribution was used with detailed knowledge of the correlation structure among the respective VGM parameters. A biased prior did not distort the results, which inspires confidence in the robustness and effectiveness of the presented method. The Bayesian framework presented in this study can be applied to a wide range of vadose zone studies and provides a blueprint for the use of prior information in inverse modelling of soil hydraulic properties at various spatial scales.

Description: Laboratory experiments to determine the preferred orientation of free-falling hexagonal prisms were performed at Reynolds numbers appropriate to falling ice crystals in the atmosphere. Hexagonal plates orient with their c axis vertical for aspect ratios 〈 0.9, whilst hexagonal columns fall with their c axis horizontal. A secondary alignment is also observed: regular hexagonal columns fall preferentially with two prism facets aligned vertically and not horizontally – the latter scenario was previously assumed to be responsible for the rare Parry arc. However, if the column is made scalene in its cross-section, it can orient such that a pair of prism facets is horizontal. This finding indicates that the development of scalene crystals may be key to the production of certain ice-crystal optical phenomena. Copyright © 2011 Royal Meteorological Society

Description: Classification of thermal waters based on their inorganic fingerprint and hydrogeothermal modelling Hydrology and Earth System Sciences Discussions, 8, 4559-4581, 2011 Author(s): I. Delgado-Outeiriño, P. Araujo-Nespereira, J. A. Cid-Fernández, J. C. Mejuto, E. Martínez-Carballo, and J. Simal-Gándara Hydrothermal features in Galicia have been used since ancient times for therapeutic purposes. A characterization of these thermal waters was carried out in order to understand their behaviour based on inorganic pattern and water-rock interaction mechanisms. In this way 15 thermal water samples were collected in the same hydrographical system. The results of the hydrogeochemistry analysis showed one main water family of bicarbonate type sodium waters, typical in the post-orogenic basins of Galicia. Principal component analysis (PCA) and partial lest squared (PLS) clustered the selected thermal waters in two groups, regarding to their chemical composition. This classification agreed with the results obtained by the use of geothermometers and the hydrogeochemical modelling. The first included thermal samples that could be in contact with surface waters and therefore, their residence time in the reservoir and their water-rock interaction would be less important than for the thermal waters of the second group.

Description: The response of Iberian rivers to the North Atlantic Oscillation Hydrology and Earth System Sciences Discussions, 8, 4459-4493, 2011 Author(s): J. Lorenzo-Lacruz, S. M. Vicente-Serrano, J. I. López-Moreno, J. C. González-Hidalgo, and E. Morán-Tejeda In this study we analyzed the influence of the North Atlantic Oscillation (NAO) on the streamflow in 187 sub-basins of the Iberian Peninsula. Monthly and one-month lagged correlations were conducted to assess the spatio-temporal extent of the NAO influence on Iberian river discharges. Analysis of the persistence of the winter NAO throughout the year was also undertaken, together with analysis of streamflow anomalies during positive and negative NAO phases. Moving-window correlation analyses were conducted to assess potential changes in the temporal evolution of the NAO influence on Iberian streamflows. The results show that the NAO has a large impact on surface water resources throughout the Iberian Peninsula during winter, and in the Atlantic watershed during autumn. We showed that water resources management and snowmelt are causing the persistent dependence of streamflows on the previous winter NAO. We found that strongly positive streamflow anomalies occurred during winter, especially in the Atlantic watershed, and provide evidence of non-stationarity and spatial variability in the NAO influence on Iberian water resources.

Description: Catchment classification: empirical analysis of hydrologic similarity based on catchment function in the eastern USA Hydrology and Earth System Sciences Discussions, 8, 4495-4534, 2011 Author(s): K. Sawicz, T. Wagener, M. Sivapalan, P. A. Troch, and G. Carrillo Hydrologic similarity between catchments, derived from their similarity in how they respond to precipitation input, is the basis for classification, for transferability, for generalization and also for understanding the potential impacts of environmental change. An important question in this context is, in how far can widely available hydrologic information (precipitation-temperature-streamflow) be used to create a first order grouping of hydrologically similar catchments? We utilize a heterogeneous dataset of 280 catchments located in the Eastern US to understand hydrologic similarity in a 6-dimensional signature space across a region with strong environmental gradients. Signatures are defined as hydrologic response characteristics that provide some insight into the hydrologic function of catchments. A Bayesian clustering scheme is used to separate the catchments into 9 classes, which are subsequently analyzed with respect to their hydrologic, as well as climatic and landscape attributes. Based on the empirical results we hypothesize the following: (1) Streamflow elasticity with respect to precipitation is modified by the soil characteristics of a catchment. (2) Spatial proximity is a good first indicator of hydrologic similarity because of the strong control climate exerts on catchment function, and because it varies slowly in space.

Description: Improving confidence in deep drainage estimates, for arid and semi-arid areas using multiple linear regression with percent clay content and rainfall Hydrology and Earth System Sciences Discussions, 8, 4535-4557, 2011 Author(s): D. L. Wohling, F. W. Leaney, and R. S. Crosbie Deep drainage estimates are required for effective management of water resources. However, field measurements are time consuming and costly so simple empirical relationships are often used. Relationships developed between clay content of the surface soil and deep drainage have been used extensively in Australia to provide regional estimates of drainage but these relationships have been poorly justified and did not include rainfall in the relationships. Here we present a rigorous appraisal of clay content of soils and rainfall as predictors of drainage using an extensive database of field observations from across Australia. This study found that annual average rainfall and the clay content of the top 2 m of the soil are statistically significant predictors of drainage. Relationships have been defined for annual, perennial and tree type vegetation as a line of best fit along with 95 % confidence intervals. This allows the uncertainty in these drainage estimates to be assessed for the first time.

Description: Catchment classification: hydrological analysis of catchment behavior through process-based modeling along a climate gradient Hydrology and Earth System Sciences Discussions, 8, 4583-4640, 2011 Author(s): G. Carrillo, P. A. Troch, M. Sivapalan, T. Wagener, C. Harman, and K. Sawicz Catchment classification is an efficient method to synthesize our understanding of how climate variability and catchment characteristics interact to define hydrological response. One way to accomplish catchment classification is to empirically relate climate and catchment characteristics to hydrologic behavior and to quantify the skill of predicting hydrologic response based on the combination of climate and catchment characteristics. Since there are important subsurface properties that cannot be readily measured, the skill of classification reflects (the lack of) the amount of cross-correlation between observable landscape features and unobservable subsurface features. The resulting empirical approach is also strongly controlled by the dataset used, and therefore lacks the power to generalize beyond the heterogeneity of characteristics found in the dataset. An alternative approach, that can partially alleviate the above-mentioned issue of observability, uses our current level of hydrological understanding, expressed in the form of a process-based model, to interrogate how climate and catchment characteristics interact to produce the observed hydrologic response. In this paper we present a general method of hydrologic analysis by means of a process-based model to support a bottom-up catchment classification system complementary to top-down classification methods. The model uses topographic, geomorphologic, soil and vegetation information at the catchment scale and conditions parameter values using readily available data on precipitation, temperature and streamflow. It is applicable to a wide range of catchments in different climate settings. We have developed a step-by-step procedure to analyze the observed hydrologic response and to assign parameter values related to specific components of the model. We applied this procedure to 12 catchments across a climate gradient east of the Rocky Mountains, USA. We show that the model is capable of reproducing the observed hydrologic behavior measured through hydrologic signatures chosen at different temporal scales. Next, we analyze the dominant time scales of catchment response and their dimensionless ratios with respect to climate and observable landscape features in an attempt to explain hydrologic partitioning. We find that only a limited number of model parameters can be related to observable landscape features. However, several climate-model time scales, and the associated dimensionless numbers, show scaling relationships with respect to the investigated hydrological signatures (runoff coefficient, baseflow index, and slope of the flow duration curve). Moreover, our analysis revealed systematic co-variation of climate, vegetation and soil related time scales along the climate gradient. If such co-variation can be shown to be robust across many catchments along different climate gradients, it opens perspective for model parameterization in ungauged catchments as well as prediction of hydrologic response in a rapidly changing environment.

Description: Insights from a joint analysis of Indian and Chinese monsoon rainfall data Hydrology and Earth System Sciences Discussions, 8, 3167-3187, 2011 Author(s): M. Zhou, F. Tian, U. Lall, and H. Hu Monsoon rainfall is of great importance for the agricultural production in both China and India. Understanding its rule and possibility of long term prediction is a challenge for research. This paper gives a joint analysis of Indian monsoon and Chinese monsoon, finds their teleconnection to Sea Surface Temperature anomaly (SSTa) and other climate indices individually and relationship in common. The results show that northern China garners less rainfall when whole Indian rainfall is below normal. Also, with cold SSTa over the Indonesia region, more rainfall would be distributed over India and South China.

Description: The within-day behaviour of 6 minute rainfall intensity in Australia Hydrology and Earth System Sciences Discussions, 8, 3189-3231, 2011 Author(s): A. W. Western, B. Anderson, L. Siriwardena, F. H. S. Chiew, A. Seed, and G. Blöschl The statistical behaviour and distribution of high-resolution (6 min) rainfall intensity within the wet part of rainy days (total rainfall depth 〉10 mm) is investigated for 42 stations across Australia. This paper compares nine theoretical distribution functions (TDFs) in representing these data. Two goodness-of-fit statistics are reported: the Root Mean Square Error (RMSE) between the fitted and observed within-day distribution; and the efficiency of prediction of the highest rainfall intensities (average intensity of the 5 highest intensity intervals). The three-parameter Generalised Pareto distribution was clearly the best performer. Good results were also obtained from Exponential, Gamma, and two-parameter Generalized Pareto distributions, each of which are two parameter functions, which may be advantageous when predicting parameter values. Results of different fitting methods are compared for different estimation techniques. The behaviour of the statistical properties of the within-day intensity distributions was also investigated and trends with latitude, Köppen climate zone (strongly related to latitude) and daily rainfall amount were identified. The latitudinal trends are likely related to a changing mix of rainfall generation mechanisms across the Australian continent.

Description: Gravity-wave-induced temperature fluctuations predicted by a mountain-wave parametrization scheme are compared with observed stratospheric temperature fluctuations from the High-Resolution Dynamics Limb Sounder (HIRDLS) instrument. The focus of the study is on the southern Andes region during the month of August 2006. The comparison reveals that while the mean amplitude of temperature fluctuations predicted by the parametrization is broadly consistent with those observed by HIRDLS, there are significant differences for individual cases. Ray-tracing calculations performed for these cases suggest that these differences are likely to be associated with horizontal propagation of mountain waves, which is not represented by current column-based parametrization schemes. The results presented in this note suggest that this is a deficiency of mountain-wave parametrization schemes that may affect models at resolutions typical of both climate and numerical weather prediction modelling. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Description: Water table fluctuation and its effects on vegetation in a semiarid environment Hydrology and Earth System Sciences Discussions, 8, 3271-3304, 2011 Author(s): L. Duan, T. Liu, X. Wang, Y. Luo, W. Wang, and X. Liu A good understanding of water table fluctuation effects on vegetation is crucial for sustaining fragile hydrology and ecology of semiarid areas such as the Horqin Sandy Land (HSL) in northern China, but such understanding is not well documented in literature. The objectives of this study were to examine spatio-temporal variations of water table and their effects on vegetation in a semiarid environment. A 9.71 km 2 area within the HSL was chosen and well-instrumented to continuously measure hydrometeorologic parameters (e.g., water table). The area comprises of meadow lands and sandy dunes as well as transitional zones in between. In addition to those measured data, this study also used Landsat TM and MODIS imageries and meteorological data at a station near the study area. The spatio-temporal variations were examined using visual plots and contour maps, while the effects on vegetation were determined by overlaying a water table depth map with a vegetation index map derived from the MODIS imageries. The results indicated that water table was mainly dependent on local topography, localized geological settings, and human activities (e.g., reclamation). At annual and monthly scales, water table was mainly a function of precipitation and potential evapotranspiration. A region within the study area where depth to water table was smaller tended to have better (i.e., more dense and productive) vegetation cover. Further, the results revealed that water table fluctuation was more sensitive for vegetations in the meadow lands than in the transitional zones, but it was least sensitive for vegetations in the sandy dunes.

Description: Comparison of catchment grouping methods for flow duration curve estimation at ungauged sites in France Hydrology and Earth System Sciences Discussions, 8, 3233-3269, 2011 Author(s): E. Sauquet and C. Catalogne The study aims at estimating flow duration curves (FDC) at ungauged sites in France and quantifying the associated uncertainties using a large dataset of 1080 FDCs. The interpolation procedure focuses here on 15 percentiles standardised by the mean annual flow, which is supposed to be known at each site. In particular, this paper discusses the relevance of different catchments grouping procedures on percentiles estimation by regional regression models. First, five parsimonious FDC parametric models were tested to approximate FDCs at gauged sites. The results show that the model based on Empirical Orthogonal Functions (EOF) expansion outperforms the other ones. In this model each FDC is interpreted as a linear combination of regional amplitude functions with weights – the parameters of the model – varying in space. Here, only one amplitude function was found sufficient to fit well most of the observed curves. Thus the considered model requires only two parameters to be estimated at ungauged locations. Second, homogeneous regions were derived according to hydrological response on one hand, and geological, climatic and topographic characteristics on the other hand. Hydrological similarity was assessed through two simple indicators: the concavity index ( IC ) that represents the shape of the standardized FDC and the seasonality ratio ( SR ) which is the ratio of summer and winter median flows. These variables were used as homogeneity criteria in three different methods for grouping catchments: (i) according to their membership in one of an a priori French classification into Hydro-Eco-Regions (HERs), (ii) by applying a regression tree clustering and (iii) by using hydrological neighbourhood obtained by canonical correlation analysis. Finally, regression models between physiographic and/or climatic variables and the two parameters of the EOF model were derived considering all the data and thereafter for each group obtained through the tested grouping techniques. Results on percentiles estimation in cross validation show a significant benefit to form homogeneous regions before developing regressions, particularly when grouping methods use hydrogeological information.

Description: Influences on flood frequency distributions in Irish river catchments Hydrology and Earth System Sciences Discussions, 8, 3305-3351, 2011 Author(s): S. Ahilan, J. J. O'Sullivan, and M. Bruen This study explores influences which result in shifts of flood frequency distributions in Irish rivers. Generalised Extreme Value (GEV) type I distributions are recommended in Ireland for estimating flood quantiles. This paper presents the findings of an investigation that identified the GEV statistical distributions that best fit the annual maximum (AM) data series extracted from 172 gauging stations of 126 rivers in Ireland. Of these 126 rivers, 25 have multiple gauging stations. Analysis of this data was undertaken to explore hydraulic and hydro-geological factors that influence flood frequency distributions and whether shifts in distributions occur in the down-river direction. The methodology involved determining the shape parameter of GEV distributions that were fitted to AM data at each site and to statistically test this shape parameter to determine whether a type I, type II or type III distribution was valid. The classification of these distributions was further supported by moment and L -moment diagrams and probability plots. Results indicated that of the 143 stations with flow records exceeding 25 yr, data for 92 was best represented by GEV type I distributions and that for another 12 and 39 stations followed type II and type III distributions respectively. The spatial, hydraulic and hydro-geological influences on flood frequency distributions were assessed by incorporating results on an Arc-GIS platform with individual layers showing karst features, flood attenuation polygons and lakes. This data reveals that type I distributions are spatially well represented throughout the country. The majority of type III distributions appear in four distinct clusters in well defined geographical areas where attenuation influences from floodplains and lakes appear to be influential. The majority of type II distributions appear to be in a single cluster in a region in the west of the country that is characterised by a karst landscape. The presence of karst in river catchments would be expected to provide additional subsurface storage and in this regard, type III distributions might be expected. The prevalence of type II distributions in this area reflects the finite nature of this storage and the effects, in extreme conditions, when the karst is saturated and further storage is no longer available. Results therefore indicate that in some instances assuming type I distributions is incorrect and may result in erroneous estimates of flood quantiles in these regions. Where actual data follows a type II distribution, flood quantiles may be underestimated and for type III distributions, overestimates may be expected.

Description: Quantifying flow and remediation zone uncertainties for partially opened wells in heterogeneous aquifers Hydrology and Earth System Sciences Discussions, 8, 3133-3166, 2011 Author(s): C.-F. Ni, C.-P. Lin, S.-G. Li, and J.-S. Chen This study presents a numerical first-order spectral model to quantify flow and remediation zone uncertainties for partially opened wells in heterogeneous aquifers. Taking advantages of spectral theories in solving unmodeled small-scale variability in hydraulic conductivity ( K ), the presented nonstationary spectral method (NSM) can efficiently estimate flow uncertainties, including hydraulic heads and Darcy velocities in r- and z profile in a cylindrical coordinate system. The velocity uncertainties associated with the particle backward tracking algorithm are then used to estimate stochastic remediation zones for scenarios with partially opened well screens. In this study the flow and remediation zone uncertainties obtained by NSM were first compared with those obtained by Monte Carlo simulations (MCS). A layered aquifer with different geometric mean of K and screen locations was then illustrated with the developed NSM. To compare NSM flow and remediation zone uncertainties with those of MCS, three different small-scale K variances and correlation lengths were considered for illustration purpose. The MCS remediation zones for different degrees of heterogeneity were presented with the uncertainty clouds obtained by 200 equally likely MCS realizations. Results of simulations reveal that the first-order NSM solutions agree well with those of MCS for partially opened wells. The flow uncertainties obtained by using NSM and MCS show identically for aquifers with small ln K variances and correlation lengths. Based on the test examples, the remediation zone uncertainties are not sensitive to the changes of small-scale ln K correlation lengths. However, the increases of remediation zone uncertainties are significant with the increases of small-scale ln K variances. The largest displacement uncertainties may have several meters of differences when the ln K variances increase from 0.1 to 1.0. Such results are also valid for the estimations of remediation zones in layered aquifers.

Description: Evaluating uncertainty estimates in hydrologic models: borrowing measures from the forecast verification community Hydrology and Earth System Sciences Discussions, 8, 3085-3131, 2011 Author(s): K. J. Franz and T. S. Hogue The hydrologic community is generally moving towards the use of probabilistic estimates of streamflow, primarily through the implementation of Ensemble Streamflow Prediction (ESP) systems, ensemble data assimilation methods, or multi-modeling platforms. However, evaluation of probabilistic outputs has not necessarily kept pace with ensemble generation. Much of the modeling community is still performing model evaluation using standard deterministic measures, such as error, correlation, or bias, typically applied to the ensemble mean or median. Probabilistic forecast verification methods have been well developed, particularly in the atmospheric sciences yet, few have been adopted for evaluating uncertainty estimates in hydrologic model simulations. In the current paper, we overview existing probabilistic forecast verification methods and apply the methods to evaluate and compare model ensembles produced from different parameter uncertainty estimation methods. The Generalized Uncertainty Likelihood Estimator (GLUE), a modified version of GLUE, and the Shuffle Complex Evolution Metropolis (SCEM) are used to generate model ensembles for the National Weather Service SACramento Soil Moisture Accounting (SAC-SMA) model for 12 forecast basins located in the Southeastern United States. We evaluate the model ensembles using relevant metrics in the following categories: distribution, correlation, accuracy, conditional statistics, and categorical statistics. We show that the probabilistic metrics are easily adapted to model simulation ensembles and provide a robust analysis of parameter uncertainty, one that is commensurate with the dimension of the ensembles themselves. Application of these methods requires no information in addition to what is already available as part of traditional model validation methodology and considers the entire ensemble or uncertainty range in the approach.

Description: A structure generator for modelling the initial sediment distribution of an artificial hydrologic catchment Hydrology and Earth System Sciences Discussions, 8, 4641-4699, 2011 Author(s): T. Maurer, A. Schneider, and H. H. Gerke Artificially-created hydrological catchments are characterized by sediment structures from technological construction processes that can potentially be important for modelling of flow and transport and for understanding initial soil and ecosystem development. The subsurface spatial structures of such catchments have not yet been sufficiently explored and described. Our objective was to develop a structure generator programme for modelling the 3-D spatial sediment distribution patterns depending on the technical earth-moving and deposition processes. For the development, the artificially-constructed hydrological catchment "Chicken Creek" located in Lower Lusatia, Germany, served as an example. The structure generator describes 3-D technological sediment distributions at two scales: (i) for a 2-D-vertical cross-section, texture and bulk density distributions are generated within individual spoil cones that result from mass dumping, particle segregation, and compaction and (ii) for the whole catchment area, the spoil cones are horizontally arranged along trajectories of mass dumping controlled by the belt stacker-machine relative to the catchment's clay layer topography. The generated 3-D texture and bulk density distributions are interpolated and visualized as a gridded 3-D-volume body using 3-D computer-aided design software. The generated subsurface sediment distribution for the Chicken Creek catchment was found to correspond to observed patterns although still without any calibration. Spatial aggregation and interpolation in the gridded volume body modified the generated distributions towards more uniform (unimodal) distributions and lower values of the standard deviations. After incorporating variations and pedotransfer approaches, generated sediment distributions can be used for deriving realizations of the 3-D hydraulic catchment structure.

Description: Skill assessment of a global hydrological model in reproducing flow extremes Hydrology and Earth System Sciences Discussions, 8, 3469-3505, 2011 Author(s): N. Candogan Yossef, L. P. H. van Beek, J. C. J. Kwadijk, and M. F. P. Bierkens As an initial step in assessing the prospect of using macro-scale hydrological models (MHMs) for hydrological forecasting, this study investigates the skill of the MHM PCR-GLOBWB in reproducing past discharge extremes on a global scale. Global terrestrial hydrology from 1958 until 2001 is simulated by forcing PCR-GLOBWB with daily meteorological data obtained by downscaling the CRU dataset to daily fields using the ERA-40 reanalysis. Simulated discharge values are compared with observed monthly streamflow records for a selection of 20 large river basins that represent all continents and a wide range of climatic zones. We assess model skill in three ways. First, the general performance of the model in reproducing hydrographs is evaluated. Second, model skill in reproducing significantly higher and lower flows than the monthly normals is assessed in terms of skill scores used for forecasts of categorical events. Third, model skill in reproducing flood and drought events is assessed by constructing binary contingency tables for floods and droughts for each basin. The results show that the model has skill in all three types of hindcasting. After bias correction the model skill in simulating hydrographs is improved considerably. For most basins it is much higher than that of the climatology. The skill in hindcasting monthly anomalies is high compared to that of an imaginary unskilled system. The model also performs better than an unskilled system in hindcasting floods and droughts, with a markedly higher skill in floods. We conclude that the prospect for using PCR-GLOBWB for monthly and seasonal hydrological forecasting is positive. Our results which we argue are representative for other similar MHMs, show that MHMs have sufficient skill for use in forecasting flow extremes.

Description: A new version of the Edwards–Slingo (ES) radiation scheme is developed using the correlated- k distribution (CKD) method. The work is conducted based on the line-by-line radiative transfer scheme GENLN2. A dataset of new ozone absorption cross-section in the ultraviolet spectrum and new oxygen collision-induced continuum data have been implemented in both the GENLN2 and new ES schemes. In order to improve the efficiency of the ES scheme, a new technique is proposed in this work to optimize the k distribution and a new method is used to deal with the gaseous overlapping absorption in a spectral band. The accuracy of the scheme is improved by replacing the scaling function used in the ES scheme with a pre-determined look-up table for consideration of the pressure and temperature dependency of the absorption oefficients. The number of spectral bands and number of absorbing species are both increased in the new scheme for better resolution of the spectral variation of absorbing species, aerosols and clouds. However, this does not increase the computational cost. Instead, it is reduced substantially compared with the previous version of the code. The treatment of transmission in the short-wave spectrum is improved by implementing the absorbing species of CH 4 , N 2 O and O 2 collision-induced continuum absorption which are not included in the ES scheme. New O 3 absorption cross-section data in the ultraviolet spectrum measured by the European Space Agency are used to generate the CKD spectral data in the short-wave spectral bands. These data have a temperature dependency and better spectral resolution. The irradiance and heating rate determined by the new scheme are tested against the same variables determined by GENLN2. It has been shown that the new scheme produces results more accurate than the ES scheme. Copyright © 2011 Royal Meteorological Society

Description: Ice cloud representation in general circulation models remains a challenging task, due to the lack of accurate observations and the complexity of microphysical processes. In this article, we evaluate the ice water content (IWC) and ice cloud fraction statistical distributions from the numerical weather prediction models of the European Centre for Medium-Range Weather Forecasts (ECMWF) and the UK Met Office, exploiting the synergy between the CloudSat radar and CALIPSO lidar. Using the last three weeks of July 2006, we analyse the global ice cloud occurrence as a function of temperature and latitude and show that the models capture the main geographical and temperature-dependent distributions, but overestimate the ice cloud occurrence in the Tropics in the temperature range from −60 °C to −20 °C and in the Antarctic for temperatures higher than −20 °C, but underestimate ice cloud occurrence at very low temperatures. A global statistical comparison of the occurrence of grid-box mean IWC at different temperatures shows that both the mean and range of IWC increases with increasing temperature. Globally, the models capture most of the IWC variability in the temperature range between −60 °C and −5 °C, and also reproduce the observed latitudinal dependencies in the IWC distribution due to different meteorological regimes. Two versions of the ECMWF model are assessed. The recent operational version with a diagnostic representation of precipitating snow and mixed-phase ice cloud fails to represent the IWC distribution in the −20 °C to 0 °C range, but a new version with prognostic variables for liquid water, ice and snow is much closer to the observed distribution. The comparison of models and observations provides a much-needed analysis of the vertical distribution of IWC across the globe, highlighting the ability of the models to reproduce much of the observed variability as well as the deficiencies where further improvements are required. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Description: Mass-flux convection schemes continue to play key roles in large-scale atmospheric modelling. Currently, however, the specification of detrainment seems a potential weakness. Using cloud-resolving model results from the European Cloud Systems (EUROCS) humidity case, we consider how detrainment adapts to its environment and formulate a minimal-complexity description of partial detrainment in convective cloud fields. This can be viewed as an approximation to the behaviour of a multi-plume scheme. The algorithm is straightforwardly implemented within the Gregory–Rowntree convection scheme as an extension of its original partial-detrainment scheme. Numerical weather prediction tests in the Met Office Unified Model show significant benefits in full-model performance. © 2011 Crown Copyright, the Met Office. Published by JohnWiley & Sons Ltd.

Description: On the colour and spin of epistemic error (and what we might do about it) Hydrology and Earth System Sciences Discussions, 8, 5355-5386, 2011 Author(s): K. Beven, P. J. Smith, and A. Wood Disinformation as a result of epistemic error is an issue in hydrological modelling. In particular the way in which the colour in model residuals resulting from epistemic errors should be expected to be non-stationary means that it is difficult to justify the spin that the structure of residuals can be properly represented by statistical likelihood functions. To do so would be to greatly overestimate the information content in a set of calibration data and increase the possibility of both Type I and Type II errors. Some principles of trying to identify periods of disinformative data prior to evaluation of a model structure of interest, are discussed. An example demonstrates the effect on the estimated parameter values of a hydrological model.

Description: Effects of freezing on soil temperature, frost propagation and moisture redistribution in peat: laboratory investigations Hydrology and Earth System Sciences Discussions, 8, 5387-5426, 2011 Author(s): R. M. Nagare, R. A. Schincariol, W. L. Quinton, and M. Hayashi The effects of freezing on soil temperature and water movement were monitored in four peat Mesocosms subjected to bidirectional freezing. Temperature gradients were applied by bringing the Mesocosm tops in contact with sub-zero air temperature while maintaining a continuously frozen layer at the bottom (proxy permafrost). Soil water movement towards the freezing front (from warmer to colder regions) was inferred from soil freezing curves and from the total water content of frozen core samples collected at the end of freezing cycle. This study illustrates how differences in initial water content influence the hydrologic functions of active layer in permafrost terrains covered with thick peat during soil freezing. A substantial amount of water, enough to raise the upper surface of frozen saturated soil within 15 cm of the soil surface at the end of freezing period, appeared to have moved upwards during freezing. Effects of temperature on soil matric potential, at least in the initial freezing period, appear to drive such movement as seen from analysis of soil freezing curves. Vapour movement from warmer to colder regions also appears to contribute in moisture movement. Frost propagation is controlled by latent heat for a long time during freezing. A simple conceptual model describing freezing of an organic active layer initially resembling a variable moisture landscape is proposed based upon the results of this study. The results of this study will help in understanding, and ultimately forecasting, the hydrologic response of wetland-dominated terrain underlain by discontinuous permafrost.

Description: We investigate the issue of variational and sequential data assimilation with nonlinear and non-smooth observation operators using a two-dimensional limited-area shallow-water equation model and its adjoint. The performance of the four-dimensional variational approach (4D-Var: two dimensions plus time) compared with that of the maximum-likelihood ensemble filter (MLEF), a hybrid ensemble/variational method, is tested in the presence of non-smooth observation operators. Following the work of Lewis & Overton and Karmitsa, we investigate minimization of the data-assimilation cost functional using the limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) quasi-Newton algorithm originally intended for smooth optimization and the limited-memory bundle method (LMBM) algorithm specifically designed to address large-scale non-smooth minimization problems. Numerical results obtained for the MLEF method show that the LMBM algorithm yields results superior to the L-BFGS method. Results for 4D-Var suggest that L-BFGS performs well when the non-smoothness is not extreme, but fails for non-smooth functions with large Lipschitz constants. The LMBM method is found to be a suitable choice for large-scale non-smooth optimization, although additional work is needed to improve its numerical stability. Finally, the results and methodologies of 4D-Var and MLEF are compared and contrasted. Copyright © 2011 Royal Meteorological Society

Description: The simultaneous presence of convective, symmetric and inertial instability in a prefrontal region of strong vertical wind shear was simulated with the non-hydrostatic model MOLOCH. Model diagnostics of absolute vorticity, pseudo-angular momentum, saturated equivalent potential temperature and vorticity reveal a sequence of events that includes ‘Δ M -adjustment’, followed by slanted ascent in symmetrically unstable regions, becoming saturated in the later stages. An idealized experiment without orography was performed to isolate the presence and role of instabilities characterizing the development. The diagnosed circulation is reminiscent of a wavenumber-two normal mode of dry symmetric instability, while moist symmetric instability is confined to a very limited region, despite the appearance of wider areas of negative moist potential vorticity. The evaluation of moist thermodynamic quantities which give proper account of condensate loading is suggested as a possible resolution of this apparent inconsistency. Copyright © 2011 Royal Meteorological Society

Description: Localization is an essential element of ensemble-based Kalman filters in large-scale systems. Two localization methods are commonly used: covariance localization and domain localization. The former applies a localizing weight to the forecast covariance matrix, while the latter splits the assimilation into local regions in which independent assimilation updates are performed. The domain localization is usually combined with observation localization, which is a weighting of the observation-error covariance matrix, resulting in a similar localization effect to that of covariance localized filters. It is shown that the use of the same localization function in covariance localization and observation localization results in distinct effective localization length-scales in the Kalman gain. In order to improve the performance of observation localization, a regulated localization scheme is introduced. Twin experiments with the Lorenz-96 model demonstrate that the regulated localization can lead to a significant reduction of estimation errors as well as increased stability of the assimilation process. Copyright © 2011 Royal Meteorological Society

Description: SCS-CN parameter determination using rainfall-runoff data in heterogeneous watersheds. The two-CN system approach Hydrology and Earth System Sciences Discussions, 8, 8963-9004, 2011 Author(s): K. X. Soulis and J. D. Valiantzas The Soil Conservation Service Curve Number (SCS-CN) approach is widely used as a simple method for predicting direct runoff volume for a given rainfall event. The CN values can be estimated by being selected from tables. However, it is more accurate to estimate the CN value from measured rainfall-runoff data (assumed available) in a watershed. Previous researchers indicated that the CN values calculated from measured rainfall-runoff data vary systematically with the rainfall depth. They suggested the determination of a single asymptotic CN value observed for very high rainfall depths to characterize the watersheds' runoff response. In this paper, the novel hypothesis that the observed correlation between the calculated CN value and the rainfall depth in a watershed reflects the effect of the inevitable presence of soil-cover complex spatial variability along watersheds is being tested. Based on this hypothesis, the simplified concept of a two-CN heterogeneous system is introduced to model the observed CN-rainfall variation by reducing the CN spatial variability into two classes. The behavior of the CN-rainfall function produced by the proposed two-CN system concept is approached theoretically, it is analyzed systematically, and it is found to be similar to the variation observed in natural watersheds. Synthetic data tests, natural watersheds examples, and detailed study of two natural experimental watersheds with known spatial heterogeneity characteristics were used to evaluate the method. The results indicate that the determination of CN values from rainfall runoff data using the proposed two-CN system approach provides reasonable accuracy and it over performs the previous original method based on the determination of a single asymptotic CN value. Although the suggested method increases the number of unknown parameters to three (instead of one), a clear physical reasoning for them is presented.

Description: In this work the dynamic behaviour of the wind in the nocturnal boundary layer is studied, with a particular focus on systematic behaviour of the near-surface wind. Recently, an extension of the well-known Blackadar model for frictionless inertial oscillations above the nocturnal boundary layer was proposed by Van de Wiel et al. , which accounts for frictional effects within the nocturnal boundary layer. It appears that the nocturnal wind velocity profile tends to perform an inertial oscillation around an equilibrium wind profile, rather than around the geostrophic wind vector (as in the Blackadar model). In the present study we propose the concept of ‘composite hodographs’ to evaluate the ideas and assumptions of the aforementioned analytical model. Composite hodographs are constructed based on a large observational dataset from the Cabauw observatory. For comparison and deeper analysis, this method is also applied to single-column model simulations that represent the same dataset. From this, it is shown that winds in the middle and upper part of the nocturnal boundary layer closely follow the dynamics predicted by the model by Van de Wiel et al. In contrast, the near-surface wind shows more complex behaviour that can be described by two different stages: (1) a decelerating phase where the wind decreases rapidly in magnitude due to enlarged stress divergence in the transition period near sunset (an aspect not included in the analytical model), and (2) a regular type of inertial oscillation, but with relatively small amplitude as compared to the oscillations in the middle and upper parts of the nocturnal boundary layer. Copyright © 2011 Royal Meteorological Society

Description: The limits of predictability of the meridional overturning circulation (MOC) and upper-ocean temperatures due to errors in ocean initial conditions and model parametrizations are investigated in an idealized configuration of an ocean general circulation model (GCM). Singular vectors (optimal perturbations) are calculated using the GCM, its tangent linear and adjoint models to determine an upper bound on the predictability of North Atlantic climate. The maximum growth time-scales of MOC and upper-ocean temperature anomalies, excited by the singular vectors, are 18.5 and 13 years respectively and in part explained by the westward propagation of upper-ocean anomalies against the mean flow. As a result of the linear interference of non-orthogonal eigenmodes of the non-normal dynamics, the ocean dynamics are found to actively participate in the significant growth of the anomalies. An initial density perturbation of merely 0.02 kg m −3 is found to lead to a 1.7 Sv MOC anomaly after 18.5 years. In addition, Northern Hemisphere upper-ocean temperature perturbations can be amplified by a factor of 2 after 13 years. The growth of upper-ocean temperature and MOC anomalies is slower and weaker when excited by the upper-ocean singular vectors than when the deep ocean is perturbed. This leads to the conclusion that predictability experiments perturbing only the atmospheric initial state may overestimate the predictability time. Interestingly, optimal MOC and upper-ocean temperature excitations are only weakly correlated, thus limiting the utility of SST observations to infer MOC variability. The excitation of anomalies in this model might have a crucial impact on the variability and predictability of Atlantic climate. The limit of predictability of the MOC is found to be different from that of the upper-ocean heat content, emphasizing that errors in ocean initial conditions will affect various measures differently and such uncertainties should be carefully considered in decadal prediction experiments. Copyright © 2011 Royal Meteorological Society

Description: A channel transmission losses model for different dryland rivers Hydrology and Earth System Sciences Discussions, 8, 8903-8962, 2011 Author(s): A. C. Costa, A. Bronstert, and J. C. de Araújo Channel transmission losses in drylands take place normally in extensive alluvial channels or streambeds underlain by fractured rocks. They can play an important role in flood prediction, groundwater recharge, freshwater supply and channel-associated ecosystems. We aim to develop a semi-distributed channel transmission losses model, a coupling of formulations which are more suitable for data-scarce dryland environments, applicable for both hydraulically disconnected losing streams and hydraulically connected losing(/gaining) streams. Hence, this approach should be able to cover a large variation in climate and hydro-geologic controls, which are typically found in dryland regions of the world. Traditionally, channel transmission losses models have been developed for site specific conditions. Our model was firstly evaluated for a losing/gaining, hydraulically connected 30 km reach of the Jaguaribe River, Ceará, Brazil, which controls a catchment area of 20 000 km 2 . Secondly, we applied it to a small losing, hydraulically disconnected 1.5 km channel reach in the Walnut Gulch Experimental Watershed (WGEW), Arizona, USA. The model based on the perceptual hydrological models of the reaches was able to predict reliably the stream flow for the both case studies. For the larger river reach, the evaluation of the hypotheses on the dominant hydrological processes was fundamental for reducing structural model uncertainties and improving the stream flow prediction, showing that both lateral stream-aquifer water fluxes and groundwater flow in the underlying alluvium parallel to the river course are necessary to predict stream flow and channel transmission losses, the former process being more relevant than the latter. The sensitivity analysis showed that even if the parameters can "potentially" produce large flow exchanges between model units in the saturated part of the modelling, large flow exchanges do not happen because they are restricted by the actual hydraulic gradient between the model units. Moreover, the saturated-part-based parameters (active in the larger river) produced much smaller variation in the sensitivity coefficient than those (active in the smaller river) which drive the unsaturated part of the channel transmission losses model.

Description: Watershed discretization based on multiple factors and its application in the Chinese Loess Plateau Hydrology and Earth System Sciences Discussions, 8, 9063-9087, 2011 Author(s): Y. Xu, B. Fu, C. He, and G. Gao The spatial discretization of watersheds is an indispensable procedure for representing landscape variations in eco-hydrological research, representing the contrast between reality and data-supported models. When discretizing a watershed, it is important to construct a scheme of a moderate number of discretized factors while adequately considering the actual eco-hydrological processes, especially in regions with unique eco-hydrological features and intense human activities. Because of their special lithological and pedologic characteristics and widespread man-made vegetation, discretization of watersheds in the Loess Plateau in Northern China is a challenge. In order to simulate the rainfall-runoff process, a watershed in the Loess Plateau, referred as Ansai, was spatially discretized into new units called land type units. These land type units were delineated under a scheme of factors including land use, vegetation condition, soil type and slope. Instead of using units delineated by overlaying land use and soil maps, the land type units were used in the Soil and Water Assessment Tool (SWAT). Curve numbers were assigned and adjusted to simulate runoff, using the US Natural Resources Conservation Service (NRCS) curve number method. The results of the runoff simulation better matched actual observations. Compared to the results that used the original units, the coefficient of determination ( R 2 ) and the Nash–Sutcliffe coefficient ( E NS ) for monthly flow simulation increased from 0.710–0.721 and 0.581–0.656 to 0.726–0.731 and 0.692–0.703, respectively. This method of delineating into land type units is an easy operation and suitable approach for eco-hydrological studies in the Chinese Loess Plateau and other similar regions. It can be further applied in soil erosion simulation and the eco-hydrological assessment of re-vegetation.

Description: Preferred jet stream positions and their link to regional circulation patterns over the winter North Atlantic/European sector are investigated to corroborate findings of multimodal behaviour of the jet positions and to analyse patterns of preferred paths and transition probabilities between jet regimes using ERA-40 data. Besides the multivariate Gaussian mixture model, hierarchical clustering and data image techniques are used for this purpose. The different approaches all yield circulation patterns that correspond to the preferred jet regimes, namely the southern, central and the northern positions associated respectively with the Greenland anticyclone or blocking, and two opposite phases of an East Atlantic-like flow pattern. Growth and decay patterns as well as preferred paths of the system trajectory are studied using the mixture model within the delay space. The analysis shows that the most preferred paths are associated with central to north and north to south jet stream transitions with a typical time-scale of about 5 days, and with life cycles of 1–2 weeks. The transition paths are found to be consistent with transition probabilities. The analysis also shows that wave breaking seems to be the dominant mechanism behind Greenland blocking. Copyright © 2011 Royal Meteorological Society

Description: The transformed Eulerian mean residual circulation is calculated from ERA-Interim for 1989–2009. Known as the Brewer–Dobson circulation, this measures the tropical upwelling of mass from troposphere to stratosphere, the mean meridional mass transport in the stratosphere and the downwelling of mass in the Extratropics. Major features of the Brewer–Dobson circulation, including the seasonal migration of the tropical upwelling toward the summer pole, are well represented. In the tropical lower stratosphere vertical velocities are less noisy than in other reanalyses, though significant tidal variations demonstrate the need for 6-hourly data. Throughout the year tropical lower stratospheric ascent rates are a minimum at the Equator and strongest in the Northern Hemisphere. In each hemisphere the maximum tropical ascent occurs during summer, whereas the strongest circulation and maximum in extratropical descent occur in the winter hemisphere. At 70 hPa the annual mean upwelling mass flux is 5.9 × 10 9 kg s −1 , with the zonal drag from resolved waves and parametrized orographic gravity wave drag (OGWD) providing 70% and 4% of the driving, respectively. Hence it is concluded that the OGWD probably underestimates the momentum deposited above 70 hPa in addition to there being an absence of drag from non-orographic gravity waves. A statistically significant trend of −5% per decade in the upwelling mass flux is considered unreliable because it is inconsistent with the negative temperature trend, assuming a mainly adiabatic temperature response at this level (70 hPa) to the changes in upwelling. Copyright © 2011 British Crown copyright, the Met Office. Published by John Wiley & Sons Ltd.

Description: Improving evapotranspiration in land surface models by using biophysical parameters derived from MSG/SEVIRI satellite Hydrology and Earth System Sciences Discussions, 8, 9113-9171, 2011 Author(s): N. Ghilain, A. Arboleda, G. Sepulcre-Cantò, O. Batelaan, J. Ardö, and F. Gellens-Meulenberghs Vegetation parameters derived from the geostationary satellite MSG/SEVIRI have been distributed at a daily frequency since 2007 over Europe, Africa and part of South America, through the LSA-SAF facility. We propose here a method to handle two new remote sensing products from LSA-SAF, leaf area index and Fractional Vegetation Cover, noted LAI and FVC respectively, for land surface models at MSG/SEVIRI scale. The developed method relies on an ordinary least-square technique and a land cover map to estimate LAI for each model plant functional types of the model spatial unit. The method is conceived to be applicable for near-real time applications at continental scale. Compared to monthly vegetation parameters from a vegetation database commonly used in numerical weather predictions (ECOCLIMAP-I), the new remote sensing products allows a better monitoring of the spatial and temporal variability of the vegetation, including inter-annual signals, and a decreased uncertainty on LAI to be input into land surface models. We assess the impact of using LSA-SAF vegetation parameters compared to ECOCLIMAP-I in the land surface model H-TESSEL at MSG/SEVIRI scale. Comparison with in-situ observations in Europe and Africa shows that the results on evapotranspiration are mostly improved, and especially in semi-arid climates. At last, the use of LSA-SAF and ECOCLIMAP-I is compared with simulations over a North-South Transect in Western Africa using LSA-SAF radiation forcing derived from remote sensing, and differences are highlighted.

Description: Single column model simulations using the UK Met Office Unified Model, as used in the Australian Community Climate Earth System Simulator, are presented for the Tropical Warm Pool–International Cloud Experiment (TWP–ICE) field study. Two formulations for the representation of clouds are compared with the extensive observations taken during the campaign, giving insight into the ability of the model to simulate tropical cloud systems. During the active monsoon phase the modelled cloud cover has a stronger dependence on relative humidity than the observations. Observed ice cloud properties during the suppressed monsoon period show that the ice water content is significantly underestimated in the simulations. The profiles of modelled ice fall speeds are faster than those observed in the levels above 12 km, implying that the observations have smaller sized particles in larger concentrations than the models. Both simulations show similar errors in the diurnal cycle of relative humidity during the active monsoon phase, suggesting that the error is less sensitive to the choice of cloud scheme and rather is driven by the convection scheme. However, during the times of suppressed convection the relative humidity error is different between the simulations, with congestus convection drying the environment too much, particularly in the prognostic cloud-scheme simulation. This result shows that the choice of cloud scheme and the way that the cloud and convection schemes interact plays a role in the temperature and moisture errors during the suppressed monsoon phase, which will impact the three-dimensional model simulations of tropical variability. Copyright © 2011 Royal Meteorological Society

Description: In-situ evaluation of internal drainage in layered soils (Tukulu, Sepane and Swartland) Hydrology and Earth System Sciences Discussions, 8, 9797-9841, 2011 Author(s): S. S. W. Mavimbela and L. D. van Rensburg The soil water release (SWC) and permeability properties of layered soils following deep infiltration depends on the structural and layering composition of the profiles diagnostic horizons. Three layered soils, the Tukulu, Sepane and Swartland soil forms, from the Free State province of South Africa, were selected for internal drainage evaluation. The soil water release curves as a function of suction ( h ) and unsaturated hydraulic conductivity ( K -coefficient) as a function of soil water content, SWC ( θ ), were characterised alongside the pedological properties of the profiles. The water hanging column in collaboration with the in-situ instantaneous profile method (IPM) was appropriate for this work. Independently, the saturated hydraulic conductivity ( K s ) was measured using double ring infiltrometers. The three soils had a generic orthic A horizon but differed remarkable with depth. A clay rich layer was found in the Tukulu and Sepane at depths of 600 to 850 mm and 300 to 900 mm, respectively. The Swartland was weakly developed with a saprolite rock found at depth of 400–700 mm. During the 1200 h drainage period, soil water loss amounted to 21, 20 and 51 mm from the respective Tukulu, Sepane and Swartland profiles. An abrupt drop in K s in conjunction with a steep K -coefficient gradient with depth was observed from the Tukulu and Sepane. Hydromorphic colours found on the clay-rich horizons suggested a wet soil water regime that implied restriction of internal drainage. It was therefore concluded that the clay rich horizons gave the Tukulu and Sepane soil types restricted internal drainage properties required for soil water storage under infield rainwater harvesting production technique. The coarseness of the Swartland promoted high drainage losses that proliferated a dry soil water regime.

Description: Promoting interdisciplinary education – the Vienna Doctoral Programme on Water Resource Systems Hydrology and Earth System Sciences Discussions, 8, 9843-9887, 2011 Author(s): G. Blöschl, G. Carr, C. Bucher, A. H. Farnleitner, H. Rechberger, W. Wagner, and M. Zessner The Vienna Doctoral Programme on Water Resource Systems (DK-WRS) is a programme that aims to educate students in interdisciplinary water science through cutting edge research at an international level. It is funded by the Austrian Science Fund and designed to run over a period of 12 yr during which 80 doctoral students are anticipated to graduate. This paper reports on our experiences of setting up and implementing the Programme. We identify three challenges: integrating the disciplines, maintaining depth in an interdisciplinary programme, and teaching subjects remote to each student's core expertise. To address these challenges we adopted a number of approaches. We use three levels of instruments to foster integration across the disciplines: joint groups (e.g. a joint study programme), joint science questions (e.g. developed in annual symposia), and joint study sites. To maintain depth we apply a system of quality control including regular feedback sessions, theses by journal publications and international study exchange. For simultaneously teaching students from civil and environmental engineering, biology, geology, chemistry, mathematics we use visually explicit teaching, learning by doing, extra mentoring and by cross relating associated subjects. Our initial assessment of the Programme shows some very positive outcomes. Joint science questions formed between students from various disciplines indicate integration is being achieved. The number of successful publications in top journals suggests that depth is maintained. Positive feedback from the students on the variety and clarity of the courses indicates the teaching strategy is working well. Our experiences have shown that implementing and running an interdisciplinary doctoral programme has its challenges and is demanding in terms of time and human resources but seeing interactions progress and watching people grow and develop their way of thinking in an interdisciplinary environment is a valuable reward.

Description: This work tackles the problem of the automated detection of the atmospheric boundary layer (BL) height h , from aerosol lidar/ceilometer observations. A new method, the Bayesian selective method (BSM), is presented. It implements a Bayesian statistical inference procedure which combines in a statistically optimal way different sources of information. Firstly, atmospheric stratification boundaries are located from discontinuities in the ceilometer backscattered signal. The BSM then identifies the discontinuity edge that has the highest probability to effectively mark the BL height. Information from the contemporaneous physical boundary layer model simulations and a climatological dataset of BL height evolution are combined in the assimilation framework to assist this choice. The BSM algorithm has been tested for 4 months of continuous ceilometer measurements collected during the BASE:ALFA project, and is shown to realistically diagnose the BL depth evolution in many different weather conditions. A standard one-dimensional processing of the ceilometer signal without the a priori support of the dynamical and climatological BL models often fails to correctly detect h , with the greatest inaccuracies occurring at night-time when residual layers can generate very strong signals, which are then classified by an automated application of the gradient or of the wavelet analysis as the most probable BL height. The BSM approach instead carries information on the low climatological probability to find elevated BL depths at night and penalizes the selection of these points. Moreover, this method is able to correctly convey information along the temporal dimension, thus filling data gaps using earlier and subsequent ceilometer information for the retrieval. Copyright © 2011 Royal Meteorological Society

Description: Retrospective optimal interpolation (ROI) is a method that is used to minimize cost functions with multiple minima without using adjoint models. We address two weaknesses associated with the cost-effective formulation of ROI and offer possible solutions. The first weakness of the cost-effective ROI formulation is that the error tolerance would become large in practical application due to computation costs. When the error tolerance is large, accuracy-saturated modes do not extract information from new incoming observations even though they are likely to be flawed. To address this problem, we modify the existing ‘reduced-resolution’ formulation by using eigen-decomposition of the background error covariance at each analysis step. We refer to the modified algorithm as a ‘reduced-rank’ algorithm. This modification allows us to deal with larger error variances while analysing the same number of control variables, because eigen-decomposition steeply reorders the error variance in descending order. As a result, when the reduced-rank algorithm is applied, the number of analysed control variables becomes smaller than when the reduced-resolution algorithm is used with the same error tolerance. The second weakness is the underestimation of the trailing mode of background error covariance that is projected onto the future observation space. This originates from errors in control variables extracted from analysis procedures. To prevent the occurrence of filter divergence due to this underestimation, we reduce the weighting of the observational increments in the analysis. By implicitly assuming that the rate of being projected onto the future observation space of the trailing eigenmodes is similar to that of the leading eigenmodes, we develop a method that inflates the observation error covariance and consequently improves analysis quality in the Lorenz 40-variable and the 960-variable model experiments. Copyright © 2011 Royal Meteorological Society

Description: El-Niño southern oscillation and rainfall erosivity in the headwater region of the Grande River Basin, Southeast Brazil Hydrology and Earth System Sciences Discussions, 8, 10707-10738, 2011 Author(s): C. R. Mello, L. D. Norton, N. Curi, S. N. M. Yanagi, and A. M. Silva Relationships between regional climate and oceanic and atmospheric anomalies are important tools in order to promote the development of models for predicting rainfall erosivity, especially in regions with substantial intra-annual variability in the rainfall regime. In this context, this work aimed to analyze the rainfall erosivity in headwaters of Grande River Basin, Southern Minas Gerais State, Brazil. This study considered the two most representative environments, the Mantiqueira Range (MR) and Plateau of Southern Minas Gerais (PSM). These areas are affected by the El Nino Southern Oscillation (ENSO) indicators Sea Surface Temperature (SST) for Niño 3.4 Region and Multivariate ENSO Index (MEI). Rainfall erosivity was calculated for individual rainfall events from January 2006 to December 2010. The analyses were conducted using the monthly data of ENSO indicators and the following rainfall variables: rainfall erosivity (EI 30 ), rainfall depth ( P ), erosive rainfall depth ( E ), number of rainfall events (NRE), number of erosive rainfall events (NEE), frequency of occurrence of an early rainfall pattern (EP), occurrence of late rainfall pattern (LP) and occurrence of intermediate rainfall patter (IP). Pearson's coefficient of correlation was used to evaluate the relationships between the rainfall variables and SST and MEI. The coefficients of correlation were significant for SST in the PSM sub-region. Correlations between the rainfall variables and negative oscillations of SST were also significant, especially in the MR sub-region, however, the Person's coefficients were lesser than those obtained for the SST positive oscillations. The correlations between the rainfall variables and MEI were also significant but lesser than the SST correlations. These results demonstrate that SST positive oscillations play a more important role in rainfall erosivity, meaning they were more influenced by El-Niño episodes. Also, these results have shown that the ENSO variables have potential to be useful for rainfall erosivity forecasting in this region.

Description: Extreme runoff response to short-duration convective rainfall in South-West Germany Hydrology and Earth System Sciences Discussions, 8, 10739-10780, 2011 Author(s): V. Ruiz-Villanueva, M. Borga, D. Zoccatelli, L. Marchi, E. Gaume, U. Ehret, and E. Zehe The 2 June 2008 flood-producing storm on the Starzel river basin in South-West Germany is examined as a prototype for organized convective systems that dominate the upper tail of the precipitation frequency distribution and are likely responsible for the flash flood peaks in this region. The availability of high-resolution rainfall estimates from radar observations and a rain gauge network, together with indirect peak discharge estimates from a detailed post-event survey, provides the opportunity to study the hydrometeorological and hydrological mechanisms associated with this extreme storm and the ensuing flood. Radar-derived rainfall, streamgauge data and indirect estimates of peak discharges are used along with a distributed hydrologic model to reconstruct hydrographs at multiple locations. The influence of storm structure, evolution and motion on the modeled flood hydrograph is examined by using the "spatial moments of catchment rainfall" (Zoccatelli et al., 2011). It is shown that downbasin storm motion had a noticeable impact on flood peak magnitude. Small runoff ratios (less than 20%) characterized the runoff response. The flood response can be reasonably well reproduced with the distributed hydrological model, using high resolution rainfall observations and model parameters calibrated at a river section which includes most of the area impacted by the storm.

Description: We develop a hydrostatic Hamiltonian particle-mesh (HPM) method for efficient long-term numerical integration of the atmosphere. In the HPM method, the hydrostatic approximation is interpreted as a holonomic constraint for the vertical position of particles. This can be viewed as defining a set of vertically buoyant horizontal meshes, with the altitude of each mesh point determined so as to satisfy the hydrostatic balance condition and with particles modelling horizontal advection between the moving meshes. We implement the method in a vertical-slice model and evaluate its performance for the simulation of idealized linear and nonlinear orographic flow in both dry and moist environments. The HPM method is able to capture the basic features of the gravity wave to a degree of accuracy comparable with that reported in the literature. The numerical solution in the moist experiment indicates that the influence of moisture on wave characteristics is represented reasonably well and the reduction of momentum flux is in good agreement with theoretical analysis. Copyright © 2011 Royal Meteorological Society

Description: This study investigates transient events of intense ocean evaporation with an amplitude exceeding 250 W m −2 , a duration of a few days, and a spatial extent of about 10ˆ6 km 2 over the eastern North Atlantic (referred to as ‘evaporation hot spots’) and their impact on southern Alpine heavy precipitation. First, moisture sources for a heavy precipitation event in the Piedmont in November 2002 are studied using a water-tagging simulation with a regional model. The results reveal three main moisture sources: land evapotranspiration, and evaporation from the Mediterranean and the North Atlantic, with the last source contributing the most. This was partly due to an evaporation hot spot that appeared along the western edge of a prominent upper-level trough about two days prior to the onset of heavy precipitation. In the hot spot area strong surface winds induced by the upper-level trough led to intense evaporation of water that was transported around the trough to the Piedmont region during subsequent days, where it contributed to the heavy precipitation. Secondly, analyses by the European Centre for Medium-Range Weather Forecasts (ECMWF) are used to investigate climatologically the potential relationship between eastern North Atlantic evaporation hot spots and southern Alpine precipitation. During a 10-year time period, 42 hot spots have been identified in the eastern North Atlantic. It is shown that they typically occur along the western flank of prominent upper-level troughs, and that the evaporating moisture is transported to Europe within one to four days. A climatological analysis of southern Alpine heavy precipitation events shows that they are frequently preceded by intense North Atlantic evaporation. Hence the climatological analysis further supports the conclusion from the Piedmont 2002 tagging experiment that intense evaporation over the North Atlantic and the subsequent moisture transport, both induced by the upper-level trough, are potential key factors for the development of southern Alpine heavy precipitation events. Copyright © 2011 Royal Meteorological Society

Description: The spectral predictability of the Met Office's Unified Model is examined using identical-twin experiments and the relative kinetic energy. In the troposphere, previous NWP results are recovered, namely the emergence of distinct regimes and the maximisation of the growth rate on synoptic scales at early times; in the middle atmosphere, the predictability decay is slower. This difference is attributed to the increase in the amplitude of gravity waves. The influence of small-scale motions is highlighted: improving their numerical representation by decreasing the timestep enhances predictability. Copyright © 2011 British Crown copyright, the MetOffice. Published by John Wiley & Sons Ltd.

Description: The effect of the relative orientation of the vertical wind shear to the surface wind on tropical cyclogenesis is explored in environments of radiative-convective equilibrium (RCE) through numerical simulation. This study serves as a companion paper to an earlier study on the thermodynamics of genesis in RCE. It is found, when the mean surface wind and shear are aligned, a negative surface wind anomaly arises from the superposition of the mean and vortex surface flows left of the shear vector. The resulting weak surface enthalpy fluxes and up-shear quasi-balanced subsidence leads to dry air being located cyclonically down-wind of the down-shear convective anomaly. Thus convection is inhibited from propagating cyclonically around the core leading to a large down-shear vortex tilt. Conversely, in a counter-aligned orientation, the negative surface wind anomaly and driest air is found right of the shear vector. Hence the driest air rotates into the down-shear flank where it moistened by shear-organized convection. Furthermore, the boundary layer is relatively moist left of shear due to the positive surface wind anomaly, therefore promoting the cyclonic propagation from down-shear and constraining the magnitude of the vortex tilt. Genesis is intimately tied to the magnitude of the tilt and is found to occur once the mid-level vortex has precessed into the up-shear flank. For smaller values of maximum tilt, vortex precession is comparatively rapid, aided by “showerhead” moistening provided by the up-shear advection of frozen condensate aloft. With the up-shear flank pre-moistened, rapid precession of the mid-level vortex, at smaller radii, leads to near saturation on the mesoscale and the onset of rapid intensification. When the magnitude of the tilt is quite large, precession is much slower and the showerhead effect is significantly reduced until just prior to the emergence of the mid-level vortex in the up-shear flank. Copyright © 2011 Royal Meteorological Society

Description: Monitoring and quantifying future climate projections of dryness and wetness extremes: SPI bias Hydrology and Earth System Sciences Discussions, 8, 10635-10677, 2011 Author(s): F. Sienz, O. Bothe, and K. Fraedrich The adequacy of the Gamma distribution (GD) for monthly precipitation totals is reconsidered. The motivation for this study is the observation that the GD fails to represent precipitation in considerable areas of global observed and simulated data. This misrepresentation may lead to erroneous estimates of the Standardised Precipitation Index (SPI), evaluations of models, and assessments of climate change. In this study, the GD is compared to the Weibull (WD), Burr Type III (BD), exponentiated Weibull (EWD) and generalised Gamma (GGD) distribution. These distributions extend the GD in terms of possible shapes (skewness and kurtosis) and the behaviour for large arguments. The comparison is based on the Akaike information criterion, which maximises information entropy, and reveals a trade-off between deviation and the numbers of parameters used. We use monthly sums of observed and simulated precipitation for 12 calendar months of the year. Assessing observed and simulated data (i) the Weibull type distributions give distinctly improved fits compared to the GD and (ii) the SPI resulting from the GD overestimates (underestimates) extreme dryness (wetness).

Description: SWAT use of gridded observations for simulating runoff – a Vietnam river basin study Hydrology and Earth System Sciences Discussions, 8, 10679-10705, 2011 Author(s): M. T. Vu, S. V. Raghavan, and S. Y. Liong Many research studies that focus on basin hydrology have used the SWAT model to simulate runoff. One common practice in calibrating the SWAT model is the application of station data rainfall to simulate runoff. But over regions lacking robust station data, there is a problem of applying the model to study the hydrological responses. For some countries and remote areas, the rainfall data availability might be a constraint due to many different reasons such as lacking of technology, war time and financial limitation that lead to difficulty in constructing the runoff data. To overcome such a limitation, this research study uses some of the available globally gridded high resolution precipitation datasets to simulate runoff. Five popular gridded observation precipitation datasets: (1) Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources (APHRODITE), (2) Tropical Rainfall Measuring Mission (TRMM), (3) Precipitation Estimation from Remote Sensing Information using Artificial Neural Network (PERSIANN), (4) Global Precipitation Climatology Project (GPCP), (5) modified Global Historical Climatology Network version 2 (GHCN2) and one reanalysis dataset National Centers for Environment Prediction/National Center for Atmospheric Research (NCEP/NCAR) are used to simulate runoff over the Dakbla River (a small tributary of the Mekong River) in Vietnam. Wherever possible, available station data are also used for comparison. Bilinear interpolation of these gridded datasets is used to input the precipitation data at the closest grid points to the station locations. Sensitivity Analysis and Auto-calibration are performed for the SWAT model. The Nash-Sutcliffe Efficiency (NSE) and Coefficient of Determination ( R 2 ) indices are used to benchmark the model performance. This entails a good understanding of the response of the hydrological model to different datasets and a quantification of the uncertainties in these datasets. Such a methodology is also useful for planning on Rainfall-runoff and even reservoir/river management both at rural and urban scales.

Description: Relative humidity (RH) measurements, as derived from wet-bulb and dry-bulb thermometers operated as a psychrometer within a thermometer screen, have limited accuracy because of natural ventilation variations. Standard RH calculations generally assume a fixed screen psychrometer coefficient, but this is too small during poor ventilation. By comparing a reference humidity probe—exposed within a screen containing a psychrometer—with wind-speed measurements under controlled conditions, a wind-speed correction for the screen psychrometer coefficient has been derived and is applicable when 2-metre wind speeds fall below 3 m s −1 . Applying this to hourly-averaged data reduced the mean moist RH bias of the psychrometer (over the reference probe) from 1.2% to 0.4%, and reduced the interquartile range of the RH differences from 2.0% to 0.8%. This correction is particularly amenable to automatic measurement systems. Copyright © 2011 Royal Meteorological Society

Description: Mass transport of contaminated soil released into surface water by landslides (Göta River, SW Sweden) Hydrology and Earth System Sciences Discussions, 8, 10589-10633, 2011 Author(s): G. Göransson, M. Larson, D. Bendz, and M. Åkesson Landslides of contaminated soil into surface water represent an overlooked exposure pathway that has not been addressed properly in existing risk analysis for landslide hazard, contaminated land, or river basin management. A landslide of contaminated soil into surface water implies an instantaneous exposure of the water to the contaminated soil, dramatically changing the prerequisites for the mobilisation and transport of pollutants. In this study, an analytical approach is taken to simulate the transport of suspended matter released in connection with landslides into rivers. Different analytical solutions to the advection-dispersion equation (ADE) were tested against the measured data from the shallow rotational, retrogressive landslide in clayey sediments that took place in 1993 on the Göta River, SW Sweden. The landslide encompassed three distinct events, namely an initial submerged slide, followed by a main slide, and a retrogressive slide. These slides generated three distinct and non-Gaussian peaks in the online turbidity recordings at the freshwater intake downstream the slide area. To our knowledge, this registration of the impact in a river of the sediment release from a landslide is one of the few of its kind in the world, and unique for Sweden considering the low frequency of landslide events, making it highly useful for evaluating how appropriate the ADE is to describe a landslide into surface water. The results yielded realistic predictions of the measured concentration variation, after proper calibration. For the three individual slides it was estimated that a total of about 0.6% (515 000 kg) of the total landslide mass went into suspension/was suspended and was transported downstream. This release corresponds to about 1 to 2% of the annual suspended sediment delivery for that river stretch. The studied landslide partly involved an industrial area and by applying the analytical solution for the transport of metals in the sediments it was found that landslides have the possibility to release a significant amount of pollutants if large contaminated areas are involved. However, further studies are needed to develop more detailed descriptions of the transport processes. There is also a need to increase the knowledge on possible environmental consequences in the near and far field, in a short and long-time perspective. Finally, the risk for the release of pollutants should not be neglected in landslide hazard and risk assessment.

Description: In early January 2011, southern China experienced another severe wintertime storm as in the winter of 2008. The storm consisted of a narrow east-west-oriented freezing rain band over central Guizhou with an extensive area of snow to the north and a rain swath to the south. This study investigates this event using conventional surface and radiosonde data as well as final (FNL) analyses data from the US National Center for Atmospheric Research/National Centers for Environmental Prediction(NCAR/NCEP). The results show that forcing by geostrophic and diabatic heating built up a positive direct secondary circulation perpendicular to the quasi-stationary front beneath 700 hPa to maintain the surface cold layer and warm layer aloft through surface cold advection and warm advection aloft. In addition, turbulence induced by strong wind shear in the middle and lower stratiform clouds with a low concentration of ice nuclei plays an important role in the growth of cloud drops and the enhancement of supercool raindrop precipitation over Guizhou. Copyright © 2011 Royal Meteorological Society

Description: Prioritization of water management under climate change and urbanization using multi-criteria decision making methods Hydrology and Earth System Sciences Discussions, 8, 9889-9925, 2011 Author(s): J.-S. Yang, E.-S. Chung, S.-U. Kim, T.-W. Kim, and Y. D. Kim This paper quantifies the transformed effectiveness of alternatives for watershed management caused by climate change and urbanization and prioritizes five options using multi-criteria decision making techniques. The climate change scenarios (A1B and A2) were obtained by using a statistical downscaling model (SDSM), and the urbanization scenario by surveying the existing urban planning. The flow and biochemical oxygen demand (BOD) concentration duration curves were derived, and the numbers of days required to satisfy the environmental flow requirement and the target BOD concentration were counted using the Hydrological Simulation Program-Fortran (HSPF) model. In addition, five feasible alternatives were prioritized by using multi-criteria decision making techniques, based on the driving force-pressure-state-impact-response (DPSIR) framework and cost component. Finally, a sensitivity analysis approach for MCDM methods was conducted to reduce the uncertainty of weights. The result indicates that the most sensitive decision criterion is cost, followed by criteria response, driving force, impact, state and pressure in that order. Since it is certain that the importance of cost component is over 0.127, use of the groundwater collected by subway stations will be the most preferred alternative in this application.

Description: Influence of feedbacks from simulated crop growth on integrated regional hydrologic simulations under climate scenarios Hydrology and Earth System Sciences Discussions, 8, 10151-10193, 2011 Author(s): P. E. V. van Walsum Climate change impact modelling of hydrologic responses is hampered by climate-dependent model parameterizations. Reducing this dependency was one of the goals of extending the regional hydrologic modelling system SIMGRO with a two-way coupling to the crop growth simulation model WOFOST. The coupling includes feedbacks to the hydrologic model in terms of the root zone depth, soil cover, leaf area index, interception storage capacity, crop height and crop factor. For investigating whether such feedbacks lead to significantly different simulation results, two versions of the model coupling were set up for a test region: one with exogenous vegetation parameters, the "static" model, and one with endogenous simulation of the crop growth, the "dynamic" model WOFOST. The used parameterization methods of the static/dynamic vegetation models ensure that for the current climate the simulated long-term average of the actual evapotranspiration is the same for both models. Simulations were made for two climate scenarios. Owing to the higher temperatures in combination with a higher CO 2 -concentration of the atmosphere, a forward time shift of the crop development is simulated in the dynamic model; the used arable land crop, potatoes, also shows a shortening of the growing season. For this crop, a significant reduction of the potential transpiration is simulated compared to the static model, in the example by 15% in a warm, dry year. In consequence, the simulated crop water stress (the unit minus the relative transpiration) is lower when the dynamic model is used; also the simulated increase of crop water stress due to climate change is lower; in the example, the simulated increase is 15 percentage points less (of 55) than when a static model is used. The static/dynamic models also simulate different absolute values of the transpiration. The difference is most pronounced for potatoes at locations with ample moisture supply; this supply can either come from storage release of a good soil or from capillary rise. With good supply of moisture, the dynamic model simulates up to 10% less actual evapotranspiration than the static one in the example. This can lead to cases where the dynamic model predicts a slight increase of the recharge in a climate scenario, where the static model predicts a decrease. The use of a dynamic model also affects the simulated demand for surface water from external sources; especially the timing is affected. The proposed modelling approach uses postulated relationships that require validation with controlled field trials. In the Netherlands there is a lack of experimental facilities for performing such validations.