High-latitude Stations Research Articles

Multi-fractal Analysis of Cosmic Rays over Mid- and High-Latitude Stations During Severe Geomagnetic Storms

Multi-fractal Analysis of Cosmic Rays over Mid- and High-Latitude Stations During Severe Geomagnetic Storms

Whistler echo trains triggered by energetic winter lightning

Lightning generated electromagnetic impulses propagating in the magnetospheric plasma disperse into whistlers – several seconds long radio wave signals with decreasing frequency. Sometimes, multiple reflections form long echo trains containing many whistlers with increasing dispersion. On January 3, 2017, two necessary prerequisites – a pronounced lightning activity and a magnetospheric plasma duct – allowed for observations of a large number of whistler echo trains by the high-latitude station in Kannuslehto, Finland. Our investigation reveals that the duct existed for nearly eight hours. We show that causative lightning sferics arrived to the duct entry from three different winter thunderstorms: a small storm at the Norwegian coast, which produced energetic lightning capable to trigger echo trains in 50% of cases, and two large storms at unexpectedly distant locations in the Mediterranean region. Our results show that intense thunderstorms can repetitively feed electromagnetic energy into a magnetospheric duct and form whistler echo trains after subionospheric propagation over distances as large as 4000 km.

Application of RiTiCE in understanding hydro-meteorological controls on ice break-up patterns in River Tornionjoki

The Arctic region experiences significant annual hydrologic events, with the spring flood and ice break-up being the most prominent. River ice break-up, in particular, poses high socioeconomic and ecological expenses, including morphological changes and damage to riverine structures. This study aims to investigate the spatiotemporal patterns of river ice in the River Tornionjoki, including the timing of ice break-up at different latitudes. We utilized observation data and remote sensing techniques to track changes in ice patterns overtime on the River Tornionjoki. The study indicates that the ice break-up in the River Tornionjoki basin typically occurs during Apr-Jun based on the reach location in different latitudes; therefore, different stations behave according to their latitudinal location. We observed significant spatial variations in ice break-up timing across the basin, with an earlier break-up in the lower latitudes compared to the upper latitudes. The average ice break-up day in lower latitude stations ranges between 200–205, while in higher latitude stations the average ice break-up day ranges between 215–228.

Occurrence of Mesospheric Frontal Structures Over the High Latitude Station, Tromsø, Norway

Occurrence of Mesospheric Frontal Structures Over the High Latitude Station, Tromsø, Norway

Influence of ionospheric resonators on daily dynamics of the first Schumann resonance spectral parameters according to data from a meridional chain of ULF magnetometers

We have experimentally studied the influence of the local ionosphere, namely the ionospheric Alfvén resonator (IAR) and the lower ionospheric resonator at altitudes 80–300 km (sub-IAR) on the amplitude and polarization of the first Schumann resonance. The study is based on spectral analysis of data from simultaneous monitoring of ULF magnetic noise components at a meridional chain of stations: high-latitude stations Barentsburg and Lovozero, mid-latitude observatory NNGU NIRFI (NL, Nizhny Novgorod Region), a low-latitude station in Israel. We have also used monitoring data from the Borok and Crete observatories. At the stations in dark conditions, significant variations were found in the spectrum of the polarization parameter ε at the frequency of the first Schumann resonance (SR). Moreover, these variations had different character at different observatories. Analysis of the daily dynamics of the parameter ε has shown that these variations are associated with the influence of local sub-IAR having different optical thickness and quality factor at these observatories. The influence of sub-IAR on polarization in the SR band was found to depend on the ratio of the boundary frequency fb (the frequency that separates the negative and positive polarization of ULF magnetic noise) to the frequency of the first SR. The IAR influence on the polarization and amplitude of magnetic fields in the frequency band of the first Schumann resonance was discovered only at the NL and Lovozero stations: a high-quality Alfvén resonator in the ionosphere above the NL station could cause the SR frequency band to change and its central frequency to shift. Analysis of low-frequency data from the observatories separated by distances of 400 km has revealed that the influence of local ionospheric resonators can lead to a difference in the amplitude characteristics of the first SR even at such distances. It has also been shown that the influence of the IAR and sub-IAR resonators on the azimuthal angle of the magnetic field vector in the frequency band of the first SR is less noticeable and can generate variations in this parameter by 10°–20°. Numerical calculations performed for the spherical waveguide model made it possible to adequately interpret the features of the daily dynamics of the parameter ε in the frequency range of the first SR.

Влияние ионосферных резонаторов на суточную динамику спектральных параметров первого шумановс-кого резонанса по данным меридиональной сети УНЧ-магнито-метров

We have experimentally studied the influence of the local ionosphere, namely the ionospheric Alfvén resonator (IAR) and the lower ionospheric resonator at altitudes 80–300 km (sub-IAR) on the amplitude and polarization of the first Schumann resonance. The study is based on spectral analysis of data from simultaneous monitoring of ULF magnetic noise components at a meridional chain of stations: high-latitude stations Barentsburg and Lovozero, mid-latitude observatory NNGU NIRFI (NL, Nizhny Novgorod Region), a low-latitude station in Israel. We have also used monitoring data from the Borok and Crete observatories. At the stations in dark conditions, significant variations were found in the spectrum of the polarization parameter ε at the frequency of the first Schumann resonance (SR). Moreover, these variations had different character at different observatories. Analysis of the daily dynamics of the parameter ε has shown that these variations are associated with the influence of local sub-IAR having different optical thickness and quality factor at these observatories. The influence of sub-IAR on polarization in the SR band was found to depend on the ratio of the boundary frequency fb (the frequency that separates the negative and positive polarization of ULF magnetic noise) to the frequency of the first SR. The IAR influence on the polarization and amplitude of magnetic fields in the frequency band of the first Schumann resonance was discovered only at the NL and Lovozero stations: a high-quality Alfvén resonator in the ionosphere above the NL station could cause the SR frequency band to change and its central frequency to shift. Analysis of low-frequency data from the observatories separated by distances of 400 km has revealed that the influence of local ionospheric resonators can lead to a difference in the amplitude characteristics of the first SR even at such distances. It has also been shown that the influence of the IAR and sub-IAR resonators on the azimuthal angle of the magnetic field vector in the frequency band of the first SR is less noticeable and can generate variations in this parameter by 10°–20°. Numerical calculations performed for the spherical waveguide model made it possible to adequately interpret the features of the daily dynamics of the parameter ε in the frequency range of the first SR.

TEC disturbances caused by CME-triggered geomagnetic storm of September 6–9, 2017

This study investigates the ionospheric response to a geomagnetic storm triggered by a Coronal Mass Ejection (CME) during 6–9 September 2017, across GPS stations located in diverse geographical regions. We analyze the changes in the magnetic field component (ΔH), the Prompt Penetration Electric Fields (PPEF), and the Total Electron Content (TEC). We find that ΔH exhibits latitude-dependent responses during the storm, with high-latitude stations experiencing more significant reductions compared to low-latitude stations. The PPEF behavior is found to be directly correlated with solar wind disturbances. Particularly during the main phase of the storm, fluctuations in PPEF were clearly associated with negative values in the Dst index. The KIRU station, located at a high latitude, shows the most pronounced PPEF effects, indicating the increased susceptibility of high-latitude regions to solar wind interactions. The time series plot of TEC, covering a full month at different stations, shows a distinct diurnal pattern driven by solar ionization. Equatorial stations such as HYDE, BOU, HON (HNLC), and DODM exhibit the highest daily TEC values. During the geomagnetic storm, TEC disturbances are evident across all stations, with significant disturbances and varying trends in TEC depletion rate observed at different locations. The TEC values differ by 5–25 TECU during the storm period, suggesting intricate ionospheric responses to geomagnetic storms at different stations. This highlights the importance of considering different geographical regions to fully understand the ionospheric dynamics related to solar activities.

Hemispheric comparison of solar flare associated cosmic noise absorption (SCNA) from high latitude stations: Maitri (70.75°S, 11.75°E) and Abisko (68.4°N, 18.9°E)

The effects of solar flares in the D-region ionosphere at two high-latitude stations: Maitri (70.75°S, 11.75°E) and Abisko (68.4°N, 18.9°E), located in different hemispheres are studied. We analyzed 37 M-class flares and 6 X-class flares of the year 2014, which occurred when either or both stations were in the sun-lit side of the Earth. Cosmic Noise Absorption (CNA) curves are obtained using the datasets of riometers located at the two stations and are analyzed for all the 43 events under study. This paper discusses: 1) relationship between CNA and flare magnitude, 2) relationship between CNA and solar zenith angle (SZA), 3) hemispheric asymmetry in the observed solar flare associate CNA (or SCNA), and 4) the effect of background ionospheric condition in the SCNA magnitude at the two high latitude stations. It is observed that the solar flare effect in SCNA strongly depends on the SZA and flare intensity. Our analysis reveals that the flare response in SCNA for the year 2014 was stronger at Abisko than at Maitri. There is an observed hemispheric asymmetry in the solar flare ionization at D-region ionosphere for the given latitude. This asymmetry can be attributed to the previously enhanced background ionospheric ionization during particle precipitation processes. This study shows the need to establish an empirical relationship between the observed CNA vs. flare intensity, SZA and latitudinal position; especially when we go higher in latitudes.

Investigating sources of variability in closing the terrestrial water balance with remote sensing

Abstract. Remote sensing (RS) data are becoming an increasingly important source of information for water resource management as they provide spatially distributed data on water availability and use. However, in order to guide appropriate use of the data, it is important to understand the impact of the uncertainties of RS data on water resource studies. Previous studies have shown that the degree of closure of the water balance from remote sensing data is highly variable across basins and that different RS products vary in their levels of accuracy depending on climatological and geographical conditions. In this paper, we analyzed the water-balance-derived runoff from global RS products for 931 catchments across the globe. We compared time series of runoff estimated through a simplified water balance equation using three precipitation (CHIRPS, GPM, and TRMM), five evapotranspiration (MODIS, SSEBop, GLEAM, CMRSET, and SEBS), and three water storage change (GRACE-CSR, GRACE-JPL, and GRACE-GFZ) RS datasets with monthly in situ discharge data for the period 2003–2016. Results were analyzed through the lens of 10 quantifiable catchment characteristics in order to investigate correlations between catchment characteristics and the quality of RS-based water balance estimates of runoff and whether specific products performed better than others under certain conditions. The median Nash–Sutcliffe efficiency (NSE) for all gauges and all product combinations was −0.02, and only 44.9 % of the time series reached a positive NSE. A positive NSE could be obtained for 73.7 % of stations with at least one product combination, while the overall best-performing product combination was positive for 58.4 % of stations. This confirms previous findings that the best-performing products cannot be globally established. When investigating the results by catchment characteristic, all combinations tended to show similar correlations between catchment characteristics and the quality of estimated runoff, with the exception of combinations using MODIS evapotranspiration, for which the correlation was frequently reversed. The combinations with the GPM precipitation product generally performed worse than the CHIRPS and TRMM data. However, this can be attributed to the fact that the GPM data are available at higher latitudes compared to the other products, where performance is generally poorer. When removing high-latitude stations, this difference was eliminated, and GPM and TRMM showed similar performance. The results show the highest positive correlation between highly seasonal rainfall and runoff NSE. On the other hand, increasing snow cover, altitude, and latitude decreased the ability of the RS products to close the water balance. The catchment's dominant climate zone was also found to be correlated with time series performance, with the tropical areas providing the highest (median NSE = 0.11) and arid areas the lowest (median NSE = −0.09) NSE values. No correlation was found between catchment area and runoff NSE. The results highlight the importance of further studies on the uncertainties of the different data products and how these interact when combining them, as well as of new approaches to using the data rather than simple water-balance-type approaches. Efforts to improve specific satellite products can also be better targeted using the results of this study.

A revised BDS Klobuchar model with a non-symmetrical processing strategy and a high-latitude amplitude constraint

The BeiDou Navigation Satellite System (BDS) broadcasts the BDS Klobuchar model (BDSK) for ionospheric delay correction in both BDS-2 regional service phase and BDS-3 global service phase. Unlike the GPS Klobuchar model, the BDSK model adopts a simple “symmetrical strategy” due to the limited ionospheric input data within China, which results in inaccurate total electron content (TEC) estimation in the southern hemisphere. Besides, the “high latitude anomaly” occurs in the BDSK model due to the lack of latitude limitation strategy, which deteriorates the accuracy of the BDSK model in the southern hemisphere further. To solve this problem, we propose a revised BDS Klobuchar model (RBDK), in which the “non-symmetrical strategy” is adopted in the southern hemisphere and an amplitude constraint is introduced at high latitudes. For comparison, the fitted BDS Klobuchar model (BDFK) that is estimated using the same ionospheric input data and constraint algorithm for the RBDK model but with the “symmetrical strategy” is also included. First, we conduct a comparative analysis of the time series of the eight parameters for the BDSK, BDFK and RBDK models in 2021. The parameters of the BDFK and RBDK models are more stable, and there will be no abnormal jumping phenomenon that occurs in the BDSK model. Then, the performance of the broadcast BDSK model, the fitted BDFK model and the proposed RBDK model during the whole year of 2021 are evaluated by comparing them with the high-precision global ionospheric map (GIM) of the International GNSS Service (IGS). In the southern hemisphere, the root mean square (RMS) values of the BDSK and BDFK are 6.5 and 5.8 TEC unit (TECU), respectively, while that of the RBDK is 4.6 TECU, showing improvements of 29% and 21%, respectively. At high latitudes in the northern hemisphere, the RMS value of the BDSK is 7.6 TECU, while they are only 3.7 TECU for the RBDK and BDFK when the amplitude constraint is introduced. Final, we select 16 IGS stations to validate the performance of the BDSK, BDFK and RBDK models in the standard point positioning (SPP) domain in 2021. The average SPP accuracy of the BDSK is 8.88 m at 6 high-latitude stations, while they are 5.48 m for the RBDK and BDFK, showing significant improvement of 38%. Excluding the high-latitude stations, the SPP accuracy of the RBDK outperforms the BDFK and BDSK by 14% and 18%, respectively, at stations in the southern hemisphere. As expected, the proposed RBDK model shows significantly better performance in the southern hemisphere and at high latitudes in contrast to the broadcast BDSK model.

Impact of geomagnetic storms on ionospheric TEC at high latitude stations: a comparative analysis of GPS observations and the IRI-2016 model

Impact of geomagnetic storms on ionospheric TEC at high latitude stations: a comparative analysis of GPS observations and the IRI-2016 model

Response of the ionosphere to equatorial electrojet at latitudes near the northern EIA crest in the East African longitudinal sector

Equatorial Electrojet (EEJ) and Equatorial Ionization Anomaly (EIA) are two distinct ionospheric phenomena which are caused by the horizontal configuration of the geomagnetic field at the equator. On the day side of the earth, these activities are controlled by a common zonal electric field, and their interactions must be researched for a number of scientific reasons. In this study, we used multiple Global Navigation Satellite System (GNSS) stations at Al Wajh, Sola village, Nama, and Sheba in 2012 to investigate the correlations between the instantaneous and integrated equatorial electrojet (EEJ) and the variations of vertical total electron content (vTEC) at the northern crest of EIA in the East Africa and Middle East longitudinal sector. The equinoctial and solstitial months have the highest and lowest values of the monthly mean EEJ and vTEC, respectively, and as a result, both parameters show a semi-annual variability. The monthly mean vTEC also demonstrates equinoctial and solstitial asymmetries, resulting in vTEC values that are higher at the September equinox and the December solstice than at the March equinox and the June solstice, respectively. These asymmetries have been also observed in the composition of [O]/[N2] ratio. On days with strong counter electrojet (CEJ), EEJ peaks are smaller and occur later than on days with weak CEJ. Additionally, when CEJ is high, EIA is suppressed and vTEC exhibits low values at all crest points. In the early hours (10:00–12:00 Lt), Nama and Sheb showed higher correlation coefficient values with significant levels than Alwj and Sola. According to the correlation analysis, the anomaly crest develops at the lower latitude stations between 10:00 and 12:00 Lt and then moves to the higher latitude stations later. Generally, the correlations are higher and more significant in the equinoxes than the solstices. The correlations between integrated EEJ (IEEJ) and vTEC were also found to be strong and significant at the lower latitude stations throughout all seasons. These findings suggest that both the instantaneous and integrated EEJ can be used as a proxy for the eastward zonal equatorial electric field and must be considered in regional and global models of total electron content of equatorial and low latitude ionosphere.

Estimating BDS-3 Satellite Differential Code Biases with the Single-Frequency Uncombined PPP Model.

Differential Code Bias (DCB) is a crucially systematic error in satellite positioning and ionospheric modeling. This study aims to estimate the BeiDou-3 global navigation satellite system (BDS-3) satellite DCBs by using the single-frequency (SF) uncombined Precise Point Positioning (PPP) model. The experiment utilized BDS-3 B1 observations collected from 25 International GNSS Service (IGS) stations located at various latitudes during March 2023. The results reveal that the accuracy of estimating B1I-B3I DCBs derived from single receiver exhibits latitude dependence. Stations in low-latitude regions show considerable variability in the root mean square (RMS) of absolute offsets for satellite DCBs estimation, covering a wide range of values. In contrast, mid- to high-latitude stations demonstrate a more consistent pattern with relatively stable RMS values. Moreover, it has been observed that the stations situated in the Northern Hemisphere display a higher level of consistency in the RMS values when compared to those in the Southern Hemisphere. When incorporating estimates from all 25 stations, the RMS of the absolute offsets in satellite DCBs estimation consistently remained below 0.8 ns. Notably, after excluding 8 low-latitude stations and utilizing data from the remaining 17 stations, the RMS of absolute offsets in satellite DCBs estimation decreased to below 0.63 ns. These enhancements underscore the importance of incorporating a sufficient number of mid- and high-latitude stations to mitigate the effects of ionospheric variability when utilizing SF observations for satellite DCBs estimation.

GNSS standard point positioning method based on spherical harmonic expansion of signal propagation path relating errors

The errors affecting standard point positioning mainly include the ionospheric delay error, tropospheric delay error, clock bias, and multipath effect. The errors are divided into two categories depending on the relation to the signal propagation path: errors related to the propagation path and errors unrelated to the propagation path. The ionospheric delay error, tropospheric delay error, and multipath effect are related to the satellite signal propagation path. Based on this characteristic, a new GPS standard point positioning algorithm based on spherical harmonic expansion (SPP-SH) is proposed, which uses spherical harmonic expansion to represent the errors related to the propagation path. We hope to continuously promote the further development of standard point positioning (SPP) by proposing this new algorithm and to further reduce the redundancy of SPP algorithms and improve the efficiency and reliability of the algorithms. The most important purpose of this manuscript is to validate the proposed SPP-SH algorithm based on GPS data provided by IGS stations. Pseudorange data of multiple continuous epochs are used to achieve a sequential sliding point positioning solution, and the truncated singular value decomposition and iteration method correcting characteristic values are used to solve the ill-conditioned equations in the SPP-SH algorithm. The SPP-SH algorithm is implemented and verified using IGS station data. The positioning accuracy of the SPP-SH algorithm for both dual-frequency data and single-frequency data is better than that of a traditional SPP algorithm. By analyzing the positioning results of the stations at different latitudes, when the spherical harmonic expansion is expanded to the same degree, the positioning accuracy in low latitude areas is poor because the station is seriously affected by factors such as the ionosphere and water vapor, while the positioning accuracy at high latitude stations is the best. The proposed SPP-SH algorithm is also realized in the kinematic scene based on the dual-frequency data collected dynamically. The kinematic data is collected by an Android phone which supports receiving dual-frequency GPS observations. To validate the reliability of the SPP-SH algorithm, two strategies are designed to validate the results, which proves that the newly proposed method can be used for positioning with low accuracy in kinematic conditions.

Total ozone trends at three northern high-latitude stations

Abstract. After the decrease of ozone-depleting substances (ODSs) as a consequence of the Montreal Protocol, it is still challenging to detect a recovery in the total column amount of ozone (total ozone) at northern high latitudes. To assess regional total ozone changes in the “ozone-recovery” period (2000–2020) at northern high latitudes, this study investigates trends from ground-based total ozone measurements at three stations in Norway (Oslo, Andøya, and Ny-Ålesund). For this purpose, we combine measurements from Brewer spectrophotometers, ground-based UV filter radiometers (GUVs), and a SAOZ (Système d'Analyse par Observation Zénithale) instrument. The Brewer measurements have been extended to work under cloudy conditions using the global irradiance (GI) technique, which is also presented in this study. We derive trends from the combined ground-based time series with the multiple linear regression model from the Long-term Ozone Trends and Uncertainties in the Stratosphere (LOTUS) project. We evaluate various predictors in the regression model and found that tropopause pressure and lower-stratospheric temperature contribute most to ozone variability at the three stations. We report significantly positive annual trends at Andøya (0.9±0.7 % per decade) and Ny-Ålesund (1.5±0.1 % per decade) and no significant annual trend at Oslo (0.1±0.5 % per decade) but significantly positive trends in autumn at all stations. Finally we found positive but insignificant trends of around 3 % per decade in March at all three stations, which may be an indication of Arctic springtime ozone recovery. Our results contribute to a better understanding of regional total ozone trends at northern high latitudes, which is essential to assess how Arctic ozone responds to changes in ODSs and to climate change.

Multi-Instrumental Observations of Midlatitude Plasma Irregularities over Eastern Asia during a Moderate Magnetic Storm on 16 July 2003

This study presents the observations of midlatitude plasma irregularities over Eastern Asia during a moderate magnetic storm on 16 July 2003. Multi-instrumental observations, including the ground-based ionosondes, the GNSS networks, and the CHAMP and ROCSAT-1 satellites, were utilized to investigate the occurrence and characteristics of midlatitude plasma irregularities. The midlatitude strong spread F (SSF) mainly occurred in the midnight–morning sector as observed by ionosondes over Japan during this storm. SSF was related to plasma depletions, which is also recorded by GNSS network in the form of the enhancement of the rate of total electron content (TEC) change index (ROTI). The possible mechanism for the generation of SSF is that the enhanced eastward electric fields, associated with the prompt penetration electric fields and disturbance dynamo electric fields, cause the uplift and latitudinal extension of equatorial plasma bubbles (EPBs) to generate the observed midlatitude SSF further. Meanwhile, plasma density increased significantly under the influence of this storm. In addition, other common type of spread F, frequency spread F (FSF), was observed over Japan on the non-storm day and/or at high latitude station WK545, which seems to be closely related to the coupling of medium-scale traveling ionospheric disturbances (MSTIDs) and sporadic E (Es) layer. The above results indicate that various types of midlatitude spread F can be produced by different physical mechanisms. It is found that SSF can significantly affect the performance of radio wave propagation compared with FSF. Our results show that space weather events have a significant influence on the day-to-day variability of the occurrence and characteristics of ionospheric F-region irregularities at midlatitudes.

Significant elevation dependent warming over the Tibetan Plateau after removing longitude and latitude factors

As the largest plateau over China and the highest plateau over the world, the Tibetan Plateau (TP hereafter) has a vast territory spanning multiple longitude and latitude zones. Warming amplification and elevation dependent warming (EDW hereafter) are the most significant characteristics of climate change patterns over the TP. However, previous studies on the EDW over the TP have not taken the differences of longitude and latitude into consideration. Based on the observational mean, maximum and minimum temperature datasets at 118 meteorological stations over the TP, this study investigates the EDW after removing the influence of longitude and latitude factors in annual, warm season and cold season during 1961–2020 using the warming components separation method. The results show that (1) The mean warming rate over the TP during 1961–2020 is 0.33 °C·decade−1, which mainly affected by latitude and elevation. The effect of latitude on the total warming accounts for more than 65% and the effect of elevation is less than 35% on an annual basis. Moreover, the latitude factor has the greatest influence of maximum temperature in cold season (72%). (2) More robust EDW is found after removing the influences of latitude factors. In particular, the maximum temperature shows change from no feature to significant EDW (0.29 to 0.35 °C·decade−1·km−1), and the most significant change occurs in cold season (0.33 °C·decade−1·km−1) due to the removed latitude factor. It shows that stronger warming means more significant EDW, and the recent warming amplification and EDW are more intense during 1991–2020 than that during 1961–1990 when latitude factor is removed. (3) The warming components separation method effectively eliminates the differences of latitudinal warming between different stations. It mainly strengthens EDW by removing the latitudinal warming at low elevation and high latitude stations over the TP. This study highlights the importance of noise removal for the study of EDW over the TP.

A regional classification of time spectral amplitudes in total electron content: Southeastern United States during solar cycle 24

To investigate the meso-scale structure (100 s of km) of the ionosphere at mid-latitudes, the spectral properties in calculated total electron content (TEC) at a cluster of GPS receivers in and around Florida are analyzed. The ionosphere does not respond exactly the same to periodic solar driving at different locations around the planet, due to the complex electrodynamic interactions of the coupled magnetosphere-ionosphere system. Therefore, at each GPS receiver in the cluster we compare spatio-temporal variations of the spectral amplitudes for diurnal, solar rotation, and seasonal oscillations. The amplitudes for these dominant oscillations are organized with respect to magnetic latitude of the receiver. A low-latitude and high-latitude station are also included to put the mid-latitude ionospheric response into a global context. The amplitudes of diurnal, seasonal, and solar rotation signals are well ordered by magnetic latitude, superposed with meso-scale deviations between stations separated by ∼100s of km. The results suggest that spatio-temporal variations of spectral amplitudes in the mid-latitude ionosphere are not dominated by a single process. This conclusion is based on our finding that at high latitude, the shape of the diurnal signal varies significantly less with solar activity compared to low- and mid-latitudes, and additionally, that the ratio of annual to semi-annual amplitudes fluctuates around 1 with time and from station-to-station only at mid-latitudes.

Statistical Study of the Ionospheric Slab Thickness at Yakutsk High-Latitude Station

The ionospheric equivalent slab thickness (EST, also named τ) is defined as the ratio of the total electron content (TEC) to the F2-layer peak electron density (NmF2), and it is a significant parameter representative of the ionosphere. This paper presents a comprehensive statistical study of the ionospheric slab thickness at Yakutsk, located at the high latitude of East Asia, using the GPS-TEC and ionosonde NmF2 data for the years 2010–2017. The results show that the τ has different diurnal and seasonal variations in high- and low-solar-activity years, and the τ is greatest in the winter, followed by the equinox, and it is smallest in the summer in both high- and low-solar-activity years, except during the noontime of low-solar-activity years. Specifically, the τ in inter of high-solar-activity year shows an approximate single peak pattern with the peak around noon, while it displays a double-peak pattern with the pre-sunrise and sunset peaks in winter of the low-solar-activity years. Moreover, the τ in the summer and equinox have smaller diurnal variations, and there are peaks with different magnitudes during the sunrise and post-sunset periods. The mainly diurnal variation of τ in different seasons of high- and low-solar-activity years can be explained within the framework of relative variation of TEC and NmF2 during the corresponding period. By defining the disturbance index (DI), which can visually assess the relationship between instantaneous values and the median, we found that the geomagnetic storm would enhance the τ at Yakutsk. An example on 7 June 2013 is also presented to analyze the physical mechanism. It should be due to the intense particle precipitation and expanded plasma convection electric field during the storm at high-latitude Yakutsk station. The results would improve the current understanding of climatological and storm-time behavior of τ at high latitudes in East Asia.

Ozone, DNA-active UV radiation, and cloud changes for the near-global mean and at high latitudes due to enhanced greenhouse gas concentrations

Abstract. This study analyses the variability and trends of ultraviolet-B (UV-B, wavelength 280–320 nm) radiation that can cause DNA damage. The variability and trends caused by climate change due to enhanced greenhouse gas (GHG) concentrations. The analysis is based on DNA-active irradiance, total ozone, total cloud cover, and surface albedo calculations with the European Centre for Medium-Range Weather Forecasts – Hamburg (ECHAM)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) chemistry–climate model (CCM) free-running simulations following the RCP 6.0 climate scenario for the period 1960–2100. The model output is evaluated with DNA-active irradiance ground-based measurements, satellite SBUV (v8.7) total-ozone measurements, and satellite MODerate-resolution Imaging Spectroradiometer (MODIS) Terra cloud cover data. The results show that the model reproduces the observed variability and change in total ozone, DNA-active irradiance, and cloud cover for the period 2000–2018 quite well according to the statistical comparisons. Between 50∘ N–50∘ S, the DNA-damaging UV radiation is expected to decrease until 2050 and to increase thereafter, as was shown previously by Eleftheratos et al. (2020). This change is associated with decreases in the model total cloud cover and negative trends in total ozone after about 2050 due to increasing GHGs. The new study confirms the previous work by adding more stations over low latitudes and mid-latitudes (13 instead of 5 stations). In addition, we include estimates from high-latitude stations with long-term measurements of UV irradiance (three stations in the northern high latitudes and four stations in the southern high latitudes greater than 55∘). In contrast to the predictions for 50∘ N–50∘ S, it is shown that DNA-active irradiance will continue to decrease after the year 2050 over high latitudes because of upward ozone trends. At latitudes poleward of 55∘ N, we estimate that DNA-active irradiance will decrease by 8.2 %±3.8 % from 2050 to 2100. Similarly, at latitudes poleward of 55∘ S, DNA-active irradiance will decrease by 4.8 % ± 2.9 % after 2050. The results for the high latitudes refer to the summer period and not to the seasons when ozone depletion occurs, i.e. in late winter and spring. The contributions of ozone, cloud, and albedo trends to the DNA-active irradiance trends are estimated and discussed.