Identification of Loriya Fault,Its Reactivation Due to 2001 Bhuj Earthquake, (Using Remote Sensing Data) and Its Bearing on the Kachchh Mainland Fault,Kachchh District,Gujarat

The 2001 Bhuj earthquake (M_w 7.7), one of the most severe earthquakes in the recent history of India, reactivated various existing active faults. It is manifested in the form of coseismic ground fissures/cracks and upheaval of land in the form of bumps. Identification and reactivation of Loriya Fault is established by 1—Geomorphic changes with the help of digital imagery (LISS III images). 2—Coseismic changes through ground checks and 3—Geophysical signatures through magnetic and gravity survey. A lineament cutting the north-western part of the Pur River alluvial fan has been revealed by satellite imagery. The streams flowing along the lineament add to the evidences of a weak plane, while the occurrence of coseismic ground fissures confirms the existence of an active fault. No dip slip movement is recorded in the trenches made across the Loriya active fault while the en-echelon pattern of ground fissures suggest strike slip movement along the fault due to 2001 earthquake.     

Gravity changes before large earthquakes in China: 1998–2005

Although it is well known that coseismic gravity changes take place during an earthquake, previous research has not yielded convincing evidence demonstrating that significant gravity changes occur before large earthquakes. Furthermore, even if we suspect that gravity changes occur before large earthquakes, we have yet to demonstrate how to consistently observe these changes for useful earthquake forecast that would bring benefits to society. We analyzed ground gravity survey data obtained in 1998, 2000, 2002, and 2005 at stations of the Crustal Movement Observation Network of China (CMONOC) and examined gravity changes be-fore the occurrence of nine large (Ms≥6.8) earthquakes that ruptured within or near mainland China and Taiwan from November 2001 to August 2008. Results from this analysis show that significant gravity changes occurred across a large region before each of these nine large earthquakes, and these changes were detected by repeated ground gravity surveys through CMONOC. Although these gravity changes were significant, more research is needed to investigate whether these gravity changes could be viewed as precursors of large earthquakes. Limitations and uncertainties in the data include sparseness of the gravity monitoring network, long time intervals between consecutive gravity surveys, inevitable measurement errors, hydrological effects on gravity, and effects of vertical crustal movements on gravity. Based on these observations, we make several recommendations about possible future direc-tions in earthquake-related research using gravity monitoring data.     

TEC variation over an equatorial and anomaly crest region in India during 2012 and 2013

Using the global positioning system (GPS) measurements, the total electron content (TEC) at station Bangalore (13.02°N, 77.57°E geographic; 04.44°N, 150.84°E geomagnetic), lying at the equatorial region, and station Lucknow (26.91°N, 80.95°E geographic; 17.96°N, 155.24°E geomagnetic), lying at equatorial ionospheric anomaly (EIA) crest region, have been estimated for the year 2012–2013. In order to evaluate the International Reference Ionosphere (IRI) model regarding simulation/modeling of ionospheric studies specially at equatorial and EIA crest regions, we have compared the TEC derived from the recent version of the IRI-2012 model and the older IRI-2007 with its three topside options, namely IRI-NeQuick (IRI-NeQ), IRI-2001, and IRI01-corr, with that of GPS-TEC over Bangalore and Lucknow. For the EIA station Lucknow, the IRI-2012 model with IRI-NeQ and IRI01-corr topside is found in good agreement with GPS-TEC during summer and equinox season, while the IRI-2012 model for all three topside options significantly overestimates the GPS-TEC during winter season. The IRI-2001 topside overestimates the GPS-TEC over both the stations during all seasons. The anomalous difference between the IRI-2012 model prediction and ground-based GPS-TEC in daytime hours during the winter season observed at Lucknow could be attributed to discrepancies in the slab thickness predicted by the model, which is more during the winter season as compared to summer and equinox. These large discrepancies in the slab thickness predicted by the IRI-2012 as well as the IRI-2007 model during the winter season have been supported by using the foF2 data from Constellation Observing System for Meteorology, Ionosphere, and Climate radio occultation-based measurements. We also observed that the discrepancies in the recent IRI-2012 model with respect to GPS-TEC are found to be slightly larger than those with the older IRI-2007 model over the EIA region Lucknow. However, over the equatorial region Bangalore, the discrepancy with the older model IRI-2007 was found to be larger than with the recent IRI-2012 model. This suggests that the performance of the IRI-2012 model is poorer than the IRI-2007 model at the EIA region while better at equatorial region, and that further improvements in the IRI-2012 models are required particularly in the low-latitude and EIA regions. The GPS-TEC showed disappearance of the winter anomaly during 2012–2013, while the IRI model failed to predict the disappearance of winter anomaly.     

中国西太平洋海域1′×1′垂线偏差模型及精度评估

联合多种测高数据和重力异常数据,设计了观测点距离和测高精度融合的定权方法,采用最小二乘方法和Vening-Meinesz公式,分别构建了西太平洋海域(0°N—40°N,105°E—145°E)1′×1′网格化垂线偏差数字模型。选取两个不同特征区域,将垂线偏差的两个数字模型和EGM2008模型三者进行相互比较分析。结果表明:卯酉分量η的均方根差大于子午分量ξ的均方根差,海底地形复杂的南海特征区域的垂线偏差均方根差大于西太平洋中部的均方根差,构建的两个垂线偏差模型总体均方根差优于1.6″。     

Morphostructural prognostication of kimberlites in parts of eastern dharwar craton: Inferences from remote sensing and gravity signatures

The varied nature of geophysical signatures encountered in the identification of blocks with favorable setting for emplacement of kimberlites, the major geological hosts of diamonds, has resulted in the development of different criteria for their prognostication. This paper attempts to bring out the association of domal morpho-structural features with kimberlite occurrences. From analysis of an IRS-ID LISS-III satellite image of the Maddur-Narayanpet-Mahboobnagar-Gadwal-Raichur region (latitudes 16° N to 17°15’ N and longitudes 77° E to 78°15’ E) in the eastern Dharwar craton, the surface morphology and structural configuration were obtained. From observed surface correlation of reported kimberlite occurrences with intersecting lineaments within regions of domal morphology, regions with such configuration were delineated as potential kimberlite bearing zones. Using this observation as a criterion, six potential kimberlite zones were delineated: near Sedam, around Maddur, north of Mahbubnagar, north of Raichur, between Raichur and Gadwal and around Gadwal. Subsequently, from reanalysis of available gravity data in a part of the Maddur area, the structural features inferred from image analysis were corroborated. The crustal configuration along a representative profile across a region associated with kimberlite emplacement was obtained to understand the subsurface pattern of its occurrence. Occurrences of kimberlites were found to be localized at intersections of lineaments in regions characterized by upwarps in the Moho and corresponding lows in the gneissic basement.     

Analysis on the global morphology of stratospheric gravity wave activity deduced from the COSMIC GPS occultation profiles

Monthly mean global morphologies of potential energy density E _p from stratospheric gravity waves are revealed by observations of COSMIC GPS radio occultation. The E _p is obtained from vertical wavelengths ranging from 2 to 10 km over cells of 1° × 2° in latitude and longitude. The computed values confirm previous results and obtain new ones. The large gravity wave E _p values found in the tropics between 25°N and 25°S could be mainly due to the strong tropical cumulus convection; July values are larger than those for January (2007). In mid and high latitudes, the most prominent features of the northern winter hemisphere are the enhanced densities above the Eurasian continent and the North Atlantic and the depressed E _p values above the North Pacific and North America for which topography, wind sources and wind filtering may be responsible. In southern winter hemisphere, large E _p values are found around 180° and 300° longitudes that are likely due to the topography of the Antarctic plateau, the Antarctic Peninsula and South America. Enhanced E _p values are found over Scandinavia. However, there is no clear evidence to show that gravity waves are localized over the Rocky Mountains, the Himalayas and the Andes. Topography and planetary wave modulations are proposed to interpret the large-scale longitudinal variations and inter-hemisphere asymmetry of the GW activity.     

Tectonic plate motion and co-seismic steps surveyed by DORIS field beacons: the New Hebrides experiment

 The New Hebrides experiment consisted of setting up a pair of DORIS beacons in remote tropical islands in the southwestern Pacific, between 1993 and 1997. Because of orbitography requirements on TOPEX/Poséidon, the beacons were only transmitting to SPOT satellites. Root-mean-square (RMS) scatters at the centimeter level on the latitude and vertical components were achieved, but 2-cm RMS scatters affected the longitude component. Nevertheless, results of relative velocity (123 mm/year N250°) are very consistent with those obtained using the global positioning system (GPS) (126 mm/yr N246°). The co-seismic step (12 mm N60°) related to the Walpole event (M _W = 7.7) is consistent with that derived from GPS (10 mm N30°) or from the centroid moment tensor (CMT) of the quake (12 mm N000°). Received: 19 November 1999 / Accepted: 17 May 2000     

Discussion on some FFT problems to determine the geoid

The aim of this investigation is to study some FFT problems related to the application of FFT to gravity field convolution integrals. And the others, such as the effect of spectral leakage, edge effects, cyclic convolution and effect of padding, are also discussed. A numerical test for these problems is made. A large area of Western China selected for the test is located between 30°N～36°N and 96°E～102°E and includes 1 858 gravity observations on land. The results show that the removal of the bias in the residual gravity anomalies is important to avoid spectral leakage. One hundred percent zero padding is highly recommended for further research of the geoid to remove cyclic convolution errors and edge effects. 1-D FFT is recommended for precise local geoid determination because it does not use kernel approximation.     

National height datum,the Gauss–Listing geoid level value w0 and its time variation 0 (Baltic Sea Level Project: epochs 1990.8, 1993.8, 1997.4)

 A methodology for precise determination of the fundamental geodetic parameter w ₀, the potential value of the Gauss–Listing geoid, as well as its time derivative ₀, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w ₀ and ₀ values (w ₀=62 636 855.75 ± 0.21 m²/s², ₀=−0.0099±0.00079 m²/s² per year, w ₀/&γmacr;=6 379 781.502 m,₀/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m²/s²) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1) the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums of countries around the Baltic Sea. Received: 14 August 2000 / Accepted: 19 June 2001     

Shallow-water loading tides in Japan from superconducting gravimetry

Gravity observations from superconducting gravimeters are used to observe loading effects from shallow-water tides on the Japanese east and west coasts. Non-linear third-diurnal and higher-frequency shallow-water tides are identified in the tide-gauge observations from these coastal areas. The most energetic constituents in the tide gauge observations are also seen in the gravity observations due to their loading effects on the deformation of the Earth. Even though the shallow-water tides at the Japanese east coast have an amplitude of only a few millimetres, they are still able to generate a loading signal at gravity sites located several hundred kilometres inland. In particular, the S₃, S₄ and S₅ solar tides occur in both gravity and tide gauge observations. It is indicated that in other shelf regions with large shallow water tides, the shallow water loading signals account for a significant signal, which should be taken into account.Acknowledgement The authors would like to thank the Hydrographic and Oceanographic Department (Japan Coast Guard), Japan Meteorological Agency and Hokkaido Development Agency for access to the tide-gauge data. Also, the Global Geodynamic Project Information System and Data Center (GGP-ISDC) is acknowledged for providing the gravity data.     

星载原子干涉技术用于地球重力场测量及其精度评估

重力梯度卫星GOCE通过搭载静电式重力梯度仪,将全球静态重力场恢复至200阶以上。目前GOCE卫星已结束寿命,亟须发展下一代更高分辨率的卫星重力梯度测量来完善200~360阶的全球静态重力场模型。原子干涉型的重力梯度测量在空间微重力环境下可获得较长的干涉时间,因此具有很高的星载测量精度,是下一代卫星重力梯度测量的候选技术之一。本文针对未来更高分辨率全球重力场测量的科学需求,提出了一种适用于空间微重力环境下的原子干涉重力梯度测量方案,其梯度测量噪声可低至0.85mE/Hz1/2。文中对不同类型的卫星重力梯度测量方案进行了重力场反演精度的对比评估,仿真结果表明,相比于现有静电式卫星重力梯度测量,原子干涉型的卫星重力梯度测量有望将重力场的恢复阶数提升至252~290阶,对应的累积大地水准面误差7~8cm,累积重力异常误差3×10-5 m/s2。     

Comparison of undulation difference accuracies using gravity anomalies and gravity disturbances

Errors are considered in the outer zone contribution to oceanic undulation differences as obtained from a set of potential coefficients complete to degree 180. It is assumed that the gravity data of the inner zone (a spherical cap), consisting of either gravity anomalies or gravity disturbances, has negligible error. This implies that error estimates of the total undulation difference are analyzed. If the potential coefficients are derived from a global field of 1°×1° mean anomalies accurate to ε_Δg=10 mgal, then for a cap radius of 10°, the undulation difference error (for separations between 100 km and 2000 km) ranges from 13 cm to 55 cm in the gravity anomaly case and from 6 cm to 36 cm in the gravity disturbance case. If ε_Δg is reduced to 1 mgal, these errors in both cases are less than 10 cm. In the absence of a spherical cap, both cases yield identical error estimates: about 68 cm if ε_Δg=1 mgal (for most separations) and ranging from 93 cm to 160 cm if ε_Δg=10 mgal. Introducing a perfect 30-degree reference field, the latter errors are reduced to about 110 cm for most separations.     

Numerical study on the mixed model in the GOCE polar gap problem

Gravity gradients acquired by the Gravity field and steady-state Ocean Circulation Explorer(GOCE) do not cover the entire earth because of its sun-synchronous orbit leaving data gaps with a radius of about 6.5° in the polar regions.Previous studies showed that the loss of data in the polar regions deteriorates the accuracy of the low order(or near zonal) coefficients of the earth gravity model,which is the so-called polar gap problem in geodesy.In order to find a stable solution for the earth gravity model from the GOCE gravity gradients,three models,i.e.the Gauss-Markov model,light constraint model and the mixed model,are compared and evaluated numerically with the gravity gradient simulated with the EGM2008.The comparison shows that the Best Linear Uniformly Unbiased Estimation(BLUUE) estimator of the mixed model can solve the polar gap problem as effectively as the light constraint model;furthermore,the mixed model is more rigorous in dealing with the supplementary information and leads to a better accuracy in determining the global geoid.     

超导重力数据检测到的2011年日本东北大地震(Mw 9.0)震前重力异常及同震重力变化

高精度的超导重力数据已广泛应用于地球动力学的研究,对大地震前重力异常和同震重力变化的探测有助于震源机制和地震预警的研究,同时高精度的同震观测数据可用于断层滑动模型的反演。本文利用日本、中国及欧洲7个超导台站2011年3月的秒采样数据研究日本东北大地震（M_w 9.0）产生的重力变化,经潮汐、气压、漂移等预处理改正后,得到改正后的重力变化。对比同时段日本岛附近发生的2861次M_b≥4级地震,分析滤波后的重力变化数据,排除非主震的影响,发现0.12 Hz≤f≤0.18 Hz频段显示了明显的震前重力异常扰动信息,所有台站在震前89 h均开始出现重力异常扰动,Medicina站的震前最大异常扰动振幅达到28×10^-8 m/s²。另外,利用球形位错理论及地震的CMT解计算得到的地表同震重力变化理论值与超导重力观测值非常接近,同震重力变化绝对值与震中距成反比,近场站点的观测值比远场站点更接近于理论计算值,前者可作为断层反演的约束条件。     

Geomorphic interpretation of the Landsat images in Tamilnadu

Landsat imagery were used to interpret the geomorphic festures in the northern part of Tamilnadu between latitudes 11° and 12°30′ N and longitudes 77° and 80° E. Different geomorphic regions were identified with the aid of ground truths. These geomorphic regions, in general, are controlled by the underlying rock types. As the different geomorphic regions have different signatures in the Landsat images, the geomorphic interpretation, which helps to bring out the regional geology is possible with the aid of ground truths.     

Precise tidal gravity recorded with superconducting gravimeters at stations Wuhan (China) and Kyoto (Japan)

 Three long series of tidal gravity observations, totalizing approximately 24 years and recorded with three superconducting gravimeters, T004, T008, and T009, at stations Wuhan (China) and Kyoto (Japan), are studied. The tidal amplitude factors and phase differences are determined precisely using Eterna and Nsv techniques. The precision of the main tidal amplitudes is at the same level of 0.01 μGal. The atmospheric gravity signals are corrected using the coefficients determined with a regression method between tidal gravity residual and station air pressure. The oceanic gravity signals are modeled based on five global oceanic models. It is found that the oceanic models developed by the analysis of measurements from Topex/Poseidon altimeters have the best fit to the superconducting gravimeter measurements, since the observed residuals and the discrepancies between the amplitude factors and the theoretical tidal models are reduced more significantly. The long-period gravity variations are dominated by the non-linear drift phenomena of the instruments, and the short-term variations in gravity are due to the background noise at the stations. Received: 20 January 2000 / Accepted: 15 September 2000     

TRIGONOMETRICAL LEVELLING ACROSS THE INLAND ICE IN NORTH GREENLAND

Abstract

As part of the scientific work of the British North Greenland Expedition (1952–1954), a programme of trigonometrical levelling was carried out from the east to the west coast of Greenland, along a line across the inland ice between latitudes 76° 40′ N., and 78° 10′ N. The primary purpose of the work was to determine accurately the heights above sea level of a series of gravity stations, the gravity measurements being made in connection with determinations of ice thickness. For meteorological purposes it was necessary to know also the altitude of the Expedition's central station, situated in latitude 78° 04′ N., longitude 38° 29′ W. The accuracy necessary for the purpose of the gravity survey was a few metres for the altitudes, while the latitude of each gravity station had to be determined with an accuracy of ± 0.1 minute.     

Improvements on existing geopotential coefficients models for precise determination of 1° global geoid and mean gravity anomalies

The concept of an idealised earth having 1° averaged heights over its land surface is introduced as a means to improve upon the existing geopotential coefficient solutions without the use of additional observed data, in order to provide more precise knowledge of the earth’s gravity field in the form of 1° global geoid and 1° mean free-air gravity anomalies especially over the mountainous regions with the visible topography condensed into the actual geoid, first by referring them to the idealised earth and then by reducing the same to the actual earth on applying appropriate corrections for the differences between the two earths.     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

A new gravity map, a new marine geoid around Japan and the detection of the Kuroshio current

About half a million marine gravity measurements over a 30^∘×30^∘ area centered on Japan have been processed and adjusted to produce a new free-air gravity map from a 5′×5′ grid. This map seems to have a better resolution than those previously published as measured by its correlation with bathymetry. The grid was used together with a high-degree and -order spherical harmonics geopotential model to compute a detailed geoid with two methods: Stokes integral and collocation. Comparisons with other available geoidal surfaces derived either from gravity or from satellite altimetry were made especially to test the ability of this new geoid at showing the sea surface topography as mapped by the Topex/Poseidon satellite. Over 2 months (6 cycles) the dynamic topography at ascending passes in the region (23^∘47^∘N and 123^∘147^∘E) was mapped to study the variability of the Kuroshio current. Received: 15 July 1994 / Accepted: 17 February 1997