Observing Global Surface Water Flood Dynamics

Flood waves moving along river systems are both a key determinant of globally important biogeochemical and ecological processes and, at particular times and particular places, a major environmental hazard. In developed countries, sophisticated observing networks and ancillary data, such as channel bathymetry and floodplain terrain, exist with which to understand and model floods. However, at global scales, satellite data currently provide the only means of undertaking such studies. At present, there is no satellite mission dedicated to observing surface water dynamics and, therefore, surface water scientists make use of a range of sensors developed for other purposes that are distinctly sub-optimal for the task in hand. Nevertheless, by careful combination of the data available from topographic mapping, oceanographic, cryospheric and geodetic satellites, progress in understanding some of the world’s major river, floodplain and wetland systems can be made. This paper reviews the surface water data sets available to hydrologists on a global scale and the recent progress made in the field. Further, the paper looks forward to the proposed NASA/CNES Surface Water Ocean Topography satellite mission that may for the first time provide an instrument that meets the needs of the hydrology community.     

Geomechanical property estimation of unconventional reservoirs using seismic data and rock physics

An extension of a previously developed rock physics model is made that quantifies the relationship between the ductile fraction of a brittle/ductile binary mixture and the isotropic seismic reflection response. By making a weak scattering (Born) approximation and plane wave (eikonal) approximation, with a subsequent ordering according to the angles of incidence, singular value decomposition analyses are performed to understand the stack weightings, number of stacks, and the type of stacks that will optimally estimate two fundamental rock physics parameters – the ductile fraction and the compaction and/or diagenesis. It is concluded that the full PP stack, i.e., sum of all PP offset traces, and the “full” PS stack, i.e., linear weighted sum of PS offset traces, are the two optimal stacks needed to estimate the two rock physics parameters. They dominate over both the second‐order amplitude variation offset “gradient” stack, which is a quadratically weighted sum of PP offset traces that is effectively the far offset traces minus the near offset traces, and the higher order fourth order PP stack (even at large angles of incidence). Using this result and model‐based Bayesian inversion, the seismic detectability of the ductile fraction (shown by others to be the important rock property for the geomechanical response of unconventional reservoir fracking) is demonstrated on a model characteristic of the Marcellus shale play.     

Characterization of complex fluvial systems using remote sensing of spatial and temporal water level variations in the Amazon,Congo, and Brahmaputra Rivers

The Surface Water and Ocean Topography (SWOT) satellite mission will provide global, space‐based estimates of water elevation, its temporal change, and its spatial slope in fluvial environments, as well as across lakes, reservoirs, wetlands, and floodplains. This paper illustrates the utility of existing remote sensing measurements of water temporal changes and spatial slope to characterize two complex fluvial environments. First, repeat‐pass interferometric SAR measurements from the Japanese Earth Resources Satellite are used to compare and contrast floodplain processes in the Amazon and Congo River basins. Measurements of temporal water level changes over the two areas reveal clearly different hydraulic processes at work. The Amazon is highly interconnected by floodplain channels, resulting in complex flow patterns. In contrast, the Congo does not show similar floodplain channels and the flow patterns are not well defined and have diffuse boundaries. During inundation, the Amazon floodplain often shows sharp hydraulic changes across floodplain channels. The Congo, however, does not show similar sharp changes during either infilling or evacuation. Second, Shuttle Radar Topography Mission measurements of water elevation are used to derive water slope over the braided Brahmaputra river system. In combination with in situ bathymetry measurements, water elevation and slope allow one to calculate discharge estimates within 2.3% accuracy. These two studies illustrate the utility of satellite‐based measurements of water elevation for characterizing complex fluvial environments, and highlight the potential of SWOT measurements for fluvial hydrology. Copyright © 2010 John Wiley & Sons, Ltd.     

GIS Methodology for Characterizing Historical Conditions of the Willamette River Flood Plain, Oregon

Recent environmental developments have stimulated an interest in conservation and restoration of the historical Willamette River flood plain, both to protect against flooding and to provide wildlife habitat. In order to best utilize scarce resources, we characterized historical and modern river channel and flood‐plain conditions to evaluate changes and help prioritize restoration sites. Using cartographic and photographic data sources, we developed a Geographic Information System (GIS) to map active channels, side channels, islands and tributaries for four separate dates, as well as riparian and flood‐plain vegetation characteristics for pre‐European settlement and modern time periods. Coverages based on flood records and other boundaries were used to partition the flood plain into spatial subsets for analysis. The GIS allowed comparisons between historical and present conditions for a variety of environmental factors. Much of the pre‐settlement channel complexity has been removed. Total channel length in 1995 was 26% less than in 1850, with almost 58% of the river's side channels disconnected from the system. In addition, we found a 72% loss of flood‐plain forest from 1850 to 1995, since it was converted to agricultural and urban land uses. Selected river and flood‐plain variables were made available for a spatial model to prioritize potential locations for flood‐plain restoration.     

Subspace inversion of electromagnetic data: application to mid-ocean-ridge exploration

Controlled-source electromagnetic (CSEM) surveys have the ability to provide tomo-graphic images of electrical conductivity within the Earth. the interpretation of such data sets has long been hampered by inadequate modelling and inversion techniques. In this paper, a subspace inversion technique is described that allows electric dipole-dipole data to be inverted for a 2-D electrical conductivity model more efficiently than with existing techniques. the subspace technique is validated by comparison with conventional inversion methods and by inverting data collected over the East Pacific Rise in 1989. A model study indicates that, with adequate data, a variety of possible mid-ocean-ridge conductivity models could be distinguished on the basis of a CSEM survey.     

The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II Λ cold dark matter simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as the structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions, the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo cannot be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.     

Backyard spectroscopy and photometry of Titan, Uranus and Neptune

We present and discuss observations of Titan (photometry and spectroscopy) and Uranus and Neptune (spectroscopy only) obtained using off-the-shelf equipment, affordable and available to many amateur astronomers and small colleges. Spectral observations compare well with published results from front-line observatories, and some evidence of seasonal change is evident in both spectra and narrowband photometry. Scattered Saturn light presents a significant problem for Titan observations, in particular, for slitless spectroscopy, and our attempts to reproduce Titan's lightcurve have so far been unsuccessful.     

The fall of the St-Robert meteorite

Abstract The St-Robert (Québec, Canada) meteorite shower occurred on 1994 June 15 at 0h02m UT accompanied by detonations audible for >200 km from the fireball endpoint. The fireball was recorded by visual observers in Vermont, New York State, New Hampshire, Québec and Ontario as well as by optical and infrared sensors in Earth-orbit. Penetration to an altitude of 36 km occurred ～60 km to the northeast of Montreal, where the bolide experienced several episodes of fragmentation. A total of 20 fragments of this H5 chondrite, comprising a total mass of 25.4 kg, were recovered in an ellipse measuring 8 × 3.5 km. One fragment of the shower partially penetrated the aluminum roof of a shed. Interpretation of the visual and satellite data suggests that the fireball traveled from south-southwest to north-northeast, with a slope from the horizontal of 55°–61°. A statistical evaluation of the likely heliocentric orbits for the body prior to collision with the Earth, coupled with theoretical modeling of the entry, suggests an entry velocity in the range of 12.7–13.3 km/s; the meteoroid had moved in a low-inclination orbit, with orbital perihelion located extremely close to the Earth's orbit. From satellite optical data, it is found that the photometric mass consumed during the largest detonation is ～1200 kg. Estimation of the amplitude of the acoustic signal detected by the most distant observer yields a source energy near 0.5 kt TNT equivalent energy, which corresponds to a mass of order 10 metric tonnes. This measure is uncertain to approximately one order of magnitude. Theoretical modeling of the entry of the object suggests a mass near 1600 kg. Cosmogenic radionuclide activities constrain the lower initial mass to be ～700 kg with an upper limit near 4000 kg. Seismic data possibly associated with the fireball suggest extremely poor coupling between the airwave and the ground. The total mass estimated to have reached the ground is ～100 kg (in material comprising >55 g fragments), while the preatmospheric mass is found to be most probably in the range of 1200–2000 kg.