The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

Modelling the ground heat flux of an urban area using remote sensing data

Summary During the Basel Urban Boundary Layer Experiment (BUBBLE) conducted in 2002, micrometeorological in-situ data were collected for different sites using a variety of instruments. This provides a unique data set for urban climate studies. Nevertheless, the spatial distribution of energy and heat fluxes can only be taken into account with remote sensing methods or numerical models. Therefore, multiple satellite images from different platforms (NOAA-AVHRR, MODIS and LANDSAT ETM+) were acquired, processed and analysed. In addition, a high resolution digital elevation model (DEM) and a 1 m resolution digital surface model (DSM) of a large part of the city of Basel was utilized. This paper focuses on the calculation and modelling of the ground (or storage) heat flux density using remotely sensed data combined with in-situ measurements using three different approaches. First, an empirical regression function was generated to estimate the storage heat flux from NDVI values second approach used the Objective Hysteresis Model (OHM) which is often used for in-situ measurements. The last method used information of the geometric parameters of urban street canyons, computed from the high resolution digital urban surface model. Modelled and measured data are found to be in agreement within ±30 Wm⁻² and result in a coefficient of determination (R²) of 0.95.     

Precursors and fluid flows in the case of the 1996, ML = 5.2 Saint-Paul-de-Fenouillet earthquake (Pyrenees,France): A complete pre-, co- and post-seismic scenario

Earthquake precursors are now regularly described but often detected only after a major or moderate seismic event. Presence and influence of fluids in the seismogenic processes are often observed at the time of earthquake studies. Even today, the understanding of the physical processes involved in the source region is a real challenge for seismic hazard assessment. Here, the aftershock sequence of the M_L = 5.2, 1996 Saint-Paul-de-Fenouillet (Eastern Pyrenees, France) earthquake is first re-examined with P-wave cross-correlations, resulting in extracting three multiplets and in determining new locations. Multiplets and spatio-temporal distribution analysis of the aftershocks allow for quantifying the hydraulic diffusivity D at a maximum value of 5 m²/s and the permeability K at 10^? 15 m² in the upper Pyrenean crust. Second, a model is established in order to explain the hydrogeochemical transient anomalies, which occurred during the 15 day-period preceding the 1996 earthquake. These anomalies consist on a temporal and spatial sequence of gas emissions in the epicentral area and on chloride and lead concentration variations in a bottled mineral water 25 km north to the main shock epicenter. The proposed model processed in a standard elastic half-space, consists of creep on a low-angle crustal normal-fault, generating volumetric strain field changes over a distance of 25 km from the epicentral area. This model is able to constrain not only the mechanisms and the locations of the geochemical anomalies, but also their timing and probable casual links to the triggering of the impending major event. Also, the active extension proposed here is compatible with seismological observations in the Pyrenees. Thus, the possibility of such creep, which can be considered as a slow-slip event, is discussed in the Pyrenean tectonic and geological context. The model is discussed and compared to previous proposed models on precursor processes of earthquakes, especially concerning the preparation zone concept. Finally, a complete seismic scenario over the period beginning 15 days before the quake and ending 5 days after is proposed and discussed.     

Rainfall regionalization on the basis of the precipitation convective features using a raingauge network and weather radar observations

The main objective of the present paper is to show a methodology for undertaking rainfall regionalization of a region taking into account the convective features of the precipitation, and useful for establishing homogeneous zones for improving the alert system. This methodology has been applied to a hydrographic region located in northeast Spain, with an area of 16000 km² and characterized by a highly contrasted topography. Information provided by meteorological radar and 5-min precipitation data for 126 automatic raingauges has been used for the period 1996–2002. The previous analysis done on the basis of the 1927–1981 rainfall rate series for the Jardí raingauge, located in Barcelona, has also been considered. To that end, the first step was to draw up a proposal for classification of the pluviometric episodes. Recourse was had for this purpose to definition of the β parameter, related with the greater or lesser convective character of the event and calculated on the basis of the rainfall intensity at the surface (Llasat, 2001) and, when data are available, on the basis of radar reflectivity. Results show that the threshold of 35 mm/h to characterize convective episodes from raingauge data can be corroborated from the radar point of view when convective precipitation is identified using 2-D algorithms with a reflectivity threshold of 43 dBZ. Once the soundness of the β parameter had been corroborated, it was applied to more than 2900 precipitation episodes recorded in the region, in order to discriminate the features of the different subregions and their time and space distribution throughout the entire series of the samples. Using this definition, 92% of the precipitation events recorded in this region, with accumulated rainfall above 35 mm, are classified as convective ones, representing 95% of the precipitation amount. Application of the β parameter combined with monthly rainfall data allows differentiation of 8 regions with different convective precipitation features.     

Analysis of mesoscale convective systems in Catalonia using meteorological radar for the period 1996–2000

The aim of this paper is to show a climatology of Mesoscale Convective Systems (MCS) in the NE of the Iberian Peninsula, on the basis of meteorological radar observations. Special attention was paid to those cases that have produced heavy rainfalls during the period 1996–2000. Identification of the MCS was undertaken using two procedures. Firstly, the precipitation structures at the lowest level were recognised by means of a 2D algorithm that distinguishes between convective and non-convective contribution. Secondly, the convective cells were identified using a 3D procedure quite similar to the SCIT (Storm Cell Identification and Tracking) algorithm that looks for the reflectivity cores in each radar volume. Finally, the convective cells (3D) were associated with the 2D structures (convective rainfall areas), in order to characterize the complete MCS. Once this methodology was presented the paper offers a proposal for classifying the precipitation systems, and particularly the MCS. 57 MCS structures were classified: 49% of them were identified as linearly well-organised systems, called TS (39%), LS (18%) and NS (43%). In addition to the classification, the following items were analysed for each MCS found: duration, season, time of day, area affected and direction of movement, and main radar parameters related with convection. The average features of those MCS show an area of about 25000 km², Z_max values of 47 dBz, an echotop of 12 km, the maximum frequency at 12 UTC and early afternoon and a displacement towards E-NE. The study was completed by analysing the field at surface, the presence of a mesoscale low near the system and the quasi-stationary features of three cases related with heavy rainfalls. Maximum rainfall (more then 200 mm in 6 h) was related with the presence of a cyclone in combination with the production of a convective train effect.     

An estimate of the influence of loading effects on tectonic velocities in the Pyrenees

Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.     

Modelling the Mixing Function to Constrain Coseismic Hydrochemical Effects: An Example from the French Pyrénées

Groundwater coseismic transient anomalies are evidenced and characterized by modelling the mixing function F characteristic of the groundwater dynamics in the Ogeu (western French Pyrénées) seismic context. Investigations of water-rock interactions at Ogeu indicate that these mineral waters from sedimentary environments result from the mixing of deep waters with evaporitic signature with surficial karstic waters. A 3-year hydrochemical monitoring of Ogeu springwater evidences that using arbitrary thresholds constituted by the mean ± 1 or 2σ, as often performed in such studies, is not a suitable approach to characterize transient anomalies. Instead, we have used a mixing function F calculated with chemical elements, which display a conservative behavior not controlled by the precipitation of a mineral phase. F is processed with seismic energy release (E_s) and effective rainfalls (R). Linear impulse responses of F to E_s and R have been calculated. Rapid responses (10 days) to rainwater inputs are evidenced, consisting in the recharge of the shallow karstic reservoir by fresh water. Complex impulse response of F to microseismic activity is also evidenced. It consists in a 2-phase hydrologic signal, with an inflow of saline water in the shallow reservoir with a response delay of 10 days, followed by an inflow of karstic water with a response delay of 70 days, the amount being higher than the saline inflow. Such a process probably results from changes in volumetric strain with subsequent microfracturation transient episodes allowing short inflow of deep salted water in the aquifer. This study demonstrates that groundwater systems in such environments are unstable systems that are highly sensitive to both rainfall inputs and microseismic activity. Impulse responses calculation of F to E_s is shown to be a powerful tool to identify transient anomalies. Similar processing is suggested to be potentially efficient to detect precursors of earthquakes when long time-series (5 years at least) are available in areas with high seismicity.     

A Mediterranean nocturnal heavy rainfall and tornadic event. Part I: Overview,damage survey and radar analysis

This study presents an analysis of a severe weather case that took place during the early morning of the 2nd of November 2008, when intense convective activity associated with a rapidly evolving low pressure system affected the southern coast of Catalonia (NE Spain). The synoptic framework was dominated by an upper level trough and an associated cold front extending from Gibraltar along the Mediterranean coast of the Iberian Peninsula to SE France, which moved north-eastward. South easterly winds in the north of the Balearic Islands and the coast of Catalonia favoured high values of 0–3 km storm relative helicity which combined with moderate MLCAPE values and high shear favoured the conditions for organized convection. A number of multicell storms and others exhibiting supercell features, as indicated by Doppler radar observations, clustered later in a mesoscale convective system, and moved north-eastwards across Catalonia. They produced ground-level strong damaging wind gusts, an F2 tornado, hail and heavy rainfall. Total lightning activity (intra-cloud and cloud to ground flashes) was also relevant, exhibiting several classical features such as a sudden increased rate before ground level severe damage, as discussed in a companion study. Remarkable surface observations of this event include 24 h precipitation accumulations exceeding 100 mm in four different observatories and 30 minute rainfall amounts up to 40 mm which caused local flash floods. As the convective system evolved northward later that day it also affected SE France causing large hail, ground level damaging wind gusts and heavy rainfall.