Variations of <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O in precipitation from the Muztagata Glacier,East Pamirs

Expeditions to Muztagata (in the eastern Pamirs) during the summer seasons of 2002 and 2003 collected precipitation samples and measured their oxygen isotopes. The δ ¹⁸O in precipitation displays a wide range, varying from −17.40‰ to +1.33‰ in June-September 2002 and from −22.31‰ to +4.59‰ in May-August 2003. The δ ¹⁸O in precipitation correlates with the initial temperature of precipitation during the observing periods. The positive correlation between δ ¹⁸O and temperature suggests that δ ¹⁸O can be used as an indicator of temperature in this region. The δ ¹⁸O values in fresh-snow samples collected from two snow events at different elevations on the Muztagata Glacier show a strong “altitude effect”, with a ratio of nearly −0.40% per 100 m from 5500 m to 7450 m.     

First <Superscript>13</Superscript>C/<Superscript>12</Superscript>C isotopic characterisation of volcanic plume CO<Subscript>2</Subscript>

We describe analytical details and uncertainty evaluation of a simple technique for the measurement of the carbon isotopic composition of CO₂ in volcanic plumes. Data collected at Solfatara and Vulcano, where plumes are fed by fumaroles which are accessible for direct sampling, were first used to validate the technique. For both volcanoes, the plume-derived carbon isotopic compositions are in good agreement with the fumarolic compositions, thus providing confidence on the method, and allowing its application at volcanoes where the volcanic component is inaccessible to direct sampling. As a notable example, we applied the same method to Mount Etna where we derived a δ¹³C of volcanic CO₂ between −0.9 ± 0.27‰ and −1.41 ± 0.27‰ (Bocca Nuova and Voragine craters). The comparison of our measurements to data reported in previous work highlights a temporal trend of systematic increase of δ¹³C values of Etna CO₂ from ~ −4‰, in the 1970’s and the 1980’s, to ~ −1‰ at the present time (2009). This shift toward more positive δ¹³C values matches a concurrent change in magma composition and an increase in the eruption frequency and energy. We discuss such variations in terms of two possible processes: magma carbonate assimilation and carbon isotopic fractionation due to magma degassing along the Etna plumbing system. Finally, our results highlight potential of systematic measurements of the carbon isotopic composition of the CO₂ emitted by volcanic plumes for a better understanding of volcanic processes and for improved surveillance of volcanic activity.     

Disturbance of the ionospheric <Emphasis Type="Italic">D</Emphasis> region by the electric current of the atmospheric-ionospheric electric circuit

The mechanism by which electron and ion densities change in the ionospheric D region due to the electric current flowing in the atmospheric-ionospheric electric circuit is studied. The current disturbance in this circuit exists over the regions of increased seismic, meteorological, and thunderstorm activity. In the framework of the model, the influence of the electron and ion transportation under the action of the electric field on the formation of a disturbance in the D region and heating of the plasma electron component by the field are considered. The calculation results show that the densities of electrons and ions can change by an order of magnitude at an increase in the current density up to ∼(10⁻⁹–10⁻⁸) A m⁻², the sign of the disturbance depending on the current direction.     

Field-aligned currents in the winter and summer hemispheres caused by IMF <Emphasis Type="Italic">By</Emphasis>

Structures controlled by the IMF By sign and season of the year have been detected based on the decomposition of field-aligned current maps constructed using magnetic field measurements on polar low-orbiting satellites. It has been indicated that field-aligned currents have identical structures, composed of the main polar circular current and the return current at the polar cap dayside boundary, at any By sign in the summer hemisphere. Two different types of structures are implemented under winter conditions depending on the By sign. For the northern winter, it is the polar circular current and the return current at the polar cap nightside boundary at By < 0; current sheets are strongly stretched along latitudes below 80° MLat, and only small part of the current is in the noon sector of the polar cap. For the summer winter, the corresponding structures are implemented at opposite By signs. The intensities of the field-aligned currents, originating as a result of the interhemispheric asymmetry and flowing along closed geomagnetic field lines near the polar cap boundary, have been estimated. The maximum of the interhemispheric current density is 0.25 μA m⁻² in the summer and 0.1 μA m⁻² in the winter; the total current is 5 × 10⁵ and 5 × 10⁴ A, respectively.     

Distribution of<Emphasis Type="Italic">L</Emphasis><Subscript>g</Subscript> coda<Emphasis Type="Italic">Q</Emphasis> in the Chinese continent and its adjacent regions

785 traces of vertical components from shallow earthquakes recorded by 10 CDSN (Chinese Digital Seismographic Network) stations and 5 GSN (Global Seismographic Network) stations were collected to study the attenuation characteristics ofL _g coda in the Chinese continent and its adjacent regions. The records were processed first using the stack spectral ratio method to obtain the average values ofQ ₀ (Q at 1Hz) and η, the frequency dependence, ofL _g coda in the ellipses corresponding to the paths. The back-projection technique was then employed to obtain the tomographic maps ofQ ₀ and η values, and the distribution of their errors. Results indicate that in the studied areaQ ₀ varies between 200 and 500. The lowest value ofQ ₀ exists in the Yun-nan-Tibetan region, while the highest value ofQ ₀ occurs in the southern edge of Siberian platform. η varies between 0.3 and 0.8. For most part of the studied area η varies inversely withQ ₀.     

Numerical simulation of the <Emphasis Type="Italic">F</Emphasis>2-layer stratification and appearance of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere: The morphology of the phenomena

This work is devoted to a numerical simulation of the equatorial ionosphere, performed using the GSM TIP model completed with a new block for calculating the electric field. It has been indicated that the usage of the wind system calculated according to the MSIS-90 model makes it possible to reproduce the electromagnetic drift velocities at the equator, the effect of the F2-layer stratification, and the appearance of the F3 layer in the equatorial ionosphere. The calculations performed using the modified GSM TIP model made it possible to detect a maximum in the electron density vertical profile at an altitude of ∼1000 km, formed by H⁺ ions, which we called the G layer. If this layer actually exists, it can be observed during sounding the low-latitude ionosphere from satellites during dark time of day.     

Dispersion Measurements of <Emphasis Type="Italic">P</Emphasis> Waves and their Implications for Mantle <Emphasis Type="Italic">Q</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript>

We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q _p ), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q _p with an interface at 660 km depth and take Q _p constant with period. Our data set is too small to solve for a possible frequency dependence of Q _p . The upper mantle Q ₁ is 476 [299–1176] and the lower mantle Q ₂ is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q _p ^–1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135.     

Dayside auroras during unusually large positive values of the IMF <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> component

The characteristics of dayside auroras during the large (16–24 nT) positive values of the IMF B _z component, observed on January 14, 1988, during the interaction between the Earth’s magnetosphere and the body of the interplanetary magnetic cloud, have been studied based on the optical observations on Heiss Island. A wide band of diffuse red luminosity with an intensity of 1–2 kilorayleigh (kR) was observed during 6 h in the interval 1030–1630 MLT at latitudes higher than 75° CGL. Rayed auroral arcs, the brightness of which in the 557.7 nm emission sharply increased to 3–7 kR in the postnoon sector immediately after the polarity reversal of the IMF B _y component from positive to negative, were continuously registered within the band. Bright auroral arcs were observed at the equatorward edge of red luminosity. It has been found out that the red auroral intensity increases and the band equatorward boundary shifts to lower latitudes with increasing solar wind dynamic pressure. However, a direct proportional dependence of the variations in the auroral features on the dynamic pressure variations has not been found. It has been concluded that the source of bright discrete auroras is located in the region of the low-latitude boundary layer (LLBL) on closed geomagnetic field lines. The estimated LLBL thickness is ∼3 R _e . It has been concluded that the intensity of the dayside red band depends on the solar wind plasma density, whereas the position of the position equatorward boundary depends on the dynamic pressure value and its variations.     

Multi-scale Tomography for Crustal <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Velocities in Southern California

Seismic tomography is a viable tool in building depth-velocity models in the presence of strong lateral velocity variations. In this study 3-D P- and S-velocity models for the crust of southern California are constrained using more than 1,000,000 P-wave first arrivals and 130,000 S-wave arrivals from local earthquakes. To cope with the uneven distribution of raypaths, a multi-scale tomography is applied with overlapping model cells of different sizes. Within the 300 × 480 × 39 km³ model volume, the smallest cell size is 10 × 10 × 3 km³. During the iterations of velocity updating, earthquake hypocenters are determined using both P and S arrivals, and full 3-D ray tracing is implemented. Except near the edges and in the lower crust, the resultant models are robust according to various tests on the effects of reference models, resolution and signal-to-noise ratio. The tomographic velocities at shallow depths correlate very well with the regional geology of southern California. In the upper crust the P-wave and S-wave models exhibit slow velocities in major sedimentary basins and fast velocities in areas of crystalline rocks. Mid-crustal low velocity zones are present under the Coso Range, San Gabriel Mountains, and a large portion of the Mojave Desert. P- and S-velocity patterns maintain their similarity in the lower crust though the models are less reliable there.     

Variability of parameters of the <Emphasis Type="Italic">F</Emphasis>2-layer maximum in the quiet midlatitude ionosphere under low solar activity: 1. Statistical properties

Results of statistical analysis of the properties of variability of F2-layer maximum parameters (critical frequency foF2 and the height hmF2) in quiet midlatitude ionosphere under low solar activity in the daytime (1000–1500 LT) and nighttime (2200–0300 LT) hours are presented on the basis of Irkutsk station data for 2007–2008. It is found that the distribution density of δfoF2 could be presented as consisting of two distinctly different normal laws of this distribution, one of which corresponds to weak (|δfoF2| < 10%) fluctuations in foF2 and the other corresponds to strong (30% > |δfoF2| > 10%) fluctuations. Weak fluctuations in foF2 to a substantial degree are related to ionospheric variability at times less of than 1–3 h and determine the δfoF2 variability in the daytime hours. Strong fluctuations in foF2 are mainly related to day-to-day variability of the ionosphere at a fixed local time, the variability increasing by approximately a factor of 3 during the transition from day to night and determining the δfoF2 variability in the nighttime hours. The distribution density of ΔhmF2 is close to the normal distribution law. An interpretation of the different character of the distribution densities of δfoF2 and ΔhmF2 is given.     

Variations of the critical frequency <Emphasis Type="Italic">foF</Emphasis>2 at middle latitudes during strong disturbance in the electric field observed in the <Emphasis Type="Italic">F</Emphasis> region over saint-santin

Morphological analysis of variations of the critical frequency foF2 in the midlatitude ionosphere at various sectors of local time is carried out on the basis of data from ground-based stations of vertical sounding of the ionosphere in the period when during use of the incoherent scatter radar at Saint-Santin an anomalously strong increase in the electric field was observed at heights of the ionospheric F region in the period of enhanced geomagnetic activity (4+ < Kp < 6−). The obtained picture of the space-time distribution of disturbances in foF2 makes it possible to assume that they could be caused by penetration to middle latitudes of the large-scale electric field of the magnetospheric convection directed westward in the nighttime and morning hours and eastward in the noon and evening sectors.     

A role of the altitude gradient of the thermospheric wind velocity in ionospheric <Emphasis Type="Italic">F</Emphasis>-region dynamics

By processing the data of vertical ionospheric sounding in Almaty for 2000–2009, we obtained the distributions of the heights of the maximum (h _m F) and bottom (h _bot F) of the F2-layer, incremental changes in its semi-thickness (δh), the characteristic time of losses (τ), and the vertical displacement velocity of the node of the thermospheric wind (V) during the transitional time of the day during nighttime increases in the electron concentration at the layer maximum. The comparison of the measured V and modeled V _m velocities showed a certain discrepancy. The influence of the altitude gradient of the meridional thermospheric wind velocity on the behaviors of h _m F, h _bot F, δh, and τ during nighttime increases in the electron concentration is studied.     

Mechanisms of stratification of the <Emphasis Type="Italic">F</Emphasis>2 layer and formation of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere

The paper is dedicated to the studies of formation mechanisms of additional layers in the equatorial ionosphere carried out using numerical simulations with use of the Global Self-Consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) modified in the part of the solution of the electric field equation in the Earth’s ionosphere. Calculations were preformed for quiet geomagnetic conditions using the MSIS-90 model for the calculation of thermospheric parameters. The obtained spatio-temporal pattern of thermospheric circulation and the variations in the dynamo electric field obtained on its basis make it possible to reproduce the stratification effect of the F2 layer and the appearance of the F3 layer in the equatorial ionosphere due to the action of the nonuniform in height zonal electric field at the geomagnetic equator. On the basis of the earlier presented results of calculations using the modified GSM TIP model, the appearance of a maximum in the vertical profile of the electron density at a height of ∼1000 km formed by H+ ions, which we called the G layer, has been predicted. Numerical simulations showed that this layer is formed by the meridional component of the thermospheric wind and is related to the formation of the nighttime midlatitude maximum at heights of the ionospheric F region.     

The 20-second regional magnitude <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript>(20R) for the Russian Far East

A new modified magnitude scale M _S (20R) is elaborated. It permits us to extend the teleseismic magnitude scale M _S (20) to the regional epicenter distances. The data set used in this study contains digital records at 12 seismic stations of 392 earthquakes that occured in the northwest Pacific Ocean in the period of 1993–2008. The new scale is based on amplitudes of surface waves of a narrow range of the periods (16–25 s) close to the period of 20 s, for distances of 80–3000 km. The digital Butterworth filter is used for processing. On the basis of the found regional features concerning distance dependence for seismic wave attenuation, all the stations of the region have been subdivided into two groups, namely, “continental” and “island-arc.” For each group of stations, its own calibration function is proposed. Individual station corrections are used to compensate for the local features.     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

Compositional characteristics of <Emphasis Type="Italic">n</Emphasis>-alkanes of the glaciers over the Tibetan Plateau and their environmental and climatic significances

We report on the concentration and compositional features of n-alkanes of natural and anthropogenic origins in the snow samples collected from the Qiyi glacier in the Qilian Mountains, the Yuzhufeng glacier in eastern Kunlun Mountains, the Xiaodongkemadi glacier in the Tanggula Mountains, and the Gurenhekou glacier in the Nyainqêntanglha Range. The results indicate a decrease in the total n-alkane concentration (T-HCs) from the northeast to the south over the Tibetan Plateau. The T-HCs in these studied areas were close to those in the Belukha and Sofiyskiy glacier, Russian Alati Mountains and the Dasuopu glacier in the Himalaya but were much higher than those in the Greenland ice sheet, suggesting that the mountain glaciers in the Asian continent may receive a higher loading of n-alkanes than the Greenland ice core. Moreover, the compositional characteristics of n-alkanes indicated that the n-alkanes in the studied areas were probably originated from the plant waxes as well as the fossil-fuel combustion exhaust, whereas the contribution from the lower organisms was small. In addition, the plant wax (C_n(wax)) and anthropogenic (non-C_n(wax)) contributions revealed that fast industrialization may have significant effects on the organic pollutant composition in glacier over the Tibetan Plateau and its circumference environment. Particularly, except for the Yuzhufeng glacier, the ΣnC₂₁ ⁻/ΣnC₂₂ ⁺ and (nC₁₅+nC₁₇+nC₁₉)/(nC₂₇+nC₂₉+nC₃₁) ratio decreased from the Qiyi glacier to the Gurenhekou glacier over the Tibetan Plateau, while the carbon preference index (CPI) values increased. These results indicate a decrease in terrigenous input while an increase in marine input from the northeast to the south over the Tibetan Plateau. These two ratios can be used as the climatic and environmental change indicators.     

Response of dayside auroras to abrupt increases in the solar wind dynamic pressure at positive and negative polarity of the IMF <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> component

The optical observations on Heiss Island (Φ′ = 75.0°) have been used to study the characteristics of auroras in the near-noon MLT sector after abrupt increases in the solar wind dynamic pressure at negative and positive polarity of the IMF B _z component. It has been found out that the 427.8 and 557.7 nm emission intensities considerably increased at B _z < 0 both equatorward of the dayside red luminosity band and within this band. The value of the emission intensities at a red luminosity maximum (I ₆₃₀₀/I ₅₅₇₇ ∼ 0.5) indicates that energetic electron precipitation is of the magnetospheric origin. At B _z > 0, fluxes of harder (E > 1 keV) precipitating electrons were superimposed on the soft spectrum of precipitating particles in the equatorial part of the red luminosity band. This red band part was hypothetically caused by the low-latitude boundary layer (LLBL) on closed lines of the geomagnetic field, the estimated thickness of which is ∼3 R _e . The 557.7 nm emission intensity increased during 3–5 min after SC/SI and was accompanied by the displacement of the red band equatorward boundary toward lower latitudes. The displacement value was ∼150–200 km when the dynamic pressure abruptly increased by a factor of 3–5. After SC/SI, the 630.0 nm emission intensity continued increasing during 16–18 min. It is assumed that the time of an increase in the red line intensity corresponds to the time of saturation of the magnetospheric boundary layers with magnetosheath particles after an abrupt increase in their density.     

High-resolution climate variability of southwest China during 57-70 ka reflected in a stalagmite δ18O record from Xinya Cave

A 26-cm-long stalagmite (XY2) from Xinya Cave in northeastern Chongqing of China has been ICP-MS ²³⁰Th/U dated, showing a depositional hiatus at 2.3 cm depth from the top. The growth of the 2.3–26 cm interval determined by four dates was between 57 ka and 70 ka, with a linear growth rate of 0.023 mm/a. We have analyzed 190 samples for δ ¹⁸O and δ ¹³C, mostly in the 2.3–26 cm part. The δ ¹⁸O and δ ¹³C values between 57 ka and 70 ka reveal decadal-to-centennial climatic variability during the glacial interval of Marine Isotope Stage 4 (MIS4), exhibiting much higher resolution than that of the published Hulu and Dongge records during this interval. Speleothem δ ¹⁸O in eastern China, including our study area can be used as a proxy of summer monsoon strength, with lighter values pointing to stronger summer monsoon and higher precipitation, and vice versa. Two decreases in the δ ¹⁸O signature of XY2 record around 59.5 and 64.5 Ka are argued to correspond to the Dansgaard-Oeschger (D-O) events 17 and 18 respectively. The Heinrich event 6 (H6) can be identified in the record as a heavy δ ¹⁸O peak around 60 ka, indicating significant weakening of the monsoon in Chongqing during the cold period. The XY2 δ ¹⁸O record shows very rapid change toward to the interstadial condition of the D-O event, but more gradual change toward to the cold stadial condition. This phenomenon found in the Greenland ice core records is rarely observed so clearly in previously published speleothem records. According to SPECMAP δ ¹⁸O record, the glacial maximum of MIS 4 was around 64.5 ka with the boundary of MIS 3/4 around 60 ka. Unlike the marine record, the speleothem record of XY2, China, exhibits much high frequency variations without an apparent glacial maximum during MIS 4. However, the timing of MIS 3/4 boundary seems to be around 60 ka when the H6 terminated, in agreement with the marine chronology. The growth period of sample XY2 during glacial times probably reflects a local karstic routing of water, rather than having climatic significance. Supported by the National Natural Science Foundation of China (Grant Nos. 40672165, 90511004, 40672202) and the Academician Special Project of Chongqing Science Committee (Grant No. 2003-7835)     

Carbon isotope ratios of C4 plants in loess areas of North China

Carbon isotope ratios (δ ¹³C) of 89 C₄ plant samples were determined from the loess area in North China. δ ¹³C values vary between −10.5‰ and −14.6‰ with a mean of −12.6‰. Along a precipitation gradient from the semi-moist area to the semiarid area, then to the arid area, the δ ¹³C values of C₄ plants show a slight decreasing trend. The δ ¹³C values of C₄ plants in the dry season are found lower than those in the wet season. These trends are opposite to those observed for C₃ species.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> and <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript> comparison for earthquakes in Greece for the period 1969 to 2007

We use 576 earthquakes of magnitude, M _w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M _w, and local magnitude, M _L, for earthquakes in Greece and the surrounding regions. We compare M _w to M _L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M _L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M _w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M _L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M _w and M _L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M _L’s. We propose scaling relations to convert the M _L values reported from the two centers to M _w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.