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: Governments, tech giants and individuals must act together if we are to clamp down on threats and abuse in the digital world

Description: From body invasion Alien -style to hallucinogenic poisons and limpet mines, fungi have astonishingly creative ways of dispatching their foes

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: Infinity is larger than large. It's only when you try to work out how much larger that you realise some infinities are larger than others

Description: Probability is one of those things we all get wrong… deeply wrong. The important thing is not to use your intuition

Description: Gravity, electromagnetism, the Higgs field – even empty space isn't empty space because of these invisible webs that permeate reality

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: The 2004 tsunami killed 227,000 people in Indonesia. Fred Pearce visits Aceh to see the living coastal defences meant to soak up the energy of future waves

Description: Adding laser tips to ordinary shoes can improve the stride and pace of people with Parkinson's disease

Description: We can now grow plants in microgravity – and crops grown in asteroid soil could sustain vast human populations off-planet

Description: Fingerprint magnets that allow you to feel electric fields, RFID chips that let you remotely control keys, tuning hearing to Wi-Fi fields: body hacking is here

Description: Got a busy calendar? Let AI assistants set your meetings for you and suggest friends for you to interact with

Description: Rumour has it that a consumer version of the virtual reality headset Oculus Rift will emerge in 2015 – and we can't wait to get one

Description: World leaders are expected to agree a global plan to halt climate change in Paris, but we're already living through a renewable energy revolution

Description: Declining snowfall in winter will leave Norwegian spruce trees at the mercy of sub-zero temperatures and insect pests

Description: Our current best snapshot of Pluto is just a handful of pixels captured by the Hubble Space Telescope – that should change when New Horizons flies by

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: There may be an economic argument for rethinking measures designed to protect health and the environment – if so, we should be free to discuss it

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: Mathematics is more fun than most of us reckon: think World of Warcraft with a pencil, says the comedian who wrote Things to Make and Do in the Fourth Dimension

Description: Australian storm-chasing photographer Nick Moir threw the kids in the car and tore off to catch this shot of squall line hitting Sydney

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: Sustained nuclear fusion on Earth always seems to be 30 years away – because our reactors have been too ambitious, say researchers designing on a small scale

Description: The root cause of the accident is still unclear, but officials say video from the cockpit shows a mechanism designed to slow the craft deployed earlier than normal

Description: Conservation of water will have to become a way of life in California as projections suggest long-term drought will become the norm

Description: The disease was thought to be rare in humans, but last year 68 per cent of people hospitalised with malaria in Malaysian Borneo were infected with it

Description: An app created by The Samaritans that flags tweets that suggest people are struggling to cope has caused a furore as it collects users data without consent

Description: Medical implants such as pacemakers can be tampered with to kill – without leaving a trace. But the race is on to find ways to spot the crime

Description: An app watches when, where and how you use your smartphone, then locks up if things suddenly change

Description: Stem cells in the breast milk of mice develop into a range of different tissues in the offspring, including the brain, liver and kidneys. The same may also be true in people

Description: Singapore's water grid monitors for leaks every millisecond. The technology is being installed in cities worldwide to conserve scarce supplies and avert floods

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: Boss of UK spy agency GCHQ censures US tech firms for facilitating terrorists with their social networks – but others say encryption is good for everyone

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: Some parents in China are sending kids hooked on the internet to rehabilitation centres where they are denied online access

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: Biohacker Andrew Kessel 's open-source drug company aims to make bespoke cancer-fighting viruses using DNA printers – paid for by a Netflix-style subscription

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: Beyoncé's tweets as sacred text, error messages that tell a tragic tale and the looping program of a failing relationship: we parse a new literary phenomenon

Description: A 23-year-old man is the first person to report a feeling of persistent déjà vu that seems to be triggered by anxiety rather than a neurological condition

Description: Just in time for Christmas, instructions for a Lego watt balance, which can fix the official value of the kilogram, have been published

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: It's easy to make a storm glass – just mix the chemicals and watch crystals form and disappear. But there might there be a scientific basis to its claims

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.