Biological and ecological characteristics ofTamarix L. and its effect on the ecological environment

Through studying biological and ecological characteristics ofTamarix L., we found wide adaptability to different ecological environment and high endurance of adversity. When planted in the edge of desert and periphery of oasis,Tamarix will act as excellent sand-fixing shrubs, improve ecological environment and play a positive role. However, introducingTamarix unreasonably will deteriorate the ecological environment and reduce biodiversity and play a negative role. Grasping the biological and ecological characteristics ofTamarix will do much to benefit further developing and utilizing of the resources

     

Water relation characteristics of<Emphasis Type="Italic">Alhagi sparsifolia</Emphasis> and consequences for a sustainable management

Water relation characteristics of the desert legumeAlhagi sparsifolia were investigated during the vegetation period from April to September 1999 in the foreland of Qira oasis at the southern fringe of the Taklamakan Desert, Xinjiang Uygur Autonomous Region of China. The seasonal variation of predawn water potentials and of diurnal water potential indicated thatAlhagi plants were well water supplied over the entire vegetation period. Decreasing values in the summer months were probably attributed to increasing temperatures and irradiation and therefore a higher evapotranspirative demand. Data from pressure-volume analysis confirmed thatAlhagi plants were not drought stressed and xylem sap flow measurements indicated thatAlhagi plants used large amounts of water during the summer months. Flood irrigation had no influence on water relations inAlhagi probably becauseAlhagi plants produced only few fine roots in the upper soil layers. The data indicate thatAlhagi sparsifolia is a drought-avoiding species that utilizes ground water by a deep roots system, which is the key characteristic to adjust the hyper-arid environment. Because growth and survival ofAlhagi depends on ground water supply, it is important that variations of ground water depth are kept to a minimum. The study will provide a theoretical basis for the restoration and management of natural vegetation around oasis in arid regions.     

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.     

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.     

A stepwise algorithm to detect significant time lags in ecological time series in terms of autocorrelation functions and ARMA model optimisation of pest population seasonal outbreaks

A central issue in population ecology is to determine the structure of negative feedback-density depend process which regulates population dynamics and seasonal fluctuations. In this work the incidence of population density dependences and seasonality was examined in fruit orchards of three closely related pest species (Adoxophyes orana, Anarsia lineatella and (Grapholita) Grapholitha molesta). Analysis included 13 moth population time series during 2003–2011. Additionally, considering that time lags and seasonality are fundamental characteristics of ecological organisation and pest management, the work aimed to introduce a step wise algorithm to detect significant population feedbacks, moth seasonality and population synchronisation of nearby locations. In the proposed procedure, each population-time series was first analysed on the basis of autocorrelation and partial autocorrelation. Moreover, assuming that each of the ecological variable, observed at successive time points, consist of a stochastic process, autoregressive moving average ARMA(p,q) models and seasonal autoregressive moving average models SARMA(p,q)x(P,Q) _S were fitted on data. The Akaike information criteria was further used by the stepwise algorithm for parameter optimization and model improvement. Model construction is accompanied by a presentation of the fitting results and a discussion of the heuristic benchmarks used to assess the forecasting performance of the models. Life cycles of populations belonging to same species appeared to synchronise by terms of their autocorrelation functions. Delayed density dependence and order was in most cases of lag:1 and 2, while lag >3 was not found more frequently as expected by chance. In A. orana and A. lineatella moth species lag = 1 delayed density dependence was significantly more frequent and in particular in nearby locations. However, the structure of the fitted models varied with respect to species and observation region. In some cases, seasonal models were considered to be more accurate in simulating moth population dynamics. Finally, to provide means in forecasting moth emergence and abundance, utile in pest management, the models were trained using 2003–2009 data sets and their forecasting performance were validated for each case using data sets of 2010–2011. In most cases, the constructed stochastic linear autoregressive models simulated the population outbreaks very well. Describing and forecasting stochastic population fluctuations is a basic tenet of theoretical and applied ecology, while detecting the relative roles of exogenous and endogenous mechanisms can partly describe the phenomenological behavior of pest population time series data and improve pest management.     

Carbon assimilation characteristics of plants in oasis-desert ecotone and their response to CO<Subscript>2</Subscript> enrichment

Six species of more than 20-year-old desert woody plants in the oasis-desert ecotone were selected for study. The results showed that: (1) in different growing seasons δ¹³C values of assimilating organ varied between -14%. and -16%. for Haloxylon ammodendron (HA),-14%. - -15%. for Calligonum mongolicum (CM) and -25%. - -28%. for Caragana korshinskii (CK), Nitraria sphaerocarpa (NS) and Hedysarum scoparium (HS). (2) The net photosynthetic rate (Pn) of HA and CM was significantly higher than those of the other species. With the decrease in Pn for the six species, their intercellular CO₂ concentration increased, but stomatal limitation value decreased under the intensive light. At the same time, the photochemical efficiency of PS II dropped to different degrees. (3) The CO₂ enrichment experiment demonstrated that, Pn of HA and CM increased to different extent under 450 μmol/mol, but their Pn reduced or approximated to the current condition under 650 μmol/mol. Under 450 μmol/mol the efficiency for solar energy utilization of CK and HS significantly reduced and under 650 μmol/mol their respiration rate exceeded photosynthesis rate. It can be concluded that HA and CM have some function of pathway for C₄ but the other three species have the function for C₃. The decline in their Pn is mainly caused by non-stomatal factors. HA, CM, CK and HS exhibited photoinhibition, which disappeared in a short time. This is a kind of positive readjustment to adapting to the desert environment. HA and CM can adapt to the high CO₂ environment, but CK and HS cannot. With the rise in atmospheric CO₂ concentration and climate warming, the latter two species in the oasis-desert ecotone may be gradually degraded or even disappear.     

Dependences of statistical characteristics of <Emphasis Type="Italic">NmE</Emphasis> on the month of the year at middle and low latitudes under daytime geomagnetically quiet conditions at low solar activity

Month-to-month changes in the statistical characteristics of the ionospheric E layer peak electron density NmE at medium and low geomagnetic latitudes under daytime geomagnetically quiet conditions are investigated. Critical frequencies of the ionospheric E layer measured by the middle latitude ionosonde Boulder and low latitude ionosondes Huancayo and Jicamarca at low solar activity from 1957 to 2015 have been used in the conducted statistical analysis. The mathematical expectation of NmE, standard deviation of NmE from the expectation of NmE, and NmE variation coefficient have been calculated for each month of the year. The months of the formation of extrema of these statistical parameters of NmE were found.     

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.     

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.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

<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.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Alcohol compounds in <Emphasis Type="Italic">Azolla imbricata</Emphasis> and potential source implication for marine sediments

This study investigated the composition of long-chain alkyl diols, triols, sec-alcohols, hydroxyl acids, and other hydroxylated compounds in Azolla imbricata and compared the organic alcohol components of Azolla filiculoides, Azolla microphylla, and South China Sea (SCS) sediments in order to investigate the possible indication of Azolla being the biological source of diols and triols in SCS sediment. Large amounts of diols, monohydroxy acids, and sec-alcohols with internal hydroxy groups at ω20 were detected in the three types of Azolla. Among these, 1,ω20-diol and ω20-hydroxy acid exhibited strong even-odd predominance distribution, whereas ω20-sec-alcohol exhibited strong odd-even predominance distribution. In addition, small amounts of diols, triols, and dihydroxy acids with internal hydroxy groups at 9, 10 or ω9, ω10 were detected, among which the chain length of C₂₉ was predominate. Compounds having similar structures as those in Azolla reflected a similar biosynthetic pathway: ω20-hydroxy acid exhibiting even-odd predominance distribution is decarboxylated to ω20-sec-alcohol exhibiting odd-even predominance distribution and converted to 1, ω20-diol with even-odd predominance distribution by acyl reduction; ω9, ω10-hydroxy acid is converted to 1,20,21(1, ω9, ω10)-triol by acyl reduction, and then converted to 9,10-diol by hydrogenation and dehydration. The alcohol components in A. imbricata were clearly not the biological source of 1,13/1,14/1,15-C_{28, 30, 32} diols and 1,3,4-C_27-29 triols in the SCS sediment. Certain marine diatoms might be the source of 1,14-C_{28, 30} diol in inshore sediment, but the biological source of diols and triols in the SCS sediment requires further investigation.     

Spatial and time variations in the <Emphasis Type="Italic">foF</Emphasis>2 (night)/<Emphasis Type="Italic">foF</Emphasis>2 (day) ratio: Specification of some effects

A series of important aspects of the time and spatial variations in the critical frequency ratio fo(night)/fo(day) considered by one of the authors [Danilov, 2007, 2008] is specified. The list of the Eastern-Hemisphere stations, for which an analysis of the above indicated ratio was performed, is completed. The available stations of the Western Hemisphere are considered. It is shown that the character of the variations in the fo(night)/fo(day) ratio is independent of limitations imposed on the Ap index of magnetic activity for the analyzed days. The consideration of the long-term variations in geomagnetic activity using the annual mean value of the Ap index does not influence the principal conclusions of the study, although in some cases changes statistical characteristics of the obtained behavior of the fo(night)/fo(day) ratio after 1980. A comparison of the results, based on the interpretation of the trends of the fo(night)/fo(day) value with the changes in the stratosphere-ionosphere coupling from the 1980s to the 1990s obtained by the authors earlier, confirms the assumption that there occurs a systematic change in the zonal wind in the upper atmosphere.     

Microbial D-amino acids and marine carbon storage

In nature, there are two conformational types of amino acids: L- and D-isomers. The L-amino acids are the predominant form and are used mainly for protein synthesis, while the D-amino acids are few in quantity but more diverse in terms of their biological functions. D-amino acids are produced by many marine microbes, which are important players in carbon and energy cycles in the ocean. As the major constituent of the marine organic carbon pool, D-amino acids can persist in the water column for a long time before being further transformed by chemical or biological processes or transported through physical processes (such as absorption and aggregation). This article reviews the microbial synthesis of D-amino acids, their physiological function and metabolism in microbes, and the contribution of D-amino acids as a carbon source to the oceanic carbon reservoir.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Modeling the distribution of <Emphasis Type="Italic">Populus euphratica</Emphasis> in the Heihe River Basin,an inland river basin in an arid region of China

Populus euphratica is a dominant tree species in riparian Tugai forests and forms a natural barrier that maintains the stability of local oases in arid inland river basins. Despite being critical information for local environmental protection and recovery, establishing the specific spatial distribution of P. euphratica has rarely been attempted via precise and reliable species distribution models in such areas. In this research, the potential geographic distribution of P. euphratica in the Heihe River Basin was simulated with MaxEnt software based on species occurrence data and 29 environmental variables. The result showed that in the Heihe River Basin, 820 km² of land primarily distributed along the banks of the lower reaches of the river is a suitable habitat for P. euphratica. We built other MaxEnt models based on different environmental variables and another eight models employing different mathematical algorithms based on the same 29 environmental variables to demonstrate the superiority of this method. MaxEnt based on 29 environmental variables performed the best among these models, as it precisely described the essential characteristics of the distribution of P. euphratica forest land. This study verified that MaxEnt can serve as an effective tool for species distribution in extremely arid regions with sufficient and reliable environmental variables. The results suggest that there may be a larger area of P. euphratica forest distribution in the study area and that ecological conservation and management of P. euphratica should prioritize suitable habitat. This research provides valuable insights for the conservation and management of degraded P. euphratica riparian forests.     

Differential physiological responses of the coastal cyanobacterium <Emphasis Type="Italic">Synechococcus</Emphasis> sp. PCC7002 to elevated <Emphasis Type="Italic">p</Emphasis>CO<Subscript>2</Subscript> at lag,exponential, and stationary growth phases

We studied the effects of expected end-of-the-century pCO₂ (1000 ppm) on the photosynthetic performance of a coastal marine cyanobacterium Synechococcus sp. PCC7002 during the lag, exponential, and stationary growth phases. Elevated pCO₂ significantly stimulated growth, and enhanced the maximum cell density during the stationary phase. Under ambient pCO₂ conditions, the lag phase lasted for 6 days, while elevated pCO₂ shortened the lag phase to two days and extended the exponential phase by four days. The elevated pCO₂ increased photosynthesis levels during the lag and exponential phases, but reduced them during the stationary phase. Moreover, the elevated pCO₂ reduced the saturated growth light (Ik) and increased the light utilization efficiency (α) during the exponential and stationary phases, and elevated the phycobilisome:chlorophyll a (Chl a) ratio. Furthermore, the elevated pCO₂ reduced the particulate organic carbon (POC):Chl a and particulate organic nitrogen (PON):Chl a ratios during the lag and stationary phases, but enhanced them during the exponential phase. Overall, Synechococcus showed differential physiological responses to elevated pCO₂ during different growth phases, thus providing insight into previous studies that focused on only the exponential phase, which may have biased the results relative to the effects of elevated pCO₂ in ecology or aquaculture.     

Solubilities of metallic sulfides and quartz in hydrothermal sulfide solutions

Thermal waters contain small amounts of dissolved sulfides which in places precipitate at or near the earth’s surface. Knowledge of the physical chemistry of hydrothermal solutions is needed at elevated temperatures and pressures. Therefore, solubilities of the epithermal minerals cinnabar, stibnite, quartz, and orpiment in aqueousNa ₂ S solutions was determined from 25–250° C, 1–2000 bars, and at severalNa ₂ S concentrations. All the minerals are appreciably soluble inNa ₂ S solutions. Pressure increase decreases solubilities of metallic sulfides but slightly increases quartz solubility. Temperature increase causes increased solubility at temperatures above 150° C, but at lower temperatures, cinnabar, orpiment and quartz show solubility decreases with increasing temperatures. Quartz and cinnabar are mutually soluble, but in the presence of stibnite only a small amount of cinnabar dissolves. The second ionization constant ofH ₂ S as calculated from the solubility data ranges from 10^?16,21 at 0°C to 10^?12,59 at 250°C. TheK ₂ ofH ₂ S is lower according to this study than any reported before, and the variation with temperature is several orders of magnitude greater than had been predicted.     

Long-term changes in the relation between <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2

The change in the dependence of the F2-layer critical frequency on its height hmF2 is considered based on two sources of initial data used earlier by the authors. It is found that the slope k of the foF2 dependence on hmF2 systematically decreases from the earlier (“etalon”) period, 1958–1980, to the later periods of 1988–2010, 1998–2010, and 1998–2014. Since the foF2 value depends on the atomic oxygen concentration in the F region much more strongly than hmF2, the found decrease in k confirms the concept of a decrease in the atomic oxygen concentration in the thermosphere with time previously formulated by the authors.