Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

Detailed Study of Time Variations in the Gutenberg–Richter <Emphasis Type="Italic">b</Emphasis>-value Based on Highly Accurate Seismic Observations at the Garm Prognostic Site,Tajikistan

The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log₁₀E^2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log₁₀E^2/3, and log₁₀N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log₁₀E^2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log₁₀N_i and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log₁₀N_i time series in the different magnitude ranges (M ≥ Mc_i) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔM_i) in the range of low magnitudes are known. It is concluded that using the parameter log₁₀N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.     

A first-order seismotectonic regionalization of Mexico for seismic hazard and risk estimation

The purpose of this work is to define a seismic regionalization of Mexico for seismic hazard and risk analyses. This seismic regionalization is based on seismic, geologic, and tectonic characteristics. To this end, a seismic catalog was compiled using the more reliable sources available. The catalog was made homogeneous in magnitude in order to avoid the differences in the way this parameter is reported by various agencies. Instead of using a linear regression to converts from m _b and M _d to M _s or M _w , using only events for which estimates of both magnitudes are available (i.e., paired data), we used the frequency-magnitude relations relying on the a and b values of the Gutenberg-Richter relation. The seismic regions are divided into three main categories: seismicity associated with the subduction process along the Pacific coast of Mexico, in-slab events within the down-going COC and RIV plates, and crustal seismicity associated to various geologic and tectonic regions. In total, 18 seismic regions were identified and delimited. For each, the a and b values of the Gutenberg-Richter relation were determined using a maximum likelihood estimation. The a and b parameters were repeatedly estimated as a function of time for each region, in order to confirm their reliability and stability. The recurrence times predicted by the resulting Gutenberg-Richter relations obtained are compared with the observed recurrence times of the larger events in each region of both historical and instrumental earthquakes.     

Study of seismicity in the NW Himalaya and adjoining regions using IMS network

The Reviewed Event Bulletin (REB) of the International Data Center (IDC) has been used in order to investigate the seismicity of the Northwest Himalaya and its neighboring region for the time period June 1999 to March 2015 within the geographical coordinates 25–40° N latitude and 65–85° E longitude. We have used a very precisely located earthquake dataset recorded by the International Monitoring System (IMS) Network containing 7,583 events with body wave magnitudes from 2.5 to 6.3. The study area has been subdivided into six regions based on the Flinn-Engdahl (F-E) seismic and geographical regionalization scheme, which was used as the region classifications of the International Data Center catalog. The examined region includes NW India, Pakistan, Nepal, Xizang, Kashmir, and Hindukush. For each region, Magnitudes of completeness (Mc) and Gutenberg-Richter (GR) recurrence parameters (a and b values) have been estimated. The Gutenberg-Richter analysis is preceded by an overview of the seismotectonics of the study area. The obtained Mc values vary from 3.5 to 3.9. The lower value of Mc was found mainly in Xizang region whereas the higher Mc threshold is evident in Pakistan region. However, the b values vary from 1.19 to 1.48. The lowest b value is recorded in Xizang region, which is mostly related to the Main Karakoram Thrust (MKT) fault, whereas the highest b values are recorded in NW India and Kashmir regions, which are mostly related to the Main Frontal Thrust (MFT) fault. The REB for the selected period has been compared to the most renowned bulletin of global seismicity, namely that issued by the National Earthquake Information Center (NEIC) of the United States Geological Survey (USGS). A study of 4,821 events recorded by USGS in the study region indicates that about 36 % of seismic events were missed and the catalog is considered as complete for events with magnitudes ≥4.0. However, both a and b values are obviously higher than those of IMS catalog. The a and b parameters in the Gutenberg-Richter magnitude–frequency relationship have been utilized to forecast the probability of future earthquakes of different magnitudes and returned periods (recurrence intervals).     

Precursor times of abnormal <Emphasis Type="Italic">b</Emphasis>-values prior to mainshocks

Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T _p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T _p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T _p , respectively, versus M are inferred from compiled data. The difference between the values of T and T _p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.     

An experimental study on variation of <Emphasis Type="Italic">b</Emphasis> value during deformation of rock samples with different structural models

A preliminary study of b value of rocks with two kinds of structural models has been made on the base of a new acoustic emission recording system. It shows that b value of the sample decreases obviously when the sample with compressive en echelon faults changes into a tensile one after interchange occurs between stress axis σ ₁ and σ ₂. A similar experiment is observed when the sample with tensile en echelon faults changes into that with a bend fault after two segments of the en echelon fault linking up. These facts indicate that the variation of b value may contain the information of the regional dominant structural model. Therefore, b-value analyses could be a new method for studying regional dominant structural models.     

Distribution of Lg coda <Emphasis Type="Italic">Q</Emphasis> in Xinjiang and its adjacent regions

In this study, we collected 1 156 broadband vertical components records at 22 digital seismic stations in Xinjiang region, Ürümqi station, and 7 stations in the adjacent regions during the period of 1999–2003. The records were firstly processed by the stacked spectral ratio method to obtain Q ₀ (Q at 1 Hz) and the frequency correlation factor η corresponding to each path. Based on the results, the distribution images of Q ₀ and η in 1°×1° grids for Xinjiang region were gained by the back-projection technique. The results indicate that Q ₀ is high (300–450) in the Tarim platform and marginal Siberian platform, while Q ₀ is low (150–250) in the southern regions as west Kunlun fold system and Songpan-Ganzi fold system. In the northern regions as Junggar fold system and Tianshan fold system, Q ₀ is also low (250–300) and η varies between 0.5 and 0.9.     

Classification of seismic events using linear discriminant function (LDF) in the Sakarya region,Turkey

The Sakarya prefecture is an interesting area with various seismicity types. This activity comes from earthquakes occurring at the North Anatolian Fault Zone and from a few quarry blast areas in the region. These quarry blast recordings produce errors in the determination of active faults and mapping of the microearthquake activity. Therefore, to recognize the tectonic activity in the region, we need to be able to discriminate between earthquakes and quarry blasts in the catalogues. In this study, a statistical analysis method (linear discriminant function) has been applied to classify seismic events occurring in the Sakarya region. We used 110 seismic events that were recorded by Sakarya University Seismic Station between 2012 and 2014. Time and frequency variant parameters, maximum S wave and maximum P wave amplitude ratio (S/P), the spectral ratio (Sr), maximum frequency (f_max), and total signal duration of the waveform were used for discrimination analyses. The maximum frequency (f_max) versus time duration of the seismic signal gives a higher classification percentage (94%) than the other discriminants. At the end of this study, 41 out of 110 events (44%) are determined as quarry blasts, and 62 (56%) are considered as earthquakes.     

Ionospheric Time-delay over Akure Using Global Positioning System Observations

Ionospheric time delay (VΔt) variability using Global Positioning System (GPS) data over Akure (7.15°N, 5.12°E), Nigeria, has been studied. The observed variability of VΔt in comparison to older results of vertical total electron content (TEC) across similar regions has shown equivalent signatures. Higher monthly mean values of VΔt (MVΔt) were observed during daytime as compared to nighttime (pre- and post-midnight) hours in all months. The highest MVΔt observed in September during daytime hours range between ~6 and ~21 ns (~1.80 and ~6.30 m) and at post-midnight, they are in the range of ~1 to ~6 ns (~0.3 to ~1.80 m). The possible mechanisms responsible for this variability were discussed. Seasonal VΔt were investigated as well.     

Influence of aquatic plants on the hydrogen isotope composition of sedimentary long-chain n-alkanes in the Lake Qinghai region,Qinghai-Tibet Plateau

The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.     

Effect of the IMF azimuthal component on the longitudinal position of the westward electrojet intensity maximum during strong geomagnetic storms

The distribution of the ionospheric currents during the geomagnetic storms of November 20–21, 2003, November 7–8, 2004, and November 9–10, 2004, depending on the IMF B _y component, has been studied based on the data from the global network of magnetic stations. It has been indicated that, during geomagnetic disturbances, the westward electrojet intensity maximum is localized in the evening sector at IMF B _y < 0 and in the morning sector at IMF B _y > 0. The region of the westward electrojet intensity maximum shifts to morning hours with increasing positive B _y values. Thus, the IMF azimuthal component forms not only the magnetospheric convection pattern during magnetic storms but is also responsible for the longitudinal position of ionospheric structures.     

Dynamics of the ionospheric electric currents and auroral luminosity boundaries during strong magnetic storms

The regularities in the southward drift of the ionospheric current centers and luminosity boundaries during strong magnetic storms of November 2003 and 2004 (with Dst ≈ ?400 and ?470 nT, respectively) are studied based on the global geomagnetic observations and TV measurements of auroras. It has been indicated that the eastward and westward electrojets in the dayside and nightside sectors simultaneously shift equatorward to minimal latitudes of Φ _min ^° ～53°–55°. It has been obtained that the Φ _min ^° latitude decreases with increasing negative values of Dst, IMF B _z component, and westward electric field strength in the solar wind. The dependence of the electrojet equatorward shift velocity (V _av) on the rate of IMF B _z variations (ΔB _z /Δt) has been determined. It is assumed that the electrojet dynamics along the meridian is caused by a change in the structure of the magnetosphere and electric fields in the solar wind and the Earth’s magnetosphere.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

Negative pressure effect on the electrical conductivity of San Carlos olivine and its implication to the electrical structure in the upper mantle

Pressure effect on the electrical conductivity of San Carlos olivine was investigated by the newly installed electrical conductivity measurement system at China University of Geosciences. Electrical conductivity of San Carlos olivine aggregates was measured up to 12 GPa and 1475 K using the Walker-type multi-anvil apparatus equipped with eight WC cubes as the second-stage anvils. The pressure generation against applied load for the experimental assemblage was examined by phase transition of Bi,quartz, forsterite under different P-T conditions. To check the data validity of this new system, electrical conductivities of the serpentinites and talc samples were measured. The results are consistent with the published data of the same samples. Electrical conductivity(σ) of the San Carlos olivine aggregates and temperature(T) satisfy the Arrhenian formula: σ=σ0exp[.(ΔE+PΔV)/kT].The pre-exponential factor(σ0), activation energy(ΔE) and activation volume(ΔV) yield value of 7.74 S/m, 0.85 eV and 0.94cm3/mol, respectively. Electrical conductivities of the San Carlos olivine aggregates decline with increasing pressure at same temperatures. The negative pressure effect can be interpreted by strain energy model of defect energy together with the lattice deformation. In addition, the electrical conductivity-depth 1-D profile of the upper mantle was constructed based on our results and some assumptions. The calculated profile is concordant with the geophysical observation at the depth of 180–350 km beneath Europe, which indicates that the upper mantle beneath Europe might be dry.     

Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes simulation

Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and physical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length z_0m and excess resistance κB^-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB^-1 at the two stations. Results indicate that z_0m has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of z_0m in the land surface model and that obtained from this study. κB^-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of z_0m and κB^-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.     

Probabilistic aftershock hazard analysis,two case studies in West and Northwest Iran

Aftershock hazard maps contain the essential information for search and rescue process, and re-occupation after a main-shock. Accordingly, the main purposes of this article are to study the aftershock decay parameters and to estimate the expected high-frequency ground motions (i.e., Peak Ground Acceleration (PGA)) for recent large earthquakes in the Iranian plateau. For this aim, the Ahar-Varzaghan doublet earthquake (August 11, 2012; M _N =6.5, M _N =6.3), and the Ilam (Murmuri) earthquake (August 18, 2014 ; M _N =6.2) have been selected. The earthquake catalogue has been collected based on the Gardner and Knopoff (Bull Seismol Soc Am 64(5), 1363-1367, 1974) temporal and spatial windowing technique. The magnitude of completeness and the seismicity parameters (a,??b) and the modified Omori law parameters (P,??K,??C) have been determined for these two earthquakes in the 14, 30, and 60 days after the mainshocks. Also, the temporal changes of parameters (a,??b,??P,??K,??C) have been studied. The aftershock hazard maps for the probability of exceedance (33%) have been computed in the time periods of 14, 30, and 60 days after the Ahar-Varzaghan and Ilam (Murmuri) earthquakes. For calculating the expected PGA of aftershocks, the regional and global ground motion prediction equations have been utilized. Amplification factor based on the site classes has also been implied in the calculation of PGA. These aftershock hazard maps show an agreement between the PGAs of large aftershocks and the forecasted PGAs. Also, the significant role of b parameter in the Ilam (Murmuri) probabilistic aftershock hazard maps has been investigated.     

Experimental Seismic Event-screening Criteria at the Prototype International Data Center

—?Experimental seismic event-screening capabilities are described, based on the difference of body-and surface-wave magnitudes (denoted as M _s :m _b ) and event depth. These capabilities have been implemented and tested at the prototype International Data Center (PIDC), based on recommendations by the IDC Technical Experts on Event Screening in June 1998. Screening scores are presented that indicate numerically the degree to which an event meets, or does not meet, the M _s :m _b and depth screening criteria. Seismic events are also categorized as onshore, offshore, or mixed, based on their 90% location error ellipses and an onshore/offshore grid with five-minute resolution, although this analysis is not used at this time to screen out events.¶Results are presented of applications to almost 42,000 events with m _b ?≥?3.5 in the PIDC Standard Event Bulletin (SEB) and to 121 underground nuclear explosions (UNE's) at the U.S. Nevada Test Site (NTS), the Semipalatinsk and Novaya Zemlya test sites in the Former Soviet Union, the Lop Nor test site in China, and the Indian, Pakistan, and French Polynesian test sites. The screening criteria appear to be quite conservative. None of the known UNE's are screened out, while about 41 percent of the presumed earthquakes in the SEB with m _b ?≥?3.5 are screened out. UNE's at the Lop Nor, Indian, and Pakistan test sites on 8 June 1996, 11 May 1998, and 28 May 1998, respectively, have among the lowest M _s :m _b scores of all events in the SEB.¶To assess the validity of the depth screening results, comparisons are presented of SEB depth solutions to those in other bulletins that are presumed to be reliable and independent. Using over 1600 events, the comparisons indicate that the SEB depth confidence intervals are consistent with or shallower than over 99.8 percent of the corresponding depth estimates in the other bulletins. Concluding remarks are provided regarding the performance of the experimental event-screening criteria, and plans for future improvements, based on recent recommendations by the IDC Technical Experts on Event Screening in May 1999.     

The spatiotemporal analysis of the minimum magnitude of completeness <Emphasis Type="Italic">Mc</Emphasis> and the Gutenberg–Richter law <Emphasis Type="Italic">b</Emphasis>-value parameter using the earthquake catalog of Greece

Spatiotemporal mapping the minimum magnitude of completeness Mc and b-value of the Gutenberg–Richter law is conducted for the earthquake catalog data of Greece. The data were recorded by the seismic network of the Institute of Geodynamics of the National Observatory of Athens (GINOA) in 1970–2010 and by the Hellenic Unified Seismic Network (HUSN) in 2011–2014. It is shown that with the beginning of the measurements at HUSN, the number of the recorded events more than quintupled. The magnitude of completeness Mc of the earthquake catalog for 1970–2010 varies within 2.7 to 3.5, whereas starting from April 2011 it decreases to 1.5–1.8 in the central part of the region and fluctuates around the average of 2.0 in the study region overall. The magnitude of completeness Mc and b-value for the catalogs of the earthquakes recorded by the old (GINOA) and new (HUSN) seismic networks are compared. It is hypothesized that the magnitude of completeness Mc may affect the b-value estimates. The spatial distribution of the b-value determined from the HUSN catalog data generally agrees with the main geotectonic features of the studied territory. It is shown that the b-value is below 1 in the zones of compression and is larger than or equal to 1 in the zones dominated by extension. The established depth dependence of the b-value is pretty much consistent with the hypothesis of a brittle–ductile transition zone existing in the Earth’s crust. It is assumed that the source depth of a strong earthquake can probably be estimated from the depth distribution of the b-value, which can be used for seismic hazard assessment.     

Correlation between foam-bubble size and drag coefficient in hurricane conditions

A recently proposed model of foam impact on the air–sea drag coefficient C _d has been employed for the estimation of the effective foam-bubble radius R _b variation with wind speed U₁₀ in hurricane conditions. The model relates C _d (U₁₀) with the effective roughness length Z _eff (U₁₀) represented as a sum of aerodynamic roughness lengths of the foam-free and foam-covered sea surfaces Z _w (U₁₀) and Z _f (U₁₀) weighted with the foam coverage coefficient α _f (U₁₀). This relation is treated for known phenomenological distributions C _d (U₁₀), Z _w (U₁₀), and α _f (U₁₀) at strong wind speeds as an inverse problem for the effective roughness parameter of foam-covered sea surface Z _f (U₁₀). The present study is aimed at the estimation of the effective roughness of the sea surface assuming that the measurement data for the effective drag coefficient are known. The effective foam-bubble size is found as a function of the wind speed.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.