Assessing the impacts of rainfall intensity and urbanization on storm runoff in an arid catchment

In the past few decades, rapid urbanization has occurred in many regions of the Kingdom of Saudi Arabia due to increasing population and urban development. Additionally, the effects of global warming on rainfall characteristics have been observed. This rapid change in urbanization and climate change has cause significant changes in the nature of land surfaces and rainfall patterns, which affect the runoff process and the amount of surface runoff during floods. This study investigated the effect of urbanization and rainfall intensity for Hafr Al-Batin watershed located in Saudi Arabia. For this purpose, a hydrologic model, HEC-HMS, was adopted to simulate the flow of different rainfall intesities and urbanization levels. Simulated results showed that for a 100-year storm, a 24-h duration, and an urbanization level of 80%, the peak flow was 213% higher than the estimated current peak and the runoff volume was 112% higher than the current runoff volume. These results show a strong linear correlation between the level of urbanization and both peak discharge and runoff volume. Furthermore, the results indicate that for short return periods, the peak flow is more sensitive to the level of urbanization compared to long periods.     

Areal extent of the 1988 flood in Bangladesh: How much did the satellite imagery show?

This paper presents the methods and results of visual interpretation of satellite imagery for estimating areal extent of the 1988 flood in Bangladesh. The main limitation of this method was that cloud-free images were not available for the peak flood period in the first week of September. Relatively cloud-free images on three selected dates in the last three weeks of September showed that areas of inundation on those dates ranged from 44000 to 60000 km² ± 5% (31 to 42% of Bangladesh). The Government of Bangladesh official estimates of the areal extent of flooding were tentative in nature, since these were based on a choropleth map of relative proportions (percentages) of flooding in different parts of the country. The official estimate of the maximum extent of flooding was 82000 km² (57% of Bangladesh), i.e. 15 to 26% larger than the area shown on satellite imagery. The actual extent of flooding was certainly larger than the estimates from satellite imagery but, perhaps, smaller than the official estimates, since a reduction of 15 to 26% of flood area in one to three weeks seemed unlikely. An alternative method of mapping flood-affected areas by using newspaper-interpreted data was attempted, but the method had limited values because of reporting bias.     

Martian water vapor: Mars Express PFS/LW observations

We present the seasonal and geographical variations of the martian water vapor monitored from the Planetary Fourier Spectrometer Long Wavelength Channel aboard the Mars Express spacecraft. Our dataset covers one martian year (end of Mars Year 26, Mars Year 27), but the seasonal coverage is far from complete. The seasonal and latitudinal behavior of the water vapor is globally consistent with previous datasets, Viking Orbiter Mars Atmospheric Water Detectors (MAWD) and Mars Global Surveyor Thermal Emission Spectrometer (MGS/TES), and with simultaneous results obtained from other Mars Express instruments, OMEGA and SPICAM. However, our absolute water columns are lower and higher by a factor of 1.5 than the values obtained by TES and SPICAM, respectively. In particular, we retrieve a Northern midsummer maximum of 60 pr-μm, lower than the 100-pr-μm observed by TES. The geographical distribution of water exhibits two local maxima at low latitudes, located over Tharsis and Arabia. Global Climate Model (GCM) simulations suggest that these local enhancements are controlled by atmospheric dynamics. During Northern spring, we observe a bulge of water vapor over the seasonal polar cap edge, consistent with the northward transport of water from the retreating seasonal cap to the permanent polar cap. In terms of vertical distribution, we find that the water volume mixing ratio over the large volcanos remains constant with the surface altitude within a factor of two. However, on the whole dataset we find that the water column, normalized to a fixed pressure, is anti-correlated with the surface pressure, indicating a vertical distribution intermediate between control by atmospheric saturation and confinement to a surface layer. This anti-correlation is not reproduced by GCM simulations of the water cycle, which do not include exchange between atmospheric and subsurface water. This situation suggests a possible role for regolith-atmosphere exchange in the martian water cycle.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part II: Mesoscale Modelling

The importance of soil moisture inputs and improved model physics in the prediction of the daytime boundary-layer structure during the Southern Great Plains Hydrology Experiment 1997 (SGP97) is investigated using the non-hydrostatic fifth-generation Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) Mesoscale Model MM5. This is Part II of a two-part study examining the relationship of surface heterogeneity to observed boundary-layer structure. Part I focuses on observations and utilizes a simple model while Part II uses observations and MM5 modelling. Soil moisture inputs tested include a lookup table based on soil type and season, output from an offline land-surface model (LSM) forced by atmospheric observations, and high-resolution ( 800 m) airborne microwave remotely sensed data. Model physics improvements are investigated by comparing an LSM directly coupled with the MM5 to a simpler force-restore method at the surface. The scale of land surface heterogeneities is compared to the scale of their effects on boundary-layer structure.The use of more detailed soil moisture fields allowed the MM5 to better represent the large-scale (hundreds of km) and small-scale (tens of km) horizontal gradients in surface-layer weather and, to a lesser degree, the atmospheric boundary-layer (ABL) height, which was evaluated against observations measured by differential absorption lidar (DIAL). The benefits of coupling an LSM to the MM5 were not readily evident in this summertime case, with the model having particular difficulty simulating the timing of maximum surface fluxes while underestimating the depth of the mixed layer.     

Feedbacks between erosion and sediment transport in experimental bedrock channels

Natural bedrock rivers flow in self‐formed channels and form diverse erosional morphologies. The parameters that collectively define channel morphology (e.g. width, slope, bed roughness, bedrock exposure, sediment size distribution) all influence river incision rates and dynamically adjust in poorly understood ways to imposed fluid and sediment fluxes. To explore the mechanics of river incision, we conducted laboratory experiments in which the complexities of natural bedrock channels were reduced to a homogenous brittle substrate (sand and cement), a single sediment size primarily transported as bedload, a single erosion mechanism (abrasion) and sediment‐starved transport conditions. We find that patterns of erosion both create and are sensitive functions of the evolving bed topography because of feedbacks between the turbulent flow field, sediment transport and bottom roughness. Abrasion only occurs where sediment impacts the bed, and so positive feedback occurs between the sediment preferentially drawn to topographic lows by gravity and the further erosion of these lows. However, the spatial focusing of erosion results in tortuous flow paths and erosional forms (inner channels, scoops, potholes), which dissipate flow energy. This energy dissipation is a negative feedback that reduces sediment transport capacity, inhibiting further incision and ultimately leading to channel morphologies adjusted to just transport the imposed sediment load. Copyright © 2007 John Wiley & Sons, Ltd.     

Large-scale three-dimensional inversion of EarthScope MT data using the integral equation method

In this paper we apply 3D inversion to MT data collected in the Northwestern United States as a part of the EarthScope project. By the end of 2009 MT data had been collected from 262 stations located throughout Oregon, Washington, Idaho, and most of Montana and Wyoming. We used data from 139 MT stations in this analysis. We developed fully parallelized rigorous 3D MT inversion software based on the integral equation method with variable background conductivity. We also implemented a receiver footprint approach which considerably reduced the computational resources needed to invert the large volumes of data covering vast areas. The data set used in the inversion was obtained through the Incorporated Research Institutions for Seismology (IRIS). The inversion domain was divided into 2.7 M cells. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region.     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.     

Stream-driven, high-density gravelly traction carpets: possible deposits in the Trabeg Conglomerate Formation, SW Ireland and some theoretical considerations of their origin

Within high-density flood flows a prominent mechanism of gravel transport and deposition is by stream-driven, high-density traction carpet (with a rheology similar to grain flow). These gravel carpets are envisaged to form the basal portion of a bipartite high-density flood flow, decoupled from an overlying sand- and silt-laden turbulent flow. Several examples already documented in the literature are reviewed and an additional case from the Lower Old Red Sandstone of southwest Ireland is presented. Two mechanisms of traction carpet initiation are discussed: by rapid entrainment of gravel into suspension on rising stage, followed by settling into the gravel traction carpet at peak and falling stage; and by overconcentration of a ‘normal’, low-density bedload. Gravel entrainment, suspension and traction carpet development are significantly easier if the flood water already carries a high concentration of sand and silt in suspension. Theoretical consideration further shows that gravelly traction carpets can be maintained in channels of relatively low gradient by the shear stress exerted by the high-density, sand-bearing turbulent flood flow above. This tangential shear stress is converted to dispersive pressure, which aids buoyancy and quasi-static grain-to-grain contacts in the support of the clasts within the gravel carpet. The carpet is thought to have a quasi-plastic rheology but behave much like a viscous fluid at high shear rates. Stream-driven gravelly traction carpets are expected to produce sheet-like units of clast- to matrix-supported conglomerate, characterized by a parallel or an a(p)a(i) clast fabric. These units may be ungraded, normally or inversely graded, depending on the rate of shear, the viscosity of the flow and the celerity of deposition.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.