Comparison of Ionospheric Parameters during Similar Geomagnetic Storms

The degree of closeness of ionospheric parameters during one magnetic storm and of the same parameters during another, similar, storm is estimated. Overall, four storms—two pairs of storms close in structure and appearance according to recording of the magnetic field Х-component—were analyzed. The examination was based on data from Sodankyla observatory (Finland). The f-graphs of the ionospheric vertical sounding, magnetometer data, and riometer data on absorption were used. The main results are as follows. The values of the critical frequencies foF2, foF1, and foE for different but similar magnetic storms differ insignificantly. In the daytime, the difference is on average 6% (from 0 to 11.1%) for all ionospheric layers. In the nighttime conditions, the difference for foF2 is 4%. The nighttime values of foEs differ on average by 20%. These estimates potentially make it possible to forecast ionospheric parameters for a particular storm.     

Gradual late stage deepening of Gega ice-dammed lake,Tsangpo gorge,southeastern Tibet,indicated by preliminary sedimentary rock magnetic properties

Gega lake, in southeastern Tibet, was formed by the blocking of Yarlung Tsangpo gorge by a glacier and is a well-known example of geomorphological damming. However, the evolution of the damming process at the site is still not understood in detail. Here, we use measurements of multiple magnetic parameters of the sediments from the Yusong (YS) 3 section, which is well-dated by optically stimulated luminescence, to provide a detailed history of the late stage of Gega dammed lake since 17.0 ka. Low-frequency field magnetic susceptibility (χ_lf) increases upwards gradually from 25 to 79?×?10^?8 m³/kg above 5.5 m, but other magnetic properties, such as frequency-dependent magnetic susceptibility (χ_fd % and χ_fd), susceptibility of anhysteretic remanent magnetization (χ_ARM), and saturation isothermal remanent magnetization (SIRM) did not show a similar degree of enhancement. The magnetic grain size indicators of χ_ARM/χ_lf, χ_ARM/SIRM, and χ_lf/SIRM all indicate a trend of increasing magnetic grain size from the 5.5 m to the top of section; however, the bulk sediment grain size decreases gradually within the same interval. The total organic carbon is very low (0.2–0.7%), and thus, it is unlikely that the sedimentary environment is sufficiently strongly reducing to lead to the dissolution of magnetic minerals. Therefore, we infer that the coarsening-upwards of the magnetic grain size, and the increasing magnetic susceptibility from 5.5 m, reflect the gradual preferential preservation of magnetic minerals, caused by the deepening of the lake. Thus, we interpret the magnetic record of the section as reflecting the gradual deepening of the late stage of Gega ice-dammed lake. Last, the abrupt disappearance of the lake may have resulted in an outburst flood.     

Seismic behavior and zoning of the sabkha soils in Jubail industrial city,Saudi Arabia

This study aimed at the micro-level seismic behavior and zoning of the saline sabkha strata in Jubail industrial area in Eastern Saudi Arabia. It encompasses the evaluation of the site-specific seismic response parameters and the liquefaction potential for various possible subsurface conditions under the probable seismic event(s). The approach to achieve the objectives of this study included the following: analysis of geologic, hydrologic, and geotechnical data of the area; performance of field and laboratory dynamic testing; and dynamic modeling and analysis of the subsurface profiles. The results of the simulation have been used to develop liquefaction potential maps and site-specific spectra of the study area, consisting of ten seismic zones under a range of probable peak horizontal ground acceleration (PHA). Results do not show significant probability of liquefaction of the loose soil layers in the study area at the maximum possible design PHA of 0.035 g; however, liquefaction is anticipated at higher PHA values. Site-specific spectral response resulted in values of S _s and S ₁ spectral accelerations to be different as compared to those suggested by local standards. The resulting seismic micro-zonation maps and the corresponding parameters are very useful for the stability analysis of the existing and planned structures in the Jubail area.     

Long-period irregular pulsations under the conditions of a quiet magnetosphere

Simultaneous observations of high-latitude long-period irregular pulsations at frequencies of 2.0–6.0 mHz (ipcl) and magnetic field disturbances in the solar wind plasma at low geomagnetic activity (Kp ~ 0) have been studied. The 1-s data on the magnetic field registration at Godhavn (GDH) high-latitude observatory and the 1-min data on the solar wind plasma and IMF parameters for 2011–2013 were used in an analysis. Ipcl (irregular pulsations continuous, long), which were observed against a background of the IMF Bz reorientation from northward to southward, have been analyzed. In this case other solar wind plasma and IMF parameters, such as velocity V, density n, solar wind dynamic pressure P = ρV² (ρ is plasma density), and strength magnitude B, were relatively stable. The effect of the IMF Bz variation rate on the ipcl spectral composition and intensity has been studied. It was established that the ipcl spectral density reaches its maximum (~10–20 min) after IMF Bz sign reversal in a predominant number of cases. It was detected that the ipcl average frequency (f) is linearly related to the IMF Bz variation rate (ΔBz/Δt). It was shown that the dependence of f on ΔBz/Δt is controlled by the α = arctan(By/Bx) angle value responsible for the MHD discontinuity type at the front boundary of magnetosphere. The results made it possible to assume that the formation of the observed ipcl spectrum, which is related to the IMF Bz reorientation, is caused by solar wind plasma turbulence, which promotes the development of current sheet instability and surface wave amplification at the magnetopause.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

Interpretation of magnetic anomalies using a genetic algorithm

A genetic algorithm (GA) is an artificial intelligence method used for optimization. We applied a GA to the inversion of magnetic anomalies over a thick dike. Inversion of nonlinear geophysical problems using a GA has advantages because it does not require model gradients or well-defined initial model parameters. The evolution process consists of selection, crossover, and mutation genetic operators that look for the best fit to the observed data and a solution consisting of plausible compact sources. The efficiency of a GA on both synthetic and real magnetic anomalies of dikes by estimating model parameters, such as depth to the top of the dike (H), the half-width of the dike (B), the distance from the origin to the reference point (D), the dip of the thick dike (δ), and the susceptibility contrast (k), has been shown. For the synthetic anomaly case, it has been considered for both noise-free and noisy magnetic data. In the real case, the vertical magnetic anomaly from the Pima copper mine in Arizona, USA, and the vertical magnetic anomaly in the Bayburt–Sar?han skarn zone in northeastern Turkey have been inverted and interpreted. We compared the estimated parameters with the results of conventional inversion methods used in previous studies. We can conclude that the GA method used in this study is a useful tool for evaluating magnetic anomalies for dike models.     

Seasonal variation of soil respiration under different land use/land cover in arid region

Through combining the soil respiration with the main environmental factors under the planting shelterbelt (Populus woodland) and the natural desert vegetation (Tamarix ramosissima+Phragmites communis community and Haloxylon ammodendron community) in the western Junngar Basin, the difference in soil respiration under different land use/land cover types and the responses of soil respiration to temperature and soil moisture were analyzed. Results showed that the rate of soil respiration increased with temperature. During the daytime, the maximum soil respiration rate occurred at 18:00 for the Populus woodland, 12:00 for T. ramosissima+Ph. communis community, and 14:00 for H. ammodendron community, while the minimum rate all occurred at 8:00. The soil respiration, with the maximum rate in June and July and then declining from August, exhibited a similar trend to the near-surface temperature from May to October. During the growing season, the mean soil respiration rates and seasonal variation differed among the land use/land cover types, and followed the order of Populus woodland > T. ramosissima+Ph. communis community > H. ammodendron community. The difference in the soil respiration rate among different land use/land cover types was significant. The soil respiration of Pouplus woodland was significantly correlated with the near-surface temperature and soil temperature at 10 cm depth (P<0.01) in an exponential manner. The soil respiration of T. ramosissima+Ph. communis and H. ammodendron communities were all linearly correlated with the near-surface temperature and soil surface temperature (P<0.01). Based on the near-surface temperature, the calculated Q₁₀ of Populus woodland, T. ramosissima+Ph. communis community and H. ammodendron community were 1.48, 1.59 and 1.63, respectively. The integrated soil respiration of the three land use/land cover types showed a significant correlation with the soil moisture at 0–5 cm, 5–15 cm and 0–15 cm depths (P<0.01). The quadratic model could best describe the relationship between soil respiration and soil moisture at 0–5 cm depth (P<0.01).     

Dependence of geomagnetic disturbances on extreme values of the solar wind <Emphasis Type="Italic">E</Emphasis><Subscript><Emphasis Type="Italic">y</Emphasis></Subscript> component

The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E _y component, dynamic pressure P _d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E _y value at E _y < 12 mV m^?1 and are inversely proportional to this value at E _y > 12 mV m^?1 for the first-type events. On the contrary, the dependence of Dst on E _y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E _y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.     

Variations in the IMF vertical component in isolated solar wind streams

The regularities of the variations in the IMF B _z component have been studied based on the data on the solar wind streams and their solar sources. Isolated solar wind streams such as magnetic clouds and shock layers before them, undisturbed heliospheric current sheets (HCSs), leading edges and bodies of high-speed streams from coronal holes (HSSs from CHs) have been considered. It has been revealed that each type of isolated streams in the interplanetary medium has it own features in the variations in the value and direction of the B _z component related to the stream immanent properties and conditions of propagation in the interplanetary plasma. The appearance of the southward B _z component is obligatory for all these streams which are, therefore, geoeffective.     

A different method for interpretation of magnetic anomalies due to 2-D dipping dikes

In this study a new method is presented to determine model parameters from magnetic anomalies caused by dipping dikes. The proposed method is applied by employing only the even component of the anomaly. First, the maximum of the even component is divided to its value at any distance x in order to obtain S1. Then, theoretical even component values are computed for the minimal depth (h) and half-width (b) values. S2 is obtained by dividing their maximum to the value computed for the same distance x. A set of S2 values is calculated by slowly increasing the half-width, and h and b for the S2 closest to S1 are determined. The same procedure is repeated by increasing the depth. The determined b values are plotted against the corresponding values of h. After repeating the process and plotting curves for different distances, it is possible to determine the actual depth and half-width values.     

Studying the Relationship between High-Latitude Geomagnetic Activity and Parameters of Interplanetary Magnetic Clouds with the Use of Artificial Neural Networks

The cause-and-effect relations of the dynamics of high-latitude geomagnetic activity (in terms of the AL index) and the type of the magnetic cloud of the solar wind are studied with the use of artificial neural networks. A recurrent neural network model has been created based on the search for the optimal physically coupled input and output parameters characterizing the action of a plasma flux belonging to a certain magnetic cloud type on the magnetosphere. It has been shown that, with IMF components as input parameters of neural networks with allowance for a 90-min prehistory, it is possible to retrieve the AL sequence with an accuracy to ~80%. The successful retrieval of the AL dynamics by the used data indicates the presence of a close nonlinear connection of the AL index with cloud parameters. The created neural network models can be applied with high efficiency to retrieve the AL index, both in periods of isolated magnetospheric substorms and in periods of the interaction between the Earth’s magnetosphere and magnetic clouds of different types. The developed model of AL index retrieval can be used to detect magnetic clouds.     

Effect of the Upper Part of the Soil Profile on the Site Response

The effect of the upper part of a layered profile on calculated seismic hazard parameters is analyzed. The behavior of the soil profile is studied as part of equivalent linear modeling with the example of six seismogeological models. Peak ground acceleration of the input signal varies from 0.1 to 0.6 g. It is shown that the character of the dependence of physicomechanical properties in the layer on the strain value leads to variations in the amplitudes of the response spectrum on the free surface and displacement of the maximum of the frequency spectrum. An increase in the thickness of the upper layer significantly influences the shape of the response spectra. Variations in the PGA values, damping coefficients, strain values with depth are calculated, as well as the velocity profiles for different soil layers with respect to the degradation curves. It is demonstrated that restricting the upper limit of the damping coefficient to 15% does not lead to noticeable changes in the shape of the response spectra on the free surface. The soil column begins to decrease the peak value of the input signal when its level exceeds 0.3 g.     

An ocean circulation model in <Emphasis Type="Italic">σ</Emphasis>S-<Emphasis Type="Italic">z</Emphasis>-<Emphasis Type="Italic">σ</Emphasis>B hybrid coordinate and its validation

A 3D, two-time-level, σS-z-σB hybrid-coordinate Marine Science and Numerical Modeling numerical ocean circulation model (HyMOM) is developed in this paper. In HyMOM, the σ coordinate is employed in the surface and bottom regions, and the z coordinate is used in the intermediate layers. This method can overcome problems with vanishing surface cells and minimize the unwanted deviation in representing bottom topography. The connection between the σ and z layers vertically includes an expanded “ghost” method and the linear interpolation. The governing equations in the σS-z-σB hybrid coordinate based on the complete Reynolds-averaged Navier-Stokes equations are derived in detail. The two-level time staggered and Eulerian forward and backward schemes, which are of second-order of accuracy, are adopted for the temporal difference in internal and external mode, respectively. The computation of the baroclinic gradient force is tested in an analytic test problem; the errors for two methods in HyMOM, which are relatively large only in the bottom layers, are obviously smaller than those in the pure σ and z models in almost all of the vertical layers. A quasi-global climatologic numerical experiment is constructed to test the simulation performance of HyMOM. With the monthly mean Levitus climatology data as reference, the HyMOM can improve the simulating accuracy compared with its pure z or σ coordinate implementation.     

Generation of superDreicer electric fields in the solar chromosphere

The electric field generation at the front of the current pulse, which originates in a coronal magnetic loop owing to the development of the Rayleigh–Taylor magnetic instability at loop footpoints, has been considered. During the τ_A ≈ l/V _A ≈ 5?25 s time (where l is the plasma plume height entering a magnetic loop as a result of the Rayleigh–Taylor instability), a disturbance related to the magnetic field tension B _?(r,t), “escapes” the instability region with the Alfvén velocity in this case. As a result, an electric current pulse I_z(z ? V _A t), at the front of which an induction magnetic field E _z, which is directed along the magnetic tube axis and can therefore accelerate particles, starts propagating along a magnetic loop with a characteristic scale of Δξ ≈ l. In the case of sufficiently large currents, when B _? ²/8π > p, an electric current pulse propagates nonlinearly, and a relatively large longitudinal electric field originates E _z ≈ 2I _z ³ V _A/c ⁴a²B_z ²l, which can be larger than the Dreicer field, depending on the electric current value.     

Many-year variations of water pollution and the state of river ecosystems in different latitudinal zones in European Russia

Many-year variations of water pollution and the state of more than two hundred of rivers (or their segments) in different latitudinal zones in European Russia are estimated. In accordance with latitudinal zonality, space and time trends in variations of water quality and the state of rivers are identified. By the level of water pollution, most river ecosystems (or their segments) were classified as polluted and very polluted, very polluted and dirty, dirty and very dirty. A tendency toward water quality improvement could be seen in the recent decade. More than 50% of the examined river ecosystems are at the transitional state of passing from natural (and/or equilibrium) to crisis state.     

A Least-squares Minimization Approach to Depth Determination from Magnetic Data

—We have developed a least-squares minimization approach to depth determination from magnetic data. By defining the anomaly value T(0) at the origin and the anomaly value T(N) at any other distance (N) on the profile, the problem of depth determination from magnetic data has been transformed into finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for a sphere, horizontal cylinder, dike, and for a geologic contact. Procedures are also formulated to estimate the effective magnetization intensity and the effective magnetization inclination. A scheme for analyzing the magnetic data has been formulated for determining the model parameters of the causative sources. The method is applied to synthetic data with and without random errors. Finally, the method is applied to two field examples from Canada and Arizona. In all cases examined, the estimated depths are found to be in goodagreement with actual values.     

Manifestations of cyclic variations in the solar magnetic field in long-term modulation of cosmic rays

The possibilities of improving the semiempirical model of cosmic ray (CR) modulation, proposed by us previously, are discussed. The following characteristics have been considered as model parameters in order to describe long-period CR variations using a unified model and to more completely reflect solar cycles in CR modulation as a complex interaction between two systems of fields (large-scale and local): the value and sign of the polar solar field, the average strength of the solar magnetic field (the B _ss integral index), partial indices (zone-even (ZE) and zone-odd (ZO) and sector-even (SE) and sector-odd (SO) indices), the tilt of the heliospheric current sheet, and the special index (F _x ) taking into account X ray flares. The role of each index in CR modulation has been revealed. When we described the long-term CR variations using many parameters and taking into account the integral index or one of four partial indices, the best results of modulation modeling during 1976–1999 were obtained for the B _ss total energetic index and SO index. A difference between the model calculations and observations increases beginning from the middle of 2000; the problem features of the CR behavior and the specific features of modeling this behavior in cycle 23 of solar activity (SA) are discussed. It is assumed that a decrease in the CR density at the last SA minimums (from cycle to cycle) can be related to a decrease in the ZO index and to a recently detected similar decrease in the vertical component of the solar dipole magnetic moment.     

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.     

Limits of transversely isotropic elastic parameters for the existence of classical Rayleigh waves

Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν _j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν _j is computed at each layer. For a given phase velocity (c), ν _j becomes complex depending on the transversely isotropic parameters. When ν _j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α _H , β _V , ξ, ? and η and find their limits beyond which ν _j becomes complex. It is seen that ν _j depends on ? and η, but does not depend on ξ. The complex ν _j occurs when ? is small and η is large. For a given c/β _V , the region of complex ν _j in a ? -η plane increases with the increase of α _H /β _V . Further, for a given α _H /β _V , the complex region of ν _j increases significantly with the decrease of c/β _V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.     

Experimental Study on Velocity Profiles with Different Roughness Elements in a Flume

The classical log law for velocity profile is applied to engineering practice. Field observations indicate that the composition of the bed materials obviously influences the shape of vertical velocity distribution. To clearly understand the roughness effect, six types of materials were laid separately at various depths for the investigation of the effects of roughness elements on the vertical velocity distribution. A down-looking 3D acoustic Doppler velocimeter was used to measure the velocity profiles. The experimental results showed that the curve characteristics of velocity profiles are strongly dependent on the roughness scale and related flow parameters. If d/R, Fr, and Re are larger than 0.15, 0.47, and 60 000, respectively, the velocity distribution may resemble an S-shape profile. The inflexion position Z*/H for a given S-shape profile was empirically deduced as Z*/H = ?0.4481d/R + 0.3225. Otherwise, the velocity profile agrees well with the logarithmic law. The findings of this study are useful in engineering practice (i.e., depth-averaged velocity and flow rate estimate). Key words: velocity profiles, acoustic Doppler velocimeter, roughness elements, logarithmic law, S-shape curve.