Sensitivity of the Sharpe and Worden gravity meters to vertical motion

Summary Using the least-squares method the constants of both gravity meters — seismometers — were derived from the amplitude responses, observed under vertical motion on a vibrating table in the 0.3–30 s range. The values of the seismogram amplitudes, above which seismic waves cause an observable displacement of the gravity meter reading beam, were determined for the SKD and Press-Ewing station seismographs.     

Ozone distribution over Mediterranean,central and southeast Europe during the International Quiet Sun Year (1964–1965)

Summary Ozone observations made during 1964 and 1965 at nine Mediterranean, central and southeast European stations (latitudes 38–52°N, longitudes 9–23°E) reveal patterns of seasonal and shorter time-variations in total ozone as well as in vertical ozone distribution. During the winter-spring season, a significant increase (20%) of ozone occurs essentially simultaneously with the spring stratospheric warming, and is noticed at all stations.—Autocorrelation coefficients show that the total ozone on any day is strongly related to the total ozone of the preceding four days in summer or one or two days in winter-spring or autumn. Changes of total ozone in southeast Europe correlate closely with those in Mediterranean Europe, and less closely with those from north central Europe.—Power spectrum analysis detects the dependence of ozone changes on processes with periods longer than 6–8 days, and indicates a significant oscillation with a period of 14–15 days, perhaps a result of the direct influence of lower stratospheric circumhemispheric circulation. — Reliable vertical ozone soundings were not available from all stations. The mean vertical profiles at Arosa, Switzerland (47°N) and Belsk, Poland (51°) are very similar. More than 60% of the variability of the total ozone is contributed by changes in ozone concentration between 10 and 24 km; less than 10% is due to variations above 33 km. Changes in ozone partial pressure at different altitudes, and relationships of those changes to total ozone, indicates that a mean vertical ozone distribution may be described adequately by considering the ozone changes in four layers: a) the troposphere, b) the lower stratosphere up to 24 km, c) a transition layer from 24 km to a variable upper border at 33–37 km, and d) the layer above 33–37 km.Part of this paper was presented at the Ozone Seminar in Potsdam, Germany, 27 September 1966.     

Seismic noise at the international gravity point in Prague-Ruzyně (Czechoslovakia)

Summary The vertical component of ground displacement was measured at the Prague - Ruzyn International Gravity Point in the frequency range of 1–300 Hz. The permanent noise, the vibrations caused by the observers during gravimetric observations and by the wind, as well as those due to normal operations at the airport, display maximum peak-to-peak amplitudes of 0.06 µm in the frequency range of 1–50 Hz; with a CG-2 gravimeter this is not detrimental to the accuracy of the observations. The taxiing of turbo-jet and jet aircraft and engine tests of aircraft generate vibrations in frequency ranges of 75–90 and 190–270 Hz. Their amplitudes, according to the results of laboratory tests published for various types of gravity meters (CG-2, GAK-PT, GVP-3, KVG), are of magnitudes which generate errors in tenths of mgl.     

Fundamental geodetic parameters of the earth's figure and the structure of the earth's gravity field derived from satellite data

Summary Using the geocentric constant GM=398 601.3 × 10 ⁹ m ³s ^–2 , the known value of the angular velocity of the Earth's rotation , Stokes' constants J _n ^(k) and S _n ^(k) upto n=21 (zonal), n=16 (tesseral and sectorial) [2], the geocentric co-ordinates and heights above sea-level of SAO satellite stations [2], the following will be derived: the potential on the geoid W_o, the scale factor for lengths R_o=GM/W_o, the radius-vector of the surface W=W_o, the parameters of the best-fitting Earth tri-axial ellipsoid, and the components of the deflections of the vertical with respect to the geocentric rotational IAG ellipsoid (Lucerne 1967), as well as to the best-fitting geocentric tri-axial ellipsoid. Some of the differences in the structure of the gravity field over the Northern and Southern Hemispheres will be given, and the mean values of gravity over the equatorial zone, determined from the dynamics of satellite orbits, on the one hand, and from terrestrial gravity data, on the other, will be compared.Presented at the Fifteenth IUGG General Assembly, Moscow, July 30 — August 14, 1971.     

Vertical velocities associated with gravity waves measured in the mesosphere and lower thermosphere with the EISCAT VHF radar

The EISCAT VHF radar (69.4°N, 19.1°E) has been used to record vertical winds at mesopause heights on a total of 31 days between June 1990 and January 1993. The data reveal a motion field dominated by quasi-monochromatic gravity waves with representative apparent periods of 30–40 min, amplitudes of up to 2.5 m s^–1 and large vertical wavelength. In some instances waves appear to be ducted. Vertical profiles of the vertical-velocity variance display a variety of forms, with little indication of systematic wave growth with height. Daily mean variance profiles evaluated for consecutive days of recording show that the general shape of the variance profiles persists over several days. The mean variance evaluated over a 10 km height range has values from 1.2 m²s^–2 to 6.5 m²s^–2 and suggests a semi-annual seasonal cycle with equinoctial minima and solsticial maxima. The mean vertical wavenumber spectrum evaluated at heights up to 86 km has a slope (spectral index) of -1.36 ± 0.2, consistent with observations at lower heights but disagreeing with the predictions of a number of saturation theories advanced to explain gravity-wave spectra. The spectral slopes evaluated for individual days have a range of values, and steeper slopes are observed in summer than in winter. The spectra also appear to be generally steeper on days with lower mean vertical-velocity variance.     

Falling sphere observations of anisotropic gravity wave motions in the upper stratosphere over Australia

A study of inertial scale gravity wave motions in the region of the atmosphere between 30 and 60 km has been undertaken, using wind and temperature data derived from rocket-borne falling sphere density experiments performed over Woomera, Australia between 1962 nad 1976. The gross features of the wave field compare favorably with those found in similar northern hemispheric studies. Wave propagation is found to be both vertically and horizontally anisotropic. A rotary spectral analysis indicates predominately upgoing wave energy, suggesting that the majority of sources of these waves lie below 30 km. A detailed statistical investigation of the waves, made using the Stokes parameters technique, reveals that phase progression is also highly directional in the horizontal, with a significant zonal component in summer, but with a strong meridional component in winter. Propagation towards the southeast is inferred in summer, with the waves possibly emanating from tropospheric sources in equatorial regions to the north of Australia. The technique also shows that, on average, the waves appear to have mean ellipse eccentricities (=f/) around 0.4–0.45. Indirect estimates of a number of important wave parameters are made. In particular,v andw flux estimates are made over several height intervals. The vertical gradient of density weighted flux implies wave-induced mean flow accelerations of the order 0.1–1 ms^–1day^–1. This suggests that dissipating gravity waves are a significant source of the momentum residuals that are encountered in studies of satellite data from this region.     

Mean equatorial gravity derived from various geopotential models

Summary Mean equatorial gravity has been computed from geopotential models GEM-10C, GEM-7, GEM-T1, GEM-T2, GEM-T3, JGM-1, JGM-2, JGM-3 and OSU91A and compared to the normal equatorial gravity, _e=978 032·699 × 10^–5 m s^–2, computed from four given parameters defining the Earth's level ellipsoid. In all models g_e>_e.     

First study of fault trench stratigraphy at Mt. Etna volcano, Southern Italy: understanding Holocene surface faulting along the Moscarello fault

Paleoseismology, the study of past earthquakes based on their geological record in the stratigraphy and landscape, is a successful newly developing field of research. The application of fault trench studies in volcanic environments is one of the youngest branches of paleoseismology. In this paper, we present the results of the first exploratory trenches excavated at Mt. Etna in Sicily, the largest European volcano. Modern surface faulting at Etna is a very well known feature, which poses significant hazard to the local community, both in terms of ground displacement of essential lifelines and ground shacking from frequent damaging earthquakes. However, while the geomorphology and the seismicity of the active fault in the Etna region consistently show very high rates of tectonic activity, the Holocene cumulative throw and slip-rates, along with the nature (coseismic vs. creeping fault slip), dimension and timing of the displacement events, are still poorly constrained. For this purpose, we selected as a sample area the Moscarello fault, one of the most outstanding segment of the Timpe system of active normal faults in the volcano’s lower eastern flank. Displaced landforms and volcanic units at the Fondo Macchia basin, in the central sector of this fault, indicate some hundreds of meters of vertical offset in the last ca. 80 kyr, with a long-term slip-rate substantially higher than 1.5–2.0 mm/yr. According to the historical sources and instrumental observations, the Moscarello fault ruptured four times in the last 150 years during shallow (H < 5 km) and moderate magnitude (M < 4.8) earthquakes. These events were associated with severe damage in a narrow epicentral area (macroseismic intensities up to the IX–X grade of the MSK scale) and extensive surface faulting (end-to-end rupture length up to 6 km, vertical offsets up to 90 cm). This clearly indicates very high modern rates of deformation along this fault. We conducted trench investigations at the Fondo Macchia site, in a point where eyewitnesses observed ca. 20 cm of coseismic vertical displacement after the April 21, 1971, M_s=3.7, earthquake. The excavated sections provided direct stratigraphic evidence for a vertical slip-rate of 1.4–2.7 mm/yr in the last ca. 6 kyr. This should be regarded as a minimum slip-rate for the central section of the fault. We explored a single scarp at a single site, while we know from recent historical observations that several parallel scarps may rupture coseismically at Fondo Macchia. Thus, the relevant deformation rate documented for the modern period might be likely extended back in the past to a time-span of some thousands of years at least. As expected, for such a volcanic environment, the activity rates of the Moscarello fault are also significantly higher than for the Apennines normal faults, typically showing slip-rates lower than 1 mm/yr. The agriculturally reworked trench hangingwall stratigraphy did not allow to recognize individual displacement events. Nevertheless, the sedimentary structures observed in the trench footwall strongly suggest that, as for the last 150–200 years of detailed historical record, fault behavior at Fondo Macchia is governed by coseismic surface displacement rather than fault creep. This research confirms that paleoseismology techniques can be effectively applied also in active volcanic environments, typically characterized by rheology and, consequently, seismicity and fault dynamics very different from those of other tectonic environments in which paleoseismology has been firstly developed and is today extensively applied.     

The mean curvature of the external equipotential surface and the vertical gravity gradient as functions of stokes's constants

Summary The mean curvature of the equipotential surface and the vertical gradient of gravity are expressed in terms of a development into a series of spherical harmonics [1, 2, 4], neglecting terms of the order of 10^–8. The curvature anomalies have been computed using the satellite data [3]. The symbols used are the same as in [5].Dedicated to 90th Birthday of Professor Frantiek Fiala     

A study of the upward continuation of gravity data from a plane surface

Summary A summation method of upward continuation of gravity data has been considered under the assumption that observations are available at regular intervals. The upward continued value has been obtained as the sum of products of individual gravity values and corresponding theoretical coefficients. Besides the usual parameter involving horizontal and vertical distances, the theoretical coefficients have been generalized to be dependent also on i) the order of a low order polynomial assumed to represent the gravity variation around a grid point and ii) the weights assigned to the gravity values at the nearest four grid points used for least-squares determination of the polynomial. Since the observations in practical cases are available over a finite area only, the effect of truncation of the area of the integration has also been discussed separately. The method has been programmed and tested on a three-dimensional model, whose true gravity effects were computed at various levels over a finite area. Upward continued values have been computed under various assumptions about the gravity field in the outside region. Comparisons of these results with the true values indicate that the truncation effect becomes increasingly important for larger values of the ratio of elevation to grid separation and/or when the gravity field is not a local one. It has also been found that the greater is the above ratio, the less important is the effect of weights on the theoretical coefficients and practically vanishes(<10 ^–4 ) when the ratio is greater than5.0.     

Regional gravity variations in Europe from superconducting gravimeters

Recent satellite missions (CHAMP, GRACE) are now returning data on the time variation of the gravity field with harmonic coefficients computed every 4 weeks. The promise is to achieve a sub-microgal accuracy that will define continental mass variations involving large-scale hydrology. With this in mind, we examine the time varying gravity field over central Europe using a limited number of high quality ground-based superconducting gravimeter stations within the Global Geodynamics Project (GGP). Our purpose is to see whether there are coherent signals between the individual stations and to compare the regional component with that predicted from models of continental hydrology. The results are encouraging. We have found, using empirical orthogonal eigenfunctions of the gravity data that a clear annual signal is present that is consistent in phase (low amplitudes in summer) and amplitude (1–3 microgal) with that determined from a large-scale model of land water in connection with global climate modeling. More work is required to define how the gravity field is related to large-scale soil moisture and other mass variations, and we have yet to compare our results to the latest satellite-derived data.     

Seismicity Variations in the Southern Aegean,Greece, Before and After the Large (M7.7) 1956 Amorgos Earthquake Due to Evolving Stress

The largest earthquake (M₀=4.9·10²⁷ dyn·cm) of the 20^th century in the territory of Greece occurred south of Amorgos Island, causing extensive destruction in the southern Aegean area. It occurred on an ENE–trending normal fault that is seated parallel to the Islands southern coastline. Changes in the rates of moderate–size earthquakes (M 5.0) that occurred before and after the Amorgos earthquake, within circular regions centered on its epicenter with radii of 100, 150 and 200 km, are investigated. The rate for moderate–size events just before the main shock appears to be considerably increased when compared to those of either preceding or subsequent periods. Further inspection reveals that more evident seismicity fluctuations are attributed to distances exceeding 100 km. These changes may be indicative of a broad region that is approaching a high stress state prior to an eventual large earthquake. Close to the main event, that is, within the 100–km radius, a remarkable quiescence period lasting about two decades before its occurrence was observed. Changes in seismicity are discussed in combination with static stress changes calculated by the application of the stress evolutionary model that takes into account the coseismic slip associated with the larger events (M 6.5) since the beginning of the 20^th century and the tectonic loading on the major faults in the study area. These larger events, as with the intermediate magnitude seismicity taking place at distances exceeding 100 km and which encircled the quiescent area observed during the last 22 years before the Amorgos earthquake, are well correlated with stress-enhanced areas in each stage of the evolutionary model.     

Combined discrete and continuous gravity observations at Mount Etna

Systematic investigation of discrete gravity measurements has continued at Mount Etna since 1986. The network now covers an area of 400 km² with about 70 stations 0.5–3 km apart. Mass redistributions occurring at depths ranging between about 8 km below sea level and a few hundred metres below the surface (magma level changes within the shallower parts of the feeding conduits) have been identified from these data. Conventional (discrete) microgravity monitoring on a network of stations furnishes only instantaneous states of the mass distribution at continuously active systems. In order to obtain information on the rate at which the volcanic processes (and thus mass transfers) occur, three stations for continuously recording gravity where installed on Mount Etna in 1998. A 16-month long sequence from one of the continuously running stations (PDN, located 2 km from the active northeast crater at the summit of Etna volcano) is presented. After removing the effects of Earth Tide and tilt, the correlation of the residual gravity sequence with simultaneous recordings of meteorological parameters acquired at the same station was analysed. Once the meteorological effects have also been removed, continuous gravity changes are within 10 μGal of gravity changes measured using conventional microgravity observations at sites very close to the continuous station. This example shows how discrete and continuous gravity observations can be used together at active volcanoes to get a fuller and more accurate picture of the spatial and temporal characteristics of volcanic processes.     

On the interpretation of vertical gravity gradients produced by magmatic intrusions

A method for the computation and interpretation of gravity and height changes and vertical gravity gradients produced by magmatic intrusions in a layered elastic–gravitational medium is presented. The methodology assumes a planar medium geometry, which consists of welded layers overlying a half-space. The medium is elastic and gravitating. The intrusion is treated as a point source and can be located at any depth inside the medium. The theoretical elastic–gravitational model allows the computation of the so-called geometric and orthometric vertical displacements as well as gravity changes of different types. The corresponding vertical gravity gradients can also be computed. We present several examples of theoretical computations and study different types of these geometric and orthometric gravity gradients and the information we can get from the application of the methodology described. The results presented show that the use of these gradients is a useful tool to obtain information on the dynamics of the injection processes, including the detection of new magma recharge. We show that using the elastic–gravitational deformation model we can explain non-linear gravity–height relationships that appear in volcanic areas. We also present the application of the methodology to Mayon volcano, Philippines, delineating the intrusion of new magma, consistent with the produced eruptions after the observation period.     

Do CRISTA Experiment/Campaign Data Represent a Typical Situation or Not?

Data from four A3 radio wave circuits in central Europe are used to analyse the representativeness of the ionization and gravity wave activity pattern in the lower ionosphere (85 – 100 km) in the CRISTA experiment interval (3 – 12 November 1994) for the given conditions. It has been found that the CRISTA experiment interval was run under conditions, which are highly representative both of October – November 1994 and autumn, low-moderate solar activity, the descending phase of solar cycle conditions, i.e. CRISTA measurements may be considered to provide values identical with, or close to climatological values (at least for central Europe, h = 85 – 100 km).     

Influence of soil consolidation and thermal expansion effects on height and gravity variations

The daily GPS height series of the Medicina station were analyzed for the period July 1996–September 2001. The station is located in the middle Po Plain on fine-grained alluvial deposits. A seasonal oscillation in the order of 18 mm (peak-to-peak amplitude) is present in the data. This crustal deformation has been modeled by including variations in the atmospheric, oceanic and hydrologic mass. The vertical positions can also be affected significantly by soil consolidation. Geotechnical parameters derived by in situ tests and laboratory analyses of the clayey soil collected at Medicina allowed the estimate of the soil settlement relevant to the seasonal oscillation of the surficial water table. Thermal expansion of the geodetic monument has to be taken into account in the case of high-precision vertical positioning. In this work models both for the soil consolidation and the thermal expansion effects are provided. The continuous gravity observations collected at Medicina by means of a superconducting gravimeter also exhibit a marked seasonal oscillation, which has been interpreted as the sum of loading and Newtonian attraction effects, as well as of the contribution due to soil consolidation. Especially the study concerning the soil consolidation effect has allowed a better insight on the seasonal vertical movements occurring at the Medicina station by providing quantitative information on soil behavior due to change of effective pressures. The results can be applied to those stations characterized by similar fine-grained soils and surficial hydrogeology.     

The effect of the Baltic Sea level on gravity at the Metsähovi station

The loading effect of the Baltic Sea is immediately recognizable in the gravity record of the superconducting gravimeter T020 in Metsähovi, Finland, by simply inspecting residual gravity together with the tide gauge record at Helsinki 30 km away. The station is 10 km from the nearest bay of the Baltic Sea and 15 km from the open sea. Sea level variations in the Baltic are non-tidal and driven at short periods primarily by wind stress, at longer periods by water exchange through the Danish straits. Locally they can have a range of 2–3 m. Loading calculations show that a uniform layer of water covering the complete Baltic Sea increases the gravity in Metsähovi by 31 nm/s² per 1 m of water, and the vertical deformation is −11 mm. The observed gravity response to the local sea level is generally less, since the variations at short periods are far from uniform areally, the same water volume just being redistributed to different places. Regression of the whole gravity record (1994-2001) on local sea level gives 50–70% of the uniform layer response, as do loading calculations using actual water distributions derived from 11 tide gauges. However, both fits are dominated by some extreme values of short duration, and parts of the gravity record with long-period variations in sea level are close to the uniform layer response. The gravity observations can be used to test corrections for other co-located geodetic observations (GPS, satellite laser ranging) which are influenced by the load effect but not sensitive enough to discriminate between models.     

Interrelationships Between Topography,Quaternary Vertical Displacement,Active Faults and Short-term Horizontal Crustal Strain in Honshu,Japan

Regional and local characteristics of active fault patterns and elevation variation throughout Honshu, Japan are characterized in terms of their fractal dimensions; this allows variation in these complex variables to be compared directly to the scalar properties of net Quaternary vertical displacement, elevation and 10- and 110-year horizontal strains. The comparisons reveal that, throughout Honshu as a whole, there is significant correlation (r=0.75) between Quaternary vertical displacement, elevation, and its fractal properties. There is poor correlation, however, of elevation and its fractal properties to horizontal crustal strain, and also between Quaternary vertical displacement and horizontal crustal strain. A slight negative correlation is observed between the fractal properties of the active fault system and horizontal crustal strain measured over 10- and 110-year time periods (–0.43 and –0.26, respectively). The correlation between the 10-year (1985–1994) and 110-year (1883–1994) area strains, 0.48, reveals the occurrence of considerable change in the distribution of regional strain over these short time frames. Local computations of the correlation between data sets made for overlapping 160 km length windows of data spaced every 20 km along analysis lines reveal internal fluctuations in the correlation between variables. The local correlation between Quaternary vertical displacement and elevation is highest through central Japan and the Kinki Triangle. There is weak negative correlation between area strain and fractal dimensions of the active fault network. The local correlation between the fractal dimensions of active faults and horizontal area strain over the recent 10-year time period averages about –0.6 through central Japan in an area that extends across the Kinki Triangle through the northern part of central Honshu and northeast across the Itoigawa Shizuoka Tectonic Line. In general, regions of greatest complexity in the active fault network are associated with persistent negative area or compressional strain. Sparsely faulted areas in general coincide with areas of positive or roughly zero area strain. The presence of negative correlation through central Japan and the Kinki Triangle area in the recent 10-year period results from a decrease of area strain within an increasingly complex active fault system that reaches maximum negative values concentrated in the Kinki Triangle during the 1985–1994 time period.     

Reduction of gravity observations with digital computer

Summary The paper gives a simple method for the reduction of Gravity observations with electronic computer. The major part of the gravity reductions is the computation of the correction for the far zones. This process was so simplified, that all far zone correction parameters depend on the topography can be determined — independent of the gravity measurement — for the whole survey area in one computer run (subprogram). The computed gravity data will be stored for plotter or other programs.

---

Zusammenfassung Die Arbeit gibt eine einfache Methode für die Reduktion gravimetrischer Meßdaten mit elektronischer Rechenmaschine. Den größten Teil der Reduktion von gravimetrischen Meßdaten machen die Korrektionen für die fernen Zonen aus. Diese Prozedur wurde so vereinfacht, daß man alle — von der Topographie abhängigen — Parameter der Fernzonenkorrektur, unabhängig von den gravimetrischen Me\daten, für das ganze Meßgebiet in einem Rechenprozeß berechnen kann. Die reduzierten gravimetrischen Meßdaten werden für Plotter oder andere Programme eingespeichert.

     

Seismic Sources on the Iberia-African Plate Boundary and their Tectonic Implications

—The plate boundary between Iberia and Africa has been studied using data on seismicity and focal mechanisms. The region has been divided into three areas: A; the Gulf of Cadiz; B, the Betics, Alboran Sea and northern Morocco; and C, Algeria. Seismicity shows a complex behavior, large shallow earthquakes (h < 30 km) occur in areas A and C and moderate shocks in area B; intermediate-depth activity (30 < h < 150 km) is located in area B; the depth earthquakes (h 650 km) are located to the south of Granada. Moment rate, slip velocity and b values have been estimated for shallow shocks, and show similar characteristics for the Gulf of Cadiz and Algeria, and quite different ones for the central region. Focal mechanisms of 80 selected shallow earthquakes (8 m_b 4) show thrust faulting in the Gulf of Cadiz and Algeria with horizontal NNW-SSE compression, and normal faulting in the Alboran Sea with E-W extension. Focal mechanisms of 26 intermediate-depth earthquakes in the Alboran Sea display vertical motions, with a predominant plane trending E-W. Solutions for very deep shocks correspond to vertical dip-slip along N-S trends. Frohlich diagrams and seismic moment tensors show different behavior in the Gulf of Cadiz, Betic-Alboran Sea and northern Morocco, and northern Algeria for shallow events. The stress pattern of intermediate-depth and very deep earthquakes has different directions: vertical extension in the NW-SE direction for intermediate depth earthquakes, and tension and pressure axes dipping about 45 ° for very deep earthquakes. Regional stress pattern may result from the collision between the African plate and Iberia, with extension and subduction of lithospheric material in the Alboran Sea at intermediate depth. The very deep seismicity may be correlated with older subduction processes.