ALBERT

Description: We examined the relationships between large wood (LW) export and precipitation patterns and intensity by analyzing the data on the annual volume of LW removed from 42 reservoirs and the daily precipitation at or near the reservoir sites. We also calculated the effective precipitation by considering the antecedent precipitation. Both daily and effective precipitation data were used as explanatory variables to explain LW export. The model selection revealed that the precipitation pattern and intensity controlling LW export varied with latitude in the Japanese archipelago. In small watersheds with narrow channel widths and low discharges, mass movements, such as landslides and debris flows, are major factors in the production and transport of LW. In this case, the effective precipitation required to initiate mass movements regulated the LW export and did not vary with the latitude. In intermediate and large watersheds with wide channel widths and high stream discharges, heavy rainfall and subsequent floods regulated buoyant depth, influencing the initiation of LW movement. In southern and central Japan, intense rainfall accompanied by typhoons or localized torrential downpours causes geomorphic disturbances, which introduce abundant pieces of LW into the channels. However, these pieces continue to be removed by repeated rainfall events. Therefore, LW export is supply-limited and potentially produces less LW accumulation. Conversely, in northern Japan, where typhoons and torrential downpours are rare, LW export is transport-limited because LW pieces recruited by bank erosion, tree mortality, and windthrow accumulate and persist on valley floors. These pieces may be easily exported by infrequent flooding.

Description: We developed a method to measure in situ the isotopic composition of liquid water with minimal supervision and, most important, with a temporal resolution of less than a minute. For this purpose a microporous hydrophobic membrane contactor (Membrana) was combined with an isotope laser spectrometer (Picarro). The contactor, originally designed for degassing liquids, was used with N2 as a carrier gas in order to transform a small fraction of liquid water to water vapor. The generated water vapor was then analyzed continuously by the Picarro analyzer. To prove the membrane's applicability, we determined the specific isotope fractionation factor for the phase change through the contactor's membrane across an extended temperature range (8°C–21°C) and with different waters of known isotopic compositions. This fractionation factor is needed to subsequently derive the liquid water isotope ratio from the measured water vapor isotope ratios. The system was tested with a soil column experiment, where the isotope values derived with the new method corresponded well (R2 = 0.998 for δ18O and R2 = 0.997 for δ2H) with those of liquid water samples taken simultaneously and analyzed with a conventional method (cavity ring-down spectroscopy). The new method supersedes taking liquid samples and employs only relatively cheap and readily available components. This makes it a relatively inexpensive, fast, user-friendly, and easily reproducible method. It can be applied in both the field and laboratory wherever a water vapor isotope analyzer can be run and whenever real-time isotope data of liquid water are required at high temporal resolution.

Description: Transit time of discharge is a hydrological characteristic used in water resource management. Previous studies have demonstrated large spatial variation in the mean transit time (MTT) of stream base flow in meso-scale catchments. Various relationships between topography and MTT have been reported. Although it is generally assumed that base flow MTT is controlled by the depth of the hydrologically active layer that recharges a stream, this hypothesis has not been tested in field studies. This study confirmed that the depth of hydrologically active soil and bedrock controls spatial variation in MTT. The study used isotopic and geochemical tracer data gathered in the 4.27 km2 Fudoji catchment, central Japan. The results, together with previously documented relationships between topography and MTT, indicate that the depth of the hydrologically active layer is sometimes, but not always, related to topography. A comprehensive understanding of the factors that control base flow production in mountainous catchments will require further study of the water flow path depths that recharge streams.

Description: This paper examines the long-term historical changes in frequency and amplitude of hydroclimatic extremes in the Blue Nile basin using data from the second half of 20th century. The temporal variability of basin-wide rainfall extremes and river flow extremes from four gauging stations was investigated under the hypothesis of no trend and no persistence in time. On the basis of a quantile anomaly analysis method, decadal variations in extreme daily, monthly, and annual quantiles were studied, and the periods of statistical significance were identified. The analysis showed that high and low river flows and rainfall depths do not vary in time in a fully random way but show a particular variation pattern. Their extremes show significant decadal variations. The 1980s had statistically significant negative anomalies in extremes in comparison with the long-term reference period of 1964–2009, while the 1960s–1970s and the 1990s–2000s had positive anomalies, although less significant. There is neither consistent increasing nor decreasing trend in rainfall and flow extremes of recent years. Therefore, anticipated trends due to global warming could not be identified. Conversely, low-flow extremes show an increasing trend during the last decade, which could be related to the effect of water regulation works at the outlet of Lake Tana. Moreover, similar patterns and statistically significant correlations were found between climatic indices representing the Pacific and Atlantic Oceans and the Blue Nile rainfall and flow extremes. Changes that occur on the Pacific Ocean appear to be a main driver for the decadal oscillations in climate and related high and low Blue Nile water availability for Ethiopia, Sudan, and Egypt.

Description: Few studies exist on infiltration processes in badlands, although infiltration and subsurface lateral flows are known to contribute to soil erosion and to control slope instability. Our investigation was carried out in a 100 m² plot located in a 0.5 ha landslide in black marls (South-East France). An artificial sprinkling was performed with an intensity of 10 mm.h -1 during 66.4 h interrupted with 8.4 h. breaks. KBr and KCl were used as tracers. A pseudo-steady state was reached after 25-35 hours and 250-350 mm of rainfall. The runoff coefficient was 40% (ratio total runoff volume/total sprinkling water amount). Pre-event water (PE) contributed to the groundwater recharge at the very beginning of the experiment but PE contribution dropped steadily while the soil was saturating. After around 200 mm cumulative rainfall, PE contribution started to rise steeply before reaching a nearly constant value. This original mechanism implies an efficient transfer process of PE. It was assumed from the description of the material structure and from hydrological evidences that PE was mainly drained from a structure porosity made of the marl's flaked nature. Total pre-event water contributions ranged from 25 to 79 % (PE contribution was over 50 % in 2/3 of the observations wells). Over the recession phase, release of pre-event water occurred from the drainage of a texture porosity. The study showed that at the plot scale, infiltration processes proved well organised despite the high heterogeneity and anisotropy of the material. It was possible to propose a general conceptual model explaining the hydrological processes over time and area. The peculiar structure of regolith originating from black marl is preserved over a large part of the weathering time, so that the material structure (type, orientation of grains, small/large pores) remains a first order control of water flow generation in Black marl soils. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In this work, the multifractal properties of hourly rainfall data recorded at a location in Southern Spain have been related to the scale properties of the corresponding Intensity-Duration-Frequency (IDF) curves. Four parametric models for the IDF curves have been fitted to the quantiles of rainfall obtained by using the Generalized Pareto frequency distribution function with the extreme data series obtained for the same place. The scaling of the rainfall intensity moments has been analysed and the empirical moments scaling exponent function has been obtained. The corresponding values of q 1 and γ 1 have been empirical and theoretically calculated and compared to some characteristics of the different IDF models. Thus, the scaling behaviour of IDF curves has been analysed and the best model has been selected. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Major-ion compositions of groundwater are employed in this study of the water-rock interactions and hydrogeochemical evolution within a carbonate aquifer system. The groundwater samples were collected from boreholes or underground tunnels in the Ordovician limestone of Yanzhou Coalfield where catastrophic groundwater inflows can be hazardous to mining and impact use of the groundwater as a water supply. The concentration of TDS ranged from 961 to 3,555 mg/l and indicates moderately to highly mineralized water. The main water-type of the middle Ordovician limestone groundwater is Ca-Mg-SO 4 , with SO 4 2- ranging from 537 mg/l to 2,297 mg/l, and average values of Ca 2+ and Mg 2+ of 455.7 mg/l and 116.6 mg/l, respectively. The water samples were supersaturated with respect to calcite and dolomite and undersaturated or saturated with respect to gypsum. Along the general flow direction, deduced from increases of TDS and Cl - , the main water-rock interactions that caused hydrogeochemical evolution of the groundwater within the aquifer were the dissolution of gypsum, the precipitation of calcite, the dissolution or precipitation of dolomite, and ion exchange. Ion exchange is the major cause for the lower mole concentration of Ca 2+ than that of SO 4 2- . The groundwater level of Ordovician aquifer is much higher than that of C-P coal-bearing aquifers, so the potential flow direction is upward and the pyrite in coal is not a possible source of sulfate, additional data on the stable sulfur and oxygen isotopic composition of the sulfate may be helpful to identify its origin. Although ion exchange probably accounts for the higher mole concentration of Na + than that of Cl - , the dissolution of aluminosilicate can not be ruled out. The data evaluation methods and results of this study could be useful in other areas to understand flow paths in aquifers and to provide information needed to identify the origin of groundwater. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Quantitative evaluation of the effect of climate variability and human activities on runoff is of great importance for water resources planning and management in term of maintaining the ecosystem integrity and sustaining the society development. In this paper, hydro-climatic data from 4 catchments (i.e., Luanhe River Catchment, Chaohe River Catchment, Hutuo River Catchment and Zhanghe River Catchment) in the Haihe River Basin from 1957–2000 were used to quantitatively attribute the hydrological response (i.e., runoff) to climate change and human activities separately. To separate the attributes, the temporal trends of annual precipitation, potential evapotranspiration (PET) and runoff during 1957–2000 were first explored by the Man-Kendall test. Despite that only Hutuo River Catchment was dominated by significant negative trend in annual precipitation, all four catchments presented significant negative trend in annual runoff varying from −0.859 mm a -1 (Chaohe River) to −1.996 mm a -1 (Zhanghe River). Change points in 1977 and 1979 are detected by precipitation-runoff double cumulative curves method and Pettitt's test for Zhanghe River and other three rivers, respectively, and are adopted to divide data set into two study periods as the pre-change period and post-change period. Three methods including hydrological model method, hydrological sensitivity analysis method and climate elasticity method were calibrated with the hydro-climatic data during the pre-change period. And then, hydrological runoff response to climate variability and human activities were quantitatively evaluated with the help of the three methods and based on the assumption that climate and human activities are the only drivers for streamflow and are independent to each other. Similar estimates of anthropogenic and climatic effects on runoff for catchments considered can be obtained from the three methods. We found that human activities were the main driving factors for the decline in annual runoff in Luanhe River Catchment, Chaohe River Catchment and Zhanghe River Catchment, accounting for over 50% of runoff reduction. However, Climate variability should be responsible for the decrease in annual runoff in the Hutuo River Catchment. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Measurements of river water quality at Haridwar, India, taken during 2005 – 2006 show that the river water does not meet the WHO and Bureau of Indian Standards criteria of drinking water quality, especially with respect to total coliform and fecal-coliform. This study investigates the removal of pathogens at a river bank filtration (RBF) site in Haridwar. Using the quality of river water and the quality of abstracted water from a nearby production well, semi-empirical models based on the concept of filtration coefficient are developed and tested for their effectiveness in removing pathogens under varying bacteriological quality of source water. A two-tier model which includes the effect of clogged layer is developed to obtain an equivalent filtration coefficient. This coefficient is found to be linearly related with natural logarithm of the concentration of pathogens in the source water. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint-source (NPS) pollution by acting as hydrologic links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field- and watershed-scale hydrology, and pollutant transport in a 38 km 2 agricultural watershed in south-central NY. Salient findings include: (i) 94% of road ditches discharged directly to natural streams, effectively doubling the drainage density, (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches, (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.

Description: We coupled the process-based NIES Integrated Catchment-based Eco-hydrology (NICE) model to an urban canopy model (UCM) and the Regional Atmospheric Modeling System (RAMS) in order to simulate the effect of urban geometry and anthropogenic exhaustion on the hydrothermal changes in the atmospheric/land and the interfacial areas of the Japanese megalopolis. The simulation was conducted with multi-scale in horizontally regional–urban–point levels, and in vertically atmosphere–surface–unsaturated–saturated layers. The model reproduced reasonably the observed hydrothermal values by using ground-truth data in various types of natural/artificial land covers. The simulated results also suggested that the latent heat flux in new water-holding pavement (consisting of porous asphalt and water-holding filler made of steel by-products based on silica compound) has a strong impact on hydrologic cycle and cooling temperature in comparison with the observed heat budget by newly incorporating the effect of water amount on the heat conductivity in the pavement. Furthermore, the model predicted the hydrothermal changes under two types of land cover scenarios to promote evaporation and to reduce air temperature against heat island phenomenon. Finally, we evaluated the relationship between the effect of groundwater use to ameliorate the heat island and the effect of infiltration on the water cycle in the catchment. These procedures to integrate the multi-scaled model simulation with political scenario based on the effective management of water resources as heat sink/source would be very powerful approaches to recovering a sound hydrologic cycle and create thermally-pleasing environments in the megalopolis. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Soil and vadose zone profiles are used as an archive of changes in groundwater recharge and water quality following changes in land use in an area of the Loess Plateau of China. A typical rain-fed loess-terrace agriculture region in Hequan, Guyuan is taken as an example and multiple tracers (chloride mass balance, stable isotopes, tritium, and water chemistry) are used to examine groundwater recharge mechanisms and to evaluate soil water chloride as an archive for recharge rate and water quality. Results show that groundwater recharge beneath natural uncultivated grassland, used as a baseline, is about 94–100 mm yr -1 and the time it takes for annual precipitation to reach water table through the thick unsaturated zone is from decades to hundreds of years (tritium free). This recharge rate is 2–3 orders of magnitude more than in the other semiarid areas with similar annual rainfall and with deep rooted vegetation and relatively high temperature. Most of the water that eventually becomes recharge originally infiltrated in the summer months. The conversion from native grassland to winter wheat has reduced groundwater recharge by 42–50% (50–55 mm yr -1 for recharge), and from winter wheat to alfalfa resulted in a significant chloride accumulation in the upper soil zone which terminated deep drainage. The paper also evaluates the time lag between potential recharge and actual recharge to aquifer and between increase in solute concentration in soil moisture and that in the aquifer following land-use change due to the deep unsaturated zone. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Shallow upland drains, grips, have been hypothesized as responsible for increased downstream flow magnitudes. Observations provide counterfactual evidence, often relating to the difficulty of inferring conclusions from statistical correlation and paired catchment comparisons; and the complexity of designing field experiments to test grip impacts at the catchment-scale. Drainage should provide drier antecedent moisture conditions, providing more storage at the start of an event; but, grips have higher flow velocities than overland flow so potentially delivering flow more rapidly to the drainage network. We develop and apply a model for assessing the impacts of grips upon flow hydrographs. The model was calibrated on the gripped case; then the gripped case was compared with the intact case by removing all grips. This comparison showed that even given parameter uncertainty, the intact case had significantly higher flood peaks and lower baseflows, mirroring field observations of the hydrological response of intact peat. The simulations suggest that this is because delivery effects may not to translate into catchment-scale impacts for three reasons. First, in our case, the proportions of flow path lengths that were hillslope were not changed significantly by gripping. Second, the structure of the grip network as compared with the structure of the drainage basin mitigated against grip-related increases in the concentration of runoff in the drainage network, although it did marginally reduce the mean timing of that concentration at the catchment outlet. Third, the effect of the latter upon downstream flow magnitudes can only be assessed by reference to the peak timing of other tributary basins, emphasizing that drain effects are both relative and scale dependent. However, given the importance of hillslope flow paths, we show that if upland drainage causes significant changes in surface roughness on hillslopes, then critical and important feedbacks may impact upon the speed of hydrological response. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Natural soil pipes are found in peatlands but little is known about their hydrological role. This paper presents the most complete set of pipe discharge data to date from a deep blanket peatland in northern England. In a 17.4-ha catchment, we identified 24 perennially-flowing and 60 ephemerally-flowing pipe outlets. Eight pipe outlets along with the catchment outlet were continuously gauged over an 18-month period. The pipes in the catchment were estimated to produce around 13.7 % of annual streamflow with individual pipes often producing large peak flows (maximum peak of 3.8 L s -1 ). Almost all pipes, whether ephemeral, perennially-flowing, shallow or deep (outlets 〉 1 m below the peat surface), showed increased discharge within a mean of 3 hours after rainfall commencement and were dominated by stormflow, indicating good connectivity between the peatland surface and the pipes. However, almost all pipes had a longer time period between hydrograph peak and return to baseflow compared to the stream (mean of 23.9 hours for pipes, 19.7 hours for stream). As a result, the proportion of streamflow produced by the pipes at any given time increased at low flows and formed the most important component of stream discharge for the lowest 10 % of flows. Thus, a small number of perennially-flowing pipes became more important to the stream system under low flow conditions and probably received water via matrix flow during periods between storms. Given the importance of pipes to streamflow in blanket peatlands, further research is required into their wider role in influencing stream-water chemistry, water temperature and fluvial carbon fluxes, as well as their role in altering local hydrochemical cycling within the peat mass itself. Enhanced piping within peatlands caused by environmental change may lead to changes in streamflow regime with larger low flows and more prolonged drainage of the peat. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The hydrology of Quebec (Canada) boreal fens is poorly documented. Many peatlands are located in watersheds with impounded rivers. In such cases, their presence influences reservoir inflows. In recent years, some fens have been subjected to an increase of their wet area, a sign that they may be evolving towards an aquatic ecosystem. This dynamic process is referred to as aqualysis. This paper presents the seasonal and monthly hydrological budgets of a small watershed including a highly aqualysed fen (James Bay region). Monitoring of precipitation ( P ), runoff ( Q ) and groundwater levels ( WL ) was conducted during the ice-free season. Three semi-empirical equations (Thornthwaite, Priestley-Taylor and Penman-Monteith) were used and compared to calculate potential evapotranspiration ( PET ). The first two equations, having fewer parameters, estimate higher PET values than the third equation. The use of pressure level gauges installed in wells, for the calculation of peatland water storage, is inconclusive. Swelling of peat, peat decomposition and plant composition could be responsible for non-negligible amounts of absorbed water, which are not entirely accounted for by well levels. The estimation of peat matrix water storage is potentially the largest source of error and the limiting factor to calculate water balances in this environment. The results show that the groundwater level and the water storage vary depending on the season and especially after a heavy rainfall. Finally, the results illustrate the complexity of water routing through the site and, thus, raise several questions to be resolved in the future. Copyright © 2012 John Wiley & Sons, Ltd.

Description: A set of three quad-pol images acquired at the L-band in interferometric repeat-pass mode by the German Aerospace Center (DLR) with the Experimental SAR (E-SAR) system, in parallel with the AgriSAR2006 campaign, has been used to provide, for the first time with airborne data, a demonstration of the retrieval of vegetation height from agricultural crops by means of polarimetric SAR interferometry (PolInSAR)-based techniques. Despite the low frequency of the data, hence providing a weak response from the vegetation volume in contrast to the ground, accurate estimates of vegetation height at field level have been obtained over winter rape and maize fields. The same procedure does not yield valid estimates for wheat, barley, and sugar beet fields due to a mismatch with the physical model employed in the inversion and to the specific crop condition at the date of acquisition. These results show the value of the information provided by both interferometry and polarimetry for some agriculture monitoring practices.

Description: The recently launched Soil Moisture and Ocean Salinity (SMOS) satellite is providing soil moisture observations at continental scales by measuring L-band microwave radiation emitted from the land surface. While its retrieval algorithms will correct for factors such as vegetation and surface roughness, it will not correct for soil salinity. This letter tests the assumption that soil salinity will have a negligible impact on L-band brightness temperature $(T_{b})$ at SMOS scales using field data; airborne $T_{b}$ observations were collected in a saline groundwater discharge area near Nilpinna Station, South Australia. At the 500-m scale, the airborne observations of $T_{b}$ could not be reproduced using the baseline algorithm of the SMOS Level 2 retrieval scheme, without accounting for soil salinity in the model. The analysis in this letter shows that soil moisture retrieval errors of at least 0.04 $hbox{m}^{3} hbox{m}^{-3}$ (i.e., the entire SMOS error budget) will occur due to salinity alone in SMOS footprints with saline coverage as low as 25% (possibly even much less). Consequently, fractional salinity coverage cannot be considered a negligible factor by microwave soil moisture satellite missions.

Description: Based on the difference in statistical distribution between the target and jamming signal in the synthetic aperture radar (SAR) image, this letter proposes a novel adaptive method for barrage jamming suppression. In this method, the covariance matrix of jamming is estimated from the SAR image. Then, the 2-D sinc function of the ideal point target is utilized as the steering vector to obtain the optimal adaptive filter. This filter can suppress the random barrage jamming effectively, thus improving the image contrast and interpretability. Additionally, this letter analyzes in detail the theoretical basis and performance of the proposed method. Finally, simulations are provided to demonstrate its effectiveness.

Description: The advent of submeter-resolution synthetic aperture radar (SAR) images from satellites such as TerraSAR-X has given a new dimension to SAR image understanding. Even though emphasis is always on discovering automatic means of target characterization, visual exploration of targets and objects is the first step in many applications. While considering the complex-valued SAR images, visual inspection of the targets in an image may provide incomplete and misleading information, as sometimes two entirely different behaving objects look quite similar in SAR images. Thus, a need was felt to develop a methodology to support visual target recognition and analysis. In this letter, we present a method which looks into the complex-valued spectrum of SAR images, thus allowing a detailed physical interpretation of the scattering behavior of objects. The presented method is a joint time–frequency analysis method based on sublook decomposition. With the presented results, we emphasize the use of complex-valued SAR images for target characterization, the use of which is primarily restricted to polarimetric and interferometric applications as of now.

Description: This letter presents a novel neighborhood-based ratio (NR) operator to produce a difference image for change detection in synthetic aperture radar (SAR) images. In order to reduce the negative influence of speckle noise on SAR images, the proposed NR operator produces a difference image by combining gray level information and spatial information of neighbor pixels. The performance comparisons of the proposed operator with a traditional ratio operator and a log-ratio operator indicate that the NR operator is superior to these traditional methods and produces better detection results.

Description: In this letter, the requirements to derive topography from a portable terrestrial radar interferometer are introduced, the instrument design and the relationship between interferometric phase shift and surface topography are explained, and two examples of topographic maps from measurements at the Rhone glacier and Grabengufer rock glacier in Switzerland are presented. In the first case, an external digital elevation model was used to assess the error of topography mapping with the portable radar interferometer and to analyze ice surface changes of the glacier in the last 14 years. We found that the height error standard deviation is about 3 m within a distance of 2 km from the sensor and observed massive thinning of the Rhone glacier. In the second case, we used the terrestrial radar interferometer in order to measure the height difference between August 2009 and March 2010 over the rock glacier as a consequence of its destabilization.

Description: We map the coseismic deformation of the 2011 Tohoku-Oki earthquake with data from three descending Envisat/ASAR tracks and six ascending ALOS/PALSAR tracks that cover most of northeastern Japan. Due to the inaccurate estimation of the satellite status, orbital ramps commonly exist in the coseismic interferograms, which resulted in inconsistency among the deformation maps released by several research groups. In this letter, calibration has been performed to accurately remove these ramps by a 2-D quadratic-phase model derived based on GPS measurements from the ARIA team at the Jet Propulsion Laboratory and Caltech. The average RMS of the interferometric synthetic aperture radar (InSAR) measurements, as compared with GPS measurements at the validation stations, has decreased from 17.8 to 7.7 cm after the orbital ramp correction is made, indicating that much more accurate InSAR measurements are achieved. The corrected coseismic deformation from the InSAR measurements is consistent with not only the GPS observations at the individual GPS stations but also with the coseismic deformation interferogram from interpolated GPS observation in the SAR viewing directions. The corrected coseismic deformation measurement results show a maximum line-of-sight displacement of up to 3.7 m from the ascending PALSAR tracks and 2.4 m from the descending ASAR tracks, respectively.

Description: Global Navigation Satellite Systems (GNSS) opportunity signals reflected at or near the Earth's surface have already shown their potential to perform retrievals of a number of geophysical parameters. Radio occultations using GNSS signals are also used for atmospheric sensing. This letter presents a GNSS technique to retrieve vegetation water content (VWC). This technique measures the received powers of the GPS signals in open sky and under the vegetation layer. From these two powers, the attenuation due to the vegetation is computed, which is related to the VWC. This letter presents the results obtained after deploying the instrument in a walnut-tree stand for 11 months.

Description: In this letter, an improved differential evolution (DE) for high-dimensional waveform inversion is proposed. In conventional evolutionary algorithms, an individual is treated as a whole, and all its variables (genes) are evaluated with a uniform fitness function. This evaluation criterion is not effective for a high-dimensional individual. Therefore, for high-dimensional waveform inversion, we incorporate the decomposition strategy of cooperative coevolution into DE to decompose the individual into some subcomponents. Another novel feature that we introduce is a local fitness function for each subcomponent, and a new mutation operator is designed to guide the mutation direction of each subcomponent with the corresponding local fitness value. Coevolution among different subcomponents is realized in the selection operation with the global fitness function. Many experiments have been carried out to evaluate the performance of this new algorithm. The results clearly show that, for high-dimensional waveform inversion, this algorithm is effective and performs better than some other methods. Finally, the new method has been applied to real seismic data.

Description: This letter proposes the kernel entropy component analysis for clustering remote sensing data. The method generates nonlinear features that reveal structure related to the Rényi entropy of the input space data set. Unlike other kernel feature-extraction methods, the top eigenvalues and eigenvectors of the kernel matrix are not necessarily chosen. Data are interestingly mapped with a distinct angular structure, which is exploited to derive a new angle-based spectral clustering algorithm based on the mapped data. An out-of-sample extension of the method is also presented to deal with test data. We focus on cloud screening from Medium Resolution Imaging Spectrometer images. Several images are considered to account for the high variability of the problem. Good results obtained show the suitability of the proposal.

Description: A major source of error for repeat-pass interferometric synthetic aperture radar is the phase delay in radio signal propagation through the atmosphere, particularly the part due to tropospheric water vapor. These effects become more significant for ScanSAR observations due to their wider coverage (e.g., 400 km $times$ 400 km for ENVISAT Advanced Synthetic Aperture Radar (ASAR) wide swath (WS) mode versus 100 km $times$ 100 km for ASAR image mode). In this letter, we demonstrate for the first time that a Medium Resolution Imaging Spectrometer water vapor correction model can significantly reduce atmospheric water vapor effects on ASAR WS interferograms, with the phase variation in non-deforming areas decreasing from 3.8 cm before correction to 0.4 cm after correction.

Description: Potential climate change effects on aspects of conjunctive management of water resources can be evaluated by linking climate models with fully integrated groundwater–surface water models. The objective of this study is to develop a modeling system that links global climate models with regional hydrologic models, using the California Central Valley as a case study. The new method is a supply and demand modeling framework that can be used to simulate and analyze potential climate change and conjunctive use. Supply-constrained and demand-driven linkages in the water system in the Central Valley are represented with the linked climate models, precipitation-runoff models, agricultural and native vegetation water use, and hydrologic flow models to demonstrate the feasibility of this method. Simulated precipitation and temperature were used from the GFDL-A2 climate change scenario through the 21st century to drive a regional water balance mountain hydrologic watershed model (MHWM) for the surrounding watersheds in combination with a regional integrated hydrologic model of the Central Valley (CVHM). Application of this method demonstrates the potential transition from predominantly surface water to groundwater supply for agriculture with secondary effects that may limit this transition of conjunctive use. The particular scenario considered includes intermittent climatic droughts in the first half of the 21st century followed by severe persistent droughts in the second half of the 21st century. These climatic droughts do not yield a valley-wide operational drought but do cause reduced surface water deliveries and increased groundwater abstractions that may cause additional land subsidence, reduced water for riparian habitat, or changes in flows at the Sacramento–San Joaquin River Delta. The method developed here can be used to explore conjunctive use adaptation options and hydrologic risk assessments in regional hydrologic systems throughout the world.

Description: An efficient approach is developed to analytically evaluate solute transport in a horizontal, divergent radial flow field with a multistep injection flow rate and an arbitrary input concentration history. By assuming a piecewise steady state flow and transforming the time domain to the cumulative injected flow domain, the concentration distribution is found to be completely determined by the total volume of injected flow and independent of specific flow rates. Thus, on the cumulative flow domain, the transport problem with a temporally varying velocity field can be transformed into a steady state flow problem. Linear convolution can then be applied on the cumulative injected flow domain to evaluate the solution for an arbitrarily time-dependent input concentration. Solutions on the regular time domain can be conveniently obtained by mapping the solution on the cumulative injected flow domain to the time domain. Furthermore, we theoretically examine the conditions for the assumption of piecewise steady state flow to be valid. On the basis of the critical time scale of the “pseudosteady state condition,” defined as when velocity changes accomplish 99% of their steady state differences, and the relative error in the mean travel time of plume front, we obtain conditions for neglecting the transitional period between two pumping steps. Such conditions include the following: (1) the duration of a pumping step, tp, must be longer than the critical time scale, tc, i.e., tp ≥ tc = 25r2S/T, where r is the radial distance, S is the storage coefficient, and T is the transmissivity, or similarly, a maximum problem domain needs to be defined for a given pumping strategy. (2) the maximum well pumping rate, qmax, should satisfy qmax ≤ πθT/25S, where θ is the effective porosity. When both conditions are satisfied, transitional periods may be neglected.

Description: Meander bends in alluvial rivers morphologically evolve toward meander cutoff with narrowing intra-meander necks, and this should steepen hydraulic gradients and intensify intra-meander hyporheic flux. This research used dye tracking and head loss measurements in a 1:500 planimetrically scaled laboratory river table to quantify the spatial and temporal intensification of intra-meander flux rates at two evolution ages. The younger meander bend, M1, had a sinuosity of 2.3, a river neck width of 0.39 cm, and 0.6% river slope, and the older meander bend, M3, had a sinuosity of 5.2, a river neck width of 0.12 cm, and 0.5% river slope. Flux into and out of the meander bend was estimated along the normalized curvilinear distance s *, with the meander neck at s * = 0.1 and s * = 0.9, the meander centroid at s * = 0.37 and s * = 0.63, and the apex at s * = 0.5. Between the meander centroid and neck we documented a 60% spatial intensification for M1 and a 90% spatial intensification for M3. Between M1 and M3 we documented a 135% temporal intensification at the neck and a 100% intensification at the centroid. Our empirical spatial and temporal intensification rates involving the M1 and the M3 scenario were 1 to 3 times lower than theoretical rates derived from a river evolution model with equivalent M1 and M3 planimetry. Over estimation by the theoretical model was attributed to exaggerated head loss caused by the model neglecting groundwater contributions to river stage. Hyporheic exchange provides critical ecosystem services and its spatial and temporal variation with meander evolution should be considered in river management. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Watershed models which combine hydrology and water quality are being widely utilized in integrated watershed management for the determination of best water management practices. In this study, the hydrology of the Lower Porsuk Stream Watershed in Turkey has been modeled with SWAT to determine optimal water management strategies. The calibration and validation process have been accomplished using data from two monitoring stations. The model has been run for the 1978–2009 period and while the 1998–2004 period has been used for calibration, the validation has spanned the whole period. The SWATCup calibration and uncertainity program has been utilized for this purpose. No significant differences have been detected among different iteration numbers in the calibration period. The monthly Nash-Sutcliffe and R 2 performance indicatiors for the upstream Esenkara station have been 0.74 and 0.88, respectively for the calibration period, and 0.87 and 0.87, respectively for the validation period. The Kıranharmanı station which is located close to the watershed outlet has shown values of 0.59 and 0.72, respectively for the calibration period, and 0.44 and 0.56, respectively for the validation period. There are uncertanities in the abstracted irrigation and groundwater quantities which have reflected in the results in the Kıranharmanı station which is more affected as it lies downstream of the irrigation areas. The effects of different irrigation practices on the flow regime have been also investigated. A scenario has been implemented in which drip irrigation wholly replaces conventional furrow and sprinkler irrigation. The scenario has shown increases in stream flows by 87 % for the whole year. The adoption of more efficient irrigation practices thus results in reducing the water stress induced by irrigation demands. With this study a modeling framework has been founded to aid water management applications in the Lower Porsuk Stream Watershed by generating scenarios for best management practices. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Many observations and studies have shown that water resources amount in the Hai River Basin decreased significantly over the last half of the twentieth century. This study attempts to attribute the observed changes in the water resources amount in the basin over a 40 year period (1961–2000) to different factors, including natural climate variability, climate change induced by anthropogenic forcing of greenhouse gas emissions (referred to as anthropogenic forcing hereafter), and local human activity. First, the temporal variation of the annual water resources amount in the basin during the past 40 years is analyzed by employing the moving-average method, the linear regression method, and the Mann-Kendall method. Second, through setting different scenarios, the effects on the water resources amount due to different factors, including natural climate variability, anthropogenic forcing, and local human activity, are obtained using the parallel climate model, the distributed hydrological model water and energy transfer processes in large river basins, and the statistical downscaling model. Third, the fingerprint-based attribution method is used to obtain the signal strengths of observed changes in water resources amount during 1961–2000 and changes in the water resources amount under different scenarios. Finally, by comparing the signal strengths, the observed changes in water resources amount in the basin can be attributed to different factors. The results indicate that natural climate variability and local human activity may be two factors responsible for the observed changes in the water resources amount during the past 40 years in the basin, with local human activity being the main factor and accounting for about 60% of the changes.

Description: Flood events can induce temporal changes in streambed elevation and particle-size composition, which may influence the bed's hydraulic properties and stream-aquifer fluxes during and after an event. This study combines a set of previously developed modeling approaches to create a synthetic flood event during which bed sediment is entrained and deposited as a function of hydraulic conditions and particle size. One simulated river reach in a state of approximate dynamic equilibrium is chosen to investigate the impacts of size-selective sediment transport on stream-aquifer interaction. Along this reach, the preferential entrainment of fine sediment during the flood's rising limb leads to overall bed coarsening, and increases in vertical hydraulic conductivity (Kbv) and downward fluxes of floodwater into the streambed. Progressively finer sediment layers are deposited during the event's falling limb, causing the redevelopment of a colmation (clogging) layer on the bed surface and a decline in overall Kbv by the event's conclusion. This reduction in Kbv leads to prolonged retention of event water in the streambed (after the reach reverts from losing to gaining river conditions) when compared with what is expected if pre-event Kbv values are used to estimate river-aquifer exchanges. This process of sequential bed coarsening and fining during a flood event provides a mechanistic explanation for the event size-and-duration threshold, inferred in some systems, that must be exceeded for significant amounts of flood recharge to occur. The major consequences of these processes—enhanced infiltration and prolonged floodwater retention—have potentially major implications for groundwater-surface water interactions, water quality, contaminant transport, and riparian biogeochemistry.

Description: Nonstationary oscillation (NSO) processes are observed in a number of hydroclimatic data series. Stochastic simulation models are useful to study the impacts of the climatic variations induced by NSO processes into hydroclimatic regimes. Reproducing NSO processes in a stochastic time series model is, however, a difficult task because of the complexity of the nonstationary behaviors. In the current study, a novel stochastic simulation technique that reproduces the NSO processes embedded in hydroclimatic data series is presented. The proposed model reproduces NSO processes by utilizing empirical mode decomposition (EMD) and nonparametric simulation techniques (i.e., k-nearest-neighbor resampling and block bootstrapping). The model was first tested with synthetic data sets from trigonometric functions and the Rössler system. The North Atlantic Oscillation (NAO) index was then examined as a real case study. This NAO index was then employed as an exogenous variable for the stochastic simulation of streamflows at the Romaine River in the province of Quebec, Canada. The results of the application to the synthetic data sets and the real-world case studies indicate that the proposed model preserves well the NSO processes along with the key statistical characteristics of the observations. It was concluded that the proposed model possesses a reasonable simulation capacity and a high potential as a stochastic model, especially for hydroclimatic data sets that embed NSO processes.

Description: Vegetation zonation and tidal hydrology are basic attributes of intertidal salt marshes, but specific links among vegetation zonation, plant water use, and spatiotemporally dynamic hydrology have eluded thorough characterization. We developed a quantitative model of an intensively studied salt marsh field site, integrating coupled 2-D surface water and 3-D groundwater flow and zonal plant water use. Comparison of model scenarios with and without heterogeneity in (1) evapotranspiration rates and rooting depths, according to mapped vegetation zonation, and (2) sediment hydraulic properties from inferred geological heterogeneity revealed the coupled importance of both sources of ecohydrological variability at the site. Complex spatial variations in root zone pressure heads, saturations, and vertical groundwater velocities emerged in the model but only when both sources of ecohydrological variability were represented together and with tidal dynamics. These regions of distinctive root zone hydraulic conditions, caused by the intersection of vegetation and sediment spatial patterns, were termed “ecohydrological zones” (EHZ). Five EHZ emerged from different combinations of sediment hydraulic properties and evapotranspiration rates, and two EHZ emerged from local topography. Simulated pressure heads and groundwater dynamics among the EHZ were validated with field data. The model and data showed that hydraulic differences between EHZ were masked shortly after a flooding tide but again became prominent during prolonged marsh exposure. We suggest that ecohydrological zones, which reflect the combined influences of topographic, sediment, and vegetation heterogeneity and do not emphasize one influence over the others, are the fundamental spatial habitat units comprising the salt marsh ecosystem.

Description: Quantification of biophysical parameters of urban trees is important for urban planning, and for assessing carbon sequestration and ecosystem services. Airborne lidar has been used extensively in recent years to estimate biophysical parameters of trees in forested ecosystems. However, similar studies are largely lacking for individual trees in urban landscapes. Prediction models to estimate biophysical parameters such as height, crown area, diameter at breast height, and biomass for over two thousand individual trees were developed using best subsets multiple linear regression for a study area in central Oklahoma, USA using point cloud distributional metrics from an Optech ALTM 2050 lidar system. A high level of accuracy was attained for estimating individual tree height (R2 = 0.89), dbh (R2 = 0.82), crown diameter (R2 = 0.90), and biomass (R2 = 0.67) using lidar-based metrics for pooled data of all tree species. More variance was explained in species-specific estimates of biomass (R2 = 0.68 for Juniperus virginiana to 0.84 for Ulmus parviflora) than in estimates from broadleaf deciduous (R2 = 0.63) and coniferous (R2 = 0.45) taxonomic groups—or the data set analysed as a whole (R2 = 0.67). The metric crown area performed particularly well for most of the species-specific biomass equations, which suggests that tree crowns should be delineated accurately, whether manually or using automatic individual tree detection algorithms, to obtain a good estimation of biomass using lidar-based metrics.

Description: In recent decades, copula functions have been applied in bivariate drought duration and severity frequency analysis. Among a number of potential copulas, Clayton has been mostly employed in drought analysis. In this research, we study the influence of the tail shape of various copula functions (i.e. Gumbel, Frank, Clayton, and Gaussian) on drought bivariate frequency analysis. The appropriateness of Clayton copula for the characterization of drought characteristics is also investigated. Drought data are extracted from standardized precipitation index (SPI) time series for four stations in Canada (La Tuque and Grande Prairie) and Iran (Anzali and Zahedan). Both duration and severity datasets are positively skewed. Different marginal distributions were first fitted to drought duration and severity data. The gamma and exponential distributions were respectively selected for drought duration and severity according to the positive skewness and Kolmogorov- Smirnov test. The results of copula modeling show that the Clayton copula function is not an appropriate choice for the employed datasets in the current study, and does not give more drought risk information than an independent model for which the duration and severity dependence is not significant. The reason is that the dependence of two variables in the upper tail of Clayton copula is very weak and similar to the independent case, while the observed data in the transformed domain of cumulative density function shows high association in the upper tail. Instead, the Frank and Gumbel copula functions show better performance than Clayton function for drought bivariate frequency analysis. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Soil hydraulic parameters were upscaled from a 30 m resolution to a 1 km resolution using a new aggregation scheme (described in the companion paper) where the scale parameter was based on the topography. When soil hydraulic parameter aggregation or upscaling schemes ignore the effect of topography, their application becomes limited at hillslope scales and beyond, where topography plays a dominant role in soil deposition and formation. Hence the new upscaling algorithm was tested at the hillslope scale (1 km) across two locations: (1) the Little Washita watershed in Oklahoma, and (2) the Walnut Creek watershed in Iowa. The watersheds were divided into pixels of 1 km resolution and the effective soil hydraulic parameters obtained for each pixel. Each pixel/domain was then simulated using the physically based HYDRUS-3-D modeling platform. In order to account for the surface (runoff/on) and subsurface fluxes between pixels, an algorithm to route infiltration-excess runoff onto downstream pixels at daily time steps and to update the soil moisture states of the downstream pixels was applied. Simulated soil moisture states were compared across scales, and the coarse scale values compared against the airborne soil moisture data products obtained during the hydrology experiment field campaign periods (SGP97 and SMEX02) for selected pixels with different topographic complexities, soil distributions, and land cover. Results from these comparisons show good correlations between simulated and observed soil moisture states across time, topographic variations, location, elevation, and land cover. Stream discharge comparisons made at two gauging stations in the Little Washita watershed also provide reasonably good results as to the suitability of the upscaling algorithm used. Based only on the topography of the domain, the new upscaling algorithm was able to provide coarse resolution values for soil hydraulic parameters which effectively captured the variations in soil moisture across the watershed domains.

Description: Hydropower accounts for about 20% of the worldwide electrical power production. In mountainous regions this ratio is significantly higher. In this study we present how future projected climatic forcing, as described in regional climate models (RCMs), will affect water resources and subsequently hydropower production in downstream hydropower plants in a glacierized alpine valley (Vispa valley, Switzerland, 778 km2). In order to estimate future runoff generation and hydropower production, we used error-corrected and downscaled climate scenarios from regional climate models (RCMs) as well as glacier retreat projections from a dynamic glacier model and coupled them to a physically based hydrological model. Furthermore, we implemented all relevant hydropower operational rules in the hydrological model to estimate future hydropower production based on the runoff projections. The uncertainty of each modeling component (climate projections, glacier retreat, and hydrological projection) and the resulting propagation of uncertainty to the projected future water availability for energy production were assessed using an analysis of variance. While the uncertainty of the projections is considerable, the consistent trends observed in all projections indicate significant changes to the current situation. The model results indicate that future melt- and rainfall-runoff will increase during spring but decline during summer. The study concludes by outlining the most relevant expected changes for hydropower operations.

Description: The effect of spatial concentration fluctuations on the reaction of two solutes, A + B ⇀ C, is considered. In the absence of fluctuations, the concentration of solutes decays as Adet = Bdet ∼ t−1. Contrary to this, experimental and numerical studies suggest that concentrations decay significantly slower. Existing theory suggests a t−d/4 scaling in the asymptotic regime (d is the dimensionality of the problem). Here we study the effect of fluctuations using the classical diffusion-reaction equation with random initial conditions. Initial concentrations of the reactants are treated as correlated random fields. We use the method of moment equations to solve the resulting stochastic diffusion-reaction equation and obtain a solution for the average concentrations that deviates from ∼t−1 to ∼t−d/4 behavior at characteristic transition time t*. We also derive analytical expressions for t* as a function of Damköhler number and the coefficient of variation of the initial concentration.

Description: Based on existing techniques in nonlinear physics that work in the Fourier domain, we develop a multivariate, wavelet-based method for the generation of synthetic discharge time series. This approach not only retains the cross-correlative structure of the original data (which makes it preferable to principal component methods that merely preserve the correlations) but also replicates the nonlinear properties of the original data. We argue that the temporal asymmetry of the typical hydrograph is the most important form of nonlinearity to preserve in the synthetic data. Using the derivative skewness as a measure of asymmetry and an example data set of 35 years of daily discharge data from 107 gauging stations in the United States, we compare two approaches that preserve the asymmetry of the original records. We generate synthetic data and then study the properties of fitting a generalized extreme value distribution to the annual maxima for a total flux time series. The synthetic series provides error bands for the fitted distribution that give a different way of assessing credible return periods. It is found that the best approach for studying extremes is to match the asymmetry of each series individually, rather than to formulate a global threshold criterion.

Description: This article focuses on household water use in Spain by analyzing the influence of a detailed set of factors. We find that, although the presence of both water-saving equipment and water-conservation habits leads to water savings, the factors that influence each are not the same. In particular, our results show that those individuals most committed to the adoption of water-saving equipment and, at the same time, less committed to water-conservation habits tend to have higher incomes.

Description: Wetlands are valuable ecosystems that provide many valuable services, yet many of these important ecosystems are at risk because of current trends in climate change. The Prairie Pothole Region (PPR) in the upper-midwest of the U.S. and south-central Canada, characterized by glacially-sculpted landscapes and abundant wetlands, is one such vulnerable region. According to regional/global climate model predictions, drought occurrence will increase in the PPR region through the 21st century and thus will probably cause the amount of water in wetlands to decline. Water surface area (WSA) of Kidder County, ND from 1984-2011 was measured by classifying TM/ETM+ images through the Modified Normalized Difference Water Index (MNDWI).We then developed a linear model based on the water surface area (WSA) of these wetlands and historical climate data, and used this to determine the wetland sensitivity to climate change and predict future wetlands WSA in the PPR. Our model based on Palmer Drought Severity Index (PDSI) of the current year (PDSI t-0 ) and two years previous (PDSI t-2 ) can explain 79% of the annual wetland WSA variance, suggesting a high sensitivity of wetlands to drought/climate change. We also predicted the PPR wetlands WSA in the 21st century under A1B scenario (a mid-carbon emission scenario) using simulated PDSI based on IPCC AR4 22-model ensemble climate. According to our prediction, the WSA of the PPR wetlands will decrease to less than half of the baseline WSA (defined as the mean wetlands WSA of the 2000s) by the mid of the 21st century, and to less than one-third by the 2080s, and will then slightly increase in the 2090s. This considerable future wetland loss caused only by climate change provides important implication to future wetland management and climate adaptation policy. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In hydrology, the storage-discharge relationship is a fundamental catchment property. Understanding what controls this relationship is at the core of catchment science. To date, there are no direct methods to measure water storage at catchment scales (10 1 -10 3  km 2 ). In this study, we use direct measurements of terrestrial water storage dynamics by means of superconducting gravimetry in a small headwater catchment of the Regen River, Germany, to derive empirical storage-discharge relationships in nested catchments of increasing scale. Our results show that the local storage measurements are strongly related to streamflow dynamics at larger scales (〉 100 km 2 ; correlation coefficient = 0.78-0.81), but at small scale no such relationship exists (~ 1 km 2 ; correlation coefficients = -0.11). The geologic setting in the region can explain both the disconnection between local water storage and headwater runoff, and the connectivity between headwater storage and streams draining larger catchment areas. More research is required to understand what controls the form of the observed storage-discharge relationships at the catchment scale. This study demonstrates that high-precision gravimetry can provide new insights into the complex relationship between state and response of hydrological systems. Copyright © 2012 John Wiley & Sons, Ltd.

Description: We present a novel fusion algorithm for electronic-reconnaissance (ER) satellite and optical imaging satellite data using coherent point set (CPS) analysis. This work is motivated by a large-scale maritime surveillance problem, where ship groups in the observations are of particular interest for tactical and strategic operations. Fusion of observations from ER satellite and optical imaging satellite is a challenging task. On the one hand, dense and continuous measurement is not available for optical imagery. On the other hand, it is difficult to extract robust features from ER measurements. Considering that the size of a ship is often less than the distance among different ships, we treat each ship as a mass point. The contributions of our work are threefold. First, multisensor data fusion is accomplished by CPS association. To the best of our knowledge, this letter is the first to investigate CPS for multimodal remote sensing data fusion. Second, a novel geometry descriptor, which encodes the topological characteristics of a point set, is presented. Third, we combine both topological features and attributive features within the framework of Dempster–Shafer theory for CPS analysis. The proposed method has been tested using different sets of simulated data and recorded data. Experimental results demonstrate the effectiveness of the proposed method.

Description: The 2008 Wenchuan Earthquake that occurred in a mountainous region of China induced massive landslides and caused numerous casualties and property losses. Analyzing the disturbances on vegetation detected from the abnormal sudden drops of the normalized difference vegetation index (NDVI) within a short period can be used for the purpose of rapid landslide identification. Although much research has confirmed the necessity of high-resolution satellite images in landslides identification, Moderate Resolution Imaging Spectroradiometry (MODIS) products still have their usefulness for high temporal resolution, as investigated by the authors. Using MODIS MOD09Q1 NDVI products at a temporal interval of 8 days during 2008, this letter presents a method that has been developed to identify landslide distribution and evolution patterns. First, to find the optimal threshold, the MODIS NDVI time series are analyzed in a training area by an iteration searching procedure. Second, the chosen threshold is used in a larger validation area. To examine the effectiveness of the proposed method, the results are compared to interpreted landslides using SPOT5 images with a spatial resolution of 2.5 m acquired before and after the main shock. An overall 75% accuracy is achieved, and better consistency is observed for landslides extending over one MODIS pixel. The proposed method has also been applied to the Wenchuan earthquake affected areas with seismic intensity IX and greater, and the similar spatial pattern of landslides distribution is obtained when compared with results by using high-resolution images and field investigation. This technique can be applied practically for rapid landslide assessment at a relatively large region after a major earthquake or other severe disturbance events.

Description: The new generation of synthetic aperture radar (SAR) sensors is providing images with very high spatial resolution, improved up to the meter scale. Such a resolution increase allows more accurate monitoring capabilities by means of interferometric approaches. The use of higher frequency enhances the sensitivity of the system even to minute changes, such as thermal dilations. This phenomenon has an impact on the interferometric products, particularly on the deformation velocity maps, if not properly handled. Man-made structures, such as steel core bridges and specific buildings, may be very sensible to thermal dilation effects. By extending the multitemporal differential interferometry SAR processing chains, in our case based on the multidimensional imaging (MDI) approach, an additional parameter related to temperature differences at acquisition instants, the thermal coefficient, can be accurately estimated. This parameter provides interesting perspectives in application to infrastructure monitoring: It brings information about the thermal behavior of the imaged objects. In this letter, we investigate the thermal response of the Musmeci bridge (Potenza, Italy), by experimenting the extended MDI approach on a real TerraSAR-X data set. Results highlight the possibility of such a technique to obtain measurements of the motion that is highly correlated with temperature, thus providing useful information about the static structure of bridges.

Description: The triplet Markov field (TMF) model is powerful in the nonstationary synthetic aperture radar (SAR) image analysis. Taking the speckle noise and the correlation of nonstationarities in two multitemporal SAR images into account, we propose a change-detection method based on the TMF model in this letter. The third field $U$ in the TMF model is redefined to describe the nonstationary textural similarity between the two images for change detection. The corresponding prior energy of $(X, U)$ is reconstructed. The adaptive weight parameter in prior energy is introduced to cope with the detection tradeoff issue. An automatic estimation of the parameter is obtained with low level of complexity. The Bayesian maximum posterior marginal criterion is utilized with the TMF model to obtain change detection. Experimental results on real SAR images validate the superiority of the proposed TMF method over the Markov random field method.

Description: A time series of COSMO-SkyMed (CSK) images is exploited for detection of seasonal snow cover in alpine areas. For the first time, a complete time series of CSK images acquired during snow fall and melt periods in winter 2010–2011 is addressed to verify the snow cover mapping capabilities of X-band radar images under different conditions (from dry to wet snow). The algorithm for snow detection is based on a multitemporal approach with the concept that free water in the snowpack attenuates the X-band synthetic aperture radar signal and wet snow can be classified by comparing images acquired under wet snow and snow-free conditions. Thresholds to make this distinction are compared across all the images to check sensitivity to different winter conditions and land-use classes. The impact of variable and fixed thresholds on the retrieved snow-covered areas is assessed. Snow maps from CSK images compared with Landsat Enhanced Thematic Mapper Plus snow maps indicate a constant underestimation in the detection of snow extent, particularly during winter season, thus showing a scarce sensitivity of X-band signals to snow in dry conditions. Probability of error maps are also calculated for each CSK snow map, thus providing information on the classification error associated to each pixel labeled as snow. The analysis of the snow line variation during spring determines good time consistency in the determination of snow maps from CSK images.

Description: Wavelet packet analysis (WPA) and gray model (GM) are investigated for nonlinear unsupervised feature extraction of hyperspectral remote sensing data in this letter. Treated as derivative series, a hyperspectral response curve of each pixel is decomposed into an approximation and various detailed compositions by WPA, and then, GM is continuously applied to find the relationship among those detailed compositions. Cluster–space representation is used for determining the optimal wavelet. New extracted features can reveal the intrinsic identities of hyperspectral data. Experimental results show the feasibility and reliability of our proposed method in terms of classification accuracy.

Description: This letter proposes a new supervised linear feature extractor for hyperspectral image classification. The criterion for feature extraction is a modified maximal relevance and minimal redundancy (MRMD), which has been used for feature selection until now. The MRMD is a function of mutual information terms, which possess higher order statistics of data; thus, it is effective for hyperspectral data with informative higher order statistics. The batch and stochastic versions of the gradient ascent are performed on the MRMD to find the optimal parameters of a linear feature extractor. Preliminary results achieve better classification performance than the traditional methods based on the first- and second-order moments of data.

Description: The high-resolution imaging of a highly squinted synthetic aperture radar remains difficult because of the severe coupling between the range and the azimuth. “Squint minimization” compensates for the range walking in the azimuth time domain, which efficiently increases the orthogonality between the range and the azimuth. However, this “squint minimization” introduces the azimuth space-variant phases, which can be compensated by the azimuth nonlinear chirp-scaling (ANCS) algorithm using large computational loads. In this letter, an azimuth overlapped subaperture algorithm (AOSA) is proposed to compensate for these phases in the Doppler frequency domain. The validity constraint of this algorithm is then analyzed. The AOSA has an advantage over ANCS in terms of the computational load and is considerably more suitable for real-time processing.

Description: The potential of polarimetric synthetic aperture radar data for the soil surface characterization of bare agricultural soils was investigated by using air- and spaceborne data acquired by Radar Aéroporté Multi-Spectral d'Etude des Signatures (RAMSES), Système Expérimental de Télédétection Hyperfréquence Imageur (SETHI), and RADARSAT-2 sensors over several study sites in France. Fully polarimetric data at ultrahigh frequency, X-, C-, L-, and P-bands were compared. The results show that the main polarimetric parameters studied (entropy, $alpha$ angle, and anisotropy) are not very sensitive to the variation of the soil surface parameters. Low correlations are observed between the polarimetric and soil parameters (moisture content and surface roughness). Thus, the polarimetric parameters are not very relevant to the characterization of the soil surface over bare agricultural areas.

Description: Airborne dual-polarization observations of sea surface normalized radar cross section (NRCS) were conducted over the North Atlantic during January–February 2011. Observations were made using the University of Massachusetts' Imaging Wind and Rain Airborne Profiler radar system installed on the National Oceanic and Atmospheric Administration's WP-3D research aircraft during several winter storm events to determine the high-wind response of the sea surface NRCS for both horizontal and vertical polarizations. During the flights, the aircraft performed several constant-roll circle maneuvers to allow collection of NRCS over a range of incidence angles. We find consistency with prior reports in the polarization ratio observed at moderate incidence angles at the winds encountered. For larger incidence angles, we observe a measurable decrease in polarization ratio with increasing wind speed.

Description: In semiarid regions, the rooting strategies employed by vegetation can be critical to its survival. Arid regions are characterized by high variability in the arrival of rainfall, and species found in these areas have adapted mechanisms to ensure the capture of this scarce resource. Vegetation roots have strong control over this partitioning, and assuming a static root profile, predetermine the manner in which this partitioning is undertaken.A coupled, dynamic vegetation and hydrologic model, tRIBS + VEGGIE, was used to explore the role of vertical root distribution on hydrologic fluxes. Point-scale simulations were carried out using two spatially and temporally invariant rooting schemes: uniform: a one-parameter model and logistic: a two-parameter model. The simulations were forced with a stochastic climate generator calibrated to weather stations and rain gauges in the semiarid Walnut Gulch Experimental Watershed (WGEW) in Arizona. A series of simulations were undertaken exploring the parameter space of both rooting schemes and the optimal root distribution for the simulation, which was defined as the root distribution with the maximum mean transpiration over a 100-yr period, and this was identified. This optimal root profile was determined for five generic soil textures and two plant-functional types (PFTs) to illustrate the role of soil texture on the partitioning of moisture at the land surface. The simulation results illustrate the strong control soil texture has on the partitioning of rainfall and consequently the depth of the optimal rooting profile. High-conductivity soils resulted in the deepest optimal rooting profile with land surface moisture fluxes dominated by transpiration. As we move toward the lower conductivity end of the soil spectrum, a shallowing of the optimal rooting profile is observed and evaporation gradually becomes the dominate flux from the land surface. This study offers a methodology through which local plant, soil, and climate can be accounted for in the parameterization of rooting profiles in semiarid regions.

Description: There has been a recent debate in the hydrological community about the relative merits of the informal generalized likelihood uncertainty estimation (GLUE) approach to uncertainty assessment in hydrological modeling versus formal probabilistic approaches. Some recent literature has suggested that the methods can give similar results in practice when properly applied. In this note, we show that the connection between formal Bayes and GLUE is not merely operational but goes deeper, with GLUE corresponding to a certain approximate Bayesian procedure even when the “generalized likelihood” is not a true likelihood. The connection we describe relates to recent approximate Bayesian computation (ABC) methods originating in genetics. ABC algorithms involve the use of a kernel function, and the generalized likelihood in GLUE can be thought of as relating to this kernel function rather than to the model likelihood. Two interpretations of GLUE emerge, one as a computational approximation to a Bayes procedure for a certain “error-free” model and the second as an exact Bayes procedure for a perturbation of that model in which the truncation of the generalized likelihood in GLUE plays a role. The intent of this study is to encourage cross-fertilization of ideas regarding GLUE and ABC in hydrologic applications. The connection we outline suggests the possibility of combining a formal likelihood with a kernel based on a generalized likelihood within the ABC framework and also allows advanced ABC computational methods to be used in GLUE applications. The model-based interpretation of GLUE may also be helpful in partially illuminating the implicit assumptions in different choices of generalized likelihood.

Description: Satellite-passive microwave remote sensing has been extensively used to estimate snow water equivalent (SWE) in northern regions. Although passive microwave sensors operate independent of solar illumination and the lower frequencies are independent of atmospheric conditions, the coarse spatial resolution introduces uncertainties to SWE retrievals due to the surface heterogeneity within individual pixels. In this article, we investigate the coupling of a thermodynamic multilayered snow model with a passive microwave emission model. Results show that the snow model itself provides poor SWE simulations when compared to field measurements from two major field campaigns. Coupling the snow and microwave emission models with successive iterations to correct the influence of snow grain size and density significantly improves SWE simulations. This method was further validated using an additional independent data set, which also showed significant improvement using the two-step iteration method compared to standalone simulations with the snow model.

Description: In this paper we present a visual analytics approach for deriving spatio-temporal patterns of collective human mobility from a vast mobile network traffic data set. More than 88 million movements between pairs of radio cells—so-called handovers—served as a proxy for more than two months of mobility within four urban test areas in Northern Italy. In contrast to previous work, our approach relies entirely on visualization and mapping techniques, implemented in several software applications. We purposefully avoid statistical or probabilistic modeling and, nonetheless, reveal characteristic and exceptional mobility patterns. The results show, for example, surprising similarities and symmetries amongst the total mobility and people flows between the test areas. Moreover, the exceptional patterns detected can be associated to real-world events such as soccer matches. We conclude that the visual analytics approach presented can shed new light on large-scale collective urban mobility behavior and thus helps to better understand the “pulse” of dynamic urban systems.

Description: Existing analytical solutions to determine aquifer response to a change in stream stage are inappropriate where an unsaturated zone exists beneath the stream, as in the case of disconnected stream-aquifer systems. A better understanding of the relationship between aquifer response and transient stream stage in disconnected systems is therefore required, as this would also aid in the field determination of the status of connection between the stream and aquifer. We use a numerical model to examine transient stream stage and the corresponding water table response. Beneath disconnected streams, the magnitude of head change in the water table level is a balance between the cumulative infiltration during a flow event and the rate at which the water can disperse laterally. Increases in wave duration, stream width, and streambed permeability result in greater infiltrated water volume and therefore a higher peak response at the water table. Conversely, higher aquifer transmissivity and aquifer hydraulic conductivity allow the water to move laterally away from the stream faster, resulting in a smaller head change below the stream. Lower unsaturated storage results in a greater and faster aquifer response because the unsaturated zone can fill more quickly. Under some combinations of parameters, the magnitude of the disconnected head response is more than seven times greater than the change in stream stage driving streambed infiltration; an effect which can never occur beneath a connected stream. The results of this sensitivity analysis are compared to field data from a river in eastern Australia to determine periods of disconnection. Where the change in aquifer head is greater than the change in stream stage, disconnection between the stream and aquifer can be determined.

Description: The ability to understand and predict the flux and fate of sediment delivered to the sea by rivers remains an outstanding scientific challenge. Approaches to this challenge are necessarily synthetic, spanning wide ranges in spatial and temporal scales. Here a conventional sediment transport theory used by engineers and sedimentologists at reach and channel scales is applied at the basin scale. Specifically, a straightforward expression proposed by Bagnold and modified accordingly predicts the observed importance of combined wetness and steepness of a source basin as a control of sediment supply to the sea. The reasonable, key assumption underlying the application of sediment transport theory in this context is that the river-mouth sites for which suspended-sediment loads are reported are alluviated, and thus characterized by transport-limited flux of sediment. This analysis also indicates the potential significance of additional, as yet poorly documented factors constraining sediment supply to the sea. These factors, some of which appear to covary systematically with climate, include river-profile concavity, river-mouth channel width and friction, and the characteristic size of sediment in transport.

Description: Key documents such as the European Water Framework Directive and the U.S. Clean Water Act state that public and stakeholder participation in water resource management is required. Participation aims to enhance resource management and involve individuals and groups in a democratic way. Evaluation of participatory programs and projects is necessary to assess whether these objectives are being achieved and to identify how participatory programs and projects can be improved. The different methods of evaluation can be classified into three groups: (i) process evaluation assesses the quality of participation process, for example, whether it is legitimate and promotes equal power between participants, (ii) intermediary outcome evaluation assesses the achievement of mainly nontangible outcomes, such as trust and communication, as well as short- to medium-term tangible outcomes, such as agreements and institutional change, and (iii) resource management outcome evaluation assesses the achievement of changes in resource management, such as water quality improvements. Process evaluation forms a major component of the literature but can rarely indicate whether a participation program improves water resource management. Resource management outcome evaluation is challenging because resource changes often emerge beyond the typical period covered by the evaluation and because changes cannot always be clearly related to participation activities. Intermediary outcome evaluation has been given less attention than process evaluation but can identify some real achievements and side benefits that emerge through participation. This review suggests that intermediary outcome evaluation should play a more important role in evaluating participation in water resource management.

Description: A planning support system for land consolidation has been developed that has, at its heart, an expert system called LandSpaCES (Land Spatial Consolidation Expert System) which contains a “design module” that generates alternative land redistributions under different scenarios and an “evaluation module” which integrates GIS with multi-criteria decision making for assessing these alternatives. This paper introduces the structural framework of the latter module which has been applied using a case study in Cyprus. Two new indices are introduced: the “parcel concentration coefficient” for measuring the dispersion of parcels; and the “landowner satisfaction rate” for predicting the acceptance of the land redistribution plan by the landowners in terms of the location of their new parcels. These two indices are used as criteria for the evaluation of the land redistribution alternatives and are transferable to any land consolidation project. Moreover, a modified version of the ratio estimation procedure, referred to as the “qualitative rating method” for assigning weights to the evaluation criteria, is presented, along with a set of non-linear value functions for standardizing the performance scores of the alternatives and incorporating expert knowledge for five evaluation criteria. The application of the module showed that it is a powerful new tool for the evaluation of alternative land redistribution plans that could be implemented in other countries after appropriate adjustments. A broader contribution has also been made to spatial planning processes, which might follow the methodology and innovations presented in this paper.

Description: The key question that is asked in this study is “how are the three independent bias components of satellite rainfall estimation, comprising hit bias, missed, and false precipitation, physically related to the estimation uncertainty of soil moisture and runoff for a physically based hydrologic model?” The study also investigated the performance of different satellite rainfall products as a function of land use and land cover (LULC) type. Using the entire Mississippi river basin as the study region and the variable infiltration capacity (VIC)-3L as the distributed hydrologic model, the study of the satellite products (CMORPH, 3B42RT, and PERSIANN-CCS) yielded two key findings. First, during the winter season, more than 40% of the rainfall total bias is dominated by missed precipitation in forest and woodland regions (southeast of Mississippi). During the summer season, 51% of the total bias is governed by the hit bias, and about 42% by the false precipitation in grassland-savanna region (western part of Mississippi basin). Second, a strong dependence is observed between hit bias and runoff error, and missed precipitation and soil moisture error. High correlation with runoff error is observed with hit bias (∼0.85), indicating the need for improving the satellite rainfall product's ability to detect rainfall more consistently for flood prediction. For soil moisture error, it is the total bias that correlated significantly (∼0.78), indicating that a satellite product needed to be minimized of total bias for long-term monitoring of watershed conditions for drought through continuous simulation.

Description: Bayesian inference is used to study the effect of precipitation and model structural uncertainty on estimates of model parameters and confidence limits of predictive variables in a conceptual rainfall-runoff model in the snow-fed Rudbäck catchment (142 ha) in southern Finland. The IHACRES model is coupled with a simple degree day model to account for snow accumulation and melt. The posterior probability distribution of the model parameters is sampled by using the Differential Evolution Adaptive Metropolis (DREAM(ZS)) algorithm and the generalized likelihood function. Precipitation uncertainty is taken into account by introducing additional latent variables that were used as multipliers for individual storm events. Results suggest that occasional snow water equivalent (SWE) observations together with daily streamflow observations do not contain enough information to simultaneously identify model parameters, precipitation uncertainty and model structural uncertainty in the Rudbäck catchment. The addition of an autoregressive component to account for model structure error and latent variables having uniform priors to account for input uncertainty lead to dubious posterior distributions of model parameters. Thus our hypothesis that informative priors for latent variables could be replaced by additional SWE data could not be confirmed. The model was found to work adequately in 1-day-ahead simulation mode, but the results were poor in the simulation batch mode. This was caused by the interaction of parameters that were used to describe different sources of uncertainty. The findings may have lessons for other cases where parameterizations are similarly high in relation to available prior information.

Description: Operational land surface models (LSMs) compute hydrologic states such as soil moisture that are needed for a range of important applications (e.g., drought, flood, and weather prediction). The uncertainty in LSM parameters is sufficiently great that several researchers have proposed conducting parameter estimation using globally available remote sensing data to identify best fit local parameter sets. However, even with in situ data at fine modeling scales, there can be significant remaining uncertainty in LSM parameters and outputs. Here, using a new uncertainty estimation subsystem of the NASA Land Information System (LIS) (described herein), a Markov chain Monte Carlo (MCMC) technique is applied to conduct Bayesian analysis for the accounting of parameter uncertainties. The Differential Evolution Markov Chain (DE-MC) MCMC algorithm was applied, for which a new parallel implementation was developed. A case study is examined that builds on previous work in which the Noah LSM was calibrated to passive (L-band) microwave remote sensing estimates of soil moisture for the Walnut Gulch Experimental Watershed. In keeping with prior related studies, the parameters subjected to the analysis were restricted to the soil hydraulic properties (SHPs). The main goal is to estimate SHPs and soil moisture simulation uncertainty before and after consideration of the remote sensing data. The prior SHP uncertainty is based on the original source of the standard SHP lookup tables for the Noah LSM. Conclusions are drawn regarding the value and viability of Bayesian analysis over alternative approaches (e.g., parameter estimation, lookup tables) and further research needs are identified.

Description: Most diversity management programs in Canada maintain that enhancing workforce diversity is of tremendous significance for business organizations in today’s competitive global urban markets. Since well-meaning diversity management initiatives have been largely ineffective thus far in dealing with workplace discrimination and racism in the Canadian workplace, this paper underscores the need to decenter the focus of diversity management from a business imperative to an antidiscrimination and social justice imperative. Within this latter perspective, the paper examines the strengths and limitations of the antiracism approach that has been implemented in various developed countries in recent years. The antiracism approach is an action-oriented strategy for institutional and systemic change that has at its core the interrogation of privilege, power disparities, and other forms of inequity within the organization. Drawing from the lessons of various initiatives that have utilized this approach, the present paper emphasizes the need for a nuanced antiracism approach in the multicultural Canadian society if diversity management is to attain its goal of greater inclusion of all individuals in informal networks and formal organizational programs.

Description: This paper examines the dynamic spatial and temporal patterns of China’s urban hierarchy from 1950 to 2005. We limit the study to mainland China’s 137 urban agglomerations of 750,000 or more population as of 2005. The paper improves upon a classic approach to measuring shifting ranks within an urban hierarchy by applying advanced spatial analysis techniques. We use a Getis-Ord , a space-time Moran scatter plot, and dynamic LISA paths to examine the regional difference and change for these urban agglomerations. Our study revealed a north-south divide in the changing urban hierarchy of China after 1990. The analysis demonstrated that the spatial and temporal shift of urban dominance in China was closely associated with policy and economic factors. The paper concludes with a discussion of the differences across the six different time periods of change from 1950 to 2005.

Description: It has long been known that surface gravity waves induce significant seepage through the porous layer found at the lake bottom. Away from coastal regions, however, the pressure signature of surface waves at the lake bottom is weak. We consider fully nonlinear internal gravity waves, whose long wavelength and slow motion implies a sustained and strong pressure perturbation even in the deep regions of the lake. We argue that internal waves can induce significant seepage through the sediment layer, in regions where surface gravity waves have negligible impact. The pressure profile at the fluid–porous layer interface is computed from the “exact” Dubreil-Jacotin-Long theory, giving a reliable profile even for large waves. This profile is used in conjunction with Darcy's law to compute the seepage within the porous region. We find that the geometric distribution of seepage is strongly controlled by both the ratio of porous media thickness to the horizontal length of the pressure perturbation, and the bottom topography, when it is present. Based on work on the interaction of internal solitary waves with the bottom boundary layer, we develop a model to account for the changes in permeability due to wave-induced instabilities in the bottom boundary layer and enhanced benthic turbulence. This turbulence acts to unplug the pores near the surface by lifting the detritus that clogs them. The resulting changes in permeability significantly enhance exchange between the free fluid and the porous medium on the downstream side of the wave.

Description: Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Description: : Exempt domestic water well provisions exist in sixteen states and allow landowners to use water wells producing relatively small amounts of water or water wells used for domestic or other specific purposes to avoid certain regulatory requirements. Exempt wells provisions apply mainly in the west, where most states use the prior appropriation rule (“first in time, first in right”) for ground water rights. Critics voice concerns over implications ranging from land use planning to water supply. This article explores existing state regulation and finds that the term “exempt well” is a misnomer. The vast majority of states impose numerous restrictions on exempt wells. Existing regulations include limitations on the quantity, location, and amount of irrigation use. Other provisions address construction standards, required filing of information and other practices. The review of regulations suggests that many steps to mitigate the impact of exempt wells are already in place.

Description: : National attention has turned to a case currently pending before the New Mexico Supreme Court. Bounds v. State of New Mexico involves a facial constitutional challenge to New Mexico's domestic well statute. The Plaintiff, Horace Bounds, claims that the statute, which requires the State Engineer to issue permits to any applicants for domestic purposes, violates the constitutional provision known as the Prior Appropriation Doctrine. The decision of the Court will likely have impacts across the country and beyond the legal realm. This decision could impact the well drilling industry, economic development, and administrative agencies. Numerous other western states have similar exempt well provisions, and the decision by the New Mexico Supreme Court could well lead to litigation or pushes for changes in legislation on other states as well.

Description: : Over the last decade, Montana's fastest growing basins have experienced a surge of reliance on a certain “permit-exempt well loophole,” which enables small ground water diversions to avoid having to go through the regulatory process that determines whether the proposed use would adversely impact senior water users and ground water supply. This loophole has generated significant conflict with senior right owners in areas where exempt wells have provided water for large-scale residential developments. As a partial solution, this article examines ground water mitigation exchanges designed to address similar issues in Washington and presents the 63rd Montana Legislature with recommendations for a Montana pilot project.

Description: : Most western states allow landowners to withdraw certain amounts of groundwater for specified purposes without obtaining a water right permit. As water demands increase in the West, some developers are using exempt wells to supply water to residential developments without acquiring the permits and water rights needed to build public water supply systems. This has led to concerns that the cumulative impact of these wells could impair senior water rights, create environmental problems, and threaten water quality in some areas. This article explores the considerations associated with mitigating exempt well impacts and argues that surgical approaches focused on specific issues and geographic areas (scalpels) are more politically and administratively feasible than broad efforts to repeal or reduce exemptions on a statewide basis (hammers). The views expressed in this article are those of the author alone and do not represent those of the Western States Water Council or its member states.

Description: : The saga over exempt wells in the western United States and Canada epitomizes a new type of water conflict – a spaghetti-western water war. The political melodrama stars local governments to serve as sheriff of water-supply planning duties. Exempt wells number in the millions, and herding the growing numbers is testing the mettle of the states and provinces responsible for the management, allocation, and protection of natural resources. The separation of laws governing ground water and surface water, coupled with changes in geography and geology within a jurisdiction, compound the administrative riddle and give rise to a broad spectrum of conflicts, from differing interpretations of hydrogeologic data, economic impacts associated with increasing the herd, to differing identities associated with the use of ground water from the exempt wells. Despite the political melodrama of exempt wells, there is room and willingness for other trails and paths to keep the herd intact. This paper describes the different breeds of conflicts associated with exempt wells and gives examples of how the mysterious stranger of collaborative decision making processes and water governance systems can ride into town and lead to successful water management and conflict resolution.

Description: : The water management practice of storing water underground for use during dry periods is a very effective technique that has been practiced for a long time. Another water management practice that has been around for a long time is collection of rainwater and storage of rainwater aboce ground. There have been significant advances more recently in our understanding of aquifer geology and technologies used for rainwater harvesting. This article combines our advance in knowledge of storing and recovery water from wells and new rainwater harvesting practices to recharge a low yielding exempt well with treated rainwater. The studies presented in this article demonstrate that exempt wells provide a valuable underground water storage option and may provide a new effective urban water management tool.

Description: Predicting long-term consequences of climate change on hydrologic processes has been limited due to the needs to accommodate the uncertainties in hydrological measurements for calibration, and to account for the uncertainties in the models that would ingest those calibrations and uncertainties in climate predictions as basis for hydrological predictions. We implemented a hierarchical Bayesian (HB) analysis to coherently admit multiple data sources and uncertainties including data inputs, parameters, and model structures to identify the potential consequences of climate change on soil moisture and streamflow at the head watersheds ranging from low to high elevations in the southern Appalachian region of the United States. We have considered climate change scenarios based on three greenhouse gas emission scenarios of the Interovernmental Panel on Climate Change: A2, A1B and B1 emission scenarios. Full predictive distributions based on hierarchical Bayesian models are capable of providing rich information and facilitating the summarization of prediction uncertainties. With predictive uncertainties taken into account, the most pronounced change in soil moisture and streamflow would occur under the A2 scenario at both low and high elevations, followed by the A1B scenario and then by the B1 scenario. Uncertainty in the change of soil moisture is less than that of streamflow for each season, especially at high elevations. A reduction of soil moisture in summer and fall, a reduction or slight increase of streamflow in summer, and an increase of streamflow in winter are predicted for all three scenarios at both low and high elevations. The hydrological predictions with quantified uncertainties from a HB model could aid more-informed water resource management in developing mitigation plans, and dealing with water security under climate change. Copyright © 2012 John Wiley & Sons, Ltd.

Description: This research proposes a combination of SWAT and MODFLOW, MD-SWAT-MODFLOW, to address the multi-aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi-aquifers so that the daily streamflow can be simulated correctly. A multi-aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in-between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river-aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD-SWAT-MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the LH-OAT sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi-aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio-temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Large areas of Amazonian evergreen forest experience seasonal droughts extending for three or more months, yet show maximum rates of photosynthesis and evapotranspiration during dry intervals. This apparent resilience is belied by disproportionate mortality of the large trees in manipulations that reduce wet season rainfall, occurring after 2–3 years of treatment. The goal of this study is to characterize the mechanisms that produce these contrasting ecosystem responses. A mechanistic model is developed based on the ecohydrological framework of TIN (Triangulated Irregular Network)-based Real Time Integrated Basin Simulator + Vegetation Generator for Interactive Evolution (tRIBS+VEGGIE). The model is used to test the roles of deep roots and soil capillary flux to provide water to the forest during the dry season. Also examined is the importance of “root niche separation,” in which roots of overstory trees extend to depth, where during the dry season they use water stored from wet season precipitation, while roots of understory trees are concentrated in shallow layers that access dry season precipitation directly. Observational data from the Tapajós National Forest, Brazil, were used as meteorological forcing and provided comprehensive observational constraints on the model. Results strongly suggest that deep roots with root niche separation adaptations explain both the observed resilience during seasonal drought and the vulnerability of canopy-dominant trees to extended deficits of wet season rainfall. These mechanisms appear to provide an adaptive strategy that enhances productivity of the largest trees in the face of their disproportionate heat loads and water demand in the dry season. A sensitivity analysis exploring how wet season rainfall affects the stability of the rainforest system is presented.

Description: Companion animals are very important to people of big cities. In Greece the families which own those animals take them for a walk in the parks and streets of their town but not always with care to not disturb other citizens. Laws and regulations for companion and stray animals are not respected in Greece, although as a result of the Athens Olympics in 2004 the Greek government started to provide funds for the collection, care, and sterilization of stray animals. This paper is a first attempt to record, through the aid of a structured questionnaire, the view of the citizens of the city of Xanthi in northern Greece regarding companion animals and what they do when animals become old. The paper also examines the extent to which the existence of stray and companion animals in parks disturbs the people who visit green areas looking for a place to get a rest or play on the grass. The majority preferred that stray animals are collected off the streets, recorded, treated from parasites, vaccinated, and sterilized and after that are available for adoption or are returned to the area they were found at.

Description: We present direct numerical simulation results for both isothermal and density-stratified turbulent flow in an open channel into and around a 120° bend with a bulk Reynolds number of 7500 and Prandtl number of 1.5. The bend is sharp with a radius-to-channel breadth ratio of 1.5. The bulk Richardson number for the stratified flow is 2.4 based on overall channel depth. The gradient Richardson number (Rig) varies between 10 and 20 at the entrance to θ ≈ 60°, where θ is the angular location. Above θ ≈ 60 − 120, Rig ≈ 1. In isothermal flow, the well-known helical flow structure is observed. In stratified conditions, the vertical variation in relative strength of the outward-directed baroclinic pressure gradient and the centrifugal acceleration leads to a more complex circulation structure. In the near bed region and immediately above the interface, the centrifugal acceleration is greater, driving flow radially inward, while just below the density interface the baroclinic pressure gradient is greater, leading to outward-directed flow. This produces a four layer circulation structure with potentially significant implications for sediment erosion and transport. Additionally, this produces a complex dynamic at the density interface where the shear orientation varies through approximately 200° over the mixing layer depth.

Description: In the majority of large river systems, flow is regulated and/or otherwise affected by operational and management activities, such as ship locking. The effect of lock operation on sediment-water oxygen fluxes was studied within a 12.9 km long impoundment at the Saar River (Germany) using eddy-correlation flux measurements. The continuous observations cover a time period of nearly 5 days and 39 individual locking events. Ship locking is associated with the generation of surges propagating back and forth through the impoundment which causes strong variations of near-bed current velocity and turbulence. These wave-induced flow variations cause variations in sediment-water oxygen fluxes. While the mean flux during time periods without lock operation was 0.5 ± 0.1 g m−2 d−1, it increased by about a factor of 2 to 1.0 ± 0.5 g m−2 d−1 within time periods with ship locking. Following the daily schedule of lock operations, fluxes are predominantly enhanced during daytime and follow a pronounced diurnal rhythm. The driving force for the increased flux is the enhancement of diffusive transport across the sediment-water interface by bottom-boundary layer turbulence and perhaps resuspension. Additional means by which the oxygen budget of the impoundment is affected by lock-induced flow variations are discussed.

Description: Existing models of spatial relations do not consider that different concepts exist on different levels in a hierarchy and in turn that the spatial relations in a given scene are a function of the specific concepts considered. One approach to determining the existence of a particular spatial relation is to compute the corresponding high level concepts explicitly using map generalization before inferring the existence of the spatial relation in question. We explore this idea through the development of a model of the spatial relation “enters” that may exist between a road and a housing estate.

Description: Model errors are inevitable in any prediction exercise. One approach that is currently gaining attention in reducing model errors is by combining multiple models to develop improved predictions. The rationale behind this approach primarily lies on the premise that optimal weights could be derived for each model so that the developed multimodel predictions will result in improved predictions. A new dynamic approach (MM-1) to combine multiple hydrological models by evaluating their performance/skill contingent on the predictor state is proposed. We combine two hydrological models, “abcd” model and variable infiltration capacity (VIC) model, to develop multimodel streamflow predictions. To quantify precisely under what conditions the multimodel combination results in improved predictions, we compare multimodel scheme MM-1 with optimal model combination scheme (MM-O) by employing them in predicting the streamflow generated from a known hydrologic model (abcd model or VIC model) with heteroscedastic error variance as well as from a hydrologic model that exhibits different structure than that of the candidate models (i.e., “abcd” model or VIC model). Results from the study show that streamflow estimated from single models performed better than multimodels under almost no measurement error. However, under increased measurement errors and model structural misspecification, both multimodel schemes (MM-1 and MM-O) consistently performed better than the single model prediction. Overall, MM-1 performs better than MM-O in predicting the monthly flow values as well as in predicting extreme monthly flows. Comparison of the weights obtained from each candidate model reveals that as measurement errors increase, MM-1 assigns weights equally for all the models, whereas MM-O assigns higher weights for always the best-performing candidate model under the calibration period. Applying the multimodel algorithms for predicting streamflows over four different sites revealed that MM-1 performs better than all single models and optimal model combination scheme, MM-O, in predicting the monthly flows as well as the flows during wetter months.

Description: The potential for riverine drinking source water to become contaminated with pathogens is related to the production and transport of fecal waste from within the local catchment area. Identifying specific relationships between land-use types and fecal contamination in riverine water provides an indication of the risk associated with land-use change and helps to target mitigation measures toward land-use types of concern. Fecal coliform (FC) data from 42 riverine sites across British Columbia (BC), Canada, were examined in relation to land-use composition (including 16 land-use types) in the local catchment area. FC concentration significantly increased in relation to anthropogenic land-use impacts but was negatively associated with undisturbed and high-elevation land types. Regression tree analysis identified that highest FC concentrations occurred in catchments characterized by more than 12.5% agricultural land and more than 1.6% urban land. Furthermore, the risk of violation of the BC partial treatment raw drinking water quality guideline for FC concentration (100 CFU 100 mL−1) increased in relation to agricultural impacts. Additional factors, such as sewage treatment discharge, low dilution in smaller streams, and higher temperatures, were associated with higher FC concentration among sites with similar levels of agricultural development. These results identify land-use types that present the greatest threat to riverine contamination, namely agricultural and urban land, and indicate the proportion of such land use associated with high contamination. Land use should be managed and source water protection should be targeted in light of these results so as to minimize the risk of surface water exposure to fecal contaminants.

Description: We report on the calibration of the one-dimensional hydrodynamic lake model DYRESM to simulate the water temperature conditions of the pre-alpine Lake Ammersee (South-east Germany) which is representative of deep and large lakes in this region. Special focus is given to the calibration in order to reproduce the correct thermal distribution and stratification including the time of onset and duration of summer stratification. To ensure the application of the model to investigate the impact of climate change on lakes, an analysis of the model sensitivity under stepwise modification of meteorological input parameters (air temperature, wind speed, precipitation, global radiation, cloud cover, vapor pressure, and tributary water temperature) was conducted. The total mean error of the calibration results is –0.23 °C, the root mean square error amounts to 1.012 °C. All characteristics of the annual stratification cycle were reproduced accurately by the model. Additionally, the simulated deviations for all applied modifications of the input parameters for the sensitivity analysis can be differentiated in the high temporal resolution of monthly values for each specific depth. The smallest applied alteration to each modified input parameter caused a maximum deviation in the simulation results of at least 0.26 °C. The most sensitive reactions of the model can be observed through modifications of the input parameters air temperature and wind speed. Hence the results show that further investigations at Lake Ammersee, such as coupling the hydrodynamic model with chemo-dynamic models to assess the impact of changing climate on biochemical conditions within lakes, can be carried out using DYRESM. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In semi-arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing-connected, transitional and losing-disconnected) was assessed at 12 sites along six Murray-Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing-disconnected conditions. A mixture of losing-disconnected, losing-connected and gaining conditions was found among the 12 sites. The losing-disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5 – 2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing-disconnected sites but never lower than 1 m below riverbed level at the losing-connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing-connected rivers by increased infiltration rates, while the decline could not be compensated at the losing-disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Abstact The distribution of groundwater inflows in a stream reach plays a major role in controlling the stream temperature, a vital component shaping the riverine ecosystem. In this study, the DTS system was installed in a small Danish lowland stream, Elverdamsåen, to assess the seasonal dynamics of groundwater inflow zones using high spatial (1 m) and temporal (3 minutes) resolution of water temperature measurements. Four simple criteria consisting of 30 min average temperature at 16:00, mean and standard deviation of diurnal temperatures, and the day-night temperature difference were applied to three DTS datasets representing stream temperature responses to the variable meteorological and hydrological conditions prevailing in summer, winter and spring. The standard deviation criterion was useful to identify groundwater discharge zones in summer and spring conditions, while the mean temperature criterion was better for the winter conditions. In total, 20 interactions were identified from the DTS datasets representing summer, 16 in winter and 19 in spring, albeit with only two interactions contributing in all three seasons. Higher baseflow to streamflow ratio, antecedent precipitation and presence of fractured clayey till in the stream reach were deemed as the vital factors causing apparent seasonal variation in the locations of upwelling zones, prompting use of DTS not only in preconceived scenarios of large diurnal temperature change but rather a long term deployment covering variable meteorological and hydrological scenarios. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Non-point source pollution is a key issue in integrated river basin management around the world, and has resulted in water contamination, aquatic ecology deterioration and eutrophication. Xin'anjiang Catchment is the key drinking water source area for Hangzhou City, China. A promising model (Soil and Water Assessment Tool: SWAT) was applied to assess the non-point source pollution and its effect on drinking water. Sensitivity analysis of model parameters was carried out using the SUFI-2 sensitivity technique. Water discharge, sediment, total nitrogen and total phosphorus loads processes from 2000 to 2010 were simulated and the spatial distributions of non-point source pollutants were evaluated at the catchment and administrative country levels. The results show that the hydrological parameters of SWAT were dominantly sensitive for non-point source pollution simulation, including CN2, RCHRG_DP, ALPHA_BF, SOL_AWC, ESCO and SOL_K, and the characteristic parameters of sub-basins (viz. HRU_SLP, SLSUBBSN). Also, water quality parameters (viz. CH_EROD, NPERCO, RSDCO and PPERCO, PHOSKD, etc.) have a significant effect on nutrients. The model performance was very satisfactory, especially for runoff, sediment, and TP simulation. The non-point source pollutant load increased from 2001 to 2010 in the whole catchment. Total nitrogen (TN) load increased from 3,428 tons (0.59 ton/km 2 ) to 7,315 tons (1.25 ton/km 2 ) and total phosphorus (TP) load increased from 299 tons (0.05 ton/km 2 ) to 867 tons (0.15 ton/km 2 ). The contribution of rice land was the largest, accounting for nearly 95%, followed by tea garden (3.56%), winter wheat (1.37%), forest (0.07%) and grassland (0.02%). Moreover, She County and Xiuning County contributed more than half of the non-point source pollutants. This study was expected to provide a method and reference for non-point source pollution quantification, and to support water quality management implementation in China. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The increased socio-economic relevance of flood risk assessment has led to the development of innovative methodologies for the hydraulic simulation of river and floodplain systems, and has promoted the development of new techniques for flood hazard and inundation mapping (e.g. Di Baldassarre et al., 2010; Vorogushyn et al., 2010). In particular, one-dimensional (1D) and two-dimensional (2D) hydraulic models have been used more and more as numerical tools (e.g. Aronica et al., 2002; Hesselink et al., 2003; Horritt et al., 2007; Pappenberger et al., 2005) as these models have proven to be able to effectively simulate river hydraulics and floodplain inundation at different levels of detail (e.g. Horritt & Bates, 2001, 2002). Flood inundation models appear also to be useful tools for the reconstruction and analysis of historical events (e.g. Di Baldassarre et al., 2009; Horritt et al., 2010), which can be very important to provide a comprehensive assessment of exposure to floods and to develop flood risk management plans as required by the recent Floods Directive 2007/60/EC (European Commission, 2007). Copyright © 2012 John Wiley & Sons, Ltd.

Description: A methodology to enhance satellite precipitation estimation using unsupervised dimensionality reduction (UDR) techniques is developed. This enhanced technique is an extension to the precipitation estimation from remotely sensed imagery using an artificial neural network (PERSIANN) and cloud classification system (CCS) method (PERSIANN-CCS) enriched using wavelet features combined with dimensionality reduction. Cloud-top brightness temperature measurements from the Geostationary Operational Environmental Satellite (GOES)-12 are used for precipitation estimation at 4 km $times$ 4 km spatial resolutions every 30 min. The study area in the continental U.S. covers parts of Louisiana, Arkansas, Kansas, Tennessee, Mississippi, and Alabama. Based on quantitative measures, root mean square error and Heidke skill score (HSS), the results show that the UDR techniques can improve the precipitation estimation accuracy. In addition, the independent component analysis is shown to have better performance than other UDR techniques; and in some cases, it achieves 10% improvement in the HSS.

Description: In the absence of model deficiencies, simulation results at the correct parameter values lead to an unbiased description of observed data with remaining deviations due to observation errors only. However, this ideal cannot be reached in the practice of environmental modeling, because the required simplified representation of the complex reality by the model and errors in model input lead to errors that are reflected in biased model output. This leads to two related problems: First, ignoring bias of output in the statistical model description leads to bias in parameter estimates, model predictions and, in particular, in the quantification of their uncertainty. Second, as there is no objective choice of how much bias to accept in which output variable, it is not possible to design an “objective” model calibration procedure. The first of these problems has been addressed by introducing a statistical (Bayesian) description of bias, the second by suggesting the use of multiobjective calibration techniques that cannot easily be used for uncertainty analysis. We merge the ideas of these two approaches by using the prior of the statistical bias description to quantify the importance of multiple calibration objectives. This leads to probabilistic inference and prediction while still taking multiple calibration objectives into account. The ideas and technical details of the suggested approach are outlined and a didactical example as well as an application to environmental data are provided to demonstrate its practical feasibility and computational efficiency.

Description: The hydraulics of step-pool streams are characterized by rapidly varied flow at the step crest, a hydraulic jump, and gradually varied flow in the pool unit of the step-pool sequence. The flow characteristics at the step crests act as the hydraulic control for the water surface profile within the upstream pool unit. Using both field and flume investigations, we demonstrate the use of weir flow concepts for assessing and categorizing the hydraulic characteristics of natural step-crests in step-pool streams. We categorize the results of our investigations in terms of the crest-clast, planform, longitudinal, and instream wood geometries of the step crests. The broad-crested weir equation can be expressed as Q = C* g0.5Wh3/2, where Q is the flowrate, C* is a dimensionless discharge coefficient, W is the crest width, g is the acceleration of gravity, and h is the upstream flow depth above the step crest. Although the flow over a natural step is generally more complex than for an engineered weir, the results of our investigations indicate that the C*-value for simulated and natural steps increases linearly as a function of the upstream head (h), with C* values ranging from 0.15 to 0.97. As a result, the application of weir flow concepts to natural steps provides means for (1) indirectly estimating flow rates; (2) characterizing the hydraulics for individual steps; (3) defining external and/or internal boundary conditions at step crests for hydraulic model simulations of natural or restored step-pool streams; and (4) estimating the upstream pressure force acting on step-crest clasts.