The electric field produced in the mantle by the dynamo in the core

A rigorous singular perturbation theory is developed to estimate the electric field E produced in the mantle M by the core dynamo when the electrical conductivity σ in M depends only on radius r, and when |r?_rln σ| ? 1 in most of M. It is assumed that σ has only one local minimum in M, either (a) at the Earth's surface ?V, or (b) at a radius b inside the mantle, or (c) at the core-mantle boundary ?K. In all three cases, the region where σ is no more than e times its minimum value constitutes a thin critical layer; in case (a), the radial electric field E_r ≈ 0 there, while in cases (b) and (c), E_r is very large there. Outside the critical layer, E_r ≈ 0 in all three cases. In no case is the tangential electric field E_S small, nearly toroidal, or nearly calculable from the magnetic vector potential A as ??_tA_S. The defects in Muth's (1979) argument which led him to contrary conclusions are identified. Benton (1979) cited Muth's work to argue that the core-fluid velocity u just below ?K can be estimated from measurements on ?V of the magnetic field B and its time derivative $\text{?}{}_{\mathrm{t}}\text{B}$. A simple model for westward drift is discussed which shows that Benton's conclusion is also wrong.In case (a), it is shown that knowledge of σ in M is unnecessary for estimating E_S on ?K with a relative error |r?_r 1nσ|^?1from measurements of E_S on ?V and knowledge of ?_tB in M (calculable from ?_tB on ?V if σ is small). Then, in case (a), u just below ?K can be estimated as $\text{?}\text{r}\text{×}\text{E}{}_{\mathrm{S}}\text{/B}{}_{\mathrm{r}}$. The method is impractical unless the contribution to E_S on ?V from ocean currents can be removed.The perturbation theory appropriate when σ in M is small is considered briefly; smallness of σ and of |r?_r ln σ|^?1 a independent questions. It is found that as σ → 0, B approaches the vacuum field in M but E does not; the explanation lies in the hydromagnetic approximation, which is certainly valid in M but fails as σ → 0. It is also found that the singular perturbation theory for |r?_r ln σ|^?1 is a useful tool in the perturbation calculations for σ when both σ and |r?_r ln σ|^?1 are small.     

Three-dimensional velocity structure around the focal area of the 2021 MS6.4 Yangbi earthquake

On May 21, 2021, an M_S6.4 earthquake occurred in Yangbi, Yunnan province, China, which exhibited typical foreshock-mainshock-aftershock characteristics. To better understand the velocity structure of the focal area and adjacent fault zones, Pg/Sg travel times at 12 seismic stations for the local earthquakes with M_L ≥ 1.5 from 2009–2019 and the Yangbi sequence in May of 2021 were used to invert the three-dimensional (3D) structures for both v_P and v_P/v_S. The obtained structure extends deeply to 15 km for area (25°N–26.5°N, 99.5°E–101°E) at a horizontal resolution of 10× 10 km, and the accuracy of the v_P velocity was verified using airgun signals excited by the Binchuan Airgun Transmitting Seismic Station (BATSS). The resulting v_P and v_P/v_S images correlate with existing fault zones and the Yangbi sequence, including: (1) The shallow velocity structure at 0 km agrees with local topography, where the Binchuan basin exhibits low-v_P and high-v_P/v_S values. From 3–15 km, v_P and v_P/v_S show variations, and the boundaries are consistent with the main faults (e.g., the Weixi-Qiaohou-Weishan, Honghe, and Chenghai faults). (2) The largest foreshock (M_S5.6), mainshock (M_S6.4), and largest aftershock (M_S5.2) occurred near the boundaries where both v_P and v_P/v_S have clear contrasts. (3) Small earthquakes are also concentrated in the transition zone between high- and low-v_P and v_P/v_S anomalies, and are biased toward low-v_P/v_S zones. (4) Boundaries in v_P and v_P/v_S are observed at 20 km west of the Weixi-Qiaohou-Weishan fault, indicating that there may exist one hidden fault.     

Influence of advection on the characteristics of turbulence over uneven surface in the oasis and the Gobi Desert

Utilizing experimental data of the atmospheric surface layer in the Gobi Oasis of Jinta in a comparative study, we demonstrate that under the condition of unstable stratification, the normalization variances of temperature in the oasis and Gobi Desert meet 1/3(z/Λ)(z/Λ)while normalization variances of both humidity and CO2 in the oasis meet(z)sz\Λ-1/3;s s z Λ the normalization variance of temperature in the oasis is large due to disturbance by advection, whereas variance of CO2 in the Gobi Desert has certain degree of deviation relative to Monin-Obukhov(M-O) scaling, and humidity variance completely deviates from variance M-O scaling. The above result indicates that under the condition of advection, humidity variance meets the relationsm sA sB D and it is determined by relative magnitude of scalar variance of advection transport. Our study reveals that, if the scalar variance of humidity or CO2 transported by advection is much larger than local scalar variance, observation value of scalar variance will deviate from M-O scaling; when scalar variance of advection transport is close to or less than local scalar variance, the observation value of scalar variance approximately meets M-O scaling.     

Possible causes of the variation in fractal dimension of the perimeter during the tropical cyclone Dan motion

We calculated the fractal dimensions Db of the perimeter of tropical cyclone(TC)Dan based on the satellite GMS-5 infrared sensor images from 1800 UTC,1 October 1999 to 1200 UTC,9 October 1999.The fractal dimensions Db were used to characterize objectively the temporal change of TC complex structure.Our results show that the change of fractal dimension during TC Dan motion can be divided into three stages.The statistically significant difference does not exist either between Dm1 and DL or between Dm3 and DL,but it exists between Dm2 and DL,where Dmi denotes the mean value of Db in i-th stage(i=1,2 and3);DL denotes Lovejoy’s fractal dimension calculated based on satellite and radar data within the size range(1–1.2×106 km2),which is used as a"normal value"of the fractal dimension of the cumulus cloud perimeter for the global tropical region.TC Dan turns to the north from the west abruptly at the end of the second stage.The emergence of the second stage with high fractal dimensions may be viewed as a possible premonition for the track turning.Our results also show that there are two kinds of processes resulting in the translation from the first stage to the second stage.One is the interaction of TC circulation and an adjacent small scale convective cloud cluster,causing to the complexity increase of a local segment of the perimeter.The other includes the fragmentation of a strong convective area within the TC inner region,the self-organization of the small strong convective cloud clusters,the emergence,development,and merger of the small scale non-convective holes,and the formation of a gap of the perimeter,causing to the complexity increase of the whole TC perimeter.     

Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

Stability of the boundary layer between the lithosphere and convecting mantle and the steady-state lithospheric geotherm

In the steady state, the convective boundary layer (CBL) (the transition from the lithosphere to the convecting mantle, the lithosphere-asthenosphere boundary) is on the verge of stability. This determines its depth, thickness, and the steady-state temperature distribution in the lithosphere. Had the mantle been homogeneous, the base of the lithosphere at the current potential temperature would lie globally at the same depth H _rh of 50 to 70 km. Actually, the regime of interaction of the mantle convection with the lithosphere is determined by the relationship between this depth and the thickness H _depl of the chemical boundary layer including the crust and the layer of the depleted rock. If the thickness of the chemical boundary layer is small H _depl < H _rh, as it is the case in the present-day oceanic mantle, the suboceanic regime is established with the mantle convection that does not reach the base of the chemical boundary layer. In this case, the top of CBL is located at depth H _rh, while the oceanic heat flow and the depth of the seafloor only depend on the potential temperature T _p and, within the areas where the crust is older than 60 to 70 Ma, are the same everywhere far from the disturbed territories (the hot points and the subduction zones). The absence of noticeable distinctions between the heat flow in the different oceanic basins suggests a global constancy of the potential temperature. If H _depl > H _rh, the subcontinental regime of the interaction of the mantle convection with the lithosphere is established. In this case, the CBL is immediately adjacent to the depleted lithosphere, its top is located at depth H _depl, and the surface heat flow q(T _p, H _depl) not only depends on the potential temperature T _p but also on the the thickness of the depleted lithosphere H _depl; it decreases with increasing H _depl and, therefore, with the age of the lithosphere. Given the potential temperature, the dependence q(T _p, H _depl) agrees well with the envelope of the results of kimberlite xenolith thermobarometry presented in the diagram of the deepest xenolith depth as a function of the heat flow. It is likely that in the lowest part of the continental lithosphere there is a zone of horizontal shear deformation, from where kimberlites entrain the strongly deformed and, at the same time, the deepest xenoliths. Besides, the azimuthal anisotropy of seismic velocities can be associated with this zone. The change in its direction with depth can be observed as the Lehmann discontinuity.     

Parameters influencing the growth and decay of the regular part of the AE index

The statistical behavior of the AE index, which describes the intensity of auroral electrojets, has been studied using hourly data. It has been shown that the AE index depends on the solar wind parameters as well as on its own value in the preceding hour. This makes it possible to describe the AE index behavior by the first-order inhomogeneous differential equation. The index growth rate is proportional mainly to the southward component of the interplanetary magnetic field (IMF). The timescale of the breakup of the currents responsible for the index is approximately 2 h. A similarity in the time variations in the AE and Dst indices has been referred to. The empirical formulas, which make it possible to restore the missing values of the IMF southward component from the available AE or Dst indices, have been proposed.     

The model of the depth variations in the properties and state of the rocks in the upper,middle, and lower continental crust of the Kola-Norwegian block,Kola Peninsula

Based on the experimental and calculated data, the model of depth variations in density, as well as the velocities of longitudinal (P) and transverse (S) waves, within the upper, middle, and lower crusts of the Kola-Norwegian block down to a depth of about 40 km is suggested. The variations in density and the velocities of the P- and S-waves are primarily caused by the changes in the mineral composition of the rocks. The relative reduction in the velocities of the P- and S-waves under the action of the increasing pressure and temperature in the depth interval from 5 to 37 km is estimated at ～2%.     

A Monte Carlo simulation of the effect of ion self-collisions on the ion velocity distribution function in the high-latitude F-region

Non-Maxwellian ion velocity distribution functions have been theoretically predicted and confirmed by observations, to occur at high latitudes. These distributions deviate from Maxwellian due to the combined effect of the E×B drift and ion-neutral collisions. The majority of previous literature, in which the effect of ion self-collisions was neglected, established a clear picture for the ion distribution under a wide range of conditions. At high altitudes and/or for solar maximum conditions, the ion-to-neutral density ratio increases and, hence, the role of ion self-collisions becomes appreciable. A Monte Carlo simulation was used to investigate the behavior of O⁺ ions that are E×B-drifting through a background of neutral O, with the effect of O⁺ (Coulomb) self-collisions included. Wide ranges of the ion-to-neutral density ratio n _i /n _n and the electrostatic field E were considered in order to investigate the change of ion behavior with solar cycle and with altitude. For low altitudes and/or solar minimum (n _i /n _n \leq10^?5), the effect of self-collisions is negligible. For higher values of n _i /n _n , the effect of self-collisions becomes significant and, hence, the non-Maxwellian features of the O⁺ distribution are reduced. For example, the parallel temperature T _i\Vert increases, the perpendicular temperature T _i⊥ decreases, the temperature anisotropy approaches unity and the toroidal features of the ion distribution function become less pronounced. Also, as E increases, the ion-neutral collision rate increases, while the ion-ion collision rate decreases. Therefore, the effect of ion self-collisions is reduced. Finally, the Monte Carlo results were compared to those that used simplified collision models in order to assess their validity. In general, the simple collision models tend to be more accurate for low E and for high n _i /n _n .     

Modelling the peak of the ionospheric E-layer

Calculations with a full time-varying model are used to study changes in the height and density of the E-layer peak, caused by known changes in the neutral atmosphere. Agreement with mean observed values of N_mE requires an increase of 10% in calculated ion densities, and an increase of 33% in the solar-maximum EUV model at λ<150 Å. At a fixed site, changes with the solar zenith angle χ agree well with the simple Chapman theory during most of the daylight hours. Simple modifications to the Chapman equations give improved accuracy near sunrise and sunset. When corrected for changes in χ, model results for summer and equinox show a decrease in the peak density N_mE at increasing latitudes. The overall change agrees well with experimental data, as summarised in the IRI model. Known changes in the neutral atmosphere also reproduce the increase in N_mE in winter, at latitudes up to 30°. The continuing increase at higher winter latitudes, in the IRI model, requires a major reduction in NO densities in winter. A suitable compromise is suggested. Equations fitted to the model results then provide a simpler and better behaved replacement for the IRI equations. Calculations at night show that known sources of ionisation, largely from starlight, can produce observed peak densities using current chemistry. There is an appreciable change with latitude, as starlight production increases in the southern hemisphere. The improbably large solar cycle change built into the IRI model, at night, cannot be reproduced and is not found in recent data. A new, simpler model is suggested. Changes in zenith angle and atmospheric composition cause the peak height (h_mE) to vary between 105 and 120 km, as a function of time, latitude, season and solar flux. These changes are approximated by simple equations that should be definitely preferable over the single, fixed height used in the IRI models.     

On the topographic coupling at the core-mantle interface

The mean tangential stresses at a corrugated interface between a solid, electrically insulating mantle and a liquid core of magnetic diffusivity λ are calculated for uniform rotation of both mantle and core at an angular velocity Ω in the presence of a corotating magnetic field B. The core and mantle are assumed to extend indefinitely in the horizontal plane. The interface has the form z = η(x, y), where z is the upward vertical distance and x, y are the zonal and latitudinal distances respectively. The function η(x, y) has a planetary horizontal length scale (i.e. of the order of the radius of the Earth) and small amplitude and vertical gradient. The liquid core flows with uniform mean zonal velocity U₀ relative to the mantle. Ω and B possess vertical and horizontal components.The vertical (poloidal) component B_p is uniform and has a value of 5 G while the horizontal (toroidal) field B_T = B_pαz, where α is a constant. When |α| ? 1, the mean horizontal stresses are found to have the same order of magnitude (10^?2 N m^?2) as those inferred from variations in the decade fluctuations in the length of the day, although the exact numerical values depend on the orientation of Ω as well as on the wavenumbers in the zonal and latitudinal directions.The influence of the steepness (as measured by α) of the toroidal field on the stresses is investigated to examine whether the constraint that the mean horizontal stresses at the core-mantle interface be of the order of 10^?2 N m^?2 might provide a selection mechanism for the behaviour of the toroidal field in the upper reaches of the outer core of the Earth. The results indicate that the restriction imposed on α is related to the value assigned to the toroidal field deep into the core. For example, if |α| ? 1 then the tangential stresses are of the right order of magnitude only if the toroidal field is comparable with the poloidal field deep in the core.     

Continuous-cyclic variations in the b-value of the earthquake frequency-magnitude distribution

Seismicity of the Earth (M ≥ 4.5) was compiled from NEIC, IRIS and ISC catalogues and used to compute b-value based on various time windows. It is found that continuous cyclic b-variations occur on both long and short time scales, the latter being of much higher value and sometimes in excess of 0.7 of the absolute b-value. These variations occur not only yearly or monthly, but also daily. Before the occurrence of large earthquakes, b-values start increasing with variable gradients that are affected by foreshocks. In some cases, the gradient is reduced to zero or to a negative value a few days before the earthquake occurrence. In general, calculated b-values attain maxima 1 day before large earthquakes and minima soon after their occurrence. Both linear regression and maximum likelihood methods give correlatable, but variable results. It is found that an expanding time window technique from a fixed starting point is more effective in the study of b-variations. The calculated b-variations for the whole Earth, its hemispheres, quadrants and the epicentral regions of some large earthquakes are of both local and regional character, which may indicate that in such cases, the geodynamic processes acting within a certain region have a much regional effect within the Earth. The b-variations have long been known to vary with a number of local and regional factors including tectonic stresses. The results reported here indicate that geotectonic stress remains the most significant factor that controls b-variations. It is found that for earthquakes with M _w ≥ 7, an increase of about 0.20 in the b-value implies a stress increase that will result in an earthquake with a magnitude one unit higher.     

Occurrence of Pyrodinium bahamense var. compressum along the southern coast of the Baja California Peninsula

As part of a continuing toxic microalgae monitoring program, 22 phytoplankton samples were collected from July to November 2010 at several sampling stations along the southern coast of the Baja California Peninsula. For the first time, the toxic dinoflagellate Pyrodinium bahamense var. compressum was found along the southeastern and southwestsern coasts of the peninsula. P. bahamense var. bahamense was first observed off San José del Cabo, which is an extension of the range of this variety. Both varieties occur as solitary cells. P. bahamense var. compressum occurred at temperatures ranging between 24.5 °C and 31 °C, whereas var. P.bahamense occurred at 28.5 °C to 29 °C, indicating its tropical and subtropical nature. Occurrence of P. bahamense var. compressum along this coastline may be related to El Niño 2009-2010.     

Influence of the Coriolis force on the formation of a seasonal thermocline

Large eddy simulation (LES) reveals that the Coriolis force plays an important role in seasonal thermocline formation. In the high-latitude ocean, a seasonal thermocline is formed at a certain depth, across which the downward transports of heat and momentum are prohibited. On the other hand, in the equatorial ocean, heat and momentum continue to propagate downward to the deeper ocean without forming a well-defined thermocline. Mechanism to clarify the latitudinal difference is suggested. The depth of a seasonal thermocline h is scaled in terms of both the Ekman length scale λ and the Monin–Obukhov length scale L, as h ??? 0.5(Lλ)^1/2, which is in contrast to the earlier suggestion as h?∝?L.     

LURR and the San Simeon M 6.5 Earthquake in 2003 and the Seismic Tendency in CA

The spatial and temporal variation of LURR (Load/Unload Respond Ratio) in California during April 2002 to June 2004 was studied in this paper. The result shows that before the San Simeon earthquake (35.7 N, 121.1 W) on Dec. 22, 2003, Y/Y _c anomalous region occurred successively near the epicenter from April 2002 to June 2002, and the maximum anomaly of Y/Y _c occurred in May, 2002. The published research work pointed out that the Y/Y _c anomaly near the San Simeon earthquake appeared from March, 2002. Compared with the five earthquake cases out of the six with M ≥ 6.5 in California during the period from 1980 to 2001, the maximum Y/Y _c and duration of Y/Y _c anomaly before this earthquake are among the normal ranges, but the time delay from the maximum anomaly time to the occurrence time of this earthquake is the longest one. The result also shows that two areas with Y/Y _c anomalies occurred from Oct. 2002 and Dec. 2002, respectively. According to statistical characteristics of the relationship between Y/Y _c anomalies and the coming earthquakes, the seismic tendency in California was discussed in this paper.     

Models of density and mechanical Q-factor distributions according to new data on the nutations and free oscillations of the Earth: 1. Ambiguity in the inverse problem

Ambiguity in the inverse problem of retrieval of the mechanical parameters of the Earth’s shell and core from the set of data on the velocities V _p and V _S , of longitudinal and transverse seismic body waves, the frequencies f _i and quality factors Q _i , of free oscillations, and the amplitudes and phases of forced nutation is considered. The numerical experiments show that the inverse problem of simultaneous retrieval of the density profile ρ in the mantle-liquid core system and the mechanical quality factor Q _μ of the mantle (if the total mass M and the total mean moment of inertia I of the Earth, and V _p and V _S are constant at all depths) has most unstable solutions. An example of depth distributions of ρ and Q _μ which are alternative to the well-known PREM model is given. In these distributions, the values of M and I and the velocities V _p and V _S at all depths for the period of oscillations T = 1 s exactly coincide with their counterparts yielded by PREM model (T = 1 s); however, the maximum deviations of the ρ and Q _μ profiles from those in the PREM model are about 3% and 40%, respectively; the mass and the moment of inertia of the liquid core are smaller than those for the PREM model by 0.75% and 0.63%, respectively. In this model, the root mean square (rms) deviations of all the measured values of f _i and Q _i from their values predicted by theory are half to third the corresponding values in the PREM model; the values of Δ for natural frequencies of the fundamental tone and overtones of radial oscillations, the fundamental tones of torsional oscillations, and the fundamental tones of spheroidal oscillations, which are measured with the highest relative accuracy, are smaller by a factor of 30, 6.6, and 2 than those in the PREM model, respectively. Such a large ambiguity in the solution of the inverse problem indicates that the current models of the depth distribution of density have relatively low accuracy, and the models of the depth distribution of the mechanical Q in the mantle are extremely unreliable. It is shown that the ambiguity in the models of depth distribution of density considerably decreases after the new data on the amplitudes and phases of the forced nutation of the Earth are taken into account. Using the same data, one may also refine by several times the recent estimates of the creep function for the lower mantle within a wide interval of periods ranging from a second to a day.     

Dielectric properties of the first 100 meters of the Moon

Reviewing 92 measurements of lunar sample dielectric constant versus density at frequencies above 100 kHz, gives the relationK′ = (1.93 ± 0.17)^p by regression analysis, where K′ is the dielectric constant of a soil or solid at a density ofpg/cm³. This formula is the geometric mean between the dielectric constant of vacuum (1) and the zero porosity dielectric constant of lunar material. Similarly, the loss tangent (D) can be described byD = [(0.00053 ± 0.00056) + (0.00025 ± 0.00009)C]p whereD is the loss tangent at densitypg/cm³ withC percent of total FeO + TiO₂ (approximately proportional to ilmenite content). Using the density versus depth relations derived from lunar surface core tubes, and from laboratory studies of lunar soil compression gives a model of the dielectric properties as a function of depth in the lunar regolith. The dielectric constant increases smoothly with depth, as a function of the soil compaction only. The loss tangent, however, is more sensitive to the ilmenite content than it is to density. Neither dielectric constant nor loss tangent varies significantly with the temperature observed in a lunar day.     

On the remote monitoring of Karenia brevis blooms of the west Florida shelf

In situ surveys (1997–2002) of Karenia brevis distribution on the west Florida shelf were used to explain spectral remote sensing reflectance, chlorophyll-a concentration, and backscattering coefficient estimates derived using SeaWiFS satellite data. Two existing approaches were tested in an attempt to differentiate K. brevis blooms from other blooms or plumes. A chlorophyll-anomaly method used operationally by the National Oceanic and Atmospheric Administration (NOAA) sometimes correctly identified K. brevis blooms but also generated false positives and false negatives. The method identified approximately 1000 km² of high chlorophyll-anomalies (>1 mg m⁻³) off southwest Florida between the 10 and 50-m isobaths nearly every day from summer to late fall. Whether these patches were K. brevis blooms or not is unknown. A second method used a backscattering:chlorophyll-a ratio to identify K. brevis patches. This method separated K. brevis from other blooms using in situ optical data, but it yielded less satisfactory results with SeaWiFS data. Spectral reflectance (R_rs) estimates for K. brevis blooms, diatom blooms, and coastal river plumes are statistically similar for many cases. Large pixel size, shallow water, and imperfect algorithms distort satellite retrievals of bio-optical parameters in patchy blooms. At present, a combination of chlorophyll-a, chlorophyll-anomaly, backscattering:chlorophyll-a ratio, RGB composites, MODIS fluorescence data, as well as time-series analysis and ancillary data such as winds, currents, and sea surface temperature can improve K. brevis bloom assessments. Progress in atmospheric correction and bio-optical inversion algorithms is required to help improve capabilities to monitor K. brevis blooms from space. Further, satellite sensors with improved radiometric capabilities and temporal/spatial resolutions are also required.