Did the 7/9/1999 M5.9 Athens Earthquake Come with a Warning?

Prior to the 7/9/1999 M_S = 5.9 Athens earthquake, regional seismicity has exhibited a power-law increase, of the form = K+A(t_c - t)ⁿ, where is estimated using an expression log = cM + d and t_c is the time of the culminating event. Such changes appeared after the 17/8/1999 M7.4 Izmit event. We quantified the performance of the power law vs. the null hypothesis of constant seismic release rates, by defining the curvature C as the ratio of the power law fit RMS/linear fit RMS, so that the smaller C is, the better the power law behaviour. By mapping C, we have established a critical radius of 110 km and observed that the region of correlated accelerating seismic release extended from the N. Aegean, through Euboea and Attica to the SW Peloponnese. A few days prior to the Athens event, min(C) was centred at the epicentral area and numerical simulation yielded t_c = 1999.676 and predicted M_S = 5.77. Seismicity rates returned to normal (quasi-constant) after the Athens event. We interpret this effect as critical point behaviour, following remote excitation of a broad area by stress redistribution due to the Izmit event which, at Athens, has triggered `premature' failure of a fault nearing its load bearing capacity. If this is correct, we have documented a case of remote earthquake triggering by another earthquake, as well as insight into the mechanisms producing it. As a corollary, we note that a large event may beget another large event in its broader region of interaction, which may be preceded by characteristic precursory seismicity changes.     

An investigation of the active tectonics in central-eastern mainland Greece with imaging and decomposition of topographic and aeromagnetic data

We report the results of a joint analysis of aeromagnetic, topographic and tectonic data in central-eastern mainland Greece. The emphasis of the analysis is placed on the detection of coherent lineations (discontinuities), collocated and correlated with faulting structures detected by geological field observation. To this effect, edge detection and image enhancement were applied to digital aeromagnetic anomaly maps and digital elevation models, comprising bidirectional differentiation, wavelet transformation (imaging) and spatial decomposition/reconstruction in the wavenumber domain. The analysis facilitated the detection of significant topographic lineaments with NNE–SSW, ENE–WSW and ESE–WNW orientations. Respectively, the aeromagnetic data exhibit two families of significant NE–SW, and one family of ESE–WNW lineaments. The major aeromagnetic and topographic lineaments coincide and have comparable width scales of the order of 2–3 km, indicating that they are produced by significant discontinuities in the upper crust. The kinematics of the NE–SW faults varies between oblique-slip and strike-slip. These faults affect Neogene to Late Quaternary deposits and have been responsible for the formation of transverse depressions and horsts. This is also corroborated by focal plane solutions from small earthquakes recorded by local networks. The nature of these structures is not yet clear. However, they have been detected by diverse methodologies, they have considerable extent and are apparently active. These attributes suggest that they may possibly be related to the propagation and diffusion of the North Anatolian and North Aegean fault systems into the Greek mainland.     

Comparative Study of Microtremor Analysis Methods

-- During a multidisciplinary microzonation pilot project in the city of Heraklion (Crete, Greece), microtremor data were collected at the top of exploratory boreholes specifically designed for the purposes of the project, over a period of 5 days, for 4 h/day at 125 Hz (continuous recordings). The data were analysed with the SSR and H/V Ratio techniques, using the standard FFT (applied to long data series) and a Multi-variate Maximum Entropy (MV-MAXENT) spectral analysis method. Both techniques, implemented with both spectral analysis methods, identify the same major resonance frequency band, albeit with different amplification levels. The MV-MAXENT however is effective in handling short data lengths while yielding high resolution spectra and addressing several shortcomings of the conventional FFT (windowing, zero padding etc.). Thus, it yields competitively similar results, with only a fraction (a few minutes) of the data required by the lower resolution (FFT) method and appears to be a powerful tool for site effect investigations. Moreover, the results of both microtremor-based techniques are consistent and remarkably similar to the results of microzonation methods that require (expensive) borehole data.     

Periodic variation ofδ 18O values from V28–239 pacific ocean deep-sea core

Spectral analysis of ¹⁸O values from V28–239 pacific ocean deep-sea core has revealed periodicities which correspond to those calculated for the eccentricity (400 and 100 Kyrs), the obliquity (41 Kyrs) and the climatic precession (23 and 19 Kyrs) as well as secondary ones spanning between 16 Kyrs to 1 million years. The methods of spectrum analysis applied were the maximum entropy, fourier and the successive approximations, where the periodicities are located and their amplitude defined.The significance and stationarity of the detected periods was examined by various tests as well as employing an evolutionary pseudosonogram.The dominant 100 Kyrs and 50 Kyrs periods are present throughout all the interval, the 30 Kyrs is at low variance during 1.2 to 2 million years interval, the precessional signal is not stationary and appears at about 600 Kyrs to 1820 Kyrs.This study shows the necessity of applying various spectral analysis techniques and several tests to extract the optimum of spectral information and also to test the stationarity of certain periodicities, especially when implied mechanisms of climatic cause and variability are involved.Deceased.     

Development of the models to estimate particulate matter from thermal infrared band of Landsat Enhanced Thematic Mapper

Particulate matter concentration and assessment of its movement pattern is crucial in air pollution studies. However, no study has been conducted to determine the PM₁₀ concentration using atmospheric correction of thermal band by temperature of nearest dark pixels group (TNDPG) of this band. For that purpose, 16 Landsat Enhanced Thematic Mapper plus ETM+ images for Sanandaj and Tehran in Iran were utilized to determine the amount of PM₁₀ concentration in the air. Thermal infrared (band 6) of all images was also used to determine the ground station temperature (GST b6) and temperature of nearest dark pixels group. Based on atmospheric correction of images using temperature retrieval from Landsat ETM+, three empirical models were established. Non-linear correlation coefficient with polynomial equation was used to analyze the correlations between particulate matter concentration and the ground station temperature for the three models. Similar analyses were also undertaken for three stations in Klang Valley, Malaysia, using 11 Landsat ETM+ images to show the effectiveness of the model in different region. The data analysis indicated a good correlation coefficient R = 0.89 and R = 0.91 between the trend of the result of temperature of nearest dark pixels group b6 ? (GST b6 ? GST) model and the trend of PM₁₀ concentration in Iran and Malaysia, respectively. This study reveals the applicability of the thermal band of Landsat TM and ETM+ to determine the PM₁₀ concentration over large areas.     

1/f noise in the UV solar spectral irradiance

The investigation of the intrinsic properties of the solar spectral irradiance as a function of the ultraviolet (UV) wavelength is attempted by exploiting rare observations performed at the Villard St. Pancrace station of the Lille University of Sciences and Technology ranging from 278 to 400 nm with a step of 0.05 nm every half an hour from nearly sunrise to sunset. To achieve this goal, the modern method of the detrended fluctuation analysis was applied on the solar spectral irradiance values versus wavelength. This analysis revealed that the solar incident flux at the top of the atmosphere and the solar spectral irradiance at the ground during two overcast sky days fluctuate with the UV wavelength exhibiting persistent long-range power-law behavior. More interestingly, the exponent of the power-law relationship between the fluctuations of the solar spectral irradiance versus UV wavelength at both the top of the atmosphere and the ground is consistently close to unity (of 1/f-type) throughout the day. This 1/f behavior has been detected in many complex dynamical systems, but despite much effort to derive a theory for its widespread occurrence in nature, it remains unexplained so far. According to the above-mentioned findings we speculate that the 1/f property of the incident solar UV flux at the top of the atmosphere could probably drive both the 1/f behavior depicted in the atmospheric components and the solar UV irradiance at the Earth's surface. The latter could influence the UV-sensitive biological ecosystems, giving rise to a 1/f-type variability in the biosphere, which has already been proven by recent observational data. We finally propose that Wien approximation could be multiplied by a 1/f function of wavelength (e.g., of the type of the fractional Brownian motion) in order to reproduce the aforementioned 1/f feature of the solar UV flux.