Intensity measures for the seismic response evaluation of buried steel pipelines under near-field pulse-like ground motions

Ground-motion Intensity Measures (IMs) are used to quantify the strength of ground motions and evaluate the response of structures. IMs act as a link between seismic demand and seismic hazard analysis and therefore, have a key role in performance-based earthquake engineering. Many studies have been carried out on the determination of suitable IMs in terms of efficiency, sufficiency and scaling robustness. The majority of these investigations focused on ordinary structures such as buildings and bridges, and only a few were about buried pipelines. In the current study, the optimal IMs for predicting the seismic demand of continuous buried steel pipelines under near-field pulse-like ground motion records is investigated. Incremental dynamic analysis is performed using twenty ground motion records. Using the results of the regression analysis, the optimality of 23 potential IMs are studied. It is concluded that specific energy density (SED) followed by \(\sqrt {VSI[{\omega _1}(PGD + RM{S_d})]} \) are the optimal IMs based on efficiency, sufficiency and scaling robustness for seismic response evaluation of buried pipelines under near-field ground motions.

     

Intensity measures for probabilistic assessment of non‐structural components acceleration demand

A fundamental issue in the framework of seismic probabilistic risk analysis is the choice of ground motion intensity measures (IMs). Based on the floor response spectrum method, the present contribution focuses on the ability of IMs to predict non‐structural components (NSCs) horizontal acceleration demand. A large panel of IMs is examined and a new IM, namely equipment relative average spectral acceleration (E‐ASA_R), is proposed for the purpose of NSCs acceleration demand prediction. The IMs efficiency and sufficiency comparisons are based on (i) the use of a large dataset of recorded earthquake ground motions; (ii) numerical analyses performed on three‐dimensional numerical models, representing actual structural wall and frame buildings; and (iii) systematic statistical analysis of the results. From the comparative study, the herein introduced E‐ASA_R shows high efficiency with respect to the estimation of maximum floor response spectra ordinates. Such efficiency is particularly remarkable in the case of structural wall buildings. Besides, the sufficiency and the simple formulation allowing the use of existing ground motion prediction models make the E‐ASA_R a promising IMs for seismic probabilistic risk assessment. Copyright © 2015 John Wiley & Sons, Ltd.     

A way of estimating the characteristic slip displacement

During the ruptures of an earthquake,the strain energy.△E,.will be transferred into,at least,three parts,i.e..the seismic radiation energy(E_s),fracture energy(E_g),and frictional energy(E_f),that is,△E = E_s + E_g + E_f.Friction,which is represented by a velocity- and state-dependent friction law by some researchers,controls the three parts.One of the main parameters of the law is the characteristic slip displacement.D_c.It is significant and necessary to evaluate the reliable value of D_c from observed and inverted seismic data.Since D_c controls the radiation efficiency.η_R = E_s/(E_s+ E_g),the value of η_r is a good constraint of estimating D_c.Integrating observed data and inverted results of source parameters from recorded seismograms.the values of E_s and E_g of an earthquake can be measured,thus leading to the value of η_R.The constraint used to estimate the reliable value of D_c will be described in this work.An example of estimates of D_c.based on the observed and inverted values of source parameters of the September 20,1999 M_S 7.6 Chi-Chi(Ji-Ji).Taiwan region,earthquake will be presented.     

Optimal intensity measures for probabilistic seismic demand modeling of extended pile-shaft-supported bridges in liquefied and laterally spreading ground

Seismic intensity measures (IMs) perform a pivotal role in probabilistic seismic demand modeling. Many studies investigated appropriate IMs for structures without considering soil liquefaction potential. In particular, optimal IMs for probabilistic seismic demand modeling of bridges in liquefied and laterally spreading ground are not comprehensively studied. In this paper, a coupled-bridge-soil-foundation model is adopted to perform an in-depth investigation of optimal IMs among 26 IMs found in the literature. Uncertainties in structural and geotechnical material properties and geometric parameters of bridges are considered in the model to produce comprehensive scenarios. Metrics such as efficiency, practicality, proficiency, sufficiency and hazard computability are assessed for different demand parameters. Moreover, an information theory based approach is adopted to evaluate the relative sufficiency among the studied IMs. Results indicate the superiority of velocity-related IMs compared to acceleration, displacement and time-related ones. In particular, Housner spectrum intensity (HI), spectral acceleration at 2.0 s (S _a-20), peak ground velocity (PGV), cumulative absolute velocity (CAV) and its modified version (CAV ₅) are the optimal IMs. Conversely, Arias intensity (I _a ) and shaking intensity rate (SIR) which are measures often used in liquefaction evaluation or related structural demand assessment demonstrate very low correlations with the demand parameters. Besides, the geometric parameters do not evidently affect the choice of optimal IMs. In addition, the information theory based sufficiency ranking of IMs shows an identical result to that with the correlation measure based on coefficient of determination (R ²). This means that R ² can be used to preliminarily assess the relative sufficiency of IMs.     

The viscosity of hydrous dacitic liquids: implications for the rheology of evolving silicic magmas

The viscosity of a series of six synthetic dacitic liquids, containing up to 5.04 wt% dissolved water, was measured above the glass transition range by parallel-plate viscometry. The temperature of the 10¹¹ Pa s isokom decreases from 1065 K for the anhydrous liquid, to 864 K and 680 K for water contents of 0.97 and 5.04 wt% H₂O. Including additional measurements at high temperatures by concentric-cylinder and falling-sphere viscometry, the viscosity (η) can be expressed as a function of temperature and water content w according to: where η is in Pa s, T is temperature in K, and w is in weight percent. Within the conditions of measurement, this parameterization reproduces the 76 viscosity data with a root-mean square deviation (RMSD) of 0.16 log units in viscosity, or 7.8 K in temperature. The measurements show that water decreases the viscosity of the dacitic liquids more than for andesitic liquids, but less than for rhyolites. At low temperatures and high water contents, andesitic liquids are more viscous than the dacitic liquids, which are in turn more viscous than rhyolitic liquids, reversing the trend seen for high temperatures and low water contents. This suggests that the relative viscosity of different melts depends on temperature and water content as much as on bulk melt composition and structure. At magmatic temperatures, rhyolites are orders of magnitude more viscous than dacites, which are slightly more viscous than andesites. During degassing, all three liquids undergo a rapid viscosity increase at low water contents, and both dacitic and andesitic liquids will degas more efficiently than rhyolitic liquids. During cooling and differentiation, changing melt chemistry, decreasing temperature and increasing crystal content all lead to increases in the viscosity of magma (melt plus crystals). Under closed system conditions, where melt water content can increase during crystallization, viscosity increases may be small. Conversely, viscosity increases are very abrupt during ascent and degassing-induced crystallization.     

A simplified analysis model for determining the seismic response of buried steel pipes at strike-slip fault crossings

The seismic response analysis of buried pipelines at fault crossings is a complex problem requiring nonlinear 3D soil-structure and large deformation analyses. Such analyses are computationally expensive and the results are hard to evaluate. Therefore, a simple numerical model is needed for engineering and design offices to determine the seismic demand of steel pipes at fault crossings. This paper presents a simplified numerical model for buried steel pipes crossing strike-slip faults and oriented perpendicular to the fault. Two pipes with different diameter to thickness (D/t) ratios and steel grades are used in the study. The proposed model permits plastic hinge formations in the pipe due to incrementally applied fault movements, allows determination of the critical length of the pipeline and measure strains developed on the tension and compression sides in the pipe. The model also considers the effect of bending as well as axial strains due to stretching.     

New magnitude scaling relations for earthquake early warning in the Alborz region,Iran

Rapid magnitude estimation relations for earthquake early warning systems in the Alborz region have been developed based on the initial first seconds of the P-wave arrival. For this purpose, a total of 717 accelerograms recorded by the Building and Housing Research Center in the Alborz region with the magnitude (M_w) range of 4.8–6.5 in the period between 1995 and 2013 were employed. Average ground motion period (\( \tau_{\text{c}} \)) and peak displacement (\( P_{\text{d}} \)) in different time windows from the P-wave arrival were calculated, and their relation with magnitude was examined. Four earthquakes that were excluded from the analysis process were used to validate the results, and the estimated magnitudes were found to be in good agreement with the observed ones. The results show that using the proposed relations for the Alborz region, earthquake magnitude could be estimated with acceptable accuracy even after 1 s of the P-wave arrival.     

Systematic retrieval of ejecta velocities and gas fluxes at Etna volcano using L-Band Doppler radar

Strombolian-type volcanic activity is characterized by a series of gas bubbles bursting at the top of a magma column and leading to the ejection of lava clots and gas emission at the surface. The quantitative analysis of physical parameters (e.g., velocity, size, and mass fluxes) controlling the emission dynamics of these volcanic products is very important for the understanding of eruption source mechanisms but remains difficult to obtain in a systematic fashion. Ground-based Doppler radar is found to be a very effective tool for measuring ejecta velocities at a high acquisition rate and close to the emission source. We present here a series of measurements carried out at Mt. Etna’s Southeast crater, using an L-band volcanological Doppler radar, during the 4 July 2001 Strombolian eruptions. Doppler radar data are supplemented by the analysis of video snapshots recorded simultaneously. We provide here a set of physical parameters systematically retrieved from 247 Strombolian explosions spanning 15 min and occurring during the paroxysm of the eruption from 21:30 to 21:45 UT. The time-average values give a maximum particle velocity of V_max^p = 94.7±24 \textm/s V_{{\max }}^p = {94}.{7}\pm {24} {\text{m/s}} , a bulk lava jet velocity of V_{\textPW - rad} = 37.6±1.9 \textm/s {V_{{{\text{PW - rad}}}}} = {37}.{6}\pm {1}.{9} {\text{m/s}} , and an initial gas velocity at the source vent of V₀^g = 118.4±36 \textm/s V_0^g = {118}.{4}\pm {36} {\text{m/s}} . The time-averaged particle diameter is found to be about D_{\textPW - rad} = 4.2±2.1 \textcm {D_{{{\text{PW - rad}}}}} = {4}.{2}\pm {2}.{1} {\text{cm}} . The volume and mass gas fluxes are estimated from time-averaged source gas velocities over the sequence duration at Q_v^g = 3 - 11 ×10³\textm³\text/s Q_v^g = {3} - {11} \times {1}{0^{{3}}}{{\text{m}}^{{3}}}{\text{/s}} and Q_m^g = 0.5 - 2 ×10³\textkg/s Q_m^g = 0.{5} - {2} \times {1}{0^{{3}}}{\text{kg/s}} , respectively.     

Average spectral acceleration as an intensity measure for collapse risk assessment

This paper investigates the performance of spectral acceleration averaged over a period range (Sa_avg) as an intensity measure (IM) for estimating the collapse risk of structures subjected to earthquake loading. The performance of Sa_avg is evaluated using the following criteria: efficiency, sufficiency, the availability or ease of developing probabilistic seismic hazard information in terms of the IM and the variability of collapse risk estimates produced by the IM. Comparisons are also made between Sa_avg and the more traditional IM: spectral acceleration at the first‐mode period of the structure (Sa(T₁)). Though most previous studies have evaluated IMs using a relatively limited set of structures, this paper considers nearly 700 moment‐resisting frame and shear wall structures of various heights to compare the efficiency and sufficiency of the IMs. The collapse risk estimates produced by Sa_avg and Sa(T₁) are also compared, and the variability of the risk estimates is evaluated when different ground motion sets are used to assess the structural response. The results of this paper suggest that Sa_avg, when computed using an appropriate period range, is generally more efficient, more likely to be sufficient and provides more stable collapse risk estimates than Sa(T₁). Copyright © 2015 John Wiley & Sons, Ltd.     

地面运动强度参数对RC桥梁地震需求的影响

在基于性能的地震工程学(PBEE)中,建立概率地震需求模型(PSDM)时需要对桥梁结构的工程需求参数(EDP)进行概率估计。其中,强地面运动参数(IM)的选择对EDP的概率估计影响很大,因此需要正确选择IM。分别采用目前最广泛使用的结构第一模态周期弹性谱加速度(5%阻尼比)Sa(T1,5%)和峰值地面加速度PGA作为IM,选择实际地震波并进行合理的调值,对一座钢筋混凝土桥墩进行IDA分析,其计算结果表明:对于不同性质EDP的概率估计值,以PGA作为IM计算所得的结果明显偏于非保守,且离散度一般也更大。说明可以针对不同性质的EDP,根据地面运动强度的大小,选择不同的IM,通过合理的调值对EDP进行概率估计,可以更加精确、高效地建立PSDM。     

Rocking amplification and strong‐motion duration

This paper characterizes the ability of natural ground motions to induce rocking demands on rigid structures. In particular, focusing on rocking blocks of different size and slenderness subjected to a large number of historic earthquake records, the study unveils the predominant importance of the strong‐motion duration to rocking amplification (ie, peak rocking response without overturning). It proposes original dimensionless intensity measures (IMs), which capture the total duration (or total impulse accordingly) of the time intervals during which the ground motion is capable of triggering rocking motion. The results show that the proposed duration‐based IMs outperform all other examined (intensity, frequency, duration, and/or energy‐based) scalar IMs in terms of both “efficiency” and “sufficiency.” Further, the pertinent probabilistic seismic demand models offer a prediction of the peak rocking demand, which is adequately “universal” and of satisfactory accuracy. Lastly, the analysis shows that an IM that “efficiently” captures rocking amplification is not necessarily an “efficient” IM for predicting rocking overturning, which is dominated by the velocity characteristics (eg, peak velocity) of the ground motion.     

Variation of Seismic Coda Wave Attenuation in the Garhwal Region, Northwestern Himalaya

Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q ₀ (Q _c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q _c ^?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.     

Fractional order intensity measures for probabilistic seismic demand modeling applied to highway bridges

Probabilistic seismic analysis of structures involves the construction of seismic demand models, often stated as probabilistic models of structural response conditioned on a seismic intensity measure. The uncertainty introduced by the model is often a result of the chosen intensity measure. This paper introduces the concept of using fractional order intensity measures (IMs) in probabilistic seismic demand analysis and uses a single frame integral concrete box‐girder bridge class and a seismically designed multispan continuous steel girder bridge class as case studies. The fractional order IMs considered include peak ground response and spectral accelerations at 0.2 and 1.0 s considering a single degree of freedom system with fractional damping, , as well as a linear single degree of freedom system with fractional response, . The study reveals the advantage of fractional order IMs relative to conventional IMs such as peak ground acceleration, peak ground velocity, or spectral acceleration at 0.2 and 1.0 s. Metrics such as efficiency, sufficiency, practicality, and proficiency are measured to assess the optimal nature of fractional order IMs. The results indicate that the proposed fractional order IMs produce significant improvements in efficiency and proficiency, whereas maintaining practicality and sufficiency, and thus providing superior demand models that can be used in probabilistic seismic demand analysis. Copyright © 2011 John Wiley & Sons, Ltd.     

The surface roughness and planetary boundary layer

Applications of the entrainment process to layers at the boundary, which meet the self similarity requirements of the logarithmic profile, have been studied. By accepting that turbulence has dominating scales related in scale length to the height above the surface, a layer structure is postulated wherein exchange is rapid enough to keep the layers internally uniform. The diffusion rate is then controlled by entrainment between layers. It has been shown that theoretical relationships derived on the basis of using a single layer of this type give quantitatively correct factors relating the turbulence, wind and shear stress for very rough surface conditions. For less rough surfaces, the surface boundary layer can be divided into several layers interacting by entrainment across each interface. This analysis leads to the following quantitatively correct formula compared to published measurements. 1 $$\begin{gathered} \frac{{\sigma _w }}{{u^* }} = \left( {\frac{2}{{9Aa}}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \left( {1 - 3^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} \frac{a}{k}\frac{{d_n }}{z}\frac{{\sigma _w }}{{u^* }}\frac{z}{L}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ = 1.28(1 - 0.945({{\sigma _w } \mathord{\left/ {\vphantom {{\sigma _w } {u^* }}} \right. \kern-\nulldelimiterspace} {u^* }})({z \mathord{\left/ {\vphantom {z L}} \right. \kern-\nulldelimiterspace} L})^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ \end{gathered} $$ where \(u^* = \left( {{\tau \mathord{\left/ {\vphantom {\tau \rho }} \right. \kern-0em} \rho }} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} \) , σ _w is the standard deviation of the vertical velocity,z is the height andL is the Obukhov scale lenght. The constantsa, A, k andd _n are the entrainment constant, the turbulence decay constant, Von Karman's constant, and the layer depth derived from the theory. Of these,a andA, are universal constants and not empirically determined for the boundary layer. Thus the turbulence needed for the plume model of convection, which resides above these layers and reaches to the inversion, is determined by the shear stress and the heat flux in the surface layers. This model applies to convection in cool air over a warm sea. The whole field is now determined except for the temperature of the air relative to the water, and the wind, which need a further parameter describing sea surface roughness. As a first stop to describing a surface where roughness elements of widely varying sizes are combined this paper shows how the surface roughness parameter,z ₀, can be calculated for an ideal case of a random distribution of vertical cylinders of the same height. To treat a water surface, with various sized waves, such an approach modified to treat the surface by the superposition of various sized roughness elements, is likely to be helpful. Such a theory is particularly desirable when such a surface is changing, as the ocean does when the wind varies. The formula, 2 $$\frac{{0.118}}{{a_s C_D }}< z_0< \frac{{0.463}}{{a_s C_D (u^* )}}$$ is the result derived here. It applies to cylinders of radius,r, and number,m, per unit boundary area, wherea _s =2rm, is the area of the roughness elements, per unit area perpendicular to the wind, per unit distance downwind. The drag coefficient of the cylinders isC _D . The smaller value ofz _o is for large Reynolds numbers where the larger scale turbulence at the surface dominates, and the drag coefficient is about constant. Here the flow between the cylinders is intermittent. When the Reynolds number is small enough then the intermittent nature of the turbulence is reduced and this results in the average velocity at each level determining the drag. In this second case the larger limit forz ₀ is more appropriate.     

The TKE dissipation rate in the northern South China Sea

The microstructure measurements taken during the summer seasons of 2009 and 2010 in the northern South China Sea (between 18°N and 22.5°N, and from the Luzon Strait to the eastern shelf of China) were used to estimate the averaged dissipation rate in the upper pycnocline 〈ε _p〉 of the deep basin and on the shelf. Linear correlation between 〈ε _p〉 and the estimates of available potential energy of internal waves, which was found for this data set, indicates an impact of energetic internal waves on spatial structure and temporal variability of 〈ε _p〉. On the shelf stations, the bottom boundary layer depth-integrated dissipation \( {\widehat{\varepsilon}}_{\mathrm{BBL}} \) reaches 17–19 mW/m², dominating the dissipation in the water column below the surface layer. In the pycnocline, the integrated dissipation \( {\widehat{\varepsilon}}_{\mathrm{p}} \) was mostly ∼10–30 % of \( {\widehat{\varepsilon}}_{\mathrm{BBL}} \). A weak dependence of bin-averaged dissipation \( \overline{\varepsilon} \) on the Richardson number was noted, according to \( \overline{\varepsilon}={\varepsilon}_0+\frac{\varepsilon_{\mathrm{m}}}{{\left(1+ Ri/R{i}_{\mathrm{cr}}\right)}^{1/2}} \), where ε ₀ + ε _m is the background value of \( \overline{\varepsilon} \) for weak stratification and Ri _cr = 0.25, pointing to the combined effects of shear instability of small-scale motions and the influence of larger-scale low frequency internal waves. The latter broadly agrees with the MacKinnon–Gregg scaling for internal-wave-induced turbulence dissipation.

     

Calibration of an M L scale for South Africa using tectonic earthquake data recorded by the South African National Seismograph Network: 2006 to 2009

A relation to determine local magnitude (M _L) based on the original Richter definition is empirically derived from synthetic Wood–Anderson seismograms recorded by the South African National Seismograph Network. In total, 263 earthquakes in the distance range 10 to 1,000 km, representing 1,681 trace amplitudes measured in nanometers from synthesized Wood–Anderson records on the vertical channel were considered to derive an attenuation relation appropriate for South Africa through multiple regression analysis. Additionally, station corrections were determined for 26 stations during the regression analysis resulting in values ranging between ?0.31 and 0.50. The most appropriate M _L scale for South Africa from this study satisfies the equation: $$ {M_{{{\bf L}}}} = {\text{lo}}{{\text{g}}_{{10}}}(A) + 1.149\;{\text{lo}}{{\text{g}}_{{10}}}(R) + 0.00063R + 2.04 - S $$ The anelastic attenuation term derived from this study indicates that ground motion attenuation is significantly different from Southern California but comparable with stable continental regions.     

A New Insight into Probabilistic Seismic Hazard Analysis for Central India

The Son-Narmada-Tapti lineament and its surroundings of Central India (CI) is the second most important tectonic regime following the converging margin along Himalayas-Myanmar-Andaman of the Indian sub-continent, which attracted several geoscientists to assess its seismic hazard potential. Our study area, a part of CI, is bounded between latitudes 18°–26°N and longitudes 73°–83°E, representing a stable part of Peninsular India. Past damaging moderate magnitude earthquakes as well as continuing microseismicity in the area provided enough data for seismological study. Our estimates based on regional Gutenberg–Richter relationship showed lower b values (i.e., between 0.68 and 0.76) from the average for the study area. The Probabilistic Seismic Hazard Analysis carried out over the area with a radius of ~300 km encircling Bhopal yielded a conspicuous relationship between earthquake return period (T) and peak ground acceleration (PGA). Analyses of T and PGA shows that PGA value at bedrock varies from 0.08 to 0.15 g for 10 % (T = 475 years) and 2 % (T = 2,475 years) probabilities exceeding 50 years, respectively. We establish the empirical relationships $ {\text{ZPA}}_{(T = 475)} = 0.1146\;[V_{\text{s}} (30)]^{ - 0.2924}, $ and $ {\text{ZPA}}_{(T = 2475)} = 0.2053\;[V_{\text{s}} (30)]^{ - 0.2426} $ between zero period acceleration (ZPA) and shear wave velocity up to a depth of 30 m [V _s (30)] for the two different return periods. These demonstrate that the ZPA values decrease with increasing shear wave velocity, suggesting a diagnostic indicator for designing the structures at a specific site of interest. The predictive designed response spectra generated at a site for periods up to 4.0 s at 10 and 2 % probability of exceedance of ground motion for 50 years can be used for designing duration dependent structures of variable vertical dimension. We infer that this concept of assimilating uniform hazard response spectra and predictive design at 10 and 2 % probability of exceedance in 50 years at 5 % damping at bedrocks of different categories may offer potential inputs for designing earthquake resistant structures of variable dimensions for the CI region under the National Earthquake Hazard Reduction Program for India.     

Ion Chemistry of the Ionosphere at E- and F-Region Altitudes: A Review

The current state of knowledge of E- and F-region ion chemistry is reviewed. Considerable attention is given to the progress in the chemistry of unexcited N₂ ⁺, O₂ ⁺, NO⁺, O⁺(⁴S), N⁺, H⁺, He⁺, Fe⁺, Mg⁺, Na⁺, Ca⁺, and K⁺ ions and electronically excited O⁺(²D), O⁺(²P), O⁺(⁴P), and $ {\text{O}}^{ + } (^{2} {\text{P}}^{*} ) $ ions. Achievements in our understanding of the role of vibrationally excited N₂ ⁺, O₂ ⁺, and NO⁺ ions in the ionosphere are discussed.