THE PROBLEM OF LONG-TERM STORAGE IN RESERVOIRS

Abstract

The geometry of a meandering stream depends strongly on the relative stream size (Q ^2/5 / g ^1/5)/D, on the valley slope, Sv, and on the charge, Q _s/Q, where Q and Q _s are the fluid and sediment discharges respectively, g is acceleration due to gravity and D is the mean sediment size. The geometry depends less strongly on the relative settling size of sediment, D/(v ^2/3 / g ^1/3), where v is the kinematic viscosity. For constant values of Q, S _v and D, the effect of increase of charge reduces the meander length, M _L, and the mean channel surface width, B, whereas meander width, M _B, bend radius, R _M, and mean channel depth, H, increase, For a constant value of (Q ^2/5 / g^1/5)/D the values of M _L, M _B, R _M and B increase with the increase of valley slope but the value of H tends to decrease.     

Lithologic control of bedrock meander dimensions in the appalachian valley and ridge province: A comment

Braun's finding that the wavelength of bedrock meanders decreases as rock hardness increases is reconciled with Hack's and Tinkler's finding that meander wavelength is greater in bedrock than alluvial channels. In populations of alluvial or bedrock channels, bank resistance controls meander wavelength via its influence on channel shape and is negatively correlated with meander wavelength. In contrast, in mixed populations of alluvial and bedrock channels, bank resistance controls meander wavelength through its effect on the recurrence interval of the channel-forming discharge and channel size and is positively correlated with meander wavelength.     

Coarse sediment supply sets the slope of bedrock channels in rapidly uplifting terrain: Field and topographic evidence from eastern Taiwan

The incision rate and steepness of bedrock channels depend on water discharge, uplift rate, substrate lithology, sediment flux, and bedload size. However, the relative role of these factors and the sensitivity of channel steepness to rapid (>1 mm yr⁻¹) uplift rates remain unclear. We conducted field and topographic analyses of fluvial bedrock channels with varying channel bed lithology and sediment source rock along the Coastal Range in eastern Taiwan, where uplift rates vary from 1.8 to 11.8 mm yr⁻¹ and precipitation is relatively consistent (1.5–2.7 m yr⁻¹), to evaluate the controls on bedrock channel steepness. We find that channel steepness is independent of rock uplift rate and annual precipitation but increases monotonically with sediment size and substrate strength. Furthermore, in reaches with uniform substrate lithology (mudstone and flysch), channel steepness systematically varies with sediment source rock but not with channel width. When applied to our data, a mechanistic incision model (saltation-abrasion model) suggests that the steepness of Coastal Range channels is set primarily by coarse-sediment supply. We also observe that larger particles are mainly composed of resistant lithologies derived from volcanic rocks and conglomerates. This result implies that hillslope bedrock properties in the source area exert a dominant control on the steepness of proximal channels through coarse-sediment production in this setting. We propose that channel steepness may be insensitive to uplift rate and flow discharge in fast-uplifting landscapes where incision processes are set by coarse sediment size and supply. Models assuming a proportionality between incision rate and basal shear stress (stream power) may not fully capture controls on fluvial channel profiles in landslide-dominated landscapes. Processes other than channel steepening, such as enhanced bedload impacts and debris-flow scour, may be required to balance rock uplift and incision in these transport-limited systems.     

Tailings fans and valley‐spur cutoffs created by hydraulic mining

Sand and gravel tailings from nineteenth century open‐pit hydraulic gold mines formed large alluvial fans at tributary con?uences in the northwestern Sierra Nevada, California. In the Bear River watershed, several of these fans were so large that they blocked main channels for decades. Some channels not only aggraded deeply, but also moved laterally and cut across the inner bends of valley spurs. Now locked in bedrock channels, these valley‐spur cutoffs impose local controls on geomorphic, hydraulic, and sedimentary processes. One cutoff has incised 25 m into bedrock over the past century (25 cm a^?1) with rapid initial incision rates of up to 50 cm a^?1 (1884–1890). Recognition of spur cutoffs in the geological record may help to identify large landslides and provide an analogue for a type of natural earth?ll dam spillway not prone to catastrophic failures. Tailing fans, valley‐spur cutoffs, and the sediment they trap are described from contemporary accounts and recent ?eld conditions in the Bear River watershed. These anthropogenic changes represent a major shift in the watershed from supply‐limited to transport‐limited sediment budgets and a change in geomorphic processes away from long‐term drainage evolution dominated by ingrown meanders. The large volumes of mining sediment stored in these landforms will be slowly released over the next millennium and could be signi?cant to contemporary ecological and public health issues due to recent ?ndings of high mercury loadings associated with hydraulic mines. Copyright © 2004 John Wiley & Sons, Ltd.     

Coda Q Estimates in the Koyna Region, India

—The coda Q, Q _c ?, have been estimated for the Koyna region of India. The coda waves of 76 seismograms from thirteen local earthquakes, recorded digitally in the region during July–August, 1996, have been analyzed for this purpose at nine central frequencies viz., 1.5, 2.0, 3.0, 4.0, 6.0, 8.0, 12.0, 16.0 and 24.0 Hz using a single backscattering model. All events with magnitude less than 3 fall in the epicentral distances less than 60 km and have focal depths which range from 0.86 to 9.43 km. For the 30 sec coda window length the estimated Q _c values vary from 81 to 261 at 1.5 Hz and 2088 to 3234 at 24 Hz, whereas the mean values of Q _c with the standard error vary from 148 ± 13.5 at 1.5 Hz to 2703 ± 38.8 at 24 Hz. Both the estimated Q _c values and their mean values exhibit the clear dependence on frequency in the region and a frequency dependence average attenuation relationship, Q _c = 96f ^1.09, has been obtained for the region, covering an approximate area of 11500 km² with the surfacial extent of about 120 km and depth of 60 km.¶Lapse time dependence of Q _c has also been studied for the region, with the coda waves analyzed at five lapse time windows from 20 to 60 sec duration with the difference of 10 sec. The frequency dependence average Q _c relationships obtained at these window lengths Q _c = 66f ^1.16 (20 sec), Q _c = 96f ^1.09 (30 sec), Q _c =131f ^1.04 (40 sec), Q _c = 148f ^1.04 (50 sec), Q _c = 182f ^1.02 (60 sec) show that the frequency dependence (exponentn) remains mostly stationary at all the lapse time window lengths, while the change in Q ₀ value is significant. Lapse time dependence of Q _c in the region is also interpreted as the function of depth.     

Coda wave attenuation in the Parecis Basin,Amazon Craton,Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos Gaúchos, Amazon craton—Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q _c) of coda waves modeled using Q_c = Q₀ f^hQ_{\rm c} =Q_{\rm 0} f^\eta , where Q ₀ is the coda quality factor at frequency of 1 Hz and η is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239–252, 2008), based in the temporal attenuation coefficient, χ(f) instead of Q(f), given by the equation c(f)=g+\fracpfQ_e \chi (f)\!=\!\gamma \!+\!\frac{\pi f}{Q_{\rm e} }, for the calculation of the geometrical attenuation (γ) and effective attenuation (Q_e^-1 )(Q_{\rm e}^{-1} ). Q _c values have been computed at central frequencies (and band) of 1.5 (1–2), 3.0 (2–4), 6.0 (4–8), 9.0 (6–12), 12 (8–16), and 18 (12–24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin Q_c = (98±12)f^(1.14±0.08)Q_{\rm c} =(98\pm 12)f^{(1.14\pm 0.08)}, for the surrounding shield Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)}, and for the whole region of Porto dos Gaúchos Q_c = (99±19)f^(1.17±0.02)Q_{\rm c} =(99\pm 19)f^{(1.17\pm 0.02)}. Using the independent frequency model, we found: for the cratonic zone, γ = 0.014 s^− 1, Q_e^-1 = 0.0001Q_{\rm e}^{-1} =0.0001, ν ≈ 1.12; for the basin zone with sediments of ~500 m, γ = 0.031 s^− 1, Q_e^-1 = 0.0003Q_{\rm e}^{-1} =0.0003, ν ≈ 1.27; and for the Parecis basin with sediments of ~1,000 m, γ = 0.047 s^− 1, Q_e^-1 = 0.0005Q_{\rm e}^{-1} =0.0005, ν ≈ 1.42. Analysis of the attenuation factor (Q _c) for different values of the geometrical spreading parameter (ν) indicated that an increase of ν generally causes an increase in Q _c, both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, Q_c = (78±23)f^(1.17±0.14)Q_{\rm c} =(78\pm 23)f^{(1.17\pm 0.14)} (in the deepest part of the basin), than in the basement, Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)} (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.     

Coda Q in the Kachchh Basin, Western India Using Aftershocks of the Bhuj Earthquake of January 26, 2001

Q_C-estimates of Kachchh Basin in western India have been obtained in a high frequency range from 1.5 to 24.0 Hz using the aftershock data of Bhuj earthquake of January 26, 2001 recorded within an epicentral distance of 80 km. The decay of coda waves of 30 sec window from 186 seismograms has been analysed in four lapse time windows, adopting the single backscattering model. The study shows that Q_c is a function of frequency and increases as frequency increases. The frequency dependent Q_c relations obtained for four lapse-time windows are: Q_c=82 f^1.17 (20–50 sec), Q_c=106 f^1.11 (30–60 sec), Q_c=126f^1.03 (40–70 sec) and Q_c=122f^1.02 (50–80 sec). These empirical relations represent the average attenuation properties of a zone covering the surface area of about 11,000, 20,000, 28,000 and 38,000 square km and a depth extent of about 60, 80, 95, 110 km, respectively. With increasing window length, the degree of frequency dependence, n, decreases marginally from 1.17 to 1.02, whereas Q₀ increases significantly from 82 to 122. At lower frequencies up to 6 Hz, Q_c⁻¹ of Kachchh Basin is in agreement with other regions of the world, whereas at higher frequencies from 12 to 24 Hz it is found to be low.     

The effects of rock uplift and rock resistance on river morphology in a subduction zone forearc,Oregon, USA

Relationships between riverbed morphology, concavity, rock type and rock uplift rate are examined to independently unravel the contribution of along-strike variations in lithology and rates of vertical deformation to the topographic relief of the Oregon coastal mountains. Lithologic control on river profile form is reflected by convexities and knickpoints in a number of longitudinal profiles and by general trends of concavity as a function of lithology. Volcanic and sedimentary rocks are the principal rock types underlying the northern Oregon Coast Ranges (between 46°30′ and 45°N) where mixed bedrock–alluvial channels dominate. Average concavity, θ, is 0·57 in this region. In the alluviated central Oregon Coast Ranges (between 45° and 44°N) values of concavity are, on average, the highest (θ = 0·82). South of 44°N, however, bedrock channels are common and θ = 0·73. Mixed bedrock–alluvial channels characterize rivers in the Klamath Mountains (from 43°N south; θ = 0·64). Rock uplift rates of ≥0·5 mm a⁻¹, mixed bedrock–alluvial channels, and concavities of 0·53–0·70 occur within the northernmost Coast Ranges and Klamath Mountains. For rivers flowing over volcanic rocks θ = 0·53, and θ = 0·72 for reaches crossing sedimentary rocks. Whereas channel type and concavity generally co-vary with lithology along much of the range, rivers between 44·5° and 43°N do not follow these trends. Concavities are generally greater than 0·70, alluvial channels are common, and river profiles lack knickpoints between 44·5° and 44°N, despite the fact that lithology is arguably invariant. Moreover, rock uplift rates in this region vary from low, ≤0·5 mm a⁻¹, to subsidence (<0 mm a⁻¹). These observations are consistent with models of transient river response to a decrease in uplift rate. Conversely, the rivers between 44° and 43°N have similar concavities and flow on the same mapped bedrock unit as the central region, but have bedrock channels and irregular longitudinal profiles, suggesting the river profiles reflect a transient response to an increase in uplift rate. If changes in rock uplift rate explain the differences in river profile form and morphology, it is unlikely that rock uplift and erosion are in steady state in the Oregon coastal mountains. Copyright © 2006 John Wiley & Sons, Ltd.     

OUTBURST FLOODS FROM GLACIER-DAMMED LAKES: THE EFFECT OF MODE OF LAKE DRAINAGE ON FLOOD MAGNITUDE

Published accounts of outburst floods from glacier-dammed lakes show that a significant number of such floods are associated not with drainage through a tunnel incised into the basal ice—the process generally assumed—but rather with ice-marginal drainage, mechanical failure of part of the ice dam, or both. Non-tunnel floods are strongly correlated with formation of an ice dam by a glacier advancing from a tributary drainage into either a main river valley or a pre-existing body of water (lake or fiord). For a given lake volume, non-tunnel floods tend to have significantly higher peak discharges than tunnel-drainage floods. Statistical analysis of data for floods associated with subglacial tunnels yields the following empirical relation between lake volume V and peak discharge Q_p : Q_p = 46V^0.66 (r² = 0.70), when Q_p is expressed in metres per second and V in millions of cubic metres. This updates the so-called Clague–Mathews relation. For non-tunnel floods, the analogous relation is Q_p = 1100V^0.44 (r² = 0.58). The latter relation is close to one found by Costa (1988) for failure of constructed earthen dams. This closeness is probably not coincidental but rather reflects similarities in modes of dam failure and lake drainage. We develop a simple physical model of the breach-widening process for non-tunnel floods, assuming that (1) the rate of breach widening is controlled by melting of the ice, (2) outflow from the lake is regulated by the hydraulic condition of critical flow where water enters the breach, and (3) the effect of lake temperature may be dealt with as done by Clarke (1982). Calculations based on the model simulate quite well outbursts from Lake George, Alaska. Dimensional analysis leads to two approximations of the form Q_p ∝ V^qf(h_i, θ₀), where q = 0.5 to 0.6, h_i is initial lake depth, θ₀ is lake temperature, and the form of f (h_i, θ₀) depends on the relative importance of viscous dissipation and the lake's thermal energy in determining the rate of breach opening. These expressions, along with the regression relations, should prove useful for assessing the probable magnitude of breach-type outburst floods.     

The effect of lithology on valley width,terrace distribution,and bedload provenance in a tectonically stable catchment with flat‐lying stratigraphy

How rock resistance or erodibility affects fluvial landforms and processes is an outstanding question in geomorphology that has recently garnered attention owing to the recognition that the erosion rates of bedrock channels largely set the pace of landscape evolution. In this work, we evaluate valley width, terrace distribution, and bedload provenance in terms of reach scale variation in lithology in the study reach and discuss the implications for landscape evolution in a catchment with relatively flat‐lying stratigraphy and very little uplift. A reach of the Buffalo National River in Arkansas was partitioned into lithologic reaches and the mechanical and chemical resistance of the main lithologies making up the catchment was measured. Valley width and the spatial distribution of terraces were compared among the different lithologic reaches. The surface grain size and provenance of coarse (2–90 mm) sediment of both modern gravel bars and older terrace deposits that make up the former bedload were measured and defined. The results demonstrate a strong impact of lithology upon valley width, terrace distribution, and bedload provenance and therefore, upon landscape evolution processes. Channel down‐cutting through different lithologies creates variable patterns of resistance across catchments and continents. Particularly in post‐tectonic and non‐tectonic landscapes, the variation in resistance that arises from the exhumation of different rocks in channel longitudinal profiles can impact local base levels, initiating responses that can be propagated through channel networks. The rate at which that response is transmitted through channels is potentially amplified and/or mitigated by differences between the resistance of channel beds and bedload sediment loads. In the study reach, variation in lithologic resistance influences the prevalence of lateral and vertical processes, thus producing a spatial pattern of terraces that reflects rock type rather than climate, regional base level change, or hydrologic variability. Copyright © 2017 John Wiley & Sons, Ltd.     

On the Causes of Frequency-Dependent Apparent Seismological Q

Variability of the Earth’s structure makes a first-order impact on attenuation measurements which often does not receive adequate attention. Geometrical spreading (GS) can be used as a simple measure of the effects of such structure. The traditional simplified GS compensation is insufficiently accurate for attenuation measurements, and the residual GS appears as biases in both Q ₀ and η parameters in the frequency-dependent attenuation law Q(f) = Q ₀ f ^η . A new interpretation approach bypassing Q(f) and using the attenuation coefficient χ(f) = γ + πf/Q _e(f) resolves this problem by directly measuring the residual GS, denoted γ, and effective attenuation, Q _e. The approach is illustrated by re-interpreting several published datasets, including nuclear-explosion and local-earthquake codas, Pn, and synthetic 50–300-s surface waves. Some of these examples were key to establishing the Q(f) concept. In all examples considered, χ(f) shows a linear dependence on the frequency, γ ≠ 0, and Q _e can be considered frequency-independent. Short-period crustal body waves are characterized by positive γ _SP values of (0.6–2.0) × 10^?2 s^?1 interpreted as related to the downward upper-crustal reflectivity. Long-period surface waves show negative γ _LP ≈ ?1.9 × 10^?5 s^?1, which could be caused by insufficient modeling accuracy at long periods. The above γ values also provide a simple explanation for the absorption band observed within the Earth. The band is interpreted as apparent and formed by levels of Q _e ≈ 1,100 within the crust decreasing to Q _e ≈ 120 within the uppermost mantle, with frequencies of its flanks corresponding to γ _LP and γ _SP. Therefore, the observed absorption band could be purely geometrical in nature, and relaxation or scattering models may not be necessary for explaining the observed apparent Q(f). Linearity of the attenuation coefficient suggests that at all periods, the attenuation of both Rayleigh and Love waves should be principally accumulated at the sub-crustal depths (~38–100 km).     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

Variation of intrinsic and scattering attenuation of seismic waves with depth in the Bam region,East-Central Iran

The attenuation properties of the lithosphere in the Bam region, East-Central Iran, have been investigated. For this purpose, 42 local earthquakes having focal depths less than 25 km have been used. The quality factor of coda waves (Q_c) has been estimated using the single back-scattering model. The quality factors Q_p, Q_d (P and direct S-waves) have been estimated using the extended-coda normalization method. Q_i and Q_s (the intrinsic and scattering attenuation parameters) have been estimated for the region. The values of Q_p, Q_d, Q_c, Q_i and Q_s show a dependence on frequency in the range of 1.5–24 Hz for the Bam region. The average frequency-dependent relationships estimated for the region are Q_p=(36±6)f^(1.03±0.06), Q_d=(59±8)f^(1.00±0.03), Q_c=(79±5)f^(1.01±0.04), Q_s=(131±4)f^(1.01±0.04) and Q_i=(104±6)f^(1.01±0.05). A comparison between Q_i and Q_s shows that intrinsic absorption is predominant over scattering.The variation of Q has also been estimated at different lapse times to observe heterogeneities variation with depth. The variation of Q with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth.The estimated Q_o values at different stations suggest a low value of Q indicating a heterogeneous and attenuative crust beneath the entire region.     

Sediment flux-driven channel geometry adjustment of bedrock and mixed gravel–bedrock rivers

Sediment supply (Q_s) is often overlooked in modelling studies of landscape evolution, despite sediment playing a key role in the physical processes that drive erosion and sedimentation in river channels. Here, we show the direct impact of the supply of coarse-grained, hard sediment on the geometry of bedrock channels from the Rangitikei River, New Zealand. Channels receiving a coarse bedload sediment supply are systematically (up to an order of magnitude) wider than channels with no bedload sediment input for a given discharge. We also present physical model experiments of a bedrock river channel with a fixed water discharge (1.5 l min⁻¹) under different Q_s (between 0 and 20 g l⁻¹) that allow the quantification of the role of sediment in setting the width and slope of channels and the distribution of shear stress within channels. The addition of bedload sediment increases the width, slope and width-to-depth ratio of the channels, and increasing sediment loads promote emerging complexity in channel morphology and shear stress distributions. Channels with low Q_s are characterized by simple in-channel morphologies with a uniform distribution of shear stress within the channel while channels with high Q_s are characterized by dynamic channels with multiple active threads and a non-uniform distribution of shear stress. We compare bedrock channel geometries from the Rangitikei and the experiments to alluvial channels and demonstrate that the behaviour is similar, with a transition from single-thread and uniform channels to multiple threads occurring when bedload sediment is present. In the experimental bedrock channels, this threshold Q_s is when the input sediment supply exceeds the transport capacity of the channel. Caution is required when using the channel geometry to reconstruct past environmental conditions or to invert for tectonic uplift rates, because multiple configurations of channel geometry can exist for a given discharge, solely due to input Q_s. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

An attenuation study using earthquakes from the 1997 Umbria-Marche sequence

We studied spatial and temporal characteristics of seismic attenuation inCentral Italy using S- and coda- waves recorded by the MarchesanSeismograph Network from earthquakes located in the epicentral area ofthe 1997 Umbria-Marche sequence. The amplitude decay of the S waveswith distance was defined calculating empirical attenuation functions at 15frequencies between 1 and 25 Hz. We analyzed separately foreshocks andaftershocks and we found the same attenuation functions, suggesting thatthe possible temporal variations could be confined in a small area. Thefrequency dependence of Q _S was approximated by the equation Q _S=18 · f ^2.0between 1 and 10 Hz. At higher frequencies (10–25 Hz), the frequencydependence of Q _s weakens, having an average value of Q _S=990. We also estimated Q from coda waves (Q _C) using the single-scattering models of Aki andChouet (1975) and Sato (1977). We found that Q _C=77 · f ^0.6, (between 2 and 20Hz) at the western side of the mountain chain, using either foreshocks oraftershocks. This relation is consistent with previous estimates of Q _Creported for the Central Apennines. For a volume sampling the Colfioritobasin, the Apennines and the Marche region we found that Q _C=55 · f ^0.8,indicating highattenuation below the mountain belt. To detect small temporal changes ofQ, we calculated spectral ratios of 5 temporal doublets located in theepicentral area and recorded at the closest station. We found temporalchanges of Q that vary from 27% to 56%, depending on the locationof the doublets. This variability suggests that the temporal change ofattenuation may depend on the spatial variation of Q and perhaps on thespatial distribution of tectonic stress in the epicentral area.     

Attenuation of Coda Waves in the Saurashtra Region, Gujarat (India)

The attenuation characteristics based on coda waves of two areas—Jamnagar and Junagarh of Saurashtra, Gujarat (India)—have been investigated in the present study. The frequency dependent relationships have been developed for both the areas using single back scattering model. The broadband waveforms of the vertical components of 33 earthquakes (Mw 1.5–3.5) recorded at six stations of the Jamnagar area, and broadband waveforms of 68 earthquakes (Mw 1.6–5) recorded at five stations of the Junagarh area have been used for the analysis. The estimated relations for the Junagarh area are: Q _c?=?(158?±?5)f^(0.99±0.04) (lapse time : 20?s), Q _c?=?(170?±?4.4)f^(0.97±0.02) (lapse time : 30?s) and Q _c?=?(229?±?6.6)f^(0.94±0.03) (lapse time : 40?s) and for the Jamnagar area are: Q _c?=?(178?±?3)f^(0.95±0.05) (lapse time : 20?s), Q _c?=?(224?±?6)f^(0.98±0.06) (lapse time : 30?s) and Q _c?=?(282?±?7)f^(0.91±0.03) (lapse time : 40?s). These are the first estimates for the areas under consideration. The Junagarh area appears to be more attenuative as compared to the Jamnagar area. The increase in Q _c values with lapse time found here for both the areas show the depth dependence of Q _c as longer lapse time windows will sample larger area. The rate of decay of attenuation (Q ^?1) with frequency for the relations obtained here is found to be comparable with those of other regions of the world though the absolute values differ. A comparison of the coda-Q estimated for the Saurashtra region with those of the nearby Kachchh region shows that the Saurashtra region is less heterogeneous. The obtained relations are expected to be useful for the estimation of source parameters of the earthquakes in the Saurashtra region of Gujarat where no such relations were available earlier. These relations are also important for the simulation of earthquake strong ground motions in the region.     

Scattering and intrinsic attenuation in Cairo metropolitan area using genetic algorithm

A total number of 46 local earthquakes (2.0≤ML≤4.0) recorded in the period 2000–2011 by the Egyptian seismographic network (ENSN) were used to estimate the total (Q_t⁻¹), intrinsic (Q_i⁻¹) and scattering attenuation (Q_sc⁻¹) in Cairo metropolitan area, Egypt. The multiple lapse time window analysis (MLTWA) under the assumption of multiple isotropic scattering with uniform distribution of scatters was firstly applied to estimate the pair of L_e⁻¹, the extinction length inverse, and B₀, the seismic albedo, in the frequency range 3–24 Hz. To take into account the effect of a depth-dependent earth model, the obtained values of B₀ and L_e⁻¹ were corrected for an earth structure characterized by a transparent upper mantle and a heterogeneous crust. The estimated values of Q_t⁻¹, Q_sc⁻¹ and Q_i⁻¹ exhibited frequency dependences. The average frequency-dependent relationships of attenuation characteristics estimated for the region are found to be: Q_t⁻¹=(0.015±0.008)f ^{(−1.02±0.02)}, Q_sc⁻¹=(0.006±0.001)f ^{(−1.01±0.02)}, and Q_i⁻¹=(0.009±0.008)f ^{(−1.03±0.02)}; showing a predominance of intrinsic absorption over scattering attenuation. This finding implies that the pore-fluid contents may have great effect on the attenuation mechanism in the upper crust where the River Nile is passing through the study area. The obtained results are comparable with those obtained in other tectonic regions.     

Geologic influences on fluvial hydrology and bedload transport in small mountainous watersheds,Northern New Mexico,U.S.A.

Discharge characteristics in six adjacent mountainous watersheds in northern New Mexico, U.S.A., vary substantially between basins underlain by different lithologies. Relatively resistant gneisses and granites underlie two basins (drainage areas: 43 and 94 km²) that have high unit discharge (0·010 to 0·14 m³s^?1 km^?2), high bankfull discharge, and sustained high discharge. Less resistant sandstones and shales underlie four basins (drainage areas: 96 to 215 km²) that have relatively low unit discharge (0·001 to 0·005 m³s^?1 km^?2), relatively low bankfull discharge, and peak discharges that are not sustained as long as those in the crystalline terrane. Analysis of snowmelt-runoff water budgets suggests that three factors control hydrologic conditions in the basins. First, area-elevation distributions appear to control the timing and amounts of water input. These distributions probably reflect the erosional resistance of the different lithologies. Second, lithology appears to control runoff production in areas having minor amounts of storage. Third, glacial deposits in headwater regions control discharge duration and timing via storage and return flow releases. The amount of return flow released by glacial deposits, however, is probably controlled by the permeability of underlying bedrock. Therefore it appears that the duration, timing, and magnitude of discharge events in the study area are controlled both directly and indirectly by lithology. Stream power and shear stress estimates derived from bankfull discharge and bed-material size data suggest that higher bedload transport rates and larger bedload particle sizes exist in streams draining crystalline rocks than in streams draining sedimentary terrane. It appears that source-area lithology, by controlling discharge production, also influences stream power, bedload transport capabilities, and therefore total amounts of bedload transport.     

Attenuation of P- and S-waves in the Chamoli Region, Himalaya, India

The attenuation properties of the crust in the Chamoli region of Himalaya have been examined by estimating the frequency-dependent relationships of quality factors for P waves (Q_α) and for S waves (Q_β) in the frequency range 1.5–24 Hz. The extended coda normalization method has been applied on the waveforms of 25 aftershocks of the 1999 Chamoli earthquake (M 6.4) recorded at five stations. The average value of Q_α is found to be varied from 68 at 1.5 Hz to 588 at 24 Hz while it varies from 126 at 1.5 Hz to 868 at 24 Hz for Q_β. The estimated frequency-dependent relations for quality factors are Q_α = (44 ± 1)f^(0.82±.04) and Q_β = (87 ± 3)f^(0.71±.03). The rate of increase of Q(f) for P and S waves in the Chamoli region is comparable with the other regions of the world. The ratio Q_β/Q_α is greater than one in the region which along with the frequency dependence of quality factors indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of attenuation relation for S wave estimated here (Q_β = 87f^0.71) with that of coda waves (Q_c = 30f^1.21) obtained by Mandal et al. (2001) for the same region shows that Q_c > Q_β for higher frequencies (>8 Hz) in the region. This indicates a possible high frequency coda enrichment which suggests that the scattering attenuation significantly influences the attenuation of S waves at frequencies >8 Hz. This observation may be further investigated using multiple scattering models. The attenuation relations for quality factors obtained here may be used for the estimation of source parameters and near-source simulation of earthquake ground motion of the earthquakes, which in turn are required for the assessment of seismic hazard in the region.