Variation characteristics and ecological risk of heavy metals in the south Yellow Sea surface sediments

Eight cruises were conducted on the south Yellow Sea (SYS) from 1998 to 2005. Variations and the potential ecological risk of heavy metals were studied using the survey data collected during October 2003. The metal content (except for As) was high in the central area where the fine grain size sediments were dominant, and low inshore area where more coarse sediments were present. This suggested that grain size was important in determining distributions of heavy metals. In some local areas, other influencing factors, such as organic content, sedimentation rate, burial efficiency and metal's existing form were discussed. The annual averages of metals showed a stable trend with appreciable fluctuations in 8 years. Using potential ecological risk index (E (RI)) to evaluate the integrated pollution effect of heavy metals, 38.7% of the investigated area was in a moderate degree of contamination, while 77.8% was under moderate ecological risk. However, no distinct correlation was found between E (RI) and plankton biomass. In conclusion, the sediment quality of SYS was good, and the ecological risk was low in general.

Effect of spatial variation on areal evapotranspiration simulation in Haibei, Tibet plateau, China

Quantification of areal evapotranspiration from remote sensing data requires the determination of surface energy balance components with support of field observations. Much attention should be given to spatial resolution sensitivity to the physics of surface heterogeneity. Using the Priestley-Taylor model, we generated evapotranspiration maps at several spatial resolutions for a heterogeneous area at Haibei, and validated the evapotranspiration maps with the flux tower data. The results suggested that the mean values for all evapotranspiration maps were quite similar but their standard deviations decreased with the coarsening of spatial resolution. When the resolution transcended about 480 m, the standard deviations drastically decreased, indicating a loss of spatial structure information of the original resolution evapotranspiration map. The absolute values of relative errors of the points for evapotranspiration maps showed a fluctuant trend as spatial resolution of input parameter data layers coarsening, and the absolute value of relative errors reached minimum when pixel size of map matched up to measuring scale of eddy covariance system. Finally, based on the analyses of the semi-variogram of the original resolution evapotranspiration map and the shapes of spatial autocorrelation indices of Moran and Geary for evapotranspiration maps at different resolutions, an appropriate resolution was suggested for the areal evapotranspiration simulation in this study area.

Local warming about 1.3 degrees C in alpine meadow has no effect on root vole (Microtus oeconomus L.) population during winter

The influence of air and soil warming on root vole (Microtus oeconomus L.) population was studied in winter period in top open chambers (OTC) (0.8-1.8 m(2) warmed by conical fiberglass material and situated in alpine meadow (3250 m) at Qinghai-Tibet Plateau, China. The OTCs were distributed on an area of 30 x 30 m of experimental warming site; another site of the same area was a control one. The root vole population was investigated on two pairs of sites in "low-grazing" and "high-grazing" (by sheep) parts of the meadow; mark-recapture method was used. The winter-season averaged air and soil temperature inside of the chambers were 1.3 degrees C higher than the temperature outside the chambers. The warming in the chambers had no statistically significant effect on root vole numbers, on average body mass of individual, and on average body mass of males and females. In conclusion, as small as 1.3 degrees C warming of soil and air introduced locally and on small (several m(2)) scale, in the alpine meadow habitat in winter period, has possibly no effect on root vole numbers and biomass.

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

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.