Regionalization of seismic hazard in Greece based on seismic sources

A semi-probabilistic approach to the seismic hazard assessment of Greece is presented. For this reason, a recent seismotectonic model for shallow and intermediate depth earthquake sources, based on historical as well as on instrumental data, was used. Different attenuation formulae were proposed for the macroseismic intensity and the strong ground motion parameters for the shallow and the intermediate focal depth shocks. The data were elaborated in terms of McGuire's computer program, which is based on the Cornell's method.A grid of equally spaced points at 20 km distance was made and the seismic hazard recurrence curves for various parameters of the seismic intensity was estimated for each point. Finally, seismic hazard maps for the area of Greece were compiled utilizing the entire range of recurrence curves. These maps depict areas of equal seismic hazard and for every area the analytical relations of the typeSI =f(Tm), whereSI is a seismic intensity parameter andTm is the mean return period, were determined.     

Development of a generic concept to analyze the accessibility of emergency facilities in critical road infrastructure for disaster scenarios: exemplary application for the 2017 wildfires in Chile and Portugal

Natural Hazards - Natural hazards such as earthquakes, floods, or wildfires pose a serious threat to road infrastructure. Especially in emergency situations, the society depends on the road...     

Dynamical response of low-latitude middle atmosphere to major sudden stratospheric warming events

The zonally averaged UK Meteorological Office (UKMO) zonal mean temperature and zonal winds for the latitudes 8.75°N and 60°N are used to investigate the low-latitude dynamical response to the high latitude sudden stratospheric warming (SSW) events that occurred during winter of the years 1998–1999, 2003–2004 and 2005–2006. The UKMO zonal mean zonal winds at 60°N show a short-term reversal to westward winds in the entire upper stratosphere and lower mesosphere and the low-latitude winds (8.75°N) show enhanced eastward flow in the upper stratosphere and strong westward flow in the lower mesosphere during the major SSW events at high latitudes. The mesosphere and lower thermosphere (MLT) zonal winds acquired by medium frequency (MF) radar at Tirunelveli (8.7°N, 77.8°E) show a change of wind direction from eastward to westward several days before the onset of SSW events and these winds decelerate and weak positive (eastward) winds prevail during the SSW events. The time variation of zonal winds over Tirunelveli is nearly similar to the one reported from high latitudes, except that the latter shows intense eastward winds during the SSW events. Besides, the comparison of daily mean meridional winds over Tirunelveli with those over Collm (52°N, 15°E) show that large equatorial winds are observed over Tirunelveli during the 2005–2006 event and over Collm during the 1998–1999 events. The variable response of MLT dynamics to different SSW events may be explained by the variability of gravity waves.     

Double origin of hydrothermal convective flux variations in the Fossa of Vulcano (Italy)

Soil-temperature measurements can provide information on the distribution of degassing fissures, their relationship to the internal structure of the volcano, and the temporal evolution of the system. At Vulcano Island (Italy), heat flux from a <3 km-deep magma body drives a hydrothermal system which extends across the main Fossa crater. This heat flux is also associated with variable magmatic gas flow. A high-density map of soil-temperatures was made in 1996 at a constant depth of 30 cm on the central and southern inner flanks of the Fossa crater. These measurements extended over an area covering about 0.04 km², across which the heat flux is predominantly associated with a shallow boiling aquifer. The map shows that hot zones relate to structures of higher permeability, mainly associated with a fissure system dating from the last eruptive cycle (1888–1890). From 1996 to January 2005, we studied the evolution of the heat flux for the high temperature part of the map, both by repeating our measurements as part of 14 visits, during which temperatures were measured at a constant depth, and using data from permanent stations which allowed soil-temperatures to be continuously measured for selected vertical profiles. These data allowed us to calculate the heat flux, and its variation, with good precision for values lower than about 100 W m⁻², which is generally the case in the study area. Above 100 W m⁻², although the heat flux value is underestimated, its variations are recorded with an error less than 10%. During the period 1996–2004, two increases in the thermal flux were recorded. The first one was related to the seismic crisis of November 1998 which opened existing or new fissures. The second, in November 2004, was probably due to magma migration, and was associated with minor seismic activity.     

Calculated optical absorption spectra of Ni2+-bearing compounds

The transition energies responsible for optical absorption spectra can be obtained by crystal-field analysis, but the transition intensities are notoriously difficult to calculate. This paper examines the basic ingredients of the calculation of optical spectrum intensities. Magnetic dipole and electric quadrupole transitions intensities are evaluated, as well as the direct d(Ni²⁺) to p(O²⁻) electric dipole transitions. All these contributions are shown to be small in the optical range, so that spectral intensities are due to the mixing of odd orbitals with the Ni²⁺ 3d ⁿ states. Received: 11 November 1997 / Revised, accepted: 6 September 1999     

Building design based on energy dissipation: a critical assessment

The basic objective of this study is the assessment of the European seismic design codes and in particular of EC2 and EC8 with respect to the recommended behaviour factor q. The assessment is performed on two reinforced concrete multi-storey buildings, having symmetrical and non-symmetrical plan view respectively, which were optimally designed under four different values of the behaviour factor. In the mathematical formulation of the optimization problem the initial construction cost is considered as the objective function to be minimized while the cross sections and steel reinforcement of the beams and the columns constitute the design variables. The provisions of Eurocodes 2 and 8 are imposed as constraints to the optimization problem. Life-cycle cost analysis, in conjunction with structural optimization, is believed to be a reliable procedure for assessing the performance of structures during their life time. The two most important findings that can be deduced are summarized as follows: (1) The proposed Eurocode behaviour factor does not lead to a more economical design with respect to the total life-cycle cost compared to other values of q (q = 1, 2). (2) The differences of the total life-cycle cost values may be substantially greater than those observed for the initial construction cost for four different q (q = 1, 2, 3, 4).     

A new isostatic model of the lithosphere and gravity field

Based on the analysis of various factors controlling isostatic gravity anomalies and geoid undulations, it is concluded that it is essential to model the lithospheric density structure as accurately as possible. Otherwise, if computed in the classical way (i.e. based on the surface topography and the simple Airy compensation scheme), isostatic anomalies mostly reflect differences of the real lithosphere structure from the simplified compensation model, and not necessarily the deviations from isostatic equilibrium. Starting with global gravity, topography and crustal density models, isostatic gravity anomalies and geoid undulations have been determined. The initial crust and upper-mantle density structure has been corrected in a least squares adjustment using gravity. To model the long-wavelength (>2000 km) features in the gravity field, the isostatic condition (i.e. equal mass for all columns above the compensation level) is applied in the adjustment to uncover the signals from the deep-Earth interior, including dynamic deformations of the Earths surface. The isostatic gravity anomalies and geoid undulations, rather than the observed fields, then represent the signals from mantle convection and deep density inhomogeneities including remnants of subducted slabs. The long-wavelength non-isostatic (i.e. the dynamic) topography was estimated to range from –0.4 to 0.5 km. For shorter wavelengths (<2000 km), the isostatic condition is not applied in the adjustment in order to obtain the non-isostatic topography due to regional deviations from classical Airy isostasy. The maximum deviations from Airy isostasy (–1.5 to 1 km) occur at currently active plate boundaries. As another result, a new global model of the lithosphere density distribution is generated. The most pronounced negative density anomalies in the upper mantle are found near large plume provinces, such as Iceland and East Africa, and in the vicinity of the mid-ocean ridge axes. Positive density anomalies in the upper mantle under the continents are not correlated with the cold and thick lithosphere of cratons, indicating a compensation mechanism due to thermal and compositional density.     

Integrated adjustment of CHAMP, GRACE, and GPS data

Various types of observations, such as space-borne Global positioning system (GPS) code and phase data, accelerometer data, K-band range and range-rate data, and ground-based satellite laser ranging data of the CHAllenging Minisatellite Payload (CHAMP) and GRAvity Climate Experiment (GRACE) satellite missions, are used together with ground-based GPS code and phase data in a rigorous adjustment to eventually solve for the ephemerides of the CHAMP, GRACE, and GPS satellites, geocenter variations, and low-degree gravity field parameters. It turns out that this integrated adjustment considerably improves the accuracy of the ephemerides for the high and low satellites, geocenter variations, and gravity field parameters, compared to the case when the adjustment is carried out stepwise or in individual satellite solutions.Acknowledgments. This study has been supported by the German Ministry of Education and Research through the Geotechnologies Programme grants 03F0333A/CHAMP and 03F0326A/GRACE.