Variation in curve numbers derived from plot runoff data for New South Wales (Australia)

The curve number method is a simple one parameter (the curve number) rainfall runoff model. While its theoretical underpinning has been questioned it remains a powerful hydrological tool in the absence of detailed data and is therefore used extensively in hydrological models. This study aims to characterize the variation in maximum retention values (S), which underlie curve numbers, for a range of agricultural treatments across a large spatial area in New South Wales (NSW), Australia. The data used for the analysis spans several decades of rainfall runoff observations. A range of different derivation methods result in variation in mean and variance of S. In particular, methods that emphasize the larger storms result in greater S and thus lower runoff. For larger spatial scales, emphasis on larger storms gives more reliable estimates of S. Systematic variation in S arises from variations in treatment, pre‐runoff soil moisture, rainfall depth, and variations in cover. On the basis of the analysis, a table of curve number values for different land uses found in NSW is presented. The resulting distributions of S and curve numbers provide guidance for rainfall runoff modelling studies in the agricultural important areas of NSW. Copyright © 2011 John Wiley & Sons, Ltd.     

Impacts into non-polar ice-rich paleodeposits on Mars: Excess ejecta craters, perched craters and pedestal craters as clues to Amazonian climate history

We compare three previously independently studied crater morphologies - excess ejecta craters, perched craters, and pedestal craters - each of which has been proposed to form from impacts into an ice-rich surface layer. Our analysis identifies the specific similarities and differences between the crater types; the commonalities provide significant evidence for a genetic relationship among the morphologies. We use new surveys of excess ejecta and perched craters in the southern hemisphere in conjunction with prior studies of all of the morphologies to create a comprehensive overview of their geographic distributions and physical characteristics. From these analyses, we conclude that excess ejecta craters and perched craters are likely to have formed from the same mechanism, with excess ejecta craters appearing fresh while perched craters have experienced post-impact modification and infilling. Impacts that led to these two morphologies overwhelmed the ice-rich layer, penetrating into the underlying martian regolith, resulting in the excavation of rock that formed the blocky ejecta necessary to armor the surface and preserve the ice-rich deposits. Pedestal craters, which tend to be smaller in diameter, have the same average deposit thickness as excess ejecta and perched craters, and form in the same geographic regions. They rarely have ejecta around their crater rims, instead exhibiting a smooth pedestal surface. We interpret this to mean that they form from impacts into the same type of ice-rich paleodeposit, but that they do not penetrate through the icy surface layer, and thus do not generate a blocky ejecta covering. Instead, a process related to the impact event appears to produce a thin, indurated surface lag deposit that serves to preserve the ice-rich material. These results provide a new basis to identify the presence of Amazonian non-polar ice-rich deposits, to map their distribution in space and time, and to assess Amazonian climate history. Specifically, the ages, distribution and physical attributes of the crater types suggest that tens to hundreds of meters of ice-rich material has been episodically emplaced at mid latitudes in both hemispheres throughout the Amazonian due to obliquity-driven climate variations. These deposits likely accumulated more frequently in the northern lowlands, resulting in a larger population of all three crater morphologies in the northern hemisphere.     

Identifying regions based on flexible user-defined constraints

The identification of regions is both a computational and conceptual challenge. Even with growing computational power, regionalization algorithms must rely on heuristic approaches in order to find solutions. Therefore, the constraints and evaluation criteria that define a region must be translated into an algorithm that can efficiently and effectively navigate the solution space to find the best solution. One limitation of many existing regionalization algorithms is a requirement that the number of regions be selected a priori. The recently introduced max-p algorithm does not have this requirement, and thus the number of regions is an output of, not an input to, the algorithm. In this paper, we extend the max-p algorithm to allow for greater flexibility in the constraints available to define a feasible region, placing the focus squarely on the multidimensional characteristics of the region. We also modify technical aspects of the algorithm to provide greater flexibility in its ability to search the solution space. Using synthetic spatial and attribute data, we are able to show the algorithm’s broad ability to identify regions in maps of varying complexity. We also conduct a large-scale computational experiment to identify parameter settings that result in the greatest solution accuracy under various scenarios. The rules of thumb identified from the experiment produce maps that correctly assign areas to their ‘true’ region with 94% average accuracy, with nearly 50% of the simulations reaching 100% accuracy.     

Water Balance Study for a Basin Integrating Remote Sensing Data and GIS

Water balance of a basin involves estimation of input precipitation, runoff, infiltration and evapotranspiration (ET). Although ET may have large variations over a big basin, it is commonly estimated using a few point measurements and this makes the estimation error prone. Satellite based remote sensing data provides few parameters for estimation of energy fluxes, at the land surface and atmosphere interaction in a distributed manner using the meteorological parameters. These parameters through surface energy balance equation have been used for the estimation of ET in this study. Various spatially distributed variables required for ET estimation; viz. NDVI, surface albedo, surface temperature etc. have been derived using remote sensing and ancillary data for Tapi basin located in western India. Beside this field data such as rainfall, air temperature, relative humidity, sunshine hours etc. have been used. For computation of runoff, Soil Conservation Services (SCS) approach has been considered. Tapi basin up to Ukai dam has been selected as the study area. Satellite data from National Oceanic and Atmospheric Administration (NOAA), Polar Orbiting Environmental Satellite, which carries the Advanced Very High Resolution Radiometer (AVHRR), have been used for preparation of various maps required for runoff and ET analysis. The results of runoff and ET have been compared with observed data for 2 years, 2002–2003 and the results have been found in good agreement with observed data.     

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.     

The strain seismograms of P- and S-waves of a local event recorded by four-gauge borehole strainmeter

At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters(FGBS) are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.