Studies on the regional feature of organic carbon in sediments off the Huanghe River Estuary waters

Organic carbon (OC) in definitely small area sediments(according to marine dimension)off the Huanghe River Estuary is investigated in order to evaluate the feature of regional difference of physical and chemical properties in marginal sea sediments. The distributions of OC in sediments with natural grain size and the relationship with the pH, Eh,Es and Fe3+/Fe2+ are discussed. In addition,OC decomposition rates in surfacial/subsurfacial sediments are estimated. OC concentrations range from 0.26% to 1.8%(wt)in the study area. Significant differences in OC content and in horizontal distribution as well as various trends in surfacial/subsurfacial sediments exhibit the feature of regional difference remarkably in marginal sea sediments. The complicated distribution of OC in surface sediments is due to the influence of bacterial activity and abundance, bioturbation of benthos and physical disturbance. The OC decomposition rate constant in surfacial/subsurfacial sediments ranges from 0.0097 to 0.076 a(-1) and the relatively high values may be mainly related to bacteria that are mainly responsible for OC mineralization;meio-and macrofauna affect OC degradation both directly, through feeding on it, and indirectly through bioturbation and at the same time coarse sediments are also disadvantageous to OC preservation. In almost all the middle and bottom sediments the contents of OC decrease with the increase of deposition depth, which indicates that mineralization of OC in the middle and bottom sediments has occurred via processes like SO42- reduction and Fe-oxide reduction.

Seasonal variations in growth and clearance rate of the Zhikong scallop Chlamys farreri suspended in the deep water of Haizhou Bay, China

The growth and survival of the Zhikong scallop Chlamys farreri suspended in deep water of Haizhou Bay were studied from July 2007 to June 2008, and the biodeposition method was used to estimate the clearance rate of C. farreri under field conditions. Results showed that the scallop grew fast during all the culture time, with the exception of summer. The condition index of the scallop increased with time and reached the highest value in spring of the second year. The survival of scallops was 60.8 +/- A 3.9% at the end of this study, mortality occurring mainly during the summer and autumn of the first year. The clearance rate fluctuated obviously with season,with the highest value in September 2007, and the lowest value in March 2008. Factors accounting for variations in growth and clearance rate of scallops are also discussed.

Spatial variations of size-fractionated chlorophyll, cyanobacteria and heterotrophic bacteria in the Central and Western Pacific

Geographic and vertical variations of size-fractionated (0.2-1 mu m, 1-10 mu m, and >10 mu m) Chlorophyll a (Chl.a) concentration, cyanobacteria abundance and heterotrophic bacteria abundance were investigated at 13 stations from 4 degrees S, 160 degrees W to 30 degrees N, 140 degrees E in November 1993. The results indicated a geographic distribution pattern of these parameters with instances of high values occurring in the equatorial region and offshore areas, and with instance of low values occurring in the oligotrophic regions where nutrients were almost undetectable. Cyanobacteria showed the highest geographic variation (ranging from 27x10(3) to 16,582x10(3) cell l(-1)), followed by Chl.a (ranging from 0.048 to 0.178 mu g l(-1)), and heterotrophic bacteria (ranging from 2.84x10(3) to 6.50 x 10(5) cell l(-1)). Positive correlations were observed between nutrients and Chl.a abundance. Correspondences of cyanobacteria and heterotrophic bacteria abundances to nutrients were less significant than that of Chl.a. The total Chl.a was accounted for 1.0-30.9%, 35.9-53.7%, and 28.1-57.3% by the >10 mu m, 1-10 mu m and 0.2-1 mu m fractions respectively. Correlation between size-fractionated Chl.a and nutrients suggest that the larger the cell size, the more nutrient-dependent growth and production of the organism. The ratio of pheophytin to chlorophyll implys that more than half of the > 10 mu m and about one third of the 1-10 mu m pigment-containing particles in the oligotrophic region were non-living fragments, while most of the 1-10 mu m fraction was living cells. In the depth profiles, cyanobacteria were distributed mainly in the surface layer, whereas heterotrophic bacteria were abundant from surface to below the euphotic zone. Chl.a peaked at the surface layer (0-20 m) in the equatorial area and at the nitracline (75-100 m) in the oligotrophic regions. Cyanobacteria were not the principle component of the picoplankton. The carbon biomass ratio of heterotroph to phytoplankton was greater than 1 in the eutrophic area and lower than 1 in oligotrophic waters.

Monitoring and simulation of water, heat,and CO2 fluxes in terrestrial ecosystems based on the APEIS-FLUX system

The Integrated Environmental Monitoring (IEM) project, part of the Asia-Pacific Environmental Innovation Strategy (APEIS) project, developed an integrated environmental monitoring system that can be used to detect, monitor, and assess environmental disasters, degradation, and their impacts in the Asia-Pacific region. The system primarily employs data from the moderate resolution imaging spectrometer (MODIS) sensor on the Earth Observation System- (EOS-) Terra/Aqua satellite,as well as those from ground observations at five sites in different ecological systems in China. From the preliminary data analysis on both annual and daily variations of water, heat and CO2 fluxes, we can confirm that this system basically has been working well. The results show that both latent flux and CO2 flux are much greater in the crop field than those in the grassland and the saline desert, whereas the sensible heat flux shows the opposite trend. Different data products from MODIS have very different correspondence, e.g. MODIS-derived land surface temperature has a close correlation with measured ones, but LAI and NPP are quite different from ground measurements, which suggests that the algorithms used to process MODIS data need to be revised by using the local dataset. We are now using the APEIS-FLUX data to develop an integrated model, which can simulate the regional water,heat, and carbon fluxes. Finally, we are expected to use this model to develop more precise high-order MODIS products in Asia-Pacific region.

Stable Carbon Isotopic Composition in Soil Organic Carbon and C3/C4 Source Variations at the Haibei Alpine Meadow

Stable carbon isotopes of organic matter originated from different soil layers (0～5 cm, 5～15 cm, 15～25 cm, 25～35 cm, 35～50 cm, 50～65 cm) were investigated in the Haibei Alpine Meadow Ecosystem Research Station of the Chinese Academy of Sciences. The preliminary results indicated that δ13C values of soil organic matter increased with increased soil depth. δ13C of soil organic carbon in 0～5 cm layer showed the lowest value, -25.09‰; while 50～65 cm soil layer possessed the lowerδ13C value, -13.87‰. Based on mass balance model of stable isotopes, it was proposed that the percentage of C4 carbon source tend to increase with increased soil depth. The preliminary study indicated that alpine meadow might have undergone a successive process from C4-dominated community to C3-dominated one. However, changing δ13C values in atmospheric CO2 overtime and different processes of soil organic carbon formation (or eluviation) might somewhat contribute to increasing δ13C values. In this case, mass balance model would underestimate C3 community and overestimate C4 community.

Diurnal and seasonal variations of UV radiation on the northern edge of the Qinghai-Tibetan Plateau

We monitored UVA, UVB, and solar radiation from August 2001 to 2003 on the northern Qinghai-Tibetan Plateau to characterize the diurnal and seasonal variations of UV radiation on the world's highest plateau. Daily UVB radiation and the ratio of UVB to total solar radiation increased significantly when the atmospheric ozone concentration decreased as estimated by the total ozone mapping spectrometer (TOMS), as well as when cloud coverage decreased. The UVB/UVA ratio also showed a significant increase when the TOMS ozone concentration decreased in the morning. The seasonal variation pattern of UVB, however, was closely correlated with solar elevation but was little affected by the seasonal pattern of the atmospheric ozone amount. Compared to observations from the central plateau, the magnitude of the UVB increase attributed to ozone depletion was smaller at the northern edge. The study suggests that the temporal variation of ground UV radiation is determined by both solar elevation and the ozone amount, but the spatial difference on the plateau is likely to be ascribed mainly to the spatial variation of the ozone amount. (c) 2007 Published by Elsevier B.V.

黄土-红粘土化学风化与晚新生代气候变化——甘肃平凉白水剖面为例

Late Cenozoic has witnessed a series of climate-environmental change which ends with a transitional shift from greenhouse to icehouse conditions. In last two decades, scientists began to employ the tectonic uplift and its weathering effect to interpret the climatic changes during the late Cenozoic. However, this endeavor has partly been restricted by the lacking of regional and global chemical weathering data. The loess-red clay deposit in the Loess Plateau may record the weathering features of the detritus material from the wide range upwind of the Loess Plateau. Therefore geochemistry of the loess-plaeosol and red clay sequences may provide insights into the regional chemical weathering regime and the connection between the chemical weathering and the late Cenozoic climate-environmental change Here we selected 319 samples from the Baishui section near the Pingliang City, Gansu Province, and analyzed them with X-ray fluorescence. Based on the result, we reconstruct the chemical weathering history of the Baishui section since 6Ma. We chose CIA as the proxy for chemical weathering intensity. The CIA ratio in soil units is higher than in adjacent loess horizons, but lower than in the red clay, in good agreement with the field observation. The CIA ratios of the Baishui section correlates well with the global ice volume fluctuations, indicating that the global cooling may contribute a lot to the chemical weathering variations in Chinese Loess Plateau. There are at least 3 million-year time scale variations that can be identified in the chemical weathering intensity curve, i.e., between 3.3 to 2.1 Ma, 1.7 to 0.9 Ma and from 0.9Ma. We think these may reflect the combined effect of the tectonic uplift and ice sheets on monsoon intensity. Other time scale variations can be also observed. In the period between 2.4 and 0.8 Ma, the CIA record display the 400,000 years cycle, which may be resulted from the Tibetan uplift during the Pliocene-early Pleistocene which have significantly amplified the monsoon response sensitivity to the orbital-scale variations in insolation. From 1.2 Ma, the 100,000 years period became intensifying, and particularly after 0.8 Ma, the earlier monsoon response at 400,000 year periodicity was overwhelmed by the ice sheet forcing at 100,000 year periodicity. These may indicate that the expansions of the Northern Hemisphere ice sheets may have crossed a threshold, which enforce the monsoon responding at the 100,000 year periodicity.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

海底金属矿床充水条件及矿坑突水的断层效应研究

The Xinli mine area of Sanshandao mine is adjacent to the Bohai Sea and its main exploitable ore deposit occurs in the undersea rock mass. The mine is the biggest undersea gold mine of China after production. The mine area faces a latent danger of water bursting, even sudden seawater inrush. There is no mature experience in undersea mining in China so far. The vein ore deposit is located in the lower wall of a fault; its possible groundwater sources mainly include bittern, Quaternary pore water and modern seawater. To ensure the safety of undersea mining, to survey the flooding conditions of the ore deposit using proper measures and study the potential seawater inrush pattern are the key technical problems. With the Xinli mine area as a case study, the engineering geological conditions of the Xinli mine area are surveyed in situ, the regional structural pattern and rock mass framework characteristics are found out, the distribution of the structural planes are modeled by a Monte Carlo method and the connectivity coefficients of rock mass structural planes are calculated. The regional hydro-geological conditions are analyzed and the in-situ hydro-geological investigation and sampling are performed in detail, the hydrochemistry and isotopes testing and groundwater dynamic monitoring are conducted, the recharge, runoff, discharge conditions are specified and the sources of flooding are distinguished. Some indices are selected from the testing results to calculate the proportion of each source in some water discharge points and in the whole water discharge of the Xinli mine area. The temporal and spatial variations of each water source of the whole ore deposit flooding are analyzed. According to the special project conditions in the Xinli mine area, the permeability coefficient tensors of the rock mass in Xinli mine area are calculated based on a fracture geometry measurement method, in terms of the connectivity and a few hydraulic testing results, a modified synthetic permeability coefficient are calculated. The hydro-geological conceptual and mathematical model are established，the water yield of mine is predicted using Visual Modflow code. The spreading law of surrounding rock mass deformation and secondary stress are studied by numerical analysis; the intrinsic mechanism of the faults slip caused by the excavation of ore deposit is analyzed. The results show that the development of surrounding rock mass deformation and secondary stress of vein ore deposit in the lower wall of a fault, is different from that in a thick-big ore deposit. The secondary stress caused by the excavation of vein ore deposit in the lower wall of a fault, is mainly distributed in the upper wall of the fault, one surface subsidence center will occur. The influences of fault on the rock mass movement, secondary stress and hydro-geological structures are analyzed; the secondary stress is blocked by the fault and the tensile stress concentration occurs in the rock mass near the fault, the original water blocking structure is destructed and the permeable structure is reconstructed, the primary structural planes begin to expand and newborn fissures occur, so the permeability of the original permeable structure is greatly enhanced, so the water bursting will probably occur. Based on this knowledge, the possible water inrush pattern and position of the Xinli mine area are predicted. Some computer programs are developed using object-oriented design method under the development platform Visual Studio.Net. These programs include a Monte Carlo simulation procedure, a joint diagrammatizing procedure, a structural planes connectivity coefficient calculating procedure, a permeability tensor calculating procedure, a water chemical formula edit and water source fixture conditions calculating procedure. A new computer mapping algorithm of joint iso-density diagram is raised. Based on the powerful spatial data management and icon functions of Geographic Information System, the pit water discharge dynamic monitoring data management information systems are established with ArcView.

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.