Temperature observations of Crater Lake,Mt Ruapehu,New Zealand,using Temperature Telemetry Bouys

The temperature of the Crater Lake of the active volcano Ruapelm has been recorded by Temperature Telemetry Buoys, to determine if lake temperature is correlated with volcanic activity. These buoys had to be specially designed to cope with the unfavourable environment of Crater Lake. A buoy contains a thermistor to measure the lake temperature, and a radio transmitter to transmit a short signal every few minutes, the interval between signals being a function of temperature. The temperature records obtained from these buoys show that the temperature near the lake surface can vary considerably within a few hours. Some of these variations appear to be caused by disturbances in convective heat transfer occurring in the lake. The occurrence of these short term temperature variations means that there is no simple relation between Crater Lake temperatures and the volcanic activity of Ruapehu. Some rapid increases in temperature followed volcanic earthquakes, but one of the biggest increases in temperature occurred just before a group of earthquakes upder the lake.     

Fault patterns by space remote sensing and the rotation of western Oregon during Cenozoic times

We used MSS-Landsat images (bands 6 and 7) as well as previous studies to interpret fault patterns in western Oregon in terms of the rotations detected by paleomagnetism. These fault zones show a great concentration between the Cascade Range and the Idaho Batholite and have several distinct trends. The north-south striking structures are attributed to the Basin and Range province, while the large NW-SE right-lateral fault zones are interpreted as resulting from an extension between the Cascades Arc and the Olympic-Wallowa lineament. The latter was a paleoplate boundary during pre-Eocene times. This extension, beginning during Late Eocene/Oligocene times and continuing through Recent times, is accompanied by a migration of the rotation pole from southeast to northwest, and by a clockwise rotation of the Coast Range—Klamath Mountains—Cascade Range block, induced by the subduction of the Farallon plate.     

The complex filling of alae crater,Kilauea Volcano,Hawaii

Since February 1969 Alae Crater, a 165-m-deep pit crater on the east rift of Kilauea Volcano, has been completely filled with about 18 million m³ of lava. The filling was episodic and complex. It involved 13 major periods of addition of lava to the crater, including spectacular lava falls as high as 100 m, and three major periods of draining of lava from the crater. Alae was nearly filled by August 3, 1969, largely drained during a violent ground-cracking event on August 4, 1969, and then filled to the low point on its rim on October 10, 1969. From August 1970 to May 1971, the crater acted as a reservoir for lava that entered through subsurface tubes leading from the vent fissure 150 m away. Another tube system drained the crater and carried lava as far as the sea, 11 km to the south. Much of the lava entered Alae by invading the lava lake beneath its crust and buoying the crust upward. This process, together with the overall complexity of the filling, results in a highly complicated lava lake that would doubtless be misinterpreted if found in the fossil record.     

Dating of paleomagnetic logs by tephra layers,Mara Lake,Canada

Summary The paleodeclination and paleoinclination logs compiled from the cores taken from Mara Lake show consistent, well defined oscillations. It was hoped that the Mt. St. Helens and Mt. Mazama tephra layers encountered in the cores would provide accurate, absolute dating control of the cores as well for the paleodeclination and paleoinclination logs. Unfortunately the dates of 3300 BP for the Mt. St. Helens tephra and 6845 ± 50 BP for the Mt. Mazama tephra are incompatible with a constant rate of sedimentation. It would appear that the Mt. Mazama date is valid and a redetermination of the Mt. St. Helens date based on the 6845 ± 50 BP date for the Mt. Mazama tephra and a constant rate of sedimentation from Mazama time to the present gives a date of 4100 BP for the Mt. St. Helens tephra. This dating control has been used to construct paleodeclination and paleoinclination logs versus time which correlate well with paleomagnetic logs from Fish Lake, Oregon and Shawnigan Lake on Vancouver Island.Presented at 2nd conference on New Trends in Geomagnetism, Castle Bechyn, Czechoslovakia, September 24–29, 1990.     

A comparison of the recent eruptions of Mt. Pelée,Martinique and Soufrière,St. Vincent

The simultaneous eruption of Mt. Pelée, Martinique and Soufrière, St. Vincent are regarded as the first recognized examples of Pelean-type and St. Vincent-type pyroclastic eruptions. Both produced nuées ardentes, the former usually laterally directed because of the presence of a dome and the latter vertically directed from an open crater. Both volcanoes have subsequently erupted for a second time this century. The 1902–05 and 1929–32 eruptions of Mt. Pelée produced andesite lava of almost identical composition and mineralogy. Both contain two generations of plagioclase, orthopyroxene, Fe-Ti oxide, corroded brown amphibole and olivine rimmed by pyroxene. In contrast, the Soufrière material is more basic in composition varying from basaltic andesite to basalt in 1902–03 and basaltic andesite in 1971–72. The Soufrière material contains two generations of plagioclase (with those of 1971–72 having additional zones of labradorite), clinopyroxene, orthopyroxene, olivine and Fe-Ti oxide. The pyroclastic deposits are strikingly different, those from the Pelean-type eruption are termed «block and ash deposits» being characterised by poorly vesicular lava blocks up to 7 m in diameter, while the St. Vincent-type eruption produced «scoria and ash deposits» containing vesicular ropey blocks or bombs no larger than 1 m in diameter. The differences in styles of eruption are attributed to differences in viscosity and mechanism of eruption of the magmas. Stratigraphic studies of Mt. Pelée reveal that the volcano has produced basaltic andesite scoria and ash deposits from St. Vincent-type eruptions. It is concluded that the recent eruptions of Pelée tapped a deep level magma during both eruptions releasing magma of similar composition, while the 1971 Soufrière magma is thought to be a remnant of the 1903 basaltic magma which remained at a high level within the volcano where it underwent enrichment in plagioclase and loss of olivine and oxide.     

Geoelectrical and geological structure of the crust in Western Slovakia

Electrical resistivity of the Earth’s crust is sensitive to a wide range of petrological and physical parameters, and it particularly clearly indicates crustal zones that have been tectonically or thermodynamically disturbed. A complex geological structure of the Alpine nappe system, remnants of older Hercynian units and Neogene block tectonics in Western Slovakia has been a target of recent magnetotelluric investigations which made a new and more precise identification of the crustal structural elements of the Western Carpathians possible. A NW-SE magnetotelluric profile, 150 km long, with 30 broad-band and 3 long-period magnetotelluric sites, was deployed, crossing the major regional tectonic elements listed from the north: Brunia (as a part of the European platform), Outer Carpathian Flysch, Klippen Belt, blocks of Penninic or Oravicum crust, Tatricum and Veporicum. Magnetotelluric models were combined with previous seismic and gravimetric results and jointly interpreted in the final integrated geological model. The magnetotelluric models of geoelectrical structures exhibit strong correlation with the geological structures of the crust in this part of the Western Carpathians. The significant resemblance in geoelectrical and crustal geological structures are highlighted in shallow resistive structures of the covering formations represented by mainly Tertiary sediments and volcanics. Also in the deeper parts of the crust highly resistive and conductive structures are shown, which reflect the original building Hercynian crust, with superposition of granitoids or granitised complexes and lower metamorphosed complexes. Another important typical feature in the construction of the Western Carpathians is the existence of young Neogene steep fault zones exhibited by conductive zones within the whole crust. The most significant fault zones separate individual blocks of the Western Carpathians and the Western Carpathians itself from the European Platform.     

Seasonal seismicity at western United States volcanic centers

We examine 20-yr data sets of seismic activity from 10 volcanic areas in the western United States for annual periodic signals (seasonality), focusing on large calderas (Long Valley caldera and Yellowstone) and stratovolcanoes (Cascade Range). We apply several statistical methods to test for seasonality in the seismic catalogs. In 4 of the 10 regions, statistically significant seasonal modulation of seismicity (> 90% probability) occurs, such that there is an increase in the monthly seismicity during a given portion of the year. In five regions, seasonal seismicity is significant in the upper 3 km of the crust. Peak seismicity occurs in the summer and autumn in Mt. St. Helens, Hebgen Lake/Madison Valley, Yellowstone Lake, and Mammoth Mountain. In the eastern south moat of Long Valley caldera (LVC) peak seismicity occurs in the winter and spring. We quantify the possible external forcing mechanisms that could modulate seasonal seismicity. Both snow unloading and groundwater recharge can generate large stress changes of > 5 kPa at seismogenic depths and may thus contribute to seasonality.     

Tectonic interpretation of the magnetic anomalies southwest of Vancouver Island

Summary A tectonic interpretation of the magnetic anomalies off the coast of California, Oregon and Washington between 40° and 52° north latitude shows that in the area surveyed the oceanic crust is cut by seven major dislocation zones which divide the region in eight areas. For five of these areas the original connection can be reconstructed. The existence of a window of young crust surrounded by older crust and of a short, isolated length of active oceanic ridge southwest of Vancouver Island as proposed recently by different authors is not confirmed.     

Generalization of the results of long-term strainmeter and GPS observations for intraplate regions

The measurements of contemporary deformations in Asia in the region of active seismic processes stretching from Northern Tien Shan through Altai and Sayan Mountains to Lake Baikal have been carried out for a few decades by strainmeters and tiltmeters. Recently, these works were expanded by the network measurements by the methods of space geodesy. Studying the contemporary displacements and deformations of the Earth’s crust is necessary for solving various problems on a wide range of spatial and time scales from exploring the contemporary deformations of the tectonic plates and present-day seismicity to estimating the man-made loads on the geological medium. The measurements of the deformations are conducted by different methods with short and long baselines, on the surface and at different depths. Based on the obtained data, it is possible to analyze the strain rates at different frequencies and to compare them with the level of seismic activity in the region. In recent decades, the seismic zone located in the south of the Siberian platforms experienced a few strong (M = 6.3–7.5) earthquakes, during which the deformation and displacement fields were captured. South of this region, the displacement fields during the earthquakes with the magnitude above 8 were recorded by the space geodetic methods in China. By excluding the coseismic effects, one can estimate the tectonic component of the displacements and deformations.     

Imbrication, flow direction and possible source areas of the pumice-flow tuffs near Bend, Oregon, U.S.A.

Imbrication, indicating flow and source direction, occurs in three Pleistocene or upper Pliocene pumice-flow tuffs exposed in a 700-km² area on the east flank of the Cascade Range near Bend, Oregon, and shows the location of previously unknown source vents of these tuffs. The imbrication is formed by inclined elongate and/or flat pumice or lithic fragments and locally by elongate plagioclase crystals. Imbrication is best developed within the lower zones of individual flow units; the pumiceous top zones also locally show imbrication directions parallel to that in the lower zones. Moreover, the areal pattern of size distribution of lithic and pumice fragments in the flows is concordant with the flow direction pattern indicated by imbrication.The upper pumice flow shows a fan-shaped pattern of flow directions indicated by imbrication which points to a western source. A possible vent, about 20 km west of Bend in the highland near Broken Top Volcano, is marked by many silicic domes and basaltic cinder cones where there is a 6–8 mgal negative Bouguer gravity anomaly. In contrast, imbrication in the middle and lower pumice flows indicates flow from a source southwest of Bend. Vents in this direction are not obvious. Possible buried vents are located about 30 km and 45 km southwest of Bend near Sitkum Butte and Lookout Mountain, respectively.     

Structural control of volcanism in the McMurdo Volcanic Group,Antarctica

Volcanoes of the McMurdo Volcanic Group occur in four volcanic provinces: Balleny, Hallett, Melbourne and Erebus. The Balleny and Hallet provinces are distributed along the Balleny Fracture Zone and Hallett Fracture respectively. Stratovolcanoes within the Melbourne province may be associated with north to northwest-trending grabens and faults in northern Victoria Land. The Erebus volcanic province is located at the intersection of the Rennick Fault and northeast trending faults along the central Transantarctic Mountains. Within the Erebus province, volcanic centres around Mt. Erebus and Mt. Discovery possess radial symmetry which may be related to radial fractures at approximately 120° to each other.     

Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the Tibetan Plateau

To investigate the relationship between velocity structure and earthquake activity on the southeastern front of the Tibetan Plateau, we make use of continuous observations of seismic ambient noise data obtained at 55 broadband stations from the regional Yunnan Seismic Network. These data are used to compute Rayleigh wave Green's Functions by cross-correlating between two stations, extracting phase velocity dispersion curves, and finally inverting to image Rayleigh wave phase velocity with periods between 5 and 34 s by ambient noise tomography. The results show significant lateral variations in crustal and uppermost mantle structures in the studied region. Phase velocity anomalies at short periods(5–12 s) are closely related to regional tectonic features such as sediment thickness and the depth of the crystalline basement. The Sichuan-Yunnan rhombic block, enclosed by the Honghe, Xiaojiang and Jianchuan faults, emerges as a large range of low-velocity anomalies at periods of 16–26 s, that inverts to high-velocity anomalies at periods of 30–34 s. The phase velocity variation in the vicinity of the Sichuan-Yunnan rhombic block suggests that the low-velocity anomaly area in the middle-lower crust may correspond to lower crustal channelized flow of the Tibetan Plateau. The spatial distribution of strong earthquakes since 1970 reveals that the Yunnan region is inhomogeneous and shows prominent characteristics of block motion. However, earthquakes mostly occur in the upper crust, with the exception of the middle-Yunnan block where earthquakes occur at the interface zone between high and low velocity as well as in the low-velocity zones, with magnitudes being generally less than 7. There are few earthquakes of magnitude 5 at the depths of 15–30 km, where gather earthquakes of magnitude 7 or higher ones which mainly occur in the interface zone between high and low velocities with others extending to the high-velocity abnormal zone.     

Electrical resistivity tomography study of Taal volcano hydrothermal system, Philippines

Taal volcano (311?m in altitude) is located in The Philippines (14°N, 121°E) and since 1572 has erupted 33 times, causing more than 2,000 casualties during the most violent eruptions. In March 2010, the shallow structures in areas where present-day surface activity takes place were investigated by DC resistivity surveys. Electrical resistivity tomography (ERT) lines were performed above the two identified hydrothermal areas located on the northern flank of the volcano and in the Main Crater, respectively. Due to rough topography, deep valleys, and dense vegetation, most measurements were collected using a remote method based on a laboratory-made equipment. This allowed retrieval of information down to a depth of 250?m. ERTs results detail the outlines of the two geothermal fields defined by previous self-potential, CO₂ soil degassing, ground temperature, and magnetic mapping (Harada et al. Japan Acad Sci 81:261–266, 2005; Zlotnicki et al. Bull Volcanol 71:29–49, 2009a, Phys Chem Earth 34:294–408, 2009b). Hydrothermal fluids originate mainly from inside the northern part of the Main Crater at a depth greater than the bottom of the Crater Lake, and flow upward to the ground surface. Furthermore, water from the Main Crater Lake infiltrates inside the surrounding geological formations. The hydrothermal fluids, outlined by gas releases and high temperatures, cross the crater rim and interact with the northern geothermal field located outside the Main Crater.     

GPS-based crustal deformations in Azerbaijan and their influence on seismicity and mud volcanism

Using Shen’s method (Shen et al., 1996), deformations of the Earth’s crust in Azerbaijan were studied based on GPS measurements. For estimating the rate of deformation, we used the field of velocity vectors for Azerbaijan, Iran, Georgia, and Armenia that were derived from GPS measurements during 1998–2012. It is established that compression is observable along the Greater Caucasus, in Gobustan, the Kura depression, Nakhchyvan Autonomous Republic, and adjacent areas of Iran. The axes of compression/contraction of the crust in the Greater Caucasus region are oriented in the S-NE direction. The maximum strain rate (approximately 200 × 10^?9 per annum) is documented in the zone of mud volcanism at the SHIK site (Shykhlar), which is marked by a sharp change in the direction of the compression axes (SW-NE). It is revealed that the deformation field also includes the zones where strain rates are very low approximating 5 × 10^?9 per annum. These zones include the Caspian-Guba and northern Gobustan areas, characterized by extensive development of mud volcanism. The extension zones are confined to the Lesser Caucasus and are revealed in the Gedabek (GEDA) and Shusha (SHOU) areas, as well as in the zone located between the DAMO and PIRM sites (Iran), where the deformation rate amounts to 100 × 10^?9 per annum. It is concluded that the predominant factor responsible for the eruption of mud volcanoes is the intensity of gas-generation processes in the earth’s interior, while deformation processes play the role of a trigger. The zone of the epicenters of strong earthquakes is correlated to the gradient zone in the crustal strain rates.     

Imaging Faults and Shear Zones Using Receiver Functions

The geometry of faults at seismogenic depths and their continuation into the ductile zone is of interest for a number of applications ranging from earthquake hazard to modes of lithospheric deformation. Teleseismic passive source imaging of faults and shear zones can be useful particularly where faults are not outlined by local seismicity. Passive seismic signatures of faults may arise from abrupt changes in lithology or foliation orientation in the upper crust, and from mylonitic shear zones at greater depths. Faults and shear zones with less than near-vertical dip lend themselves to detection with teleseismic mode-converted waves (receiver functions) provided that they have either a contrast in isotropic shear velocity (V _s), or a contrast in orientation or strength of anisotropic compressional velocity (V _p). We introduce a detection method for faults and shear zones based on receiver functions. We use synthetic seismograms to demonstrate common features of dipping isotropic interfaces and contrasts in dipping foliation that allows determination of their strike and depth without making further assumptions about the model. We proceed with two applications. We first image a Laramide thrust fault in the western U.S. (the Wind River thrust fault) as a steeply dipping isotropic velocity contrast in the middle crust near the surface trace of the fault; further downdip and across the range, where basin geometry suggests the fault may sole into a subhorizontal shear zone, we identify a candidate shear zone signal from midcrustal depths. The second application is the use of microstructural data from exhumed ductile shear zones in Scotland and in the western Canadian Shield to predict the character of seismic signatures of present-day deep crustal shear zones. Realistic anisotropy in observed shear fabrics generates a signal in receiver functions that is comparable in amplitude to first-order features like the Moho. Observables that can be robustly constrained without significant tradeoffs are foliation strike and the depth of the foliation contrast. We find that an anisotropy of only a few percent in the shear zone is sufficient to generate a strong signal, but that the shear zone width is required to be >2 km for typical frequencies used in receiver function analysis to avoid destructive interference due to the signals from the boundaries of the shear zone.     

Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

The Beni-Ilmane (North-Central Algeria) Earthquake Sequence of May 2010

Starting on 14 May 2010 and lasting several months, the village of Beni-Ilmane (Msila District, North-Central Algeria) and its surroundings were struck by an important seismic crisis marked by three successive moderate shocks (5.0 ≤ M _d ≤ 5.2). This sequence of events caused severe damage in the Beni-Ilmane village and in the epicentral area. The poor quality of masonry construction and the cumulative effects of the large number of aftershock events played a key role in the destruction. To follow this earthquake sequence, 11 temporary seismic stations, in addition to the permanent stations of the Algerian seismic network, were deployed in the region. A representative set of well located aftershocks in the period of maximum activity (lasting 18 days) were selected. The horizontal distribution of the aftershocks shows two main earthquake clusters located near Beni-Ilmane village, one cluster oriented E–W and the other oriented NNE–SSW, crossing the first cluster at its eastern tip. The aftershocks distribution suggests that the three main shocks ruptured two distinct and adjacent fault segments of about 8 km length. The focal mechanisms of the first and third events, located in the NNE–SSW cluster, show near-vertical left-lateral strike-slip fault planes. In the second cluster, oriented E–W, focal mechanisms show a high-angle reverse fault. A field survey, initiated immediately after the first main shock, identified surface fissures generated by the three largest events in the sequence. The fissures, concentrated in a narrow area at the western termination of the NE–SW Jebel Choukchot anticline (location of Beni-Ilmane village), showed several orientations which were mainly related to gravity instabilities. The 2010 Beni-Ilmane earthquake sequence, located in the Bibans–Hodna Mountains transition zone, demonstrates that the Tellian Atlas–High Plateaus border region is an active seismic zone marked by moderate and possibly strong earthquakes; thus, a reevaluation of the seismic hazard in the region is needed.     

Crystallization of tholeiitic basalt in Alae Lava Lake,Hawaii

The eruption of Kilauea Volcano August 21–23, 1963, left 600,000 cubic meters of basaltic lava in a lava lake as much as 15 meters deep in Alae pit crater. Field studies of the lake began August 27 and include repeated core drilling, measurements of temperature in the crust and melt, and precise level surveys of the lake surface. The last interstitial melt in the lake solidified late in September 1964; by mid August 1965 the maximum temperature was 690°C at a depth of 11.5 meters. Pumice air-quenched from about 1140°C contains only 5 percent crystals — clinopyroxene, cuhedral olivine (Fo ₈₀), and a trace of plagioclase, (An ₇₀). Drill cores taken from the zone of crystallization in the lake show that olivine continued crystallizing to about 1070°C; below that it reacts with the melt, becoming corroded and mantled by pyroxene and plagioclase. Below 1070°C, pyroxene and plagioclase crystallized at a constant ratio. Ilmenite first appeared at about 1070°C and was joined by magnetite at about 1050°C; both increased rapidly in abundance to 1000°C. Apatite first appeared as minute needles in interstitial glass at 1000°C. Both the abundance and index of refraction of glass quenched from melt decreased nearly linearly with falling temperature. At 1070°C the quenched lava contains about 65 percent dark-brown glass with an index of 1.61; at 980°C it contains about 8 percent colorless glass with an index of 1.49. Below 980°C, the percentage of glass remained constant. Progressive crystallization forced exsolution of gases from the melt fraction; these formed vesicles and angular pores, causing expansion of the crystallizing lava and lifting the surface of the central part of the lake an average of 19.5 cm. The solidified basalt underwent pneumatolitic alteration, including deposition of cristobalite at 800°C, reddish alteration of olivine at 700°C, tarnishing of ilmenite at 550°C, deposition of anhydrite at 250°C, and deposition of native sulfur at 100°C. Ferric-ferrous ratios suggest that oxidation with maximum intensity between 550°C and 610°C moved downward in the crust as it cooled; this was followed by reduction at a temperature of about 100°C. The crystallized basalt is a homogeneous fine-grained rock containing on the average 48.3 percent by volume intergranular pyroxene (augite > pigeonite), 34.2 percent plagioclase laths (An_{60 70}), 7.9 percent interstitial glass, 6.9 percent opaques (ilmenite > magnetite), 2.7 percent olivine (Fo_{70 80}), and a trace of apatite. Chemical analyses of 18 samples, ranging from initially quenched pumice to lava cored more than a year after the eruption from the center and from near the base of the lake, show little variation from silica-saturated tholeiitic basalt containing 50.4 percent SiO₂, 2.4 percent Na₂O, and 0.54 percent K₂O. Apparently there was no significant crystal settling and no appreciable vapor-phase transport of these components during the year of crystallization. However, seven samples of interstitial liquid that had been filter-pressed into gash fractures and drill holes from partly crystalline mush near the base of the crust show large differences from the bulk composition of the solidified crust—lower MgO, CaO, and Al₂O₃; and higher total iron, TiO₂, Na₂O, K₂O, P₂O₅, and F, and, in most samples, SiO₂. The minor elements Ba, Ga, Li, Y, and Yb and possibly Cu tend to be enriched in the filter-pressed liquids, and Cr and possibly Ni tend to be depleted.     

Isostatic status in North China (1)

Following Airy and Pratt principles, five kinds of local-compensation models are analysed and a rapid 3-D gravity formula is utilized to calculate isostatic anomalies for 66 models with different parameters. It is noted that isostatic gravity maps appear nearly identical in their main patterns and features. The optimum compensation model in North China is one of modified Airy models in which the different density distribution in the surface, upper crust and lower crust is taken into account and the standard crustal thickness is about 50km. The position of the complete compensation interface is located in the upper mantle. The North China platform as a whole is under sub-isostatic equilibrium status with an isostatic anomaly of about 18·10^?5 m/s² on an average. The distribution of isostatic gravity anomaly shows an obvious blockwise pattern. Most positive anomaly areas occur over the eastern part, the Jiao-Liao Block, Mt. Yan block and northern margin of the Hebei-Shandong block, whereas a negative area occurs in the Shanxi graben. The comparison of models indicates that the Moho discontinuity is not suitable to be taken as a compensation interface, and the compensation effects in Airy model are better than that in Pratt model, which is consistent with the feature of dominant layered structure and less lateral inhomogeneity in crust. Some results about composite compensation, the basic characteristics of isostatic anomaly and deep stucture will be published later in the second part of this paper.