Spatial and temporal variation of gravity field in the capital region

Ｓｐａｔｉａｌａｎｄｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆ　ｇｒａｖｉｔｙ　ｆｉｅｌｄｉｎｔｈｅｃａｐｉｔａｌｒｅｇｉｏｎＣｈａｎｇ－ＣａｉＨＵＡ；（华昌才）ＹｏｎｇＧＵＯ；（果勇）Ｄｕａｎ－ＦａＬＩＵ；（刘瑞法）ＧａｎｇＸＩＡＯ；（肖钢），Ｊ．Ｔ．...     

Comparative study of temporal variations in the earth’s gravity field using GRACE gravity models in the regions of three recent giant earthquakes

Comparative analysis of coseismic and postseismic variations of the Earth’s gravity field is carried for the regions of three giant earthquakes (Andaman-Sumatra, December 26, 2004, magnitude M _w = 9.1; Maule-Chile, February 27, 2010, M _w = 8.8, and Tohoku-Oki, March 11, 2011, M _w = 9.0) with the use of GRACE satellite data. Within the resolution of GRACE models, the coseismic changes of gravity caused by these seismic events manifest themselves by large negative anomalies located in the rear of the subduction zone. The real data are compared with the synthetic anomalies calculated from the rupture surface models based on different kinds of ground measurements. It is shown that the difference between the gravity anomalies corresponding to different rupture surface models exceeds the uncertainties of the GRACE data. There-fore, the coseismic gravity anomalies are at least suitable for rejecting part of the models that are equivalent in the ground data. Within the first few months after the Andaman-Sumatra earthquake, a positive gravity anomaly started to grow above the deep trench. This anomaly rapidly captured the area of the back-arc basin and largely compensated the negative coseismic anomaly. The processes of viscoelastic stress relaxation do not fully allow for these rapid changes of gravity. According to the calculations, even with a sufficiently low viscosity of the upper mantle, relaxation only covers about a half of the observed change of the field. In order to explain the remaining temporal variations, we suggested the process of downdip propagation of the coseismic rupture surface. The feasibility of such a process was supported by numerical simulations. The sum of the gravity anomalies caused by this process and the anomaly generated by the processes of viscoelastic relaxation accounts well for the observed changes of the gravity field in the region of the earthquake. The similar postseismic changes of gravity were also detected for the region of the Tohoku-Oki earthquake. Just as in the case discussed above, this earthquake was also followed by a rapid growth of a positive postseismic anomaly, which partially counterbalanced the negative coseismic anomaly. The time variations of the gravity field in the region of the Maule-Chile earthquake differ from the pattern of changes observed in the island arcs described above. The postseismic gravity variations are in this case concentrated in a narrower band above the deep trench and shelf, and they do not spread over the continental territory, where the negative coseismic anomaly is located. These discrepancies reflect the difference in the geodynamical settings of the studied earthquakes.     

On the gradient of curvature of the plumblines of the Earth’s normal gravity field and its isocurvature lines

This paper presents an approach to determine the gradient of curvature of the normal plumblines at a point P above the ellipsoid and introduces a new geometrical object which is the isocurvature line. The assumed facts are the coordinates of the point P and the formula for the normal gravity potential U. For the determination of the gradient of the normal plumbline curvature k at the point P we define a small circle on the meridian plane of P whose center is at the point P. The circle has the radius of one meter and interior D. In this circle we construct a curvature replacement function to approximate the curvature function k. This replacement function is a quotient of polynomials hence it is easy to find its partial derivatives at the point P. For the construction of replacement function we make the assumption that in the interior of the circle D the first order partial derivatives of U behave linearly and the second order partial derivatives have constant values which equal their value at the point P. Then we set the gradient of the curvature function to be equal with the gradient of the aforementioned replacement function at P. An isocurvature line of the normal gravity field passing through a point P is a curve such that the value of the function of the plumblines’ curvature k is constant and equals k(P). We give a formula to find the direction of the isocurvature line on the meridian plane and we prove that there are infinitely many isocurvature lines passing through the point P and they all lie on a special surface, the isocurvature surface.     

An approach to refined mapping of the anomalous gravity field in the Earth’s polar caps

This work describes the different sets of instruments and methodic approaches for testing the models of gravity anomalies by repeated airborne gravimetric surveys in the polar cap regions of the Earth. The survey design including the specifications for flying the survey profiles and the arrangement of the base stations in polar areas is described, and the necessary modifications of the airborne gravity metering complexes for high-latitude measurements are suggested.     

Far-zone gravity field contributions corrected for the effect of topography by means of molodensky’s truncation coefficients

A spectral representation of the topographic corrections to gravity field quantities is formulated in terms of spherical height functions. When computing the far-zone contributions to the topographic corrections, various types of the truncation coefficients are applied to a spectral representation of Newton’s integral. In this study we utilise Molodensky’s truncation coefficients in deriving the expressions for the far-zone contributions to topographic corrections. The expressions for computing the far-zone gravity field contributions corrected for the effect of topography are then obtained by combining the expressions for the far-zone contributions to the gravity field quantities and to the respective topographic corrections, both expressed in terms of Molodensky’s truncation coefficients. The numerical examples of the far-zone contributions to the topographic corrections and to the topography-corrected gravity field quantities are given over the study area situated in the Canadian Rocky Mountains with adjacent planes. Coefficients of the global elevation and geopotential models are used as the input data.     

Intensity of geomagnetic field in the Precambrian and evolution of the Earth’s deep interior

Reliable data on the paleointensity of the geomagnetic field can become an important source of information both about the mechanisms of generation of the field at present and in the past, and about the internal structure of the Earth, especially the structure and evolution of its core. Unfortunately, the reliability of these data remains a serious problem of paleomagnetic research because of the limitations of experimental methods, and the complexity and diversity of rocks and their magnetic carriers. This is true even for relatively “young” Phanerozoic rocks, but investigation of Precambrian rocks is associated with many additional difficulties. As a consequence, our current knowledge of paleointensity, especially in the Precambrian period, is still very limited. The data limitations do not preclude attempts to use the currently available paleointensity results to analyze the evolution and characteristics of the Earth’s internal structure, such as the age of the Earth’s solid inner core or thermal conductivity in the liquid core. However, such attempts require considerable caution in handling data. In particular, it has now been reliably established that some results on the Precambrian paleointensity overestimate the true paleofield strength. When the paleointensity overestimates are excluded from consideration, the range of the field strength changes in the Precambrian does not exceed the range of its variation in the Phanerozoic. This result calls into question recent assertions that the Earth’s inner core formed in the Mesoproterozoic, about 1.3 billion years ago, triggering a statistically significant increase in the long-term average field strength. Instead, our analysis has shown that the quantity and quality of the currently available data on the Precambrian paleointensity are insufficient to estimate the age of the solid inner core and, therefore, cannot be useful for solving the problem of the thermal conductivity of the Earth’s core. The data are consistent with very young or very “old” inner core ages and, correspondingly, with high or low values of core thermal conductivity.     

Large-scale aseismic creep in the areas of the strong earthquakes revealed from the GRACE data on the time variations of the Earth’s gravity field

The refinement of the accuracy and resolution of the monthly global gravity field models from the GRACE satellite mission, together with the accumulation of more than a decade-long series of these models, enabled us to reveal the processes that occur in the regions of large (M_w≥8) earthquakes that have not been studied previously. The previous research into the time variations of the gravity field in the regions of the giant earthquakes, such as the seismic catastrophes in Sumatra (2004) and Chile (2010), and the Tohoku mega earthquake in Japan (2011), covered the coseismic gravity jump followed by the long postseismic changes reaching almost the same amplitude. The coseismic gravity jumps resulting from the lower-magnitude events are almost unnoticeable. However, we have established a long steady growth of gravity anomalies after a number of such earthquakes. For instance, in the regions of the subduction earthquakes, the growth of the positive gravity anomaly above the oceanic trench was revealed after two events with magnitudes M_w=8.5 in the Sumatra region (the Nias earthquake of March 2005 and the Bengkulu event of September 2007 near the southern termination of Sumatra Island), after the earthquake with M_w=8.5 on Hokkaido in September 2007, a doublet Simushir earthquake with the magnitudes M_w = 8.3 and 8.1 in the Kuriles in November 2006 and January 2007, and after the earthquake off the Samoa Island in September 2009 (M_w=8.1). The steady changes in the gravity field have also been recorded after the earthquake in the Sichuan region (May 2008, M_w = 8.0) and after the doublet event with magnitudes 8.6 and 8.2, which occurred in the Wharton Basin of the Indian Ocean on April 11, 2012. The detailed analysis of the growth of the positive anomaly in gravity after the Simushir earthquake of November 2006 is presented. The growth started a few months after the event synchronously with the seismic activation on the downdip extension of the coseismically ruptured fault plane zone. The data demonstrating the increasing depth of the aftershocks since March 2007 and the approximately simultaneous change in the direction and average velocity of the horizontal surface displacements at the sites of the regional GPS network indicate that this earthquake induced postseismic displacements in a huge area extending to depths below 100 km. The total displacement since the beginning of the growth of the gravity anomaly up to July 2012 is estimated at 3.0 m in the upper part of the plate’s contact and 1.5 m in the lower part up to a depth of 100 km. With allowance for the size of the region captured by the deformations, the released total energy is equivalent to the earthquake with the magnitude M_w = 8.5. In our opinion, the growth of the gravity anomaly in these regions indicates a large-scale aseismic creep over the areas much more extensive than the source zone of the earthquake. These processes have not been previously revealed by the ground-based techniques. Hence, the time series of the GRACE gravity models are an important source of the new data about the locations and evolution of the locked segments of the subduction zones and their seismic potential.     

Far-zone contributions to the gravity field quantities by means of Molodensky’s truncation coefficients

To reduce the numerical complexity of inverse solutions to large systems of discretised integral equations in gravimetric geoid/quasigeoid modelling, the surface domain of Green’s integrals is subdivided into the near-zone and far-zone integration sub-domains. The inversion is performed for the near zone using regional detailed gravity data. The farzone contributions to the gravity field quantities are estimated from an available global geopotential model using techniques for a spherical harmonic analysis of the gravity field. For computing the far-zone contributions by means of Green’s integrals, truncation coefficients are applied. Different forms of truncation coefficients have been derived depending on a type of integrals in solving various geodetic boundary-value problems. In this study, we utilise Molodensky’s truncation coefficients to Green’s integrals for computing the far-zone contributions to the disturbing potential, the gravity disturbance, and the gravity anomaly. We also demonstrate that Molodensky’s truncation coefficients can be uniformly applied to all types of Green’s integrals used in solving the boundaryvalue problems. The numerical example of the far-zone contributions to the gravity field quantities is given over the area of study which comprises the Canadian Rocky Mountains. The coefficients of a global geopotential model and a detailed digital terrain model are used as input data.     

Prospects for the Development and Application of the Earthquake Monitoring Network

With the rapid development of the economy in China, the seismic network has been changing rapidly, in that the capability of instruments, technological systems and network density are approaching those of developed countries and a large quantity of observation data has been accumulated. How to apply these resources to economic construction and public safety has become an important issue worth studying. In order to improve earthquake prediction and earthquake emergency response, it is suggested in this paper that extracting valuable precursor information, improving earthquake rapid reporting ability and extending rapid intensity reporting function are key issues. Integrating network resources, building unified standards and a multifunction seismic monitoring network are preconditions of establishing a public safety service platform and earthquake observation resources will contribute significantly to the fields of engineering, ocean, meteorology, and environmental protection. Thus, the future directions of the development of the seismic network are exploring monitoring resources, enhancing independent innovation, constructing a technological platform and enlarging the service field.     

The contribution of seismic data in microzonation studies for downtown L��Aquila

After the 2009 April 6th Mw 6.3 L??Aquila earthquake (Central Italy) the Italian Civil Defense Department promoted the microzoning study in the ten zones in the epicentral area that suffered major damage. In this paper we present the activities and the results concerning a temporary seismic network installed in the historical L??Aquila city center indicated as ??macroarea 1?? in the microzoning project. Seismic data were collected to investigate the amplification effects in the city and to support the microzoning activities in verifying both geological profiles and 1D numerical modeling of the seismic response of the city. The conventional spectral approaches using both microtremor and earthquake data allowed to determine the fundamental resonance frequencies and the amplification factors within the city respectively. The spatial variability of these quantities can be related to the geological and geomorphologic characteristics of the investigated area. A comparison between the network data and the data recorded by the two strong motion instruments installed in the city was also made. This allows verifying the relative response of the accelerometric stations that recorded in the city the major events of the sequence.     

Computation and analysis of the geomagnetic field model in China and its adjacent area for 2003

Based on the geomagnetic data at 135 stations and 35 observatories in China in 2003, the Taylor polynomial model and the spherical cap harmonic model in China and its adjacent area for 2003 were established. In the model calculation, the truncation order of the model and the influences of the boundary restriction on the model calculation were carefully analyzed. The results show that the geomagnetic data used are precise and reliable, and the selection of the truncation order is reasonable. The Taylor polynomial model and the spherical cap harmonic model in China and its adjacent area established in this paper are consistent very well.     

Incentives for field hydrology and data sharing: collaboration and compensation: reply to “A need for incentivizing field hydrology,especially in an era of open data”*

We thank Allen and Berghuijs for continuing the discussion on field hydrology and data sharing and discuss two incentives to promote data collection and sharing in hydrological sciences: a collaborative attitude and additional funding to make data publicly available.     

Harmonic sources of the main geomagnetic field in the Earth’s core

The large-scale harmonic magnetic-convective sources of the main geomagnetic field in the Earth’s core have been determined for the first time. The determination is based on a complete system of eigenfunctions of the magnetic diffusion equation in a homogeneously conducting sphere, which is surrounded by an insulator. The sources of the main geomagnetic field observed, which is responsible for the distribution of the electric currents generating this field in the core, are expressed in terms of large-scale eigenfunctions. In this case, the dipole sources are directly related to the observed geomagnetic dipole, whereas the quadrupole sources are related to the quadrupole, etc. The time variations in the obtained sources are responsible for individual spatiotemporal features in the generation or suppression of each Gaussian component of the observed geomagnetic field. When the commonly accepted observational international geomagnetic reference field (IGRF) models were used to partially reveal these time variations, it became possible to specify the estimate of the Earth’s core conductivity and determine the minimum period that can separate us from the commencement of further inversion or excursion.     

Inversion of gravity and topography data for the crust thickness of China and its adjacent region

Introduction Since the middle of the century, gravitational isostasy has been a fundamental hypothesis for inverting the gravity data to find the crust thickness. Geophysicists have done a lot of researches on using gravity data to investigate the depth of Moho discontinuity. Since 1980, the International Lithosphere Program emphasized the importance of investigating the Moho depth variation. Thereafter a lot of results have been published in the world (Braitenberg et al, 2000; Kaban et al,…     

The effect of wind on the gravity wave propagation in the Earth’s ionosphere

The effect of ionospheric wind on the gravity wave propagation is studied. These waves arise in the ionosphere due to intensification of their sources near the Earth’s surface during enhanced seismic activity. The influence of the wind on these waves is connected with the Ampere’s force that produces the ion-drag force acting on the atmosphere. This results in the occurrence of the discrete wave spectrum the maximum of which increases in proportion to the numbers of the natural scale. Furthermore, these waves are amplified during propagation from the source region in the direction perpendicular to the wind direction. These peculiarities of the gravity waves can be used for monitoring of seismic activity based on the ionosphere sounding.     

Variation of the intensity of the Earth’s magnetic field in Portugal in the 1st millennium BC

The study of magnetization of the ceramic material from 21 archeological monuments of Portugal (the Evora province), dated archeologically from the Bronze Age to the end of the Iron Age has been carried out. For the purpose of more detailed timing of the material from the monuments the method of ceramic age dating on the basis of its porosity has been used. In order to take into account the distorted factors in the determination of the parameters of the ancient geomagnetic field with the aim of the maximal approximation to the actual values the diagnostic features of magnetite weathering have been considered and the level of weathering of the magnetic fraction in the ceramics from archeological monuments has been determined. The data of geomagnetic field-strength variation in the time interval of the 12th century BC to the beginning of the Common Era have been obtained. The field-strength at this time interval varied in the range of 60–90 micro Tesla with the maximal values in the 9th, 8th, and the second half of the 5th to the beginning of the 4th century BC. In addition, the timing of the ceramic material from the urns of the megalithic complex Monte de Tera of the Evora province has been clarified.     

Estimation of the temperature in the Earth’s interior from surface measurements of the electromagnetic field

The possibility of contactless remote estimation of the temperature in the Earth’s interior from surface magnetotelluric (MT) measurements is examined. The neuronet analysis of MT and temperature measurements in the Bishkek geodynamic research area (the Northern Tien Shan) showed that a contactless electromagnetic geothermometer can in principle be realized. An optimal method including MT measurements and treatment of available thermograms is developed. The method minimizes uncertainties of the remote temperature estimation. The use of six to eight thermograms for calibration of electromagnetic data is shown to provide a 12% relative error of prediction, and a priori geological information available for the region under study can reduce this error. Areas of practical application of a contactless electromagnetic geothermometer are outlined.