2011年日本东北地区太平洋近海地震对亚洲板块及韩国大地控制网的影响分析

本研究分析了2011年3月11日发生的M_w9.0日本东北地区太平洋近海地震对亚洲地区和韩国国内GPS卫星常年跟踪站的位移影响.为此,利用了日本东北地区太平洋近海地震发生前后两周(2011年3月4日到3月18日)的GPS站点数据,包括震中附近地区(韩国,中国,中国台湾地区,日本和俄罗斯)55个GPS卫星常年跟踪站和284个IGS 全球跟踪站,并采用GAMIT/GLOBK软件进行处理和平差,估算出所有GPS站点的同震形变.结果显示,日本东北地区太平洋近海地震引起的同震形变影响在亚洲地区比较明显,包括日本和附近国家,距离震中2702 km的中国武汉(WUHN)站也观测到同震形变.为精确分析日本东北地区太平洋近海地震对韩国国家大地控制网的影响,通过GAMIT/GLOBK软件计算出韩国GPS卫星常年跟踪站之间的基线长度变形,并分析出弹性变形量.结果表明:大部分GPS站点均向震中方向膨胀,且向震中的垂直方向收缩.由日本东北地区太平洋近海地震导致的最大剪应变达到韩国国家大地控制网年均变形率的约7倍,对韩国的地壳产生14.5~57.7 mm的水平位移,并导致韩国国家大地控制网产生弹性变形.因此,在不及时更新维护韩国国家大地控制网的情况下,GPS测量成果将会发生最大20 mm的位置误差.     

Relation of surface movements in West Bohemia to earthquake swarms

GPS observations in the Western Bohemia/Vogtland earthquake swarm region revealed indications of horizontal displacements of low amplitude, and no clear long-term trend in 1993–2007. On the other hand, in 1998–2001 there was relatively significant active movement along NNE-SSW oriented line that we called the “Cheb-Kraslice GPS Boundary” (ChKB), identical with an important limitation of earthquake activity. The most impressive were dextral (right-lateral) movements in the 1998–1999 period followed by reverse sinistral (left-lateral) movements in 1999–2000 that correlate with prevailing motion defined by fault plane solutions of the Autumn 2000 earthquake swarm. Before the February 2004 micro-swarm, two points located on opposite sides of the Mariánské Lázně fault showed extension in the order of about 7 mm in the same NNE-SSW direction of ChKB. The new NOKO permanent GPS station in Novy Kostel showed the peak-to-peak vertical changes up to 10 mm before and during the February 2007 micro-swarm. Annual precise levelling campaigns in the local network around Novy Kostel revealed regular vertical displacements during the 1994, 1997 and 2000 earthquake swarms. The points around the Novy Kostel seismological station showed uplift during the active periods, including the micro-swarm February 2004. However, no such indication was observed on levelling points in the period of the February 2007 swarm. Long-term vertical displacements depend on the same direction NNE-SSW (ChKB) as the GPS displacements. Both geodetic techniques have revealed oscillating displacements, GPS horizontal, and levelling vertical, rather than any long-term trends in the study period 1993–2007. The displacements exhibited significant spatial and temporal relation to tectonic activity (earthquake swarms) including their coincidence with the seismologically determined sense of motion along the fault plane during earthquakes.     

Approximate recovery of coseismic deformation from Taiwan strong-motion records

Since 1990, digital strong-motion accelerographs and global positioning system (GPS) instruments have been widely deployed in the Taiwan region (Shin et al. 2003; Yu et al. 2001). The 1999 Chi-Chi, Mw 7.6 earthquake and the 2003 Chengkung, Mw 6.8 earthquake were well recorded by both digital accelerographs and GPS instruments. These data offer a good opportunity to determine coseismic displacements from strong-motion records and to compare the results with those derived from GPS measurements. As noted by Boore (2001), a double integration of the acceleration data often leads to unreasonable results, and baseline corrections are therefore required in most cases before the integration. Based on the works of Iwan et al. (1985) and Boore (2001), we developed an improved method for baseline correction and validated it using an extensive set of data from shake-table tests of a known “step” displacement on 249 accelerographs. Our baseline correction method recovered about 97% of the actual displacement from the shake-table data. We then applied this baseline correction method to compute coseismic displacements from the strong-motion data of the Chi-Chi and Chengkung earthquakes. Our results agree favorably with the coseismic displacements determined by the GPS measurements at nearby sites. The ratio of seismic to geodetic displacement varies from 0.78 to 1.41, with an average of about 1.05.     

Non-seismological indications of recent tectonic activity in the West Bohemia earthquake swarm region

The West Bohemia seismoactive region is characterised by repeated occurrence of earthquake swarms with the strongest event of magnitude M=4.5 recorded during the 1985/1986 swarm. In the 1990s, complex observations of some geodynamic parameters have been introduced. Special attention was paid to geodetic measurements, using the techniques of precise levelling and GPS. In the most active focal zone of the region, vertical movements were distributed with an organised pattern during swarm periods but randomly during quiet periods. The GPS data have not shown any systematic trend of horizontal displacements yet. Consistent compression-like indications were obtained between 1998 and 1999 campaigns, followed by a reverse displacement in the subsequent 2000 swarm period. Temporal gravity changes were analysed from a set of data recorded twice a year since 1993. During the last three swarm periods in December 1994, January 1997 and Autumn 2000, significant temporal variations were observed in the most active focal zone of the region. On the contrary, during seismically quiet periods only weak signals of gravity changes were recorded. The variations probably reflected strain evolution in the upper crust with possible effect of fluid migration and mass displacement on gravity. The groundwater level has been recorded in 2 wells. The water level changes differ significantly in their characteristics. Water level record of the first well provides a smooth curve with a noticeable minimum coinciding exactly with the December 1994 swarm and with some changes before the January 1997 and after the Autumn 2000 swarms. The water level record of the other well reflects significant earth tide and air pressure effects. After a complex analysis, sharp residual signals indicated the first strong, macroseismically recorded, events of the Autumn 2000 swarm.     

GPS-derived Deformation of the Central Andes Including the 1995 Antofagasta Mw = 8.0 Earthquake

—In order to study both the interplate seismic loading cycle and the distribution of intraplate deformation of the Andes, a 215 site GPS network covering Chile and the western part of Argentina was selected, monumented and observed in 1993 and 1994. A dense part of the network in northern Chile and northwest Argentina, comprising some 70 sites, was re-observed after two years in October/November, 1995. The M _w = 8.0 Antofagasta (North Chile) earthquake of 30th July, 1995 took place between the two observations. The city of Antofagasta shifted 80 cm westwards by this event and the displacement still reached 10 cm at locations 300 km from the trench. Three different deformation processes have been considered for modeling the measured displacements (1) interseismic accumulation of elastic strain due to subduction coupling, (2) coseismic strain release during the Antofagasta earthquake and (3) crustal shortening in the Sub-Andes.¶Eastward displacement of the sites to the north and to the south of the area affected by the earthquake is due to the interseismic accumulation of elastic deformation. Assuming a uniform slip model of interseismic coupling, the observed displacements at the coast require a fully locked subduction interface and a depth of seismic coupling of 50 km. The geodetically derived fault plane parameters of the Antofagasta earthquake are consistent with results derived from wave-form modeling of seismolog ical data. The coseismic slip predicted by the variable slip model reaches values of 3.2 m in the dip-slip and 1.4 m in the strike-slip directions. The derived rake is 66°. Our geodetic results suggest that the oblique Nazca–South American plate convergence is accommodated by oblique earthquake slip with no slip partitioning. The observed displacements in the back-arc indicate a present-day crustal shortening rate of 3–4 mm/year which is significantly slower than the average of 10 mm/year experienced during the evolution of the Andean plateau.     

Simultaneous estimation of secular and episodic crustal motion via geodetic observations

Repeated high quality geodetic observations allow the estimation of the free surface velocity field for a region. Usually, a yearly secular rate is estimated, while the possible episodic motion (seismic slip) due to an earthquake is evaluated via inversion of the geodetic data. The episodic motion influence is subtracted from the total field. The actual region and the respective geodetic observations affected by the seismic event are often assessed by rather vague criteria.This paper deals with an attempt to estimate simultaneously the secular and the episodic two-dimensional crustal motion for a region by means of repeated GPS observations of a geodetic network, carried out for a number of epochs.The model is applied to the Gulf of Corinth available GPS data, spanning a time interval of more than a decade, and the results are discussed.     

Analysis of Coseismic Fault Slip Models of the 2012 Indian Ocean Earthquake: Importance of GPS Data for Crustal Deformation Studies

Based on continuous GPS data, we analyze coseismic deformation due to the 2012 Indian Ocean earthquake. We use the available coseismic slip models of the 2012 earthquake, derived from geodetic and/or seismic waveform inversion, to calculate the coseismic displacements in the Andaman-Nicobar, Sumatra and Java. In our analysis, we employ a spherical, layered model of the Earth and we find that Java Island experienced coseismic displacements up to 8 mm, as also observed by our GPS network. Compared to coseismic offsets measured from GPS data, a coseismic slip model derived from multiple observations produced better results than a model based on a single type of observation.     

Plate boundary deformation and continuing deflation of the Askja volcano,North Iceland,determined with GPS, 1987–1993

GPS geodetic measurements were conducted around the Askja central volcano located at the divergent plate boundary in north Iceland in 1987, 1990, 1992 and 1993. The accuracy of the 1987 and 1990 measurements is in the range of 10 mm for horizontal components; the accuracy of the 1992 and 1993 measurements is about 4 mm in the horizontal plane. Regional deformation in the Askja region is dominated by extension. Points located outside a 30–45 km wide plate boundary deformation zone indicate a displacement of 2.4±0.5 cm/a in the direction N 99°E±12° of the Eurasian plate relative to the North American plate in the period 1987–1990. Within the plate boundary deformation zone extensional strain accumulates at a rate of 0.8 strain/a. Displacement of control points next to Askja (>7 km from the caldera center) in the periods 1990–1993 and 1992–1993 show deflation and contraction towards the caldera. These results are in accordance with the results obtained by other geodetic methods in the area, which indicate that the deflation at Askja occurs in response to a pressure decrease at about 2.8 km depth, located close to the center of the main Askja caldera. A Mogi point source was fixed at this location and the GPS data used to solve for the source strength. A central subsidence of 11±2.5 cm in the period 1990–1993 is indicated, and 5.5±1.5 cm in the period 1992–1993. The maximum tensional strain rate, according to the point source model, occurs at a horizontal distance of 2.5–6 km from the source, at the same location as the main caldera boundary. Discrepancies between the observed displacements and predicted displacements from the Mogi model near the Askja caldera can be attributed to the regional eastwest extension that occurs at Askja.     

Horizontal crustal movement before the great Wenchuan earthquake obtained from GPS observations in the regional network

The great Wenchuan earthquake of M8.0 on May 12, 2008, occurred in an area with dense GPS observation stations in the regional network of the Crustal Movement Observation Network of China (CMONOC). Non-continuous observations were carried out at the 1 000 GPS stations of the regional network in 1999, 2001, 2004 and 2007. The horizontal displacements at GPS stations in the regional network before the Wenchuan earthquake show that the main driving tectonic force of the earthquake was the northward pushing of the Indian plate, added at the same time by the pushing of plates on the east and south. In comparison to the displacements in other regions, the horizontal displacements near and around the seismic area is characterized by diverging eastward displacements, that is, the stations to the north of the epicenter moved in the ENE direction while those to the south of epicenter moved in ESE direction with smaller displacements at stations near the epicenter. The accuracy of the estimated strain results is briefly discussed. In order to obtain the anomalous information before the earthquake, the methods of both best fits by trend surface and statistics have been used in the study for finding the future epicentral area from the strain accumulations in the regional network observed from 1999 to 2007 before the Wenchuan earthquake. Besides the epicentral area of the western Kunlun mountain pass earthquake of M8.1 in 2001, the results of best fits by trend surfaces of the strain accumulations from 1999 to 2007 in the regional network show that the Wenchuan earthquake occurred at the eastern fringe of a large area with relatively large accumulations of the first shear strains and also at the northeastern fringe of a smaller area with significant accumulated areal compressions. The statistics of the accumulations of the strain components demonstrates that they also showed anomalous distribution patterns in this area and its neighborhood with increasing accumulations of both shear strains and areal compressions.     

Recent crustal movement studies at seiyal area, northwest of aswan Lake, Egypt

A local geodetic network of 20 points was established in 1988 around the Seiyal fault for monitoring lateral movements. Two levellinglines, crossing the fault, were prepared for monitoring vertical movements. The initial horizontal geodetic measurements were carried out in March 1989. The measurements were repeated twice a year. The initial vertical geodetic measurements (levelling) were carried out in November 1989 and repeated in November 1990.Four horizontal geodetic measurements and two levelling measurements were analysed. The results revealed remarkable horizontal and vertical changes during the period of measurement.     

Determination of recent crustal deformations based on precise GPS measurements in the Vogtland earthquake area

This paper presents methods and results of recent geodetic investigations in the Saxon part of the NW Bohemia/Vogtland earthquake area. The technique of precise GPS measurements and the intercomparison of successive observation epochs are described in detail. Horizontal surface deformations of the earth's crust were determined with an accuracy level of 1 mm over a range of some 10⁴ m. The results of five GPS campaigns carried out from 1994 to 2001 reveal significant relative displacements of up to 5 mm between successive observation epochs. Furthermore, there are clear indications that the deformation process is not linear with time but correlated with the seismic activity.     

Geodetic and Structural Research in La Palma, Canary Islands, Spain: 1992–2007 Results

We review the results of the geodetic and structural studies carried out on La Palma Island using geodetic and geophysical data during the period 1992–2007. InSAR and GPS observation techniques were applied to study the existence of deformation on the island and gravity observations were carried out for structural studies. Gravity data were inverted using a nonlinear three-dimensional gravity inversion approach to obtain the geometry of the anomalous bodies constructed in a random growth process with respect to an exponentially stratified background. The main structural feature is a large central body (under the Caldera de Taburiente) with high density, which was interpreted as the Pliocene-age uplifted seamount and a relatively dense intrusive plutonic complex/magma body. The Cumbre Vieja series is characterized by elongated minima distributed according to the rift structure. InSAR results show a clear subsidence located on the Teneguía volcano, where the last eruption took place in 1971. A thermal source is the most probable origin for this deformation. A GPS network composed of 26 stations covering the total island surface was set up. Vertical displacements determined comparing the GPS coordinates obtained in 2007 with coordinates determined in 1994 are consistent with the InSAR results obtained in the southern part of the island. This is not the case for the northern part. From the comparison of 2006 and 2007 coordinates it is clear that more time is needed to obtain significant displacements, but observed trends are also consistent with InSAR results. All the observed significant displacements are in stations located outside of the large high-density central body.     

Present state of geodynamic properties of Kalabsha area, nortwest of Aswan Lake, Egypt

Aswan Lake is the second largest man-made lake in the world. Its filling started 1964 and reached the maximum water level in 1978. An earthquake of magnitude 5.5 took place in 1981 along the most active fault near the lake (Kalabsha fault). This earthquake was follwed by a tremendous number of smaller events that continue till now.Seismicity and the underground water table around the lake are monitored continuously through a radio-telemetered network. A local geodetic network was established around parts of the active faults in the northwestern part of the High Dam Lake, for monitoring vertical and lateral movements. The Kalabsha local geodetic network (the first one) was established around an active part of the Kalabsha fault in 1983. Precise geodetic measurements have benn carried out twice a year since 1984.On the basis of the repeated geodetic measurements, seismicity of the area and geophysical as well as geological data, the present state of the geodynamical properties of the Kalabsha area is studied.Remarkable horizontal movements were detected; they are correlated with the seismicity of the area and are attributed to the differential loading by the lake. The Kalabsha fault is a right-lateral strike-slip motion on an E-W plane. The magnitude of the movements detected along the fault is variable for the different epochs of measurements and is correlated with both seismicity and water loading in the lake.     

用多种数据构建2008年汶川特大地震同震位移场

本文主要以GPS、精密水准观测和卫星SAR遥感图像分析2008年汶川特大地震同震位移特征.GPS数据包括:(1)四川盆地和川西高原地区各类国家等级GPS网点复测;(2)沿破裂带国家天文大地网GPS复测.前者推算的同震位移测定精度优于2 cm,后者6~8 cm.SAR遥感资料包括:(1)ALOS 卫星升轨相位干涉图像,精度优于8 cm;(2)ALOS和ENVISAT卫星影像合成的三维位移图,精度优于0.5 m.同震位移场显示,断层下盘(四川盆地)变形总体呈扇形集中指向震中,断层上盘(龙门山)变形总体上呈逆时针旋转态势,最大的实测水平位移5.5 m.汶川、理县、茂县等地测站位移指向破裂带方向,而平武、青川等地测站逐渐转变为平行,乃至远离破裂带方向,与汶川地震逆冲兼走滑的破裂特征一致.断层上盘大幅隆升,下盘靠近断层的区域以下沉为主,远场表现为幅度很小的隆升,垂直升降区域间,有一条与龙泉山断裂带平行的升降过渡带,调节龙泉断层的应力状态.用实测变形场检验多个地震波破裂模型表明,近场(距离断层50 km) 模型形变准确度可达40~50 cm, 远场精度优于5 cm.     

Verification of the concept of seismoionospheric coupling under quiet heliogeomagnetic conditions,using the Wenchuan (China) earthquake of May 12, 2008, as an example

The variability of the ionosphere during April–May, 2008, has been analyzed in detail in order to reveal anomalous variations related to seismic activity, initiated by the strongest Wenchuan earthquake (M = 7.9) in the Sichuan province on May 12, 2008. Information about the total electron content (TEC) from the network of GPS receivers in the earthquake region, the global IONEX TEC maps, and the reconstructed vertical profiles of electron density according to the data of GPS receivers were used as a data source. The spatial and time localization of the observed anomalies, their morphological features, and the absence of geomagnetic disturbances during the observation period undoubtedly demonstrate that the observed variations were caused by seismic activity.     

A note on the L’Aquila earthquake of 6 April 2009: Permanent ground displacements obtained from strong-motion accelerograms

On 6 April 2009, 01:32 GMT, an M_w 6.3 earthquake hit the Abruzzi region of central Italy causing widespread damage in the City of L’Aquila and its nearby villages. The mainshock of this earthquake was recorded by 57 digital strong-motion instruments, four of which are located on the hanging wall of the Paganica Fault near L’Aquila. These stations are no more than 6 km from the epicentre. We use accelerometric data from these four stations to estimate permanent ground displacements caused by the mainshock. Our numerical results reveal south-east and downwards directed permanent co-seismic displacements which are in fair agreement with the outcomes of GPS and InSAR measurements reported in preliminary Istituto Nazionale di Geofisica e Vulcanologia (INGV) reports.     

Vertical Displacements in the Nový Kostel Seismoactive Area

Since 1994 annual campaigns of precise levelling have been performed on the network established in the surroundings of Nový Kostel. The network covers the most dynamic part of the West Bohemia seismoactive region with a total of 70 points in 1999. The measurements are connected to a reference point of the national levelling system and to our GPS-gravity network. It was proved by the error and confidence level of the data that small displacements of three and more millimetres can be recognised. The data analysis showed that all points could be divided into several groups according to temporal changes of height. Mutual comparison of these changes enabled to detect both general and very local short-term movements. It appeared that, in general, the involved part of the mountain block is subsiding relative to the Cheb basin. This is in contradiction to the recent uplift of the Kruné hory Mts. The correlation between vertical displacements and earthquake swarms was examined with the conclusion that during a swarm period the movements have special homogeneous pattern, contrary to inter-periods. The division line of different displacements (a fault zone) for the swarms 1994 and 1997 was determined.     

Horizontal crustal movement in Chinese mainland from 1999 to 2001

The paper introduces the horizontal crustal movement obtained from GPS observations in the regional networks (including the basic network and the fiducial network) of the Crustal Movement Observation Network of China (CMONOC) carried out in 1999 and 2001. This paper is characterized by the acquisition of the horizontal displacement velocities during the period from 1999 to 2001 at the observation stations in the regional networks with datum definition of a group of stable stations with small mutual displacements in east China. Based on the most detailed map of horizontal crustal movement in Chinese mainland, the division of blocks, their displacements and deformations are studied. An approach to analysis of the intensity of the horizontal crustal deformation is proposed. The general characteristics of the recent horizontal crustal movement in Chinese mainland and that before the Kunlunshan earthquake of M=8.1 on November 14, 2001 are analyzed. Foundation item: The National Development and Programming Project for Key Basic Research (95-13-03-07).     

DORIS as a potential part of a Global Geodetic Observing System

We have processed all available DORIS data from all available satellites, except Jason-1 over the past 10 years (from January 1993 to April 2003). Weekly solutions have been produced for stations positions coordinates, geocenter motion and scale factor stability. We present here accuracy presently achievable for all types of potential geodetic products. Typically weekly stations positions can be derived with a repeatability of 1.0–1.5 cm using data from 5 satellites simultaneously, showing the significant improvement in precision that has been gained recently using the additional new DORIS satellites. As an example, we show how such new results can detect displacement from large magnitude earthquakes, such as the 2003 Denali fault earthquake in Alaska. Displacements of −5 cm in latitude and +2 cm in longitude were easily detected using the DORIS data and are confirmed by recent GPS determination. The terrestrial reference frame was also well be monitored with DORIS during this 10-year period. Other geodetic products, such as tropospheric corrections for atmospheric studies are also analyzed. Finally, we discuss here the possible advantages and weaknesses of the DORIS system as additional geodetic tool, in conjunction with the already existing GPS, VLBI and SLR services, to participate in an Global Geodetic Observing System (GGOS).     

Comparison of geodetic and seismic strain rates in Greece by using a uniform processing approach to campaign GPS measurements over the interval 1994–2000

In this study we rigorously combine 18 old campaign GPS data sets from Greece covering the period 1994–2000. Although the majority of these old datasets have been analyzed and reported previously, it has not been possible to combine them into a single velocity field and apply strain analysis. Here a uniform, final coordinate solution is given by reprocessing 43 global, long-running International GNSS Service (IGS) sites together with 280 local sites. The 221 daily SINEX files are then combined in a least squares approach and the geodetic horizontal velocity field in ITRF2000 and Europe-fixed reference frame is derived. Two methods are used to compute the geodetic strain rates: (i) discrete estimates within contiguous polygons, and (ii) a continuous curvature surface fitted to the velocity field. The seismic hazard potential can be determined by comparing the geodetic and seismic strain rates. The published 300 year earthquake catalogue best describes the major active tectonic features at the scale of geodetic strain determination. The geodetic strain appears larger than the seismic strain for the majority of the region, suggesting that accumulated strain has not yet been released by earthquakes. The geodetic field is consistent with the detailed constraints implied by the observed orientations of faulting as these are given in the 300-year catalogue. We have shown that with the GPS dataset used in this work and following this processing scheme reasonable results can be obtained comparable with more recent studies, CGPS data and by recent earthquake activity.