ALBERT

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

Description: Extending periodic eddy covariance latent heat fluxes through tree sapflow measurements to estimate long-term total evaporation in a peat swamp forest Hydrology and Earth System Sciences Discussions, 11, 13607-13661, 2014 Author(s): A. D. Clulow, C. S. Everson, M. G. Mengistu, J. S. Price, A. Nickless, and G. P. W. Jewitt A combination of measurement and modelling was used to find a pragmatic solution to estimate the annual total evaporation (ET) from the rare and indigenous Nkazana Peat Swamp Forest (PSF) on the east coast of Southern Africa to improve the water balance estimates within the area. Total evaporation was measured during three window periods (between seven and nine days each) using an eddy covariance (EC) system on a telescopic mast above the forest canopy. Sapflow of an understory and an emergent tree was measured using a low maintenance heat pulse velocity system for an entire hydrological year (October 2009 to September 2010). An empirical model was derived, describing the relationship between the observed ET of the Nkazana PSF measured during two of the window periods ( R 2 = 0.92 and 0.90) which, overlapped with sapflow measurements, thereby providing hourly estimates of predicted ET of the Nkazana PSF for a year, totalling 1125 mm (while rainfall was 650 mm). In building the empirical model, it was found that including the understory tree sapflow provided no benefit to the model performance. In addition, the observed emergent tree sapflow relationship with observed ET between the two field campaigns was consistent and could be represented by a single empirical model ( R 2 = 0.90; RMSE = 0.08 mm). During the window periods of EC measurement, no single meteorological variable was found to describe the Nkazana PSF ET satisfactorily. However, in terms of evaporation models, the hourly FAO56 Penman–Monteith equation best described the observed ET from EC during the August 2009 ( R 2 = 0.75), November 2009 ( R 2 = 0.85) and March 2010 ( R 2 = 0.76) field campaigns, compared to the Priestley–Taylor model ( R 2 = 0.54, 0.74 and 0.62 during the respective field campaigns). From the empirical model of ET and the FAO56 Penman–Monteith equation, a monthly crop factor ( K c ) was derived for the Nkazana PSF providing a method of estimating long-term swamp forest ET from meteorological data. The monthly crop factor indicated two distinct periods. From February to May, it was between 1.2 and 1.4 compared with June to January, when the crop factor was 0.8 to 1.0. The derived monthly K c values were verified as accurate (to one significant digit) using historical data measured at the same site, also using EC, from a~previous study. The measurements provided insights into the microclimate within a subtropical peat swamp forest and the contrasting sapflow of emergent and understory trees. They showed that expensive, high maintenance equipment can be used during manageable window periods in conjunction with low maintenance systems, dedicated to individual trees, to derive a model to estimate long-term ET over remote heterogeneous forests. In addition, the contrast in ET and rainfall emphasises the reliance of the Nkazana PSF on groundwater.

Description: Long-term variability of storm surge frequency in the Venice Lagoon: an update thanks to eighteenth century sea level observations Natural Hazards and Earth System Sciences Discussions, 2, 7465-7486, 2014 Author(s): F. Raicich Sea level observations made in the Venice Lagoon between 1751 and 1792 have been recovered, consisting of two time series of daily data on high and low waters at Venice and Chioggia. From comparisons with modern observations the quality of the 18th century data appears to be good enough to allow a useful analysis. A composite time series of daily mean sea level is obtained by merging the 18th century data and 1872–2004 observations at Venice Punta della Salute. The absence of reliable information on vertical references prevents the connection of the two 18th century time series with each other and with modern observations. However, daily sea level anomalies relative to the mean sea level enable to recognize storm surge events, that appear to occur more frequently in the second half of the 18th century than in the late 19th and 20th centuries, particularly during the 1751–1769 period. The record-breaking storm surge of 4 November 1966 turns out to be a remarkable event also in comparison with the events extracted from the 18th century time series. Further work is required to fill the gap between the old and modern observations.

Description: Ground-penetrating radar observations for estimating the vertical displacement of rotational landslides Natural Hazards and Earth System Sciences Discussions, 2, 7487-7506, 2014 Author(s): C. Lissak, O. Maquaire, J.-P. Malet, F. Lavigne, C. Virmoux, C. Gomez, and R. Davidson The objective of this paper is to demonstrate the applicability of Ground Penetrating Radar (GPR) for monitoring the displacement of slow-moving landslides. GPR data is used to estimate the vertical movement of rotational slides in combination with other surveying techniques. The experimental site is located along the Normandy coast (North East France) here several rotational landslides are continuously affected by a seasonal kinematic pattern (low displacement rates of 0.01 to 0.10 m yr −1 ) and periodically by major acceleration events (high displacement of 1.0 to 7.0 m per event).

Description: Continental Portuguese Territory Flood Susceptibility Index – contribution for a vulnerability index Natural Hazards and Earth System Sciences Discussions, 2, 7521-7552, 2014 Author(s): R. Jacinto, N. Grosso, E. Reis, L. Dias, F. D. Santos, and P. Garrett This work defines a national flood susceptibility index for the Portuguese continental territory, by proposing the aggregation of different variables which represent natural conditions for permeability, runoff and accumulation. This index is part of the national vulnerability index developed in the scope of Flood Maps in Climate Change Scenarios (CIRAC) project, supported by the Portuguese Association of Insurers (APS). This approach expands on previous works by trying to bridge the gap between different floods mechanisms (e.g. progressive and flash floods) occurring at different spatial scales in the Portuguese territory through: (a) selecting homogeneously processed datasets, (b) aggregating their values to better translate the spatially continuous and cumulative influence in floods at multiple spatial scales. Results show a good ability to capture, in the higher susceptibility classes, different flood types: progressive floods and flash floods. Lower values are usually related to: mountainous areas, low water accumulation potential and more permeable soils. Validation with independent flood datasets confirmed these index characteristics, although some overestimation can be seen in the southern region of Alentejo where, due to a dense hydrographic network and an overall low slope, floods are not as frequent as a result of lower precipitation mean values. Future work will focus on: (i) including extreme precipitation datasets to represent the triggering factor, (ii) improving representation of smaller and stepper basins, (iii) optimizing variable weight definition process, (iii) developing more robust independent flood validation datasets.

Description: Dye tracing for investigating flow and transport properties of hydrocarbon-polluted Rabots glaciär, Kebnekaise, Sweden Hydrology and Earth System Sciences Discussions, 11, 13711-13744, 2014 Author(s): C. C. Clason, C. Coch, J. Jarsjö, K. Brugger, P. Jansson, and G. Rosqvist Over 11 000 L of hydrocarbon pollution was deposited on the surface of Rabots glaciär on the Kebnekaise Massif, northern Sweden, following the crash of a Royal Norwegian Air Force aircraft in March 2012. An environmental monitoring programme was subsequently commissioned, including water, snow and ice sampling. The scientific programme further included a series of dye tracing experiments during the 2013 melt season, conducted to investigate flow pathways for pollutants through the glacier hydrological system, and to gain new insight to the internal hydrological system of Rabots glaciär. Results of dye tracing reveal a degree of homogeneity in the topology of the drainage system throughout July and August, with an increase in efficiency as the season progresses, as reflected by decreasing temporary storage and dispersivity. Early onset of melting likely led to formation of an efficient, discrete drainage system early in the melt season, subject to decreasing sinuosity and braiding as the season progressed. Analysis of turbidity-discharge hysteresis further supports the formation of discrete, efficient drainage, with clockwise diurnal hysteresis suggesting easy mobilisation of readily-available sediments in channels. Dye injection immediately downstream of the pollution source zone revealed prolonged storage of dye followed by fast, efficient release. Twinned with a low dye recovery, and supported by sporadic detection of hydrocarbons in the proglacial river, we suggest that meltwater, and thus pollutants in solution, may be released periodically from this zone of the glacier hydrological system. The here identified dynamics of dye storage, dispersion and breakthrough indicate that the ultimate fate and permanence of pollutants in the glacier system is likely to be governed by storage of pollutants in the firn layer and ice mass, or within the internal hydrological system, where it may refreeze. This shows that future studies on the fate of hydrocarbons in pristine, glaciated mountain environments should address the extent to which pollutants in solution act like water molecules or whether they are more susceptible to, for example, refreezing into the surrounding ice, becoming stuck in micro-fractures and pore spaces, or sorption onto subglacial sediments.

Description: Air temperature variability over three glaciers in the Ortles-Cevedale (Italian Alps): effects of glacier disintegration, intercomparison of calculation methods, and impacts on mass balance modeling The Cryosphere Discussions, 8, 6147-6192, 2014 Author(s): L. Carturan, F. Cazorzi, F. De Blasi, and G. Dalla Fontana Glacier mass balance models rely on accurate spatial calculation of input data, in particular air temperature. Lower temperatures (the so-called glacier cooling effect), and lower temperature variability (the so-called glacier damping effect) generally occur over glaciers, compared to ambient conditions. These effects, which depend on the geometric characteristics of glaciers and display a high spatial and temporal variability, have been mostly investigated on medium- to large-size glaciers so far, while observations on smaller ice bodies are scarce. Using a dataset from 8 on-glacier and 4 off-glacier weather stations, collected in summer 2010 and 2011, we analyzed the air temperature variability and wind regime over three different glaciers in the Ortles-Cevedale. The magnitude of the cooling effect and the occurrence of katabatic boundary layer (KBL) processes showed remarkable differences among the three ice bodies, suggesting the likely existence of important reinforcing mechanisms during glacier decay and disintegration. None of the methods proposed in the literature for calculating on-glacier temperature from off-glacier data fully reproduced our observations. Among them, the more physically-based procedure of Greuell and Böhm (1998) provided the best overall results where the KBL prevail, but it was not effective elsewhere (i.e. on smaller ice bodies and close to the glacier margins). The accuracy of air temperature estimations strongly impacted the results from a mass balance model which was applied to the three investigated glaciers. Most importantly, even small temperature deviations caused distortions in parameter calibration, thus compromising the model generalizability.

Description: Land surface temperature (LST) retrieval is a key issue in infrared quantitative remote sensing. In this paper, a split window algorithm is proposed to estimate LST with daytime data in two mid-infrared channels (channel 66 (3.746~4.084 μm) and channel 68 (4.418~4.785 μm)) from Airborne Hyperspectral Scanner (AHS). The estimation is conducted after eliminating reflected direct solar radiance with the aid of water vapor content (WVC), the view zenith angle (VZA), and the solar zenith angle (SZA). The results demonstrate that the LST can be well estimated with a root mean square error (RMSE) less than 1.0 K. Furthermore, an error analysis for the proposed method is also performed in terms of the uncertainty of LSE and WVC, as well as the Noise Equivalent Difference Temperature (NEΔT). The results show that the LST errors caused by a LSE uncertainty of 0.01, a NEΔT of 0.33 K, and a WVC uncertainty of 10% are 0.4~2.8 K, 0.6 K, and 0.2 K, respectively. Finally, the proposed method is applied to the AHS data of 4 July 2008. The results show that the differences between the estimated and the ground measured LST for water, bare soil and vegetation areas are approximately 0.7 K, 0.9 K and 2.3K, respectively.

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Inversion method for initial tsunami waveform reconstruction Natural Hazards and Earth System Sciences Discussions, 2, 7735-7772, 2014 Author(s): V. V. Voronin, T. A. Voronina, and V. A. Tcheverda This paper deals with the application of r-solution method to recover the initial tsunami waveform in a tsunami source area by remote water-level measurements. Wave propagation is considered within the scope of a linear shallow-water theory. An ill-posed inverse problem is regularized by means of least square inversion using a truncated SVD approach. The properties of obtained solution are determined to a large extent by the properties of an inverse operator, which were numerically investigated. The method presented allows one to control instability of the numerical solution and to obtain an acceptable result in spite of ill-posedness of the problem. It is shown that the accuracy of tsunami source reconstruction strongly depends on the signal-to-noise ratio, the azimuthal coverage of recording stations with respect to the source area and bathymetric features along the wave path. The numerical experiments were carried out with synthetic data and various computational domains including a real bathymetry. The method proposed allows us to make a preliminary prediction of the efficiency of the inversion with a given set of the recording stations and to find out the most informative part of the existing observation system. This essential property of the method can prove to be useful in designing a monitoring system for tsunamis.

Description: Light Detection and Ranging (LiDAR) is an active sensor that can effectively acquire a large number of three-dimensional (3-D) points. LiDAR systems can be equipped on different platforms for different applications, but to integrate the data, point cloud registration is needed to improve geometric consistency. The registration of airborne and terrestrial mobile LiDAR is a challenging task because the point densities and scanning directions differ. We proposed a scheme for the registration of airborne and terrestrial mobile LiDAR using the least squares 3-D surface registration technique to minimize the surfaces between two datasets. To analyze the effect of point density in registration, the simulation data simulated different conditions and estimated the theoretical errors. The test data were the point clouds of the airborne LiDAR system (ALS) and the mobile LiDAR system (MLS), which were acquired by Optech ALTM 3070 and Lynx, respectively. The resulting simulation analysis indicated that the accuracy of registration improved as the density increased. For the test dataset, the registration error of mobile LiDAR between different trajectories improved from 40 cm to 4 cm, and the registration error between ALS and MLS improved from 84 cm to 4 cm. These results indicate that the proposed methods can obtain 5 cm accuracy between ALS and MLS.

Description: Nighttime satellite imagery from the Defense Meteorological Satellite Programs’ Operational Linescan System (DMSP OLS) is being used for myriad applications including population mapping, characterizing economic activity, disaggregate estimation of CO2 emissions, wildfire monitoring, and more. Here we present a method for in situ radiance calibration of the DMSP OLS using a ground based light source as an active target. We found that the wattage of light used by our active target strongly correlates with the signal measured by the DMSP OLS. This approach can be used to enhance our ability to make intertemporal and intersatellite comparisons of DMSP OLS imagery. We recommend exploring the possibility of establishing a permanent active target for the calibration of nocturnal imaging systems.

Description: For near real-time water applications, the Moderate Resolution Imaging Spectroradiometers (MODIS) on Terra and Aqua are currently the only satellite instruments that can provide well-calibrated top-of-atmosphere (TOA) radiance data over the global aquatic environments. However, TOA radiance data in the MODIS ocean bands over turbid atmosphere in east China often saturate, leaving only four land bands to use. In this study, an approach based on Empirical Orthogonal Function (EOF) analysis has been developed and validated to estimate chlorophyll a concentrations (Chla, μg/L) in surface waters of Taihu Lake, the third largest freshwater lake in China. The EOF approach analyzed the spectral variance of normalized Rayleigh-corrected reflectance (Rrc) data at 469, 555, 645, and 859 nm, and subsequently related that variance to Chla using 28 concurrent MODIS and field measurements. This empirical algorithm was then validated using another 30 independent concurrent MODIS and field measurements. Image analysis and radiative transfer simulations indicated that the algorithm appeared to be tolerant to aerosol perturbations, with unbiased RMS uncertainties of 〈80% for Chla ranging between 3 and 100 μg/L. Application of the algorithm to a total of 853 MODIS images between 2000 and 2013 under cloud-free conditions revealed spatial distribution patterns and seasonal changes that are consistent to previous findings based on floating algae mats. The current study can provide additional quantitative estimates of Chla that can be assimilated in an existing forecast model, which showed improved performance over the use of a previous Chla algorithm. However, the empirical nature, relatively large uncertainties, and limited number of spectral bands all point to the need of further improvement in data availability and accuracy with future satellite sensors.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Strategy of valid 14 C dates choice in syngenetic permafrost The Cryosphere Discussions, 8, 5589-5621, 2014 Author(s): Y. K. Vasil'chuk and A. C. Vasil'chuk The main problem of radiocarbon dating within permafrost is the uncertain reliability of the 14 C dates. Syngenetic sediments contain allochthonous organic deposit that originated at a distance from its present position. Due to the very good preservation of organic materials in permafrost conditions and numerous re-burials of the fossils from ancient deposits into younger ones the dates could be both younger and older than the true age of dated material. The strategy for the most authentic radiocarbon date selection for dating of syncryogenic sediments is considered taking into account the fluvial origin of the syngenetic sediments. The re-deposition of organic material is discussed in terms of cyclic syncryogenic sedimentation and also the possible re-deposition of organic material in subaerial-subaqueous conditions. The advantages and the complications of dating organic micro-inclusions from ice wedges by the accelerator mass spectrometry (AMS) method are discussed applying to true age of dated material search. Radiocarbon dates of different organic materials from the same samples are compared. The younger age of the yedoma from cross-sections of Duvanny Yar in Kolyma River and Mamontova Khayata in the mouth of Lena River is substantiated due to the principle of the choice of the youngest 14 C date from the set.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Simulating the Antarctic ice sheet in the Late-Pliocene warm period: PLISMIP-ANT, an ice-sheet model intercomparison project The Cryosphere Discussions, 8, 5539-5588, 2014 Author(s): B. de Boer, A. M. Dolan, J. Bernales, E. Gasson, H. Goelzer, N. R. Golledge, J. Sutter, P. Huybrechts, G. Lohmann, I. Rogozhina, A. Abe-Ouchi, F. Saito, and R. S. W. van de Wal In the context of future climate change, understanding the nature and behaviour of ice sheets during warm intervals in Earth history is of fundamental importance. The Late-Pliocene warm period (also known as the PRISM interval: 3.264 to 3.025 million years before present) can serve as a potential analogue for projected future climates. Although Pliocene ice locations and extents are still poorly constrained, a significant contribution to sea-level rise should be expected from both the Greenland ice sheet and the West and East Antarctic ice sheets based on palaeo sea-level reconstructions. Here, we present results from simulations of the Antarctic ice sheet by means of an international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP-ANT). For the experiments, ice-sheet models including the shallow ice and shelf approximations have been used to simulate the complete Antarctic domain (including grounded and floating ice). We compare the performance of six existing numerical ice-sheet models in simulating modern control and Pliocene ice sheets by a suite of four sensitivity experiments. Ice-sheet model forcing fields are taken from the HadCM3 atmosphere–ocean climate model runs for the pre-industrial and the Pliocene. We include an overview of the different ice-sheet models used and how specific model configurations influence the resulting Pliocene Antarctic ice sheet. The six ice-sheet models simulate a comparable present-day ice sheet, although the models are setup with their own parameter settings. For the Pliocene simulations using the Bedmap1 bedrock topography, some models show a small retreat of the East Antarctic ice sheet, which is thought to have happened during the Pliocene for the Wilkes and Aurora basins. This can be ascribed to either the surface mass balance, as the HadCM3 Pliocene climate shows a significant increase over the Wilkes and Aurora basin, or the initial bedrock topography. For the latter, our simulations with the recently published Bedmap2 bedrock topography indicate a significantly larger contribution to Pliocene sea-level rise from the East Antarctic ice sheet for all six models relative to the simulations with Bedmap1.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

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

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

Description: Building change detection is useful for land management, disaster assessment, illegal building identification, urban growth monitoring, and geographic information database updating. This study proposes an automatic method that applies object-based analysis to multi-temporal point cloud data to detect building changes. The aim of this building change detection method is to identify areas that have changed and to obtain from-to information. In this method, the data are first preprocessed to generate two sets of digital surface models (DSMs), digital elevation models, and normalized DSMs from registered old and new point cloud data. Thereafter, on the basis of differential DSM, candidates for changed building objects are identified from the points in the smooth areas by using a connected component analysis technique. The random sample consensus fitting algorithm is then used to distinguish the true changed buildings from trees. The changed building objects are classified as “newly built”, “taller”, “demolished” or “lower” by using rule-based analysis. Finally, a test data set consisting of many buildings of different types in an 8.5 km2 area is selected for the experiment. In the test data set, the method correctly detects 97.8% of buildings larger than 50 m2. The accuracy of the method is 91.2%. Furthermore, to decrease the workload of subsequent manual checking of the result, the confidence index for each changed object is computed on the basis of object features.

Description: In this study, we test and demonstrate the utility of disturbance and recovery information derived from annual Landsat time series to predict current forest vertical structure (as compared to the more common approaches, that consider a sample of airborne Lidar and single-date Landsat derived variables). Mean Canopy Height (MCH) was estimated separately using single date, time series, and the combination of single date and time series variables in multiple regression and random forest (RF) models. The combination of single date and time series variables, which integrate disturbance history over the entire time series, overall provided better MCH prediction than using either of the two sets of variables separately. In general, the RF models resulted in improved performance in all estimates over those using multiple regression. The lowest validation error was obtained using Landsat time series variables in a RF model (R2 = 0.75 and RMSE = 2.81 m). Combining single date and time series data was more effective when the RF model was used (opposed to multiple regression). The RMSE for RF mean canopy height prediction was reduced by 13.5% when combining the two sets of variables as compared to the 3.6% RMSE decline presented by multiple regression. This study demonstrates the value of airborne Lidar and long term Landsat observations to generate estimates of forest canopy height using the random forest algorithm.

Description: Rice is the most important food security crop in Asia. Information on its seasonal extent forms part of the national accounting of many Asian countries. Synthetic Aperture Radar (SAR) imagery is highly suitable for detecting lowland rice, especially in tropical and subtropical regions, where pervasive cloud cover in the rainy seasons precludes the use of optical imagery. Here, we present a simple, robust, rule-based classification for mapping rice area with regularly acquired, multi-temporal, X-band, HH-polarized SAR imagery and site-specific parameters for classification. The rules for rice detection are based on the well-studied temporal signature of rice from SAR backscatter and its relationship with crop stages. We also present a procedure for estimating the parameters based on “temporal feature descriptors” that concisely characterize the key information in the rice signatures in monitored field locations within each site. We demonstrate the robustness of the approach on a very large dataset. A total of 127 images across 13 footprints in six countries in Asia were obtained between October 2012, and April 2014, covering 4.78 m ha. More than 1900 in-season site visits were conducted across 228 monitoring locations in the footprints for classification purposes, and more than 1300 field observations were made for accuracy assessment. Some 1.6 m ha of rice were mapped with classification accuracies from 85% to 95% based on the parameters that were closely related to the observed temporal feature descriptors derived for each site. The 13 sites capture much of the diversity in water management, crop establishment and maturity in South and Southeast Asia. The study demonstrates the feasibility of rice detection at the national scale using multi-temporal SAR imagery with robust classification methods and parameters that are based on the knowledge of the temporal dynamics of the rice crop. We highlight the need for the development of an open-access library of temporal signatures, further investigation into temporal feature descriptors and better ancillary data to reduce the risk of misclassification with surfaces that have temporal backscatter dynamics similar to those of rice. We conclude with observations on the need to define appropriate SAR acquisition plans to support policies and decisions related to food security.

Description: Snow mass decrease in the Northern Hemisphere (1979/80–2010/11) The Cryosphere Discussions, 8, 5623-5644, 2014 Author(s): Z. Li, J. Liu, L. Huang, N. Wang, B. Tian, J. Zhou, Q. Chen, and P. Zhang Snow cover has a key effect on climate change and hydrological cycling, as well as water supply to a sixth of the world's population across the Northern Hemisphere. However, reliable data on trends in snow cover in the Northern Hemisphere is lacking. Snow water equivalent (SWE) is a common measure of the amount of equivalent water of the snow pack. Here we verify the accuracy of three existing global SWE products and merge the most accurate aspects of them to generate a new SWE product covering the last 32 years (1979/80–2010/11). Using this new SWE product, we show that there has been a significant decreasing trend in the total mass of snow in the Northern Hemisphere. The most notable changes in total snow mass are −16.45 ± 6.68 and −13.55 ± 7.80 Gt year −1 in January and February, respectively. These are followed by March and December, which have trends of −12.58 ± 6.88 and −10.70 ± 5.62 Gt year −1 , respectively, from 1979/80 to 2010/11. During the same period, the temperature in the study area raised 0.17 °C decade −1 , which is thought to be the main reason of SWE decline.

Description: Factors controlling Slope Environmental Lapse Rate (SELR) of temperature in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya The Cryosphere Discussions, 8, 5645-5686, 2014 Author(s): R. J. Thayyen and A. P. Dimri Moisture, temperature and precipitation interplay forced through the orographic processes sustains the Himalayan cryospheric system. However, factors controlling the Slope Environmental Lapse Rate (SELR) of temperature along the higher Himalayan mountain slopes across various glacio-hydrologic regimes remain as a key knowledge gap. Present study dwells on the orographic processes driving the moisture–temperature interplay in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya. Systematic data collection at three altitudes between 2540 and 3763 m a.s.l. in the Garhwal Himalaya (hereafter called monsoon regime) and between 3500 and 5600 m a.s.l. in the Ladakh Himalaya (herefater called cold-arid regime) revealed moistrue control on temperature distribution at temporal and spatial scales. Observed daily SELR of temperature ranges between 9.0 to 1.9 °C km −1 and 17.0 to 2.8 °C km −1 in the monsoon and cold-arid regimes respectively highlighting strong regional variability. Moisture influx to the region, either from Indian summer monsoon (ISM) or from Indian winter monsoon (IWM) forced lowering of SELR. This phenophena of "monsoon lowering" of SELR is due to the release latent heat of condensation from orographically focred lifted air parcel. Seasonal response of SELR in the monsoon regime is found to be closly linked with the variations in the local lifting condensation levels (LCL). Contrary to this, cold-arid system is characterised by the extremely high values of daily SELR upto 17 °C km −1 signifying the extremely arid conditions prevailing in summer. Distinctly lower SELR devoid of monsoon lowering at higher altitude sections of monsoon and cold-arid regimes suggests sustained wetter high altitude regimes. We have proposed a SELR model for both glacio-hydrological regimes demostrating with two sections each using a derivative of the Clausius–Clapeyron relationship by deriving monthly SELR indices. It has been proposed that the manifestations of presence or absence of moisture is the single most important factor determining the temperature distribution along the higher Himalayan slopes driven by the orographic forcings. This work also suggests that the arbitary use of temperature lapse rate to extrapolate temperature to the higher Himalaya is extremely untenable.

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

Description: Sampling frequency trade-offs in the assessment of mean transit times of tropical montane catchment waters under semi-steady-state conditions Hydrology and Earth System Sciences Discussions, 11, 12443-12488, 2014 Author(s): E. Timbe, D. Windhorst, R. Celleri, L. Timbe, P. Crespo, H.-G. Frede, J. Feyen, and L. Breuer Stream and soil waters were collected on a weekly basis in a tropical montane cloud forest catchment for two years and analyzed for stable water isotopes in order to infer transit time distribution functions and to define the mean transit times. Depending on the water type (stream or soil water), lumped distribution functions such as Exponential-Piston flow, Linear-Piston flow and Gamma models using temporal isotopic variations of precipitation event samples as input, were fitted. Samples were aggregated to daily, weekly, biweekly, monthly and bimonthly time scales in order to check the sensitivity of temporal sampling on model predictions. The study reveals that the effect of decreasing sampling frequency depends on the water type. For soil waters with transit times in the order of weeks to months, there was a clear trend of over prediction. In contrast, the trend of prediction for stream waters, with a dampened isotopic signal and mean transit times in the order of 2 to 4 years, was less clear and depending on the type of model used. The trade-off to coarse data resolutions could potentially lead to misleading conclusions on how water actually moves through the catchment, while at the same time predictions can reach better fitting efficiencies, lesser uncertainties, errors and biases. For both water types an optimal sampling frequency seems to be one or at most two weeks. The results of our analyses provide information for the planning (in particular in terms of cost-benefit and time requirements) of future fieldwork in similar Andean or other catchments.

Description: Global insurance markets are vast and diverse, and may offer many opportunities for remote sensing. To date, however, few operational applications of remote sensing for insurance exist. Papers claiming potential application of remote sensing typically stress the technical possibilities, without considering its contribution to customer value for the insured or to the profitability of the insurance industry. Based on a systematic search of available literature, this review investigates the potential and actual support of remote sensing to the insurance industry. The review reveals that research on remote sensing in classical claim-based insurance described in the literature revolve around crop damage and flood and fire risk assessment. Surprisingly, the use of remote sensing in claim-based insurance appears to be instigated by government rather than the insurance industry. In contrast, insurance companies are offering various index insurance products that are based on remote sensing. For example, remotely sensed index insurance for rangelands and livestock are operational, while various applications in crop index insurance are being considered or under development. The paper discusses these differences and concludes that there is particular scope for application of remote sensing by the insurance industry in index insurance because (1) indices can be constructed that correlate well with what is insured; (2) these indices can be delivered at low cost; and (3) it opens up new markets that are not served by claim-based insurance. The paper finally suggests that limited adoption of remote sensing in insurance results from a lack of mutual understanding and calls for greater cooperation between the insurance industry and the remote sensing community.

Description: Vegetation cover can affect the lithological mapping capability of space- and airborne instruments because it obscures the spectral signatures of the underlying geological substrate. Despite being widely accepted as a hindrance, few studies have explicitly demonstrated the impact vegetation can have on remote lithological mapping. Accordingly, this study comprehensively elucidates the impact of vegetation on the lithological mapping capability of airborne multispectral data in the Troodos region, Cyprus. Synthetic spectral mixtures were first used to quantify the potential impact vegetation cover might have on spectral recognition and remote mapping of different rock types. The modeled effects of green grass were apparent in the spectra of low albedo lithologies for 30%–40% fractional cover, compared to just 20% for dry grass cover. Lichen was found to obscure the spectra for 30%–50% cover, depending on the spectral contrast between bare rock and lichen cover. The subsequent impact of vegetation on the remote mapping capability is elucidated by considering the outcomes of three airborne multispectral lithological classifications alongside the spectral mixing analysis and field observations. Vegetation abundance was found to be the primary control on the inability to classify large proportions of pixels in the imagery. Matched Filtering outperformed direct spectral matching algorithms owing to its ability to partially unmix pixel spectra with vegetation abundance above the modeled limits. This study highlights that despite the limited spectral sampling and resolution of the sensor and dense, ubiquitous vegetation cover, useful lithological information can be extracted using an appropriate algorithm. Furthermore, the findings of this case study provide a useful insight to the potential capabilities and challenges faced when utilizing comparable sensors (e.g., Landsat 8, Sentinel-2, WorldView-3) to map similar types of terrain.

Description: The origin of large areas dominated by pristine open vegetation that is in sharp contrast with surrounding dense forest within the Amazonian lowland has generally been related to past arid climates, but this is still an issue open for debate. In this paper, we characterize a large open vegetation patch over a residual megafan located in the northern Amazonia. The main goal was to investigate the relationship between this paleolandform and vegetation classes mapped based on the integration of optical and SAR data using the decision tree. Our remote sensing dataset includes PALSAR and TM/Landsat images. Five classes were identified: rainforest; flooded forest; wooded open vegetation; grassy-shrubby open vegetation; and water body. The output map resulting from the integration of PALSAR and TM/Landsat images showed an overall accuracy of 94%. Narrow, elongated and sinuous belts of forest within the open vegetation areas progressively bifurcate into others revealing paleochannels arranged into distributary pattern. Such characteristics, integrated with pre-existing geological information, led us to propose that the distribution of vegetation classes highlight a morphology attributed to a Quaternary megafan developed previous to the modern fluvial tributary system. The characterization of such megafan is important for reconstructing landscape changes associated with the evolution of the Amazon drainage basin.

Description: Flash Flood Guidance consists of indices that estimate the amount of rain of a certain duration that is needed over a given small basin in order to cause minor flooding. Backwater catchment inundation from swollen rivers or regional groundwater inputs are not significant over the spatial and temporal scales for the majority of upland flash flood prone basins, as such, these effects are not considered. However, some lowland areas and flat terrain near large rivers experience standing water long after local precipitation has ceased. NASA is producing an experimental product from the MODIS that detects standing water. These observations were assimilated into the hydrologic model in order to more accurately represent soil moisture conditions within basins, from sources of water from outside of the basin. Based on the upper soil water content, relations are used to derive an error estimate for the modeled soil saturation fraction; whereby, the soil saturation fraction model state can be updated given the availability of satellite observed inundation. Model error estimates were used in a Monte Carlo ensemble forecast of soil water and flash flood potential. Numerical experiments with six months of data (July 2011–December 2011) showed that MODIS inundation data, when assimilated to correct soil moisture estimates, increased the likelihood that bankfull flow would occur, over non-assimilated modeling, at catchment outlets for approximately 44% of basin-days during the study time period. While this is a much more realistic representation of conditions, no actual events occurred allowing for validation during the time period.

Description: For a foam-covered sea surface, it is difficult to retrieve sea surface salinity (SSS) with L-band brightness temperature (1.4 GHz) because of the effect of a foam layer with wind speeds stronger than 7 m/s, especially at low sea surface temperature (SST). With foam-controlled experiments, emissivities of a foam-covered water surface at low SST (−1.4 °C to 1.7 °C) are measured for varying SSS, foam thickness, incidence angle, and polarization. Furthermore, a theoretical model of emissivity is introduced by combining wave approach theory with the effective medium approximation method. Good agreement is obtained upon comparing theoretical emissivities with those of experiments. The results indicate that foam parameters have a strong influence on increasing emissivity of a foam-covered water surface. Increments of experimental emissivities caused by foam thickness of 1 cm increase from about 0.014 to 0.131 for horizontal polarization and 0.022 to 0.150 for vertical polarization with SSS increase and SST decrease. Contributions of the interface between the foam layer and water surface to the foam layer emissivity increments are discussed for frequencies between 1 and 37 GHz.

Description: A model study of Abrahamsenbreen, a surging glacier in northern Spitsbergen The Cryosphere Discussions, 8, 5687-5726, 2014 Author(s): J. Oerlemans and W. J. J. van Pelt The climate sensitivity of Abrahamsenbreen, a 20 km long surge-type glacier in northern Spitsbergen, is studied with a simple glacier model. A scheme to describe the surges is included, which makes it possible to account for the effect of surges on the total mass budget of the glacier. A climate reconstruction back to AD 1300, based on ice-core data from Lomonosovfonna and climate records from Longyearbyen, is used to drive the model. The model is calibrated by requesting that it produces the correct Little Ice Age maximum glacier length and simulates the observed magnitude of the 1978-surge. Abrahamsenbreen is strongly out of balance with the current climate. If climatic conditions will remain as they were for the period 1989–2010, the glacier will ultimately shrink to a length of about 4 km (but this will take hundreds of years). For a climate change scenario involving a 2 m yr −1 rise of the equilibrium line from now onwards, we predict that in the year 2100 Abrahamsenbreen will be about 12 km long. The main effect of a surge is to lower the mean surface elevation and to increase the ablation area, thereby causing a negative perturbation of the mass budget. We found that the occurrence of surges leads to a somewhat stronger retreat of the glacier in a warming climate. Because of the very small bed slope, Abrahamsenbreen is sensitive to small perturbations in the equilibrium-line altitude E . For a decrease of E of only 160 m, the glacier would steadily grow into the Woodfjorddalen until after 2000 years it would reach the Woodfjord and calving could slow down the advance.

Description: In contrast to ground-based solar-induced chlorophyll fluorescence (Fs) detection, the influence of atmospheric radiation transfer is the major difficulty in Fs retrieval from space. In this study, we first simulated top-of-atmosphere (TOA) radiance using FluorMODgui3.1 and MODTRAN5 code. Based on the simulated dataset, we analyzed the sensitivities of five potential Fs retrieval bands (Hα, K I, Fe, O2-A, and O2-B) to different atmospheric transfer parameters, including atmosphere profile, aerosol optical depth (AOD550), vertical water vapor column (H2O), vertical ozone column (O3), solar zenith angle (SZA), view zenith angle (VZA), relative azimuth angle (RAA) and elevation. The results demonstrate that the Hα, O2-A and O2-B bands are the most sensitive to these atmospheric parameters. However, only the O2-A and O2-B bands were found to be sensitive to the imaging geometric parameters. When the spectral resolution was sufficient, the K I and Fe bands proved to have the best potential for space-based Fs retrieval given the current available accuracies of atmospheric products, while the O2-A band was shown to perform better at lower spectral resolutions. The band sensitivity analysis presented here will be useful for band selection and atmospheric correction for space-based Fs retrieval.

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

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

Description: Quantitative high-resolution observations of soil water dynamics in a complicated architecture with time-lapse Ground-Penetrating Radar Hydrology and Earth System Sciences Discussions, 11, 12365-12404, 2014 Author(s): P. Klenk, S. Jaumann, and K. Roth High-resolution time-lapse Ground-Penetrating Radar (GPR) observations of advancing and retreating water tables can yield a wealth of information about near-surface water content dynamics. In this study, we present and analyze a series of imbibition, drainage and infiltration experiments which have been carried out at our artificial ASSESS test site and observed with surface based GPR. The test site features a complicated but known subsurface architecture constructed with three different kinds of sand. It allows studying soil water dynamics with GPR under a wide range of different conditions. Here, we assess in particular (i) the accurate determination of soil water dynamics averaged over the whole vertical extent by evaluating the bottom reflection and (ii) the feasibility of monitoring the dynamic shape of the capillary fringe reflection. The phenomenology of the GPR response of a dynamically changing capillary fringe is developed from a soil physical point of view. We then explain experimentally observed phenomena based on numerical simulations of both the water content dynamics and the expected GPR response.

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

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

Description: Simulating long-term past changes in the balance between water demand and availability and assessing their main drivers at the river basin management scale Hydrology and Earth System Sciences Discussions, 11, 12315-12364, 2014 Author(s): J. Fabre, D. Ruelland, A. Dezetter, and B. Grouillet The aim of this study was to assess the balance between water demand and availability and its spatial and temporal variability from 1971 to 2009 in the Herault (2500 km 2 , France) and the Ebro (85 000 km 2 , Spain) catchments. Natural streamflow was evaluated using a conceptual hydrological model. The regulation of river flow was accounted for through a widely applicable demand-driven reservoir management model applied to the largest dam in the Herault basin and to 11 major dams in the Ebro basin. Urban water demand was estimated from population and monthly unit water consumption data. Water demand for irrigation was computed from irrigated area, crop and soil data, and climatic forcing. Finally, a series of indicators comparing water supply and water demand at strategic resource and demand nodes were computed at a 10 day time step. Variations in water stress in each catchment over the past 40 years were successfully modeled, taking into account climatic and anthropogenic pressures and changes in water management strategies over time. Observed changes in discharge were explained by separating human and hydro-climatic pressures on water resources: respectively 20 and 3% of the decrease in the Ebro and the Herault discharges were linked to human-induced changes. Although key areas of the Herault basin were shown to be highly sensitive to hydro-climatic variability, the balance between water uses and availability in the Ebro basin appears to be more critical, owing to high agricultural pressure on water resources. The proposed modeling framework is currently being used to assess water stress under climatic and socio-economic prospective scenarios. Further research will investigate the effectiveness of adaptation policies aimed at maintaining the balance between water use and availability.

Description: Time-series analysis of the long-term hydrologic impacts of afforestation in the Águeda watershed of North-Central Portugal Hydrology and Earth System Sciences Discussions, 11, 12223-12256, 2014 Author(s): D. Hawtree, J. P. Nunes, J. J. Keizer, R. Jacinto, J. Santos, M. E. Rial-Rivas, A.-K. Boulet, F. Tavares-Wahren, and K.-H. Feger The north-central region of Portugal has undergone significant afforestation of the species Pinus pinaster and Eucalyptus globulus since the early 1900s; however, the long-term hydrologic impacts of this land cover change are not fully understood. To contribute to a better understanding of the potential hydrologic impacts of this land cover change, this study examines the temporal trends in 7 years of data from the Águeda watershed (part of the Vouga Basin) over the period of 1936 to 2010. Meteorological and hydrological records were analysed using a combined Thiel–Sen/Mann–Kendall trend testing approach, to assess the magnitude and significance of patterns in the observed data. These trend tests indicated that there had been no significant reduction in streamflow yield over either the entire test period, or during sub-record periods, despite the large-scale afforestation which had taken place. This lack of change is attributed to both the characteristics of the watershed and the nature of the land cover change. By contrast, a number of significant trends were found for baseflow index, which showed positive trends in the early data record (primarily during Pinus pinaster afforestation), followed by a reversal to negative trends later in the data record (primarily during Eucalyptus globulus afforestation). These changes are attributed to vegetation impacts on streamflow generating processes, both due to the species differences and to alterations in soil properties (i.e. promoting water repellency of the topsoil). These results highlight the importance of considering both vegetation types/dynamics and watershed characteristic when assessing hydrologic impacts, in particular with respect to soil properties.

Description: Reply to D. L. Peters' comment on "Streamflow input to Lake Athabasca, Canada" by Rasouli et al. (2013) Hydrology and Earth System Sciences Discussions, 11, 12257-12270, 2014 Author(s): K. Rasouli, M. A. Hernández-Henríquez, and S. J. Déry This paper provides a reply to a comment from Peters (2014) on our recent effort focused on evaluating changes in streamflow input to Lake Athabasca, Canada. Lake Athabasca experienced a 21.2% decline in streamflow input between 1960 and 2010 that has led to a marked decline in its water levels in recent decades. A reassessment of trends in naturalized Lake Athabasca water levels shows insignificant changes from our previous findings reported in Rasouli et al. (2013), and hence our previous conclusions remain unchanged. The reply closes with recommendations for future research to minimize uncertainties in historical assessments of trends in Lake Athabasca water levels and to better project its future water levels driven by climate change and anthropogenic activities in the Athabasca Lake Basin.

Description: Is sinuosity a function of slope and bankfull discharge? – A case study of the meandering rivers in the Pannonian Basin Hydrology and Earth System Sciences Discussions, 11, 12271-12290, 2014 Author(s): J. Petrovszki, G. Timár, and G. Molnár Pre-regulation channel sinuosities of the meandering rivers of the Pannonian Basin are analysed in order to define a mathematical model to estimate the influence of the bankfull discharge and the channel slope on them. As a primary database, data triplets of slope, discharge and sinuosity values were extracted from historical and modern datasets and pre-regulation historical topographic maps. Channel slope values were systematically modified to estimate figures valid before the river regulation works. The bankfull discharges were estimated from the average discharges using a robust yet complex method. The "classical" graphs of Leopold and Wolman (1957), Ackers and Charlton (1970b) and Schumm and Khan (1972) were compiled to a set up a theoretical surface, whose parameters are estimated by the real values of the above database, containing characteristics of the Pannonian Basin rivers. As a result it occurred that there is a two-dimensional function of the bankfull discharges, which provides a good estimation of the most probable sinuosity values of the rivers with the given slope and discharge characteristics. The average RMS error of this estimation is around 15% on this dataset and believed to be the effect of the non-analysed changes in the sediment discharge and size distribution.

Description: Technical Note: Field experiences using UV/VIS sensors for high-resolution monitoring of nitrate in groundwater Hydrology and Earth System Sciences Discussions, 11, 12291-12314, 2014 Author(s): M. Huebsch, F. Grimmeisen, M. Zemann, O. Fenton, K. G. Richards, P. Jordan, A. Sawarieh, P. Blum, and N. Goldscheider Two different in-situ spectrophotometers are compared that were used in the field to determine nitrate-nitrogen (NO 3 -N) concentrations at two distinct spring discharge sites. One sensor was a double wavelength spectrophotometer (DWS) and the other a multiple wavelength spectrophotometer (MWS). The objective of the study was to review the hardware options, determine ease of calibration, accuracy, influence of additional substances and to assess positive and negative aspects of the two sensors as well as troubleshooting and trade-offs. Both sensors are sufficient to monitor highly time-resolved NO 3 -N concentrations in emergent groundwater. However, the chosen path length of the sensors had a significant influence on the sensitivity and the range of detectable NO 3 -N. The accuracy of the calculated NO 3 -N concentrations of the sensors can be affected, if the content of additional substances such as turbidity, organic matter, nitrite or hydrogen carbonate significantly varies after the sensors have been calibrated to a particular water matrix. The MWS offers more possibilities for calibration and error detection, but requires more expertise compared with the DWS.

Description: A catalog of high-impact windstorms in Switzerland since 1859 Natural Hazards and Earth System Science, 14, 2867-2882, 2014 Author(s): P. Stucki, S. Brönnimann, O. Martius, C. Welker, M. Imhof, N. von Wattenwyl, and N. Philipp In recent decades, extremely hazardous windstorms have caused enormous losses to buildings, infrastructure and forests in Switzerland. This has increased societal and scientific interest in the intensity and frequency of historical high-impact storms. However, high-resolution wind data and damage statistics mostly span recent decades only. For this study, we collected quantitative (e.g., volumes of windfall timber, losses relating to buildings) and descriptive (e.g., forestry or insurance reports) information on the impact of historical windstorms. To define windstorm severity, normalized and declustered quantitative data were processed by extreme value statistics. Descriptive information was classified using a conceptual guideline. Validation with independent damage information, as well as comparison with wind measurements and a reanalysis, indicates that the most hazardous winter storms are captured, while too few moderate windstorms are detected. Strong storms in the wind measurements and reanalysis are thus added to the catalog. The final catalog encompasses approximately 240 high-impact windstorms in Switzerland since 1859. It features three robust severity classes and contains eight extreme windstorms. Evidence of high winter storm activity in the early and late 20th century compared to the mid-20th century in both damage and wind data indicates a co-variability of hazard and related damage on decadal timescales.

Description: Generalized drought assessment in Dongliao river basin based on water resources system Natural Hazards and Earth System Sciences Discussions, 2, 6703-6746, 2014 Author(s): B. S. Weng, D. H. Yan, H. Wang, T. L. Qin, and J. Yin Drought is firstly a resource issue, and with its development it transforms into a disaster issue. The occurrences of drought events usually feature determinacy and randomness. Drought issue has become one of the major factors to affect sustainable economic and social development. In this paper, we propose the generalized drought assessment index (GDAI) based on water resources system for assessing drought events. The GDAI considers water supply and water demand using a distributed hydrological model. We demonstrate the use of the proposed index in the Dongliao river basin (DRB) in the northeast China. The results simulated by the GDAI are then compared to observed drought disaster records in DRB. As second, the temporal distribution of drought events and the spatial distribution of drought frequency from the GDAI are compared with the traditional approach (i.e. the SPI, the PDSI, and the RWD). Then, generalized drought times (GDT), generalized drought duration (GDD), and generalized drought severity (GDS) were calculated by theory of runs. Application of the GDT, the GDD, and the GDS of various drought levels (i.e. mild drought, moderate drought, severe drought, and extreme drought) to the period 1960–2010 shows that the centers of gravity of them are all distributed in the middle reached of DRB, and change with time. The proposed methodology helps water managers in water-stressed regions to quantify the impact of drought, consequently, to make decisions regarding coping with drought issue.

Description: A GIS-based model to estimate flood consequences and the degree of accessibility and operability of strategic emergency response structures in urban areas Natural Hazards and Earth System Science, 14, 2847-2865, 2014 Author(s): R. Albano, A. Sole, J. Adamowski, and L. Mancusi Efficient decision-making regarding flood risk reduction has become a priority for authorities and stakeholders in many European countries. Risk analysis methods and techniques are a useful tool for evaluating costs and benefits of possible interventions. Within this context, a methodology to estimate flood consequences was developed in this paper that is based on GIS, and integrated with a model that estimates the degree of accessibility and operability of strategic emergency response structures in an urban area. The majority of the currently available approaches do not properly analyse road network connections and dependencies within systems, and as such a loss of roads could cause significant damages and problems to emergency services in cases of flooding. The proposed model is unique in that it provides a maximum-impact estimation of flood consequences on the basis of the operability of the strategic emergency structures in an urban area, their accessibility, and connection within the urban system of a city (i.e. connection between aid centres and buildings at risk), in the emergency phase. The results of a case study in the Puglia region in southern Italy are described to illustrate the practical applications of this newly proposed approach. The main advantage of the proposed approach is that it allows for defining a hierarchy between different infrastructure in the urban area through the identification of particular components whose operation and efficiency are critical for emergency management. This information can be used by decision-makers to prioritize risk reduction interventions in flood emergencies in urban areas, given limited financial resources.

Description: Modeling rapid mass movements using the shallow water equations Natural Hazards and Earth System Sciences Discussions, 2, 6775-6809, 2014 Author(s): S. Hergarten and J. Robl We propose a new method to model rapid mass movements on complex topography using the shallow water equations in Cartesian coordinates. These equations are the widely used standard approximation for the flow of water in rivers and shallow lakes, but the main prerequisite for their application – an almost horizontal fluid table – is in general not satisfied for avalanches and debris flows in steep terrain. Therefore, we have developed appropriate correction terms for large topographic gradients. In this study we present the mathematical formulation of these correction terms and their implementation in the open source flow solver GERRIS. This novel approach is evaluated by simulating avalanches on synthetic and finally natural topographies and the widely used Voellmy flow resistance law. The results are tested against analytical solutions and the commercial avalanche model RAMMS. The overall results are in excellent agreement with the reference system RAMMS, and the deviations between the different models are far below the uncertainties in the determination of the relevant fluid parameters and involved avalanche volumes in reality. As this code is freely available and open source, it can be easily extended by additional fluid models or source areas, making this model suitable for simulating several types of rapid mass movements. It therefore provides a valuable tool assisting regional scale natural hazard studies.

Description: Magnetotelluric investigation in the High Agri Valley (southern Apennine, Italy) Natural Hazards and Earth System Sciences Discussions, 2, 6747-6773, 2014 Author(s): M. Balasco, A. Giocoli, S. Piscitelli, G. Romano, A. Siniscalchi, T. A. Stabile, and S. Tripladi In this paper we present the result of a Magnetotelluric (MT) investigation carried out across the High Agri Valley (HAV), southern Italy. Several MT soundings were carried out in order to obtain a ~15 km long 2-D resistivity model with an investigation depth of ~10 km. The main aim was to provide valuable data on the geological and structural setting of the HAV. The MT model was compared with pre-existing geological, geophysical and seismic data. The MT model can be schematized as a superposition of three stack lateral varying layers with different thickness and resistivity values: a surficial low-medium resistivity layer, associated to the Quaternary deposits and to the allochthonous units, and a deeper high resistivity layer, related to the Apulia Platform, separated by a thin layer connected to the mélange zone and to the Pliocene terrigenous marine deposits. Sharp lateral resistivity variations are interpreted as faults that, on the basis of accurate focal mechanism computations, display normal-faulting kinematics.

Description: In this paper we explore the performances and the opportunities provided by the European satellite Sentinel-1 for water resource management applications in low-income countries. The analysis is supported by a synthetic aperture radar (SAR) simulator, which allowed the quantification of the expected characteristics of Sentinel-1 products in three applications: interferometric digital elevation models (DEMs) generation, land cover mapping and estimation of water volumes retained by small reservoirs. The obtained results quantitatively show that Sentinel-1 data characteristics are fully suitable for most of the application already explored in the recent SAR literature.

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

Description: Wind waves in the Black Sea: results of a hindcast study Natural Hazards and Earth System Science, 14, 2883-2897, 2014 Author(s): V. S. Arkhipkin, F. N. Gippius, K. P. Koltermann, and G. V. Surkova In this study we describe the wind wave fields in the Black Sea. The general aims of the work were the estimation of statistical wave parameters and the assessment of interannual and seasonal wave parameter variability. The domain of this study was the entire Black Sea. Wave parameters were calculated by means of the SWAN wave model on a 5 × 5 km rectangular grid. Initial conditions (wind speed and direction) for the period between 1949 and 2010 were derived from the NCEP/NCAR reanalysis. According to our calculations the average significant wave height on the Black Sea does not exceed 0.7 m. Areas of most significant heavy sea are the southwestern and the northeastern parts of the sea as expressed in the spatial distribution of significant wave heights, wave lengths and periods. Besides, long-term annual variations of wave parameters were estimated. Thus, linear trends of the annual total duration of storms and of their quantity are nearly stable over the hindcast period. However, an intensification of storm activity is observed in the 1960s–1970s.

Description: Estimation of three-dimensional crustal movements in the 2011 Tohoku-oki, Japan earthquake from TerraSAR-X intensity images Natural Hazards and Earth System Sciences Discussions, 2, 7383-7408, 2014 Author(s): W. Liu, F. Yamazaki, M. Matsuoka, T. Nonaka, and T. Sasagawa The Tohoku-Oki earthquake on 11 March 2011 caused significant widespread crustal movements. In a previous study, we proposed a method for capturing two-dimensional (2-D) surface displacements from a pair of pre- and post-event TerraSAR-X (TSX) intensity images. However, it is difficult to detect three-dimensional (3-D) displacements from one pair of TSX images. In this study, three pairs of pre- and post-event TSX images taken on different paths were used to estimate 3-D crustal movements. The relationship between the actual 3-D displacements and the converted 2-D movements in the SAR images was derived based on the observation model of a SAR sensor. The 3-D movements were then calculated from three sets of detected 2-D movements that occurred within a short time period. Compared with GPS observations, the proposed method was found to be capable of detecting the 3-D crustal movements with sub-pixel accuracy.

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

Description: Complex networks, streamflow, and hydrometric monitoring system design Hydrology and Earth System Sciences Discussions, 11, 13663-13710, 2014 Author(s): M. Halverson and S. Fleming Network theory is applied to an array of streamflow gauges located in the Coast Mountains of British Columbia and Yukon, Canada. The goal of the analysis is to assess whether insights from this branch of mathematical graph theory can be meaningfully applied to hydrometric data, and more specifically, whether it may help guide decisions concerning stream gauge placement so that the full complexity of the regional hydrology is efficiently captured. The streamflow data, when represented as a complex network, has a global clustering coefficient and average shortest path length consistent with small-world networks, which are a class of stable and efficient networks common in nature, but the results did not clearly suggest a scale-free network. Stability helps ensure that the network is robust to the loss of nodes; in the context of a streamflow network, stability is interpreted as insensitivity to station removal at random. Community structure is also evident in the streamflow network. A community detection algorithm identified 10 separate communities, each of which appears to be defined by the combination of its median seasonal flow regime (pluvial, nival, hybrid, or glacial, which in this region in turn mainly reflects basin elevation) and geographic proximity to other communities (reflecting shared or different daily meteorological forcing). Betweenness analyses additionally suggest a handful of key stations which serve as bridges between communities and might therefore be highly valued. We propose that an idealized sampling network should sample high-betweenness stations, as well as small-membership communities which are by definition rare or undersampled relative to other communities, while retaining some degree of redundancy to maintain network robustness.

Description: We processed X-band COSMO-SkyMed 3-m resolution StripMap HIMAGE time series (March 2011–June 2013) with the Stanford Method for Persistent Scatterers (StaMPS), to retrieve an updated picture of the condition and structural health of the historic centre of Rome, Italy, and neighbouring quarters. Taking advantage of an average target density of over 2800 PS/km2, we analysed the spatial distribution of more than 310,000 radar targets against: (1) land cover; (2) the location of archaeological ruins and restoration activities; and (3) the size, orientation and morphology of historical buildings. Radar interpretation was addressed from the perspective of conservators, and the deformation estimates were correlated to local geohazards and triggering factors of structural collapse. In the context of overall stability, deformation was identified at the single-monument scale, e.g., for the Roman cistern and exedra in the Oppian Hill. Comparative assessment against InSAR processing of C-band imagery (1992–2010) published in the literature confirms the persistence of ground motions affecting monuments and subsidence in southern residential quarters adjacent to the Tiber River, due to the consolidation of compressible deposits. Vertical velocity estimated from COSMO-SkyMed PS exceeds −7.0 mm/y in areas of recent urbanization.

Description: The purpose of this study was to examine how different polarimetric parameters and an object-based approach influence the classification results of various land use/land cover types using fully polarimetric ALOS PALSAR data over coastal wetlands in Yancheng, China. To verify the efficiency of the proposed method, five other classifications (the Wishart supervised classification, the proposed method without polarimetric parameters, the proposed method without an object-based analysis, the proposed method without textural and geometric information and the proposed method using the nearest-neighbor classifier) were applied for comparison. The results indicated that some polarimetric parameters, such as Shannon entropy, Krogager_Kd, Alpha, HAAlpha_T11, VanZyl3_Vol, Derd, Barnes2_T33, polarization fraction, Barnes1_T33, Neuman_delta_mod and entropy, greatly improved the classification results. The shape index was a useful feature in distinguishing fish ponds and rivers. The distance to the sea can be regarded as an important factor in reducing the confusion between herbaceous wetland vegetation and grasslands. Furthermore, the decision tree algorithm increased the overall accuracy by 6.8% compared with the nearest neighbor classifier. This research demonstrated that different polarimetric parameters and the object-based approach significantly improved the performance of land cover classification in coastal wetlands using ALOS PALSAR data.

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

Description: Amalgamation in landslide maps: effects and automatic detection Natural Hazards and Earth System Sciences Discussions, 2, 7651-7678, 2014 Author(s): O. Marc and N. Hovius Inventories of individually delineated landslides are a key to understanding landslide physics and mitigating their impact. They permit assessment of area-frequency distributions and landslide volumes, and testing of statistical correlations between landslides and physical parameters such as topographic gradient or seismic strong motion. Amalgamation, i.e. the mapping of several adjacent landslides as a single polygon, can lead to potentially severe distortion of the statistics of these inventories. This problem can be especially severe in datasets produced by automated mapping. We present 5 inventories of earthquake-induced landslides mapped with different materials and techniques and affected by varying degrees of amalgamation. Errors on the total landslide volume and power-law exponent of the area-frequency distribution, resulting from amalgamation, may be up to 200 and 50%, respectively. We present an algorithm based on image and DEM analysis, for automatic identification of amalgamated polygons. On a set of about 2000 polygons larger than 1000 m 2 , tracing landslides triggered by the 1994 Northridge earthquake, the algorithm performs well, with only 2.7–3.6% wrongly amalgamated landslides missed and 3.9–4.8% correct polygons wrongly identified as amalgams. This algorithm can be used broadly to check landslide inventories and allow faster correction by automating the identification of amalgamation.

Description: Due to an error on our part, the pagination and the doi number of this manuscript [1] were missing during updating. The correct pagination is 9600–9652, and the doi number is 10.3390/rs6109600. Here is the correct version:[...]

Description: Factors affecting flood insurance purchase in residential properties in Johor, Malaysia Natural Hazards and Earth System Science, 14, 3297-3310, 2014 Author(s): U. G. Aliagha, T. E. Jin, W. W. Choong, M. Nadzri Jaafar, and H. M. Ali High-impact floods have become a virtually annual experience in Malaysia, yet flood insurance has remained a grossly neglected part of comprehensive integrated flood risk management. Using discriminant analysis, this study seeks to identify the demand-side variables that best predict flood insurance purchase and risk aversion between two groups of residential homeowners in three districts of Johor State, Malaysia: those who purchased flood insurance and those who did not. Our results revealed an overall 34% purchase rate, with Kota Tinggi district having the highest (44%) and thus the highest degree of flood risk aversion. The Wilks' lambda F test for equality of group means, standardised discriminant function coefficients, structure correlation, and canonical correlation has clearly shown that there are strong significant attribute differences between the two groups of homeowners, based on the measures of objective flood risk exposure, subjective risk perception, and socio-economic cum demographic variables. However, the measures of subjective risk perception were found to be more predictive of flood insurance purchase and flood risk aversion.

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

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

Description: Monitoring and modelling of soil–plant interactions: the joint use of ERT, sap flow and Eddy Covariance data to characterize the volume of an orange tree root zone Hydrology and Earth System Sciences Discussions, 11, 13353-13384, 2014 Author(s): G. Cassiani, J. Boaga, D. Vanella, M. T. Perri, and S. Consoli Mass and energy exchanges between soil, plants and atmosphere control a number of key environmental processes involving hydrology, biota and climate. The understanding of these exchanges also play a critical role for practical purposes e.g. in precision agriculture. In this paper we present a methodology based on coupling innovative data collection and models in order to obtain quantitative estimates of the key parameters of such complex flow system. In particular we propose the use of hydro-geophysical monitoring via 4-D Electrical Resistivity Tomography (ERT) in conjunction with measurements of plant transpiration via sap flow and evapotranspiration from Eddy Covariance (EC). This abundance of data is fed to a spatially distributed soil model in order to characterize the distribution of active roots. We conducted experiments in an orange orchard in Eastern Sicily (Italy), characterized by the typical Mediterranean semi-arid climate. The subsoil dynamics, particularly influenced by irrigation and root uptake, were characterized mainly by the ERT setup, consisting of 48 buried electrodes on 4 instrumented micro boreholes (about 1.2 m deep) placed at the corners of a square (about 1.3 m in side) surrounding the orange tree, plus 24 mini-electrodes on the surface spaced 0.1 m on a square grid. During the monitoring, we collected repeated ERT and TDR soil moisture measurements, soil water samples, sap flow measurements from the orange tree and EC data. We conducted a laboratory calibration of the soil electrical properties as a function of moisture content and pore water electrical conductivity. Irrigation, precipitation, sap flow and ET data are available allowing knowledge of the system's long term forcing conditions on the system. This information was used to calibrate a 1-D Richards' equation model representing the dynamics of the volume monitored via 3-D ERT. Information on the soil hydraulic properties was collected from laboratory and field experiments. The successful results of the calibrated modeling exercise allow the quantification of the soil volume interested by root water uptake. This volume is much smaller (with a surface area less than 2 m 2 , and about 40 cm thickness) than expected and assumed in the design of classical drip irrigation schemes that prove to be losing at least half of the irrigated water that is not uptaken by the plants.

Description: Moraine-dammed lake failures in Patagonia and assessment of outburst susceptibility in the Baker Basin Natural Hazards and Earth System Science, 14, 3243-3259, 2014 Author(s): P. Iribarren Anacona, K.P. Norton, and A. Mackintosh Glacier retreat since the Little Ice Age has resulted in the development or expansion of hundreds of glacial lakes in Patagonia. Some of these lakes have produced large (≥ 10 6 m 3 ) Glacial Lake Outburst Floods (GLOFs) damaging inhabited areas. GLOF hazard studies in Patagonia have been mainly based on the analysis of short-term series (≤ 50 years) of flood data and until now no attempt has been made to identify the relative susceptibility of lakes to failure. Power schemes and associated infrastructure are planned for Patagonian basins that have historically been affected by GLOFs, and we now require a thorough understanding of the characteristics of dangerous lakes in order to assist with hazard assessment and planning. In this paper, the conditioning factors of 16 outbursts from moraine-dammed lakes in Patagonia were analysed. These data were used to develop a classification scheme designed to assess outburst susceptibility, based on image classification techniques, flow routine algorithms and the Analytical Hierarchy Process. This scheme was applied to the Baker Basin, Chile, where at least seven moraine-dammed lakes have failed in historic time. We identified 386 moraine-dammed lakes in the Baker Basin of which 28 were classified with high or very high outburst susceptibility. Commonly, lakes with high outburst susceptibility are in contact with glaciers and have moderate (〉 8°) to steep (〉 15°) dam outlet slopes, akin to failed lakes in Patagonia. The proposed classification scheme is suitable for first-order GLOF hazard assessments in this region. However, rapidly changing glaciers in Patagonia make detailed analysis and monitoring of hazardous lakes and glaciated areas upstream from inhabited areas or critical infrastructure necessary, in order to better prepare for hazards emerging from an evolving cryosphere.

Description: Seismology of the Oso-Steelhead landslide Natural Hazards and Earth System Sciences Discussions, 2, 7309-7327, 2014 Author(s): C. Hibert, C. P. Stark, and G. Ekström We carry out a combined analysis of the short- and long-period seismic signals generated by the devastating Oso-Steelhead landslide that occurred on 22 March 2014. The seismic records show that the Oso-Steelhead landslide was not a single slope failure, but a succession of multiple failures distinguished by two major collapses that occurred approximately three minutes apart. The first generated long-period surface waves that were recorded at several proximal stations. We invert these long-period signals for the forces acting at the source, and obtain estimates of the first failure runout and kinematics, as well as its mass after calibration against the mass-center displacement estimated from remote-sensing imagery. Short-period analysis of both events suggests that the source dynamics of the second are more complex than the first. No distinct long-period surface waves were recorded for the second failure, which prevents inversion for its source parameters. However, by comparing the seismic energy of the short-period waves generated by both events we are able to estimate the volume of the second. Our analysis suggests that the volume of the second failure is about 15–30% of the total landslide volume, which is in agreement with ground observations.

Description: In volcanic areas, where it can be difficult to perform direct surveys, digital photogrammetry techniques are rarely adopted for routine volcano monitoring. Nevertheless, they have remarkable potentialities for observing active volcanic features (e.g., fissures, lava flows) and the connected deformation processes. The ability to obtain accurate quantitative data of definite accuracy in short time spans makes digital photogrammetry a suitable method for controlling the evolution of rapidly changing large-area volcanic phenomena. The systematic acquisition of airborne photogrammetric datasets can be adopted for implementing a more effective procedure aimed at long-term volcano monitoring and hazard assessment. In addition, during the volcanic crisis, the frequent acquisition of oblique digital images from helicopter allows for quasi-real-time monitoring to support mitigation actions by civil protection. These images are commonly used to update existing maps through a photo-interpretation approach that provide data of unknown accuracy. This work presents a scientific tool (Orthoview) that implements a straightforward photogrammetric approach to generate digital orthophotos from single-view oblique images provided that at least four Ground Control Points (GCP) and current Digital Elevation Models (DEM) are available. The influence of the view geometry, of sparse and not-signalized GCP and DEM inaccuracies is analyzed for evaluating the performance of the developed tool in comparison with other remote sensing techniques. Results obtained with datasets from Etna and Stromboli volcanoes demonstrate that 2D features measured on the produced orthophotos can reach sub-meter-level accuracy.

Description: When object-based analysis is applied to very high-resolution imagery, pixels within the segments reveal large spectral inhomogeneity; their distribution can be considered complex rather than normal. When normality is violated, the classification methods that rely on the assumption of normally distributed data are not as successful or accurate. It is hard to detect normality violations in small samples. The segmentation process produces segments that vary highly in size; samples can be very big or very small. This paper investigates whether the complexity within the segment can be addressed using multiple random sampling of segment pixels and multiple calculations of similarity measures. In order to analyze the effect sampling has on classification results, statistics and probability value equations of non-parametric two-sample Kolmogorov-Smirnov test and parametric Student’s t-test are selected as similarity measures in the classification process. The performance of both classifiers was assessed on a WorldView-2 image for four land cover classes (roads, buildings, grass and trees) and compared to two commonly used object-based classifiers—k-Nearest Neighbor (k-NN) and Support Vector Machine (SVM). Both proposed classifiers showed a slight improvement in the overall classification accuracies and produced more accurate classification maps when compared to the ground truth image.

Description: Comparison of aeolian snow transport events and snow mass fluxes between observations and simulations made by the regional climate model MAR in Adélie Land, East Antarctica The Cryosphere Discussions, 8, 6007-6032, 2014 Author(s): A. Trouvilliez, H. Gallée, F. Naaim-Bouvet, C. Genthon, C. Amory, V. Favier, C. Agosta, L. Piard, and H. Bellot The regional climate model MAR including a coupled snow pack/aeolian snow transport parameterisation is compared with aeolian snow mass fluxes at a fine spatial resolution (5 km horizontally and 2 m vertically) and at a fine temporal resolution (30 min) over 1 month in Antarctica. Numerous feedbacks are taken into account in the MAR including the drag partitioning caused by the roughness elements. Wind speed is correctly simulated with a positive value of the Nash test (0.60 and 0.37) but the wind speeds above 10 m s −1 are underestimated. The aeolian snow transport events are correctly reproduced with a good temporal resolution except for the aeolian snow transport events with a particles' maximum height below 1 m. The simulated threshold friction velocity, calculated without snowfall, is overestimated. The simulated aeolian snow mass fluxes between 0 to 2 m have the same variations but are underestimated compared to the second-generation FlowCapt values and so is the simulated relative humidity at 2 m. This underestimation is not entirely due to the underestimation of the simulated wind speed. The MAR underestimates the aeolian snow quantity that pass through the first two meters by a factor ten compared to the second-generation FlowCapt value (13 990 kg m −1 and 151 509 kg m −1 respectively). It will conduct the MAR, with this parametrisation, to underestimate the effect of the aeolian snow transport on the Antarctic surface mass balance.

Description: Influence of weak layer heterogeneity and slab properties on slab tensile failure propensity and avalanche release area The Cryosphere Discussions, 8, 6033-6057, 2014 Author(s): J. Gaume, G. Chambon, N. Eckert, M. Naaim, and J. Schweizer Dry-snow slab avalanches are generally caused by a sequence of fracture processes including failure initiation in a weak snow layer underlying a cohesive slab followed by crack propagation within the weak layer (WL) and tensile fracture through the slab. During past decades, theoretical and experimental work has gradually improved our knowledge of the fracture process in snow. However, our limited understanding of crack propagation and fracture arrest propensity prevents the evaluation of avalanche release sizes and thus impedes hazard assessment. To address this issue, slab tensile failure propensity is examined using a mechanically-based statistical model of the slab–WL system based on the finite element method. This model accounts for WL heterogeneity, stress redistribution by elasticity of the slab and the slab possible tensile failure. Two types of avalanche release are distinguished in the simulations: (1) full-slope release if the heterogeneity is not sufficient to stop crack propagation and to trigger a tensile failure within the slab, (2) partial-slope release if fracture arrest and slab tensile failure occurs due to the WL heterogeneity. The probability of these two release types is presented as a function of the characteristics of WL heterogeneity and of the slab. One of the main outcomes is that, for realistic values of the parameters, the tensile failure propensity is mainly influenced by slab properties. Hard and thick snow slabs are more prone to wide-scale crack propagation and thus lead to larger avalanches (full-slope release). In this case, the avalanche size is mainly influenced by topographical and morphological features such as rocks, trees, slope curvature and the spatial variability of the snow depth as it is often claimed in the literature.

Description: Epic landslide erosion from mountain roads in Yunnan, China – challenges for sustainable development Natural Hazards and Earth System Science, 14, 3093-3104, 2014 Author(s): R. C. Sidle, M. Ghestem, and A. Stokes Expanding systems of mountain roads in developing countries have significantly increased the risk of landslides and sedimentation, and have created vulnerabilities for residents and aquatic resources. We measured landslide erosion along seven road segments in steep terrain in the upper Salween River basin, Yunnan, China and estimated sediment delivery to channels. Landslide erosion rates along the roads ranged from 2780 to 48 235 Mg ha −1 yr −1 , the upper end of this range being the highest rate ever reported along mountain roads. The two roads with the highest landslide erosion (FG1 = 12 966 Mg ha −1 yr −1 ; DXD = 48 235 Mg ha −1 yr −1 ) had some of the highest sediment delivery rates to channels (about 80 and 86%, respectively). Overall, 3 times more landslides occurred along cut slopes compared to fill slopes, but fill slope failures had a combined mass 〉 1.3 times that of cut slope failures. Many small landslides occurred along road cuts, but these were often trapped on the road surface. Given the magnitude of the landslide problem and the lack of attention to this issue, a more sustainable approach for mountain road development is outlined based on an analysis of landslide susceptibility and how thresholds for landslide trigger mechanisms would be modified by road location and different construction techniques.

Description: Evaluation of forest fire models on a large observation database Natural Hazards and Earth System Science, 14, 3077-3091, 2014 Author(s): J. B. Filippi, V. Mallet, and B. Nader This paper presents the evaluation of several fire propagation models using a large set of observed fires. The observation base is composed of 80 Mediterranean fire cases of different sizes, which come with the limited information available in an operational context (burned surface and approximative ignition point). Simulations for all cases are carried out with four different front velocity models. The results are compared with several error scoring methods applied to each of the 320 simulations. All tasks are performed in a fully automated manner, with simulations run as first guesses with no tuning for any of the models or cases. This approach leads to a wide range of simulation performance, including some of the bad simulation results to be expected in an operational context. Disregarding the quality of the input data, it is found that the models can be ranked based on their performance and that the most complex models outperform the more empirical ones. Data and source codes used for this paper are freely available to the community.

Description: Monitoring and prediction in Early Warning Systems (EWS) for rapid mass movements Natural Hazards and Earth System Sciences Discussions, 2, 7149-7179, 2014 Author(s): M. Stähli, M. Sättele, C. Huggel, B. W. McArdell, P. Lehmann, A. Van Herwijnen, A. Berne, M. Schleiss, A. Ferrari, A. Kos, D. Or, and S. M. Springman Rapid mass movements (RMM) pose a substantial risk to people and infrastructure. Reliable and cost-efficient measures have to be taken to reduce this risk. One of these measures includes establishing and advancing the State of Practice in the application of Early Warning Systems (EWS). EWS have been developed during the past decades and are rapidly increasing. In this document, we focus on the technical part of EWS, i.e. the prediction and timely recognition of imminent hazards, as well as on monitoring slopes at risk and released mass movements. Recent innovations in assessing spatial precipitation, as well as monitoring and modelling precursors, the triggering and deformation of RMM offer new opportunities for next-generation EWS. However, technical advancement can only be transferred into more reliable, operational EWS with an intense dialog between scientists, engineers and those in charge of warning. To this end, further experience with new comprehensive prototype systems jointly operated by scientists and practitioners will be essential.

Description: A new approach to mapping permafrost and change incorporating uncertainties in ground conditions and climate projections The Cryosphere, 8, 2177-2194, 2014 Author(s): Y. Zhang, I. Olthof, R. Fraser, and S. A. Wolfe Spatially detailed information on permafrost distribution and change with climate is important for land use planning, infrastructure development, and environmental assessments. However, the required soil and surficial geology maps in the North are coarse, and projected climate scenarios vary widely. Considering these uncertainties, we propose a new approach to mapping permafrost distribution and change by integrating remote sensing data, field measurements, and a process-based model. Land cover types from satellite imagery are used to capture the general land conditions and to improve the resolution of existing permafrost maps. For each land cover type, field observations are used to estimate the probabilities of different ground conditions. A process-based model is used to quantify the evolution of permafrost for each ground condition under three representative climate scenarios (low, medium, and high warming). From the model results, the probability of permafrost occurrence and the most likely permafrost conditions are determined. We apply this approach at 20 m resolution to a large area in Northwest Territories, Canada. Mapped permafrost conditions are in agreement with field observations and other studies. The data requirements, model robustness, and computation time are reasonable, and this approach may serve as a practical means to mapping permafrost and changes at high resolution in other regions.

Description: Results from a full coupling of the HIRHAM regional climate model and the MIKE SHE hydrological model for a Danish catchment Hydrology and Earth System Sciences, 18, 4733-4749, 2014 Author(s): M. A. D. Larsen, J. C. Refsgaard, M. Drews, M. B. Butts, K. H. Jensen, J. H. Christensen, and O. B. Christensen A major challenge in the emerging research field of coupling of existing regional climate models (RCMs) and hydrology/land-surface models is the computational interaction between the models. Here we present results from a full two-way coupling of the HIRHAM RCM over a 4000 km × 2800 km domain at 11 km resolution and the combined MIKE SHE-SWET hydrology and land-surface models over the 2500 km 2 Skjern River catchment. A total of 26 one-year runs were performed to assess the influence of the data transfer interval (DTI) between the two models and the internal HIRHAM model variability of 10 variables. DTI frequencies between 12 and 120 min were assessed, where the computational overhead was found to increase substantially with increasing exchange frequency. In terms of hourly and daily performance statistics the coupled model simulations performed less accurately than the uncoupled simulations, whereas for longer-term cumulative precipitation the opposite was found, especially for more frequent DTI rates. Four of six output variables from HIRHAM, precipitation, relative humidity, wind speed and air temperature, showed statistically significant improvements in root-mean-square error (RMSE) by reducing the DTI. For these four variables, the HIRHAM RMSE variability corresponded to approximately half of the influence from the DTI frequency and the variability resulted in a large spread in simulated precipitation. Conversely, DTI was found to have only a limited impact on the energy fluxes and discharge simulated by MIKE SHE.

Description: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice The Cryosphere, 8, 2219-2233, 2014 Author(s): S. Arndt and M. Nicolaus Arctic sea ice has not only decreased in volume during the last decades, but has also changed in its physical properties towards a thinner and more seasonal ice cover. These changes strongly impact the energy budget, and might affect the ice-associated ecosystems. In this study, we quantify solar shortwave fluxes through sea ice for the entire Arctic during all seasons. To focus on sea-ice-related processes, we exclude fluxes through open water, scaling linearly with sea ice concentration. We present a new parameterization of light transmittance through sea ice for all seasons as a function of variable sea ice properties. The maximum monthly mean solar heat flux under the ice of 30 × 10 5 Jm −2 occurs in June, enough heat to melt 0.3 m of sea ice. Furthermore, our results suggest that 96% of the annual solar heat input through sea ice occurs during only a 4-month period from May to August. Applying the new parameterization to remote sensing and reanalysis data from 1979 to 2011, we find an increase in transmitted light of 1.5% yr −1 for all regions. This corresponds to an increase in potential sea ice bottom melt of 63% over the 33-year study period. Sensitivity studies reveal that the results depend strongly on the timing of melt onset and the correct classification of ice types. Assuming 2 weeks earlier melt onset, the annual transmitted solar radiation to the upper ocean increases by 20%. Continuing the observed transition from a mixed multi-year/first-year sea ice cover to a seasonal ice cover results in an increase in light transmittance by an additional 18%.

Description: Estimating forest canopy height from large-footprint satellite LiDAR waveforms is challenging given the complex interaction between LiDAR waveforms, terrain, and vegetation, especially in dense tropical and equatorial forests. In this study, canopy height in French Guiana was estimated using multiple linear regression models and the Random Forest technique (RF). This analysis was either based on LiDAR waveform metrics extracted from the GLAS (Geoscience Laser Altimeter System) spaceborne LiDAR data and terrain information derived from the SRTM (Shuttle Radar Topography Mission) DEM (Digital Elevation Model) or on Principal Component Analysis (PCA) of GLAS waveforms. Results show that the best statistical model for estimating forest height based on waveform metrics and digital elevation data is a linear regression of waveform extent, trailing edge extent, and terrain index (RMSE of 3.7 m). For the PCA based models, better canopy height estimation results were observed using a regression model that incorporated both the first 13 principal components (PCs) and the waveform extent (RMSE = 3.8 m). Random Forest regressions revealed that the best configuration for canopy height estimation used all the following metrics: waveform extent, leading edge, trailing edge, and terrain index (RMSE = 3.4 m). Waveform extent was the variable that best explained canopy height, with an importance factor almost three times higher than those for the other three metrics (leading edge, trailing edge, and terrain index). Furthermore, the Random Forest regression incorporating the first 13 PCs and the waveform extent had a slightly-improved canopy height estimation in comparison to the linear model, with an RMSE of 3.6 m. In conclusion, multiple linear regressions and RF regressions provided canopy height estimations with similar precision using either LiDAR metrics or PCs. However, a regression model (linear regression or RF) based on the PCA of waveform samples with waveform extent information is an interesting alternative for canopy height estimation as it does not require several metrics that are difficult to derive from GLAS waveforms in dense forests, such as those in French Guiana.

Description: Quality segment generation is a well-known challenge and research objective within Geographic Object-based Image Analysis (GEOBIA). Although methodological avenues within GEOBIA are diverse, segmentation commonly plays a central role in most approaches, influencing and being influenced by surrounding processes. A general approach using supervised quality measures, specifically user provided reference segments, suggest casting the parameters of a given segmentation algorithm as a multidimensional search problem. In such a sample supervised segment generation approach, spatial metrics observing the user provided reference segments may drive the search process. The search is commonly performed by metaheuristics. A novel sample supervised segment generation approach is presented in this work, where the spectral content of provided reference segments is queried. A one-class classification process using spectral information from inside the provided reference segments is used to generate a probability image, which in turn is employed to direct a hybridization of the original input imagery. Segmentation is performed on such a hybrid image. These processes are adjustable, interdependent and form a part of the search problem. Results are presented detailing the performances of four method variants compared to the generic sample supervised segment generation approach, under various conditions in terms of resultant segment quality, required computing time and search process characteristics. Multiple metrics, metaheuristics and segmentation algorithms are tested with this approach. Using the spectral data contained within user provided reference segments to tailor the output generally improves the results in the investigated problem contexts, but at the expense of additional required computing time.

Description: Using HyMap data, mineral identification and mineral mapping were conducted on the basis of the spectral absorption index (SAI) and other spectral absorption features in a study area in Tudun, eastern Tien Shan. Alteration minerals, such as calcite, alumina-rich (Al-rich) muscovite, epidote, and antigorite, were explored, and their relative abundance was depicted. A cross-validation was performed, and it showed a high degree of consistency between the imagery results and the results of previous literature. To further validate the mineral mapping from HyMap data, a field survey was carried out and rock samples were collected for quantitative analysis using a Por Infrared Mineral Analyzer (PIMA) and the software affiliated with it. Minerals were discriminated, and their relative abundance was calculated from the spectra. Although we found that the absorption band-depth and SAI agreed well with each other and with the relative abundance of mineral alterations, the spectral absorption band-depth provided a better representation. Finally, ore prospecting of the study area was presented, and we found the distribution and close spatial relationships among the minerals extracted using the HyMap data. In the northern and northwestern part of the Gold-mine area, there was a mineralized muscovite alteration showing a sheet or block distribution. In the Copper-mine area, Al-poor muscovite with a sheet distribution was distributed in the north and northeast region, and Al-rich muscovite showed a block distribution enclosed by the distribution area of Al-poor muscovite. These all showed good ore prospects for the study area.

Description: Testing gridded land precipitation data and precipitation and runoff reanalyses (1982–2010) between 45° S and 45° N with Normalized Difference Vegetation Index data Hydrology and Earth System Sciences Discussions, 11, 13175-13205, 2014 Author(s): S. O. Los The realistic simulation of key components of the land-surface hydrological cycle – precipitation, runoff, evaporation and transpiration – in general circulation models of the atmosphere is crucial to assess adverse weather impacts on environment and society. Here, gridded precipitation data from observations and precipitation and runoff fields from reanalyses were tested with satellite-derived global vegetation index data for 1982–2010 and latitudes between 45° S and 45° N. Data were obtained from the Climate Research Unit (CRU), the Global Precipitation Climatology Project (GPCP) and Tropical Rainfall Monitoring Mission (TRMM; analysed for 1998–2010 only) and (precipitation and runoff) reanalyses were obtained from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), the European Centre for Medium-Range Weather Forecasts (ECMWF) and the NASA Global Modelling and Assimilation Office (GMAO). Annual land-surface precipitation was converted to annual potential vegetation net primary productivity (NPP) and was compared to mean annual Normalized Difference Vegetation Index data measured by the Advanced Very High Resolution Radiometer (1982–1999) and MODIS (2001–2010). The effect of spatial resolution on the agreement between NPP and NDVI was investigated as well. The CRU and TRMM derived NPP agreed most closely with the NDVI data. The GPCP data showed weaker spatial agreement, largely because of their lower spatial resolution, but similar temporal agreement. MERRA Land and ERA Interim precipitation reanalyses showed similar spatial agreement as the GPCP data and good temporal agreement in semi-arid regions of the Americas, Asia, Australia and southern Africa. The NCEP/NCAR reanalysis showed the lowest spatial agreement which could only in part be explained by its lower spatial resolution. No reanalysis showed realistic interannual precipitation variations for northern tropical Africa. Inclusion of runoff in the NPP prediction resulted only in (marginally) better agreement for the MERRA Land reanalysis and worse agreement for the NCEP/NCAR and ERA Interim reanalyses.

Description: A global dataset of the extent of irrigated land from 1900 to 2005 Hydrology and Earth System Sciences Discussions, 11, 13207-13258, 2014 Author(s): S. Siebert, M. Kummu, M. Porkka, P. Döll, N. Ramankutty, and B. R. Scanlon Irrigation intensifies land use by increasing crop yield but also impacts water resources. It affects water and energy balances and consequently the microclimate in irrigated regions. Therefore, knowledge of the extent of irrigated land is important for hydrological and crop modelling, global change research, and assessments of resource use and management. Information on the historical evolution of irrigated lands is limited. The new global Historical Irrigation Dataset (HID) provides estimates of the temporal development of the area equipped for irrigation (AEI) between 1900 and 2005 at 5 arc-minute resolution. We collected subnational irrigation statistics from various sources and found that the global extent of AEI increased from 63 million ha (Mha) in 1900 to 112 Mha in 1950 and 306 Mha in 2005. We developed eight gridded versions of time series of AEI by combining subnational irrigation statistics with different data sets on the historical extent of cropland and pasture. Different rules were applied to maximize consistency of the gridded products to subnational irrigation statistics or to historical cropland and pasture data sets. The HID reflects very well the spatial patterns of irrigated land in the western United States as shown on historical maps. Mean aridity on irrigated land increased and river discharge decreased from 1900–1950 whereas aridity decreased from 1950–2005. The dataset and its documentation are made available in an open data repository at https://mygeohub.org/publications/8 (doi: 10.13019/M2MW2G ).

Description: Debris flows in the Eastern Italian Alps: seasonality and atmospheric circulation patterns Natural Hazards and Earth System Sciences Discussions, 2, 7197-7224, 2014 Author(s): E. I. Nikolopoulos, M. Borga, F. Marra, S. Crema, and L. Marchi The work examines the seasonality and large-scale atmospheric circulation patterns of debris flows in the Trentino-Alto Adige region (Eastern Italian Alps). Analysis is based on classification algorithms applied on a uniquely dense archive of debris flows and hourly rain gauge precipitation series covering the period 2000–2009. Results highlight the seasonal and synoptic forcing patterns linked to debris flows in the study area. Summer and fall season account for 92% of the debris flows in the record, while atmospheric circulation characterized by Zonal West, Mixed and Meridional South, Southeast patterns account for 80%. Both seasonal and circulation patterns exhibit geographical preference. In the case of seasonality, there is a strong north–south separation of summer–fall dominance while spatial distribution of dominant circulation patterns exhibits clustering, with both Zonal West and Mixed prevailing in the northwest and central east part of the region, while the southern part relates to Meridional South, Southeast pattern. Seasonal and synoptic pattern dependence is pronounced also on the debris flow triggering rainfall properties. Examination of rainfall intensity–duration thresholds derived for different data classes (according to season and synoptic pattern) revealed a distinct variability in estimated thresholds. These findings imply a certain control on debris-flow events and can therefore be used to improve existing alert systems.

Description: Climatic signals from 76 shallow firn cores in Dronning Maud Land, East Antarctica The Cryosphere Discussions, 8, 5961-6005, 2014 Author(s): S. Altnau, E. Schlosser, E. Isaksson, and D. Divine The spatial and temporal distribution of surface mass balance (SMB) and δ 18 O were investigated in the first comprehensive study of a set of 76 firn cores retrieved by various expeditions during the past three decades in Dronning Maud Land, East Antarctica. The large number of cores was used to calculate stacked records of SMB and δ 18 O, which considerably increased the signal-to-noise ratio compared to earlier studies and facilitated the detection of climatic signals. Considerable differences between cores from the interior plateau and the coastal cores were found. The δ 18 O of both the plateau and the ice shelf cores exhibit a slight positive trend over the second half of the 20th century. In the corresponding period, the SMB has a negative trend in the ice shelf cores, but increases on the plateau. Comparison with meteorological data from Neumayer Station revealed that for the ice shelf regions atmospheric dynamic effects are more important than thermodynamics, while on the plateau, the temporal variations of SMB and δ 18 O occur mostly in parallel, thus can be explained by thermodynamic effects. The Southern Annular Mode (SAM) exhibits a positive trend since the mid-1960s, which is assumed to lead to a cooling of East Antarctica. This is not confirmed by the firn core data in our data set. Changes in the atmospheric circulation that result in a changed seasonal distribution of precipitation/accumulation could partly explain the observed features in the ice shelf cores.

Description: Post-LIA glacier changes along a latitudinal transect in the Central Italian Alps The Cryosphere, 8, 2235-2252, 2014 Author(s): R. Scotti, F. Brardinoni, and G. B. Crosta The variability of glacier response to atmospheric temperature rise in different topo-climatic settings is still a matter of debate. To address this question in the Central Italian Alps, we compile a post-LIA (Little Ice Age) multitemporal glacier inventory (1860–1954–1990–2003–2007) along a latitudinal transect that originates north of the continental divide in the Livigno Mountains and extends south through the Disgrazia and Orobie ranges, encompassing continental-to-maritime climatic settings. In these sub-regions, we examine the area change of 111 glaciers. Overall, the total glacierized area has declined from 34.1 to 10.1 km 2 , with a substantial increase in the number of small glaciers due to fragmentation. The average annual decrease (AAD) in glacier area has risen by about 1 order of magnitude from 1860–1990 (Livigno: 0.45; Orobie: 0.42; and Disgrazia: 0.39 % a −1 ) to 1990–2007 (Livigno: 3.08; Orobie: 2.44; and Disgrazia: 2.27 % a −1 ). This ranking changes when considering glaciers smaller than 0.5 km 2 only (i.e., we remove the confounding caused by large glaciers in Disgrazia), so that post-1990 AAD follows the latitudinal gradient and Orobie glaciers stand out (Livigno: 4.07; Disgrazia: 3.57; and Orobie: 2.47 % a −1 ). More recent (2007–2013) field-based mass balances in three selected small glaciers confirm post-1990 trends showing the consistently highest retreat in continental Livigno and minimal area loss in maritime Orobie, with Disgrazia displaying transitional behavior. We argue that the recent resilience of glaciers in Orobie is a consequence of their decoupling from synoptic atmospheric temperature trends, a decoupling that arises from the combination of local topographic configuration (i.e., deep, north-facing cirques) and high winter precipitation, which ensures high snow-avalanche supply, as well as high summer shading and sheltering. Our hypothesis is further supported by the lack of correlations between glacier change and glacier attributes in Orobie, as well as by the higher variability in ELA,sub〉0 positioning, post-LIA glacier change, and interannual mass balances, as we move southward along the transect.

Description: High-resolution land surface modeling utilizing remote sensing parameters and the Noah UCM: a case study in the Los Angeles Basin Hydrology and Earth System Sciences, 18, 4791-4806, 2014 Author(s): P. Vahmani and T. S. Hogue In the current work we investigate the utility of remote-sensing-based surface parameters in the Noah UCM (urban canopy model) over a highly developed urban area. Landsat and fused Landsat–MODIS data are utilized to generate high-resolution (30 m) monthly spatial maps of green vegetation fraction (GVF), impervious surface area (ISA), albedo, leaf area index (LAI), and emissivity in the Los Angeles metropolitan area. The gridded remotely sensed parameter data sets are directly substituted for the land-use/lookup-table-based values in the Noah-UCM modeling framework. Model performance in reproducing ET (evapotranspiration) and LST (land surface temperature) fields is evaluated utilizing Landsat-based LST and ET estimates from CIMIS (California Irrigation Management Information System) stations as well as in situ measurements. Our assessment shows that the large deviations between the spatial distributions and seasonal fluctuations of the default and measured parameter sets lead to significant errors in the model predictions of monthly ET fields (RMSE = 22.06 mm month −1 ). Results indicate that implemented satellite-derived parameter maps, particularly GVF, enhance the capability of the Noah UCM to reproduce observed ET patterns over vegetated areas in the urban domains (RMSE = 11.77 mm month −1 ). GVF plays the most significant role in reproducing the observed ET fields, likely due to the interaction with other parameters in the model. Our analysis also shows that remotely sensed GVF and ISA improve the model's capability to predict the LST differences between fully vegetated pixels and highly developed areas.

Description: Given the advances in remotely sensed imagery and associated technologies, several global land cover maps have been produced in recent times including IGBP DISCover, UMD Land Cover, Global Land Cover 2000 and GlobCover 2009. However, the utility of these maps for specific applications has often been hampered due to considerable amounts of uncertainties and inconsistencies. A thorough review of these global land cover projects including evaluating the sources of error and uncertainty is prudent and enlightening. Therefore, this paper describes our work in which we compared, summarized and conducted an uncertainty analysis of the four global land cover mapping projects using an error budget approach. The results showed that the classification scheme and the validation methodology had the highest error contribution and implementation priority. A comparison of the classification schemes showed that there are many inconsistencies between the definitions of the map classes. This is especially true for the mixed type classes for which thresholds vary for the attributes/discriminators used in the classification process. Examination of these four global mapping projects provided quite a few important lessons for the future global mapping projects including the need for clear and uniform definitions of the classification scheme and an efficient, practical, and valid design of the accuracy assessment.

Description: Comprehensive evaluation of water resources security in the Yellow River basin based on a Fuzzy Multi-Attribute Decision Analysis Approach Hydrology and Earth System Sciences Discussions, 11, 371-410, 2014 Author(s): K. K. Liu, C. H. Li, Y. P. Cai, M. Xu, and X. H. Xia In this paper, a Fuzzy Multi-Attribute Decision Analysis Approach (FMADAA) was adopted in water resources security evaluation for the nine provinces in the Yellow River basin in 2006. A numerical approximation system and a modified left-right scoring approach were adopted to cope with the uncertainties in the acquired information. Four multi-attribute decision making methods were implemented in the evaluation model for impact evaluation, including simple weighted addition (SWA), weighted product (WP), cooperative game theory (CGT) and technique for order preference by similarity to ideal solution (TOPSIS) which could be used for helping rank the water resources security in those nine provinces as well as the criteria alternatives. Moreover, several aggregation methods including average ranking procedure, borda and copeland methods were used to integrate the ranking results. The ranking results showed that the water resources security of the entire basin is in critical, insecurity and absolute insecurity state, especially in Shanxi, Inner Mongolia and Ningxia provinces in which water resources were lower than the average quantity in China. Hence, future planning of the Yellow River basin should mainly focus on the improvement of water eco-environment status in the provinces above.

Description: Impact of the Hoa Binh Dam (Vietnam) on water and sediment budgets in the Red River basin and delta Hydrology and Earth System Sciences Discussions, 11, 333-370, 2014 Author(s): D. V. Vu, S. Ouillon, D. T. Tran, and V. C. La The Hoa Binh Dam, located on a tributary of the Red River in Vietnam, has a capacity of 9.45 × 10 9 m 3 and was commissioned in December 1988. Although being important for flood prevention, electricity production, and irrigation in northern Vietnam, the Hoa Binh Dam has also highly influenced the suspended sediment distribution in the lower Red River basin, in the delta and in the coastal zone. Its impact was analysed from 50 yr dataset of water discharge and suspended sediment concentration (1960–2010) and the distribution of water and sediment across the nine mouths of the delta was calculated using the MIKE 11 numerical model before and after the dam settlement. Although water discharge at the delta inlet decreased by only 8.8%, the yearly suspended sediment flux dropped, on average from 119 to 43 × 10 6 t yr −1 at Son Tay near Hanoi, and from 85 to 35 × 10 6 t yr −1 in the river mouths. Water regulation has led to decreased water discharge in the wet season and increased water discharge in the dry season. Suspended sediment discharge proportionally increased in northern and southern estuaries and decreased through the main and central Ba Lat mouth. Tidal pumping, which causes a net sediment flux from the coast to the estuary at low discharge, is high in the northern delta, as a consequence of the high tidal range (up to 4.5 m in spring tide; diurnal tide). The shifts in the dynamic and characteristics of the turbidity maximum zone in the Cam-Bach Dang estuary are probably the cause of the enhanced sediment deposition in the Haiphong harbor. Along the coast, the reduced sedimentation rates are coincident with the lower sediment delivery that has been observed since the impoundment of the Hoa Binh Dam.

Description: Theory of the generalized chloride mass balance method for recharge estimation in groundwater basins characterised by point and diffuse recharge Hydrology and Earth System Sciences Discussions, 11, 307-332, 2014 Author(s): N. Somaratne and K. R. J. Smettem Application of the conventional chloride mass balance (CMB) method to point recharge dominant groundwater basins can substantially under-estimate long-term average annual recharge by not accounting for the effects of localized surface water inputs. This is because the conventional CMB method ignores the duality of infiltration and recharge found in karstic systems, where point recharge can be a contributing factor. When point recharge is present in groundwater basins, recharge estimation is unsuccessful using the conventional CMB method with, either unsaturated zone chloride or groundwater chloride. In this paper we describe a generalized CMB that can be applied to groundwater basins with point recharge. Results from this generalized CMB are shown to be comparable with long-term recharge estimates obtained using the watertable fluctuation method, groundwater flow modelling and Darcy flow calculations. The generalized CMB method provides an alternative, reliable long-term recharge estimation method for groundwater basins characterised by both point and diffuse recharge.

Description: Detecting and monitoring forest degradation in the tropics has implications for various fields of interest (biodiversity, emission calculations, self-sustenance of indigenous communities, timber exploitation). However, remote-sensing-based detection of forest degradation is difficult, as these subtle degradation signals are not easy to detect in the first place and quickly lost over time due to fast re-vegetation. To overcome these shortcomings, a time series analysis has been developed to map and monitor forest degradation over a longer period of time, with frequent updates based on Landsat data. This time series approach helps to reduce both the commission and the omission errors compared to, e.g., bi- or tri-temporal assessments. The approach involves a series of pre-processing steps, such as geometric and radiometric adjustments, followed by spectral mixture analysis and classification of spectral curves. The resulting pixel-based classification is then aggregated to degradation areas. The method was developed on a study site in Cameroon and applied to a second site in Central African Republic. For both areas, the results were finally evaluated against visual interpretation of very high-resolution optical imagery. Results show overall accuracies in both study sites above 85% for mapping degradation areas with the presented methods.

Description: Radar sensor networks, including bi- and multi-static radars, provide several operational advantages, like reduced vulnerability, good system flexibility and an increased radar cross-section. However, radar-to-radar interference suppression is a major problem in distributed radar sensor networks. In this paper, we present a cross-matched filtering-based radar-to-radar interference suppression algorithm. This algorithm first uses an iterative filtering algorithm to suppress the radar-to-radar interferences and, then, separately matched filtering for each radar. Besides the detailed algorithm derivation, extensive numerical simulation examples are performed with the down-chirp and up-chirp waveforms, partially overlapped or inverse chirp rate linearly frequency modulation (LFM) waveforms and orthogonal frequency division multiplexing (ODFM) chirp diverse waveforms. The effectiveness of the algorithm is verified by the simulation results.

Description: Brief Communication: Light-absorbing impurities can reduce the density of melting snow The Cryosphere Discussions, 8, 259-271, 2014 Author(s): O. Meinander, A. Kontu, A. Virkkula, A. Arola, L. Backman, P. Dagsson-Waldhauserová, O. Järvinen, T. Manninen, J. Svensson, G. de Leeuw, and M. Leppäranta Climatic effects of Black Carbon (BC) deposition on snow have been proposed to result from reduced snow albedo and increased melt due to light-absorbing particles. In this study, we hypothesize that BC may decrease the liquid water retention capacity of melting snow, and present our first data, where both the snow density and elemental carbon content were measured. In our experiments, artificially added light-absorbing impurities decreased the density of seasonally melting natural snow. We also suggest three possible processes that might lead to the lower snow density.

Description: An algorithm based on a Bayesian network classifier was adapted to produce 10-day burned area (BA) maps from the Long Term Data Record Version 3 (LTDR) at a spatial resolution of 0.05° (~5 km) for the North American boreal region from 2001 to 2011. The modified algorithm used the Brightness Temperature channel from the Moderate Resolution Imaging Spectroradiometer (MODIS) band 31 T31 (11.03 μm) instead of the Advanced Very High Resolution Radiometer (AVHRR) band T3 (3.75 μm). The accuracy of the BA-LTDR, the Collection 5.1 MODIS Burned Area (MCD45A1), the MODIS Collection 5.1 Direct Broadcast Monthly Burned Area (MCD64A1) and the Burned Area GEOLAND-2 (BA GEOLAND-2) products was assessed using reference data from the Alaska Fire Service (AFS) and the Canadian Forest Service National Fire Database (CFSNFD). The linear regression analysis of the burned area percentages of the MCD64A1 product using 40 km × 40 km grids versus the reference data for the years from 2001 to 2011 showed an agreement of R2 = 0.84 and a slope = 0.76, while the BA-LTDR showed an agreement of R2 = 0.75 and a slope = 0.69. These results represent an improvement over the MCD45A1 product, which showed an agreement of R2 = 0.67 and a slope = 0.42. The MCD64A1, BA-LTDR and MCD45A1 products underestimated the total burned area in the study region, whereas the BA GEOLAND-2 product overestimated it by approximately five-fold, with an agreement of R2 = 0.05. Despite MCD64A1 showing the best overall results, the BA-LTDR product proved to be an alternative for mapping burned areas in the North American boreal forest region compared with the other global BA products, even those with higher spatial/spectral resolution.

Description: This paper presents a new spaceborne synthetic aperture radar (SAR) burst mode named “Extended Terrain Observation by Progressive Scans (ETOPS)” for wide swath imaged coverage. This scheme extends the imaging performance of the conventional Terrain Observation by Progressive Scans (TOPS) mode with a very limited azimuth beam steering capability. Compared with the TOPS mode with the same azimuth beam steering range for the same swath width, a finer azimuth resolution could be obtained. With the same system parameters, examples of four burst SAR imaging modes named ScanSAR, TOPS, inverse TOPS (ITOPS) and ETOPS are given, and their corresponding system performances are analyzed and compared. Simulation results show that the proposed ETOPS mode could obtain a better high-resolution wide-swath imaging performance under the same conditions.

Description: Hard paths, soft paths or no paths? Cross-cultural perceptions of water solutions Hydrology and Earth System Sciences, 18, 109-120, 2014 Author(s): A. Wutich, A. C. White, D. D. White, K. L. Larson, A. Brewis, and C. Roberts In this study, we examine how development status and water scarcity shape people's perceptions of "hard path" and "soft path" water solutions. Based on ethnographic research conducted in four semi-rural/peri-urban sites (in Bolivia, Fiji, New Zealand, and the US), we use content analysis to conduct statistical and thematic comparisons of interview data. Our results indicate clear differences associated with development status and, to a lesser extent, water scarcity. People in the two less developed sites were more likely to suggest hard path solutions, less likely to suggest soft path solutions, and more likely to see no path to solutions than people in the more developed sites. Thematically, people in the two less developed sites envisioned solutions that involve small-scale water infrastructure and decentralized, community-based solutions, while people in the more developed sites envisioned solutions that involve large-scale infrastructure and centralized, regulatory water solutions. People in the two water-scarce sites were less likely to suggest soft path solutions and more likely to see no path to solutions (but no more likely to suggest hard path solutions) than people in the water-rich sites. Thematically, people in the two water-rich sites seemed to perceive a wider array of unrealized potential soft path solutions than those in the water-scarce sites. On balance, our findings are encouraging in that they indicate that people are receptive to soft path solutions in a range of sites, even those with limited financial or water resources. Our research points to the need for more studies that investigate the social feasibility of soft path water solutions, particularly in sites with significant financial and natural resource constraints.

Description: Goal-oriented networks and capacity building for natural hazards – examples in the Dresden region Natural Hazards and Earth System Science, 14, 81-94, 2014 Author(s): G. Hutter Networks and networking are important for building social capacities for natural hazards. However, up to now, it has been an open question which types of networks contribute to capacity building under certain circumstances. The paper focuses on the type of a goal-oriented network. The distinction between goal orientation and goal directedness is used to show the following: goal directedness of networks to build capacities for natural hazards involves intensive and continuous processes of "sensemaking" (Weick, 1995) to specify the network goal. This process of specifying an initial goal statement is important in small and large networks at the regional level. The governance form of a lead organization network facilitates goal specification. The paper illustrates these findings through evidence from two case studies conducted in the Dresden region in Germany.

Description: Mountain chains of Central Asia host a large number of glaciated areas that provide critical water supplies to the semi-arid populated foothills and lowlands of this region. Spatio-temporal variations of glacier flows are a key indicator of the impact of climate change on water resources as the glaciers react sensitively to climate. Satellite remote sensing using optical imagery is an efficient method for studying ice-velocity fields on mountain glaciers. In this study, temporal and spatial changes in surface velocity associated with the Inylchek glacier in Kyrgyzstan are investigated. We present a detailed map for the kinematics of the Inylchek glacier obtained by cross-correlation analysis of Landsat images, acquired between 2000 and 2011, and a set of ASTER images covering the time period between 2001 and 2007. Our results indicate a high-velocity region in the elevated part of the glacier, moving up to a rate of about 0.5 m/day. Time series analysis of optical data reveals some annual variations in the mean surface velocity of the Inylchek during 2000–2011. In particular, our findings suggest an opposite trend between periods of the northward glacial flow in Proletarskyi and Zvezdochka glacier, and the rate of westward motion observed for the main stream of the Inylchek.

Description: A decade (2002–2012) of supraglacial lake volume estimates across Russell Glacier, West Greenland The Cryosphere, 8, 107-121, 2014 Author(s): A. A. W. Fitzpatrick, A. L. Hubbard, J. E. Box, D. J. Quincey, D. van As, A. P. B. Mikkelsen, S. H. Doyle, C. F. Dow, B. Hasholt, and G. A. Jones Supraglacial lakes represent an ephemeral storage buffer for meltwater runoff and lead to significant, yet short-lived, episodes of ice-flow acceleration by decanting large meltwater and energy fluxes into the ice sheet's hydrological system. Here, a methodology for calculating lake volume is used to quantify storage and drainage across Russell Glacier, West Greenland, between 2002 and 2012. Using 502 MODIS scenes, water volume at ~200 seasonally occurring lakes was derived using a depth–reflectance relationship, which was independently calibrated and field validated against lake bathymetry. The inland expansion of lakes is strongly correlated with air temperature: during the record melt years of 2010 and 2012, lakes formed and drained earlier, attaining their maximum volume 38 and 20 days earlier than the 11 yr mean, as well as occupying a greater area and forming at higher elevations (〉 1800 m) than previously. Despite occupying under 2% of the study area, lakes delay the transmission of up to 7–13% of the bulk meltwater discharged. Although the results are subject to an observational bias caused by periods of cloud cover, we estimate that across Russell Glacier, 28% of supraglacial lakes drain rapidly ( 〈 4 days). Clustering of such events in space and time suggests a synoptic trigger mechanism. Further, we find no evidence to support a unifying critical size or depth-dependent drainage threshold.

Description: Fabric measurement along the NEEM ice core, Greenland, and comparison with GRIP and NGRIP ice cores The Cryosphere Discussions, 8, 307-335, 2014 Author(s): M. Montagnat, N. Azuma, D. Dahl-Jensen, J. Eichler, S. Fujita, F. Gillet-Chaulet, S. Kipfstuhl, D. Samyn, A. Svensson, and I. Weikusat Fabric (distribution of crystallographic orientations) profile along the full NEEM ice core, Greenland, is presented in this work. Data were measured in the field by an Automatic Ice Texture Analyzer every 10 m, from 33 m down to 2461 m depth. The fabric evolves from a slightly anisotropic fabric at the top, toward a strong single maximum at about 2300 m, which is typical of a deformation pattern mostly driven by uniaxial compression and simple shearing. A sharp increase in the fabric strengthening is observed at the Holocene to Wisconsin climatic transition. A similar strengthening, toward an anisotropic single maximum-type fabric, has been observed in several ice cores from Greenland and Antarctica, and can be attributed to a positive feedback between changes in ice viscosity at the climatic transition, and the impact of a shear component of stress. Centimeter scale abrupt texture (fabric and microstructure) variations are observed in the bottom part of the core. Their positions are in good agreement with the folding hypothesis used for a climatic reconstruction by Dahl-Jensen et al. (2013). Comparison is made to two others ice cores drilled along the same ridge; the GRIP ice core drilled at the summit of the ice sheet, and the NorthGRIP ice core, drilled 325 km to the NNW of the summit along the ridge, and 365 km upstream from NEEM. The fabric profile clearly reflects the increase in shear deformation when moving NW along the ridge from GRIP to NorthGRIP and NEEM. The difference in fabric profiles between NEEM and NorthGRIP also evidences a stronger lateral extension associated with a sharper ridge at NorthGRIP.