Two attempts of study of seismic source from teleseismic data by simulated annealing non-linear inversion

Two seismic source studies usingteleseismic data are performed by Simulated Annealing(SA), a non-linear inversion method. The Very FastSimulated Annealing (VFSA) algorithm is used and onlyteleseismic data are inverted. We have designed a fastand efficient way of multiple direct problemevaluation, which is based on pre-calculating theelementary Green's function. During the process we setthe values of the inversion control parameters(temperature, number of iterations) and modified thecooling schedule. In the current version, theinversion seeks for the point source mechanism, thedepth of the source, the scalar moment and the sourcetime function (STF). The method is applied to twoearthquakes: 18 Nov. 1992 in Greece, M 5.9 and 14 Sep.1995 in Mexico, M7.3. The calculation is performed ona simple 1D model of the structure. For the firstearthquake the inversion recovered the solution fairlywell; for the second the solution was less acceptable.However, we do not consider this fact to be a failureof the method, but a consequence of an inadequatemodel of the source and of the medium structure. Acasual attempt of reliability determination was alsoperformed; the obtained values of errors arereasonable, except for a few cases when the methodfailed.     

Tracer measurements in groyne fields for the quantification of mean hydraulic residence times and of the exchange with the stream

Tracer experiments were carried out in groyne fields (GF) of the river Elbe near Havelberg (Germany) in order to estimate the hydraulic connectivity with the river channel. The characteristic times of the five groyne fields, which were estimated from the exponentially declining tracer curves in 43 runs, ranged between 15 min and 69 min and did not correlate with the water level. Methodological investigations show that single point injection and two measurements (in the outflowing water and in the dominant region) are sufficient to provide robust in‐situ tracer curves. Using simplified mathematical simulations with connected stirred tanks, the conditions are investigated for the development of breaks in tracer curves and for the occurrence of significant errors in the estimation of intrinsic residence times. It was shown that an initial uniform dye distribution is not mandatory for the estimation during steady states. In special cases, point injections are more advisable. Moreover, the mean hydraulic residence time was found to be not equivalent to the estimated characteristic time. In fact, it is mostly overestimated by tracer experiments. The degree of overestimation depends on mixing and the volumetric proportions between the different parts of the GF and can be calculated from measured dye concentration differences. For example, an overestimation of 32% was calculated for a groyne field with a commonly found circulation flow pattern.     

Assessment of gully erosion rates through interviews and measurements: a case study from northern Ethiopia

Gullying has been widespread in the Ethiopian Highlands during the 20th century. It threatens the soil resource, lowers crop yields in intergully areas through enhanced drainage and desiccation, and aggravates flooding and reservoir siltation. Knowing the age and rates of gully development during the last few decades will help explain the reasons for current land degradation. In the absence of historical written or photographic documentation, the AGERTIM method (Assessment of Gully Erosion Rates Through Interviews and Measurements) has been developed. It comprises measurements of contemporary gully volumes, monitoring of gully evolution over several years and semi‐structured interview techniques. Gully erosion rates in the Dogu'a Tembien District, Tigray, Ethiopia, were estimated in three representative case‐study areas. In Dingilet, gullying started around 1965 after gradual environmental changes (removal of vegetation from cropland in the catchment and eucalyptus plantation in the valley bottom); rill‐like incisions grew into a gully, which increased rapidly in the drier period between 1977 and 1990. The estimated evolution of the total gully volume in the other areas show patterns similar to those of the Dingilet gully. Average gully erosion rate over the last 50 years is 6·2 t ha^?1 a^?1. Since 1995, no new gullies have developed in the study area. Area‐specific short‐term gully erosion rates are now on average 1·1 t ha^?1 a^?1. The successful application of the AGERTIM method requires an understanding of the geomorphology of the study area and an integration of the researchers with the rural society. It reveals that rapid gully development in the study area is some 50 years old and is mainly caused by human‐induced environmental degradation. Under the present‐day conditions of ‘normal’ rain and catchment‐wide soil and water conservation, gully erosion rates are decreasing. Copyright © 2006 John Wiley & Sons, Ltd.     

Glacier development and topographic context

This paper analyses the topographic context of the remaining glaciated areas in the Maladeta Massif (Central Spanish Pyrenees). These ice‐covered surfaces have been incorporated into a geographic information system (GIS) in an attempt at correlating the presence of ice with a range of topographic variables obtained from a digital elevation model. The use of generalized additive models and binary regression tree models enabled us (i) to quantify the spatial variability in the distribution of glaciers attributable to characteristics of the local terrain, (ii) to investigate the interaction between the variables that account for the ice cover distribution and (iii) to map the probability of glacier development. Our results show that although the development of glaciers depends on regional climate conditions, the topographic context is of paramount importance in determining the location, extent, shape and recent evolution of each glacial body. Thus, the joint effect of altitude, exposure to incoming solar radiation, slope and mean curvature is able to explain more than 70 per cent of the observed variance. Copyright © 2006 John Wiley & Sons, Ltd.     

Development and validation of a metallic haunch seismic retrofit solution for existing under‐designed RC frame buildings

The feasibility and efficiency of a seismic retrofit solution for existing reinforced concrete frame systems, designed before the introduction of modern seismic‐oriented design codes in the mid 1970s, is conceptually presented and experimentally investigated. A diagonal metallic haunch system is introduced at the beam–column connections to protect the joint panel zone from extensive damage and brittle shear mechanisms, while inverting the hierarchy of strength within the beam–column subassemblies and forming a plastic hinge in the beam. A complete step‐by‐step design procedure is suggested for the proposed retrofit strategy to achieve the desired reversal of strength hierarchy. Analytical formulations of the internal force flow at the beam–column‐joint level are derived for the retrofitted joints. The study is particularly focused on exterior beam–column joints, since it is recognized that they are the most vulnerable, due to their lack of a reliable joint shear transfer mechanism. Results from an experimental program carried out to validate the concept and the design procedure are also presented. The program consisted of quasi‐static cyclic tests on four exterior, ? scaled, beam–column joint subassemblies, typical of pre‐1970 construction practice using plain round bars with end‐hooks, with limited joint transverse reinforcement and detailed without capacity design considerations. The first (control specimen) emulated the as‐built connection while the three others incorporated the proposed retrofitted configurations. The experimental results demonstrated the effectiveness of the proposed solution for upgrading non‐seismically designed RC frames and also confirmed the applicability of the proposed design procedure and of the analytical derivations. Copyright © 2006 John Wiley & Sons, Ltd.     

Seismic quantification of the explosions that destroyed the dome of Volcán de Colima, Mexico, in July–August 2003

In July–August 2003, the andesitic lava dome at Volcán de Colima, México, was destroyed by a sequence of explosions that replaced the 2×10⁶ m³ dome with a crater 200 m across and 30 m deep. The two strongest explosions occurred on July 17 and August 28. The initial low-frequency impulses that they produced, which were recorded on broadband seismic records, allowed an estimation of the counter forces of the initiating process as being equal to 0.3×10¹¹ N and 1×10¹¹ N for the July and August events, respectively. The seismic characteristics follow the Nishimura-Hamaguchi scaling law for volcanic explosions, reflecting self-similarity in the processes initiating explosive events. The results also show that counter forces can discriminate between the sizes of explosive eruptions that are assigned the same magnitude by conventional methods of classification such as the Volcanic Explosivity Index. The increasing use of broadband seismometers may therefore provide the basis for using counter forces to determine the magnitude of explosive eruptions.     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Diffuse Emission of CO2 from Showa-Shinzan, Hokkaido, Japan: A Sign of Volcanic Dome Degassing

Two soil CO₂ efflux surveys were carried out in September 1999 and June 2002 to study the spatial distribution of diffuse CO₂ degassing and estimate the total CO₂ output from Showa-Shinzan volcanic dome, Japan. Seventy-six and 81 measurements of CO₂ efflux were performed in 1999 and 2002, respectively, covering most of Showa-Shinzan volcano. Soil CO₂ efflux data showed a wide range of values up to 552 g m^-2 d^-1. Carbon isotope signatures of the soil CO₂ ranged from -0.9‰ to -30.9‰, suggesting a mixing between different carbon reservoirs. Most of the study area showed CO₂ efflux background values during the 1999 and 2002 surveys (B = 8.2 and 4.4 g m^-2 d^-1, respectively). The spatial distribution of CO₂ efflux anomalies for both surveys showed a good correlation with the soil temperature, indicating a similar origin for the extensive soil degassing generated by condensation processes and fluids discharged by the fumarolic system of Showa-Shinzan. The total diffuse CO₂ output of Showa-Shinzan was estimated to be about 14.0–15.6 t d^-1 of CO₂ for an area of 0.53 km².     

Explicit formula for the geoid-quasigeoid separation

The explicit formula for the geoid-to-quasigeoid correction is derived in this paper. On comparing the geoidal height and height anomaly, this correction is found to be a function of the mean value of gravity disturbance along the plumbline within the topography. To evaluate the mean gravity disturbance, the gravity field of the Earth is decomposed into components generated by masses within the geoid, topography and atmosphere. Newton’s integration is then used for the computation of topography-and atmosphere-generated components of the mean gravity, while the combined solution for the downward continuation of gravity anomalies and Stokes’ boundary-value problem is utilized in computing the component of mean gravity disturbance generated by mass irregularities within the geoid. On application of this explicit formulism a theoretical accuracy of a few millimetres can be achieved in evaluation of the geoid-to-quasigeoid correction. However, the real accuracy could be lower due to deficiencies within the numerical methods and to errors within the input data (digital terrain and density models and gravity observations).