Probability distribution of random wave-current forces

Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave-wave and wave-current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan-Pierson-Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. (c) 2006 Elsevier Ltd. All rights reserved.

高阶广义屏波动方程叠前深度偏移

Attaining sufficient accuracy and efficiency of generalized screen propagator and improving the quality of input gathers are often problems　of wave equation presack depth migration,　in this paper,a high order formula of generalized screen propagator　for one-way wave equation is proposed by using the asymptotic expansion of single-square-root operator. Based on the formula，a new generalized screen propagator is developed ,which is composed of split-step Fourier propagator and high order correction terms，the new generalized screen propagator not only improving calculation precision without sharply increasing the quantity of computation，facilitates the suitability of generalized screen propagator to the media with strong lateral velocity variation.　As wave-equation prestack depth migration is sensitive to the quality of input gathers, which greatly affect the output，and the available seismic data processing system has inability to obtain traveltimes corresponding to the multiple arrivals, to estimate of great residual statics， to merge seismic datum from different projects and to design inverse Q filter, we establish　difference equations with an embodiment of Huygens’s principle for obtaining traveltimes corresponding to the multiple arrivals,bring forward a time variable matching filter for seismic datum merging by using the fast algorithm called Mallat tree for wavelet transformations， put forward a method for estimation of residual statics by applying the optimum model parameters estimated by iterative inversion with three organized algorithm,i.e,the CMP intertrace cross-correlation algorithm,the Laplacian image edge extraction algorithm,and the DFP algorithm, and present phase-shift inverse Q filter based on Futterman’s amplitude and phase-velocity dispersion formula and wave field extrapolation theory. All of their numerical and real data calculating results shows that our theory and method are practical and efficient. Key words: prestack depth migration, generalized screen propagator, residual statics，inverse Q filter ,traveltime，3D seismic datum mergence

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.

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

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.

Reduction Approaches for Vibration Control of Repetitive Structures

The reduction approaches are presented for vibration control of symmetric, cyclic periodic and linking structures. The condensation of generalized coordinates, the locations of sensors and actuators, and the relation between system inputs and control forces are assumed to be set in a symmetric way so that the control system posses the same repetition as the structure considered. By employing proper transformations of condensed generalized coordinates and the system inputs, the vibration control of an entire system can be implemented by carrying out the control of a number of sub-structures, and thus the dimension of the control problem can be significantly reduced.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

Statistics of Turbulence in Benthal Boundary Layers

This paper deals with turbulence behavior inbenthalboundarylayers by means of large eddy simulation (LES). The flow is modeled by moving an infinite plate in an otherwise quiescent water with an oscillatory and a steady velocity components. The oscillatory one aims to simulate wave effect on the flow. A number of large-scale turbulence databases have been established, based on which we have obtained turbulencestatisticsof the boundarylayers, such as Reynolds stress, turbulence intensity, skewness and flatness ofturbulence, and temporal and spatial scales of turbulent bursts, etc. Particular attention is paid to the dependences of those statistics on two nondimensional parameters, namely the Reynolds number and the current-wave velocity ratio defined as the steady current velocity over the oscillatory velocity amplitude. It is found that the Reynolds stress and turbulence intensity profile differently from phase to phase, and exhibit two types of distributions in an oscillatory cycle. One is monotonic occurring during the time when current and wave-induced components are in the same direction, and the other inflectional occurring during the time when current and wave-induced components are in opposite directions. Current component makes an asymmetrical time series of Reynolds stress, as well as turbulence intensity, although the mean velocity series is symmetrical as a sine/cosine function. The skewness and flatness variations suggest that the turbulence distribution is not a normal function but approaches to a normal one with the increasing of Reynolds number and the current-wave velocity ratio as well. As for turbulent bursting, the dimensionless period and the mean area of all bursts per unit bed area tend to increase with Reynolds number and current-wave velocity ratio, rather than being constant as in steady channel flows.

DIFFERENTIAL GEOMETRICAL METHODS IN THE STUDY OF OPTICAL-TRANSMISSION (SCALAR THEORY) .1. STATIC TRANSMISSION CASE

Static optical transmission is restudied by postulation of the optical path as the proper element in a three-dimensional Riemannian manifold (no torsion); this postulation can be applied to describe the light-medium interactive system. On the basis of the postulation, the behaviors of light transmitting through the medium with refractive index n are investigated, the investigation covering the realms of both geometrical optics and wave optics. The wave equation of light in static transmission is studied modally, the postulation being employed to derive the exact form of the optical field equation in a medium (in which the light is viewed as a single-component field). Correspondingly, the relationships concerning the conservation of optical fluid and the dynamic properties are given, and some simple applications of the theories mentioned are presented.

基于广义模型和分类回归树的中国常见树种地理分布模拟研究

研究植被、物种分布与环境的关系一直是生态学中的重点。长期以来，在全球变化与陆地生态系统的研究中，主要研究重点是对大尺度植被分布的模拟和预测，并建立了大量的气候－植被分布关系模型。而对于物种潜在分布的模拟和预测，国内外相关的研究较少。近年来，随着统计技术和地理信息系统的发展，用于预测物种分布的统计模型技术得到了迅速的发展。统计模型技术已被广泛应用于生物地理分布、植物群落、生物多样性、气候变化影响评估等方面。 本论文基于当前在物种分布研究中应用广泛的广义线性模型、广义加法模型及分类回归树3种统计模型技术，对我国常见树种的地理分布进行模拟分析，并比较不同模型模拟精度的优劣，将模拟精度较高的模型应用于预测未来气候情景下我国几种主要树种的未来潜在地理分布。 基于建立的广义线性模型（GLM）、二次项逐步回归广义线性模型（SGLM）、广义加法模型（GAM）和分类回归树（CART）4个模型对我国20种常见树种地理分布进行模拟，结果表明，4个模型均有较高的模拟精度。GAM的模拟精度最高;添加二次项并进行逐步回归有效的提高了GLM的模拟精度;CART是一种基于规则的模型技术，模拟结果比GLM稍好，比GAM略差。 对不同树种的模拟分析表明，4个模型对于主要分布在暖温带落叶阔叶林区域的油松、辽东栎分布的模拟结果较差;GLM对分布在温带针阔混交林中红松、蒙古栎、胡桃楸和糠椴的模拟结果不太理想;4个模型对分布在中国亚热带常绿阔叶林区域的树种均表现出较高的模拟精度;对广布种也表现出很高的模拟精度。 结合地理信息系统，以地图形式将青冈、油松的模拟结果表示出来。结果表明：地理信息系统直观的反映出了模型模拟结果差异。4个模型均能很好模拟青冈的分布，且模拟结果接近;而对油松分布模拟结果4个模型均不甚理想，以GLM最差。这些结果与模型模拟评估结果相吻合。 在未来气候变化情景下，基于4个模型模拟结果优劣，以我国三种主要造林树种马尾松、油松、红松和两种常见树种青冈、蒙古栎为研究对象，分析其未来变化趋势。结果表明，未来气候变化情景下，对于马尾松而言，4个模型均预测马尾松在基本保持原有分布的基础上，其未来潜在分布区域均有所扩大，且有向西和向北扩展的趋势;对于油松而言，基于GLM、SGLM和GAM3个模型，油松的未来潜在分布除有北移的趋势外，其分布区还将向东北和西南两个方向扩展;对于红松而言，基于SGLM、GAM和CART3个模型的预测结果较为接近，即红松的未来潜在分布区域将有所减少;对蒙古栎而言，4个模型预测蒙古栎未来分布均将向西扩展;对青冈而言，4个模型预测青冈能基本保持其原有分布区，并向西和向北扩展，其中CART预测结果还表明，青冈在广东南部及广西南部的分布区域将消失。

Proposal for the direct optical detection of pure spin currents in semiconductors

We suggest a different practical scheme for the direct detection of pure spin current by using the two-color Faraday rotation of optical quantum interference process (QUIP) in a semiconductor system. We demonstrate theoretically that the Faraday rotation of QUIP depends sensitively on the spin orientation and wave vector of the carriers, and can be tuned by the relative phase and the polarization direction of the omega and 2 omega laser beams. By adjusting these parameters, the magnitude and direction of the spin current can be detected.

Longitudinal optic phonon-plasmon coupling in delta-doped metamorphic InAlAs/InGaAs high-electron-mobility transistor structures on GaAs substrates

InAlAs/InGaAs metamorphic high-electron-mobility transistor structures with different spacer layers on GaAs substrates are characterized by Raman measurements. The influence of In0.52Al0.48As spacer thickness on longitudinal optic phonon-plasmon coupling is investigated. It is found that the intensity of GaAs-like longitudinal optic phonon, which couples with collective intersubband transitions of two-dimensional electron gas, is strongly affected by the different subband energy spacings, subband electron concentrations, and wave function distributions, which are determined by different spacer thicknesses. (C) 2001 American Institute of Physics.

ELECTRONIC-STRUCTURES OF QUANTUM WIRES FORMED BY LATERAL STRAIN

The electronic structures of quantum wires formed by lateral strain are studied in the framework of the effective-mass envelope-function method. The hole energy levels, wave functions, and optical transition matrix elements are calculated for the real quantum-wire structure, and the results are compared with experiment. It is found that one-dimensional confinement effects exist for both electronic and hole states related to the n (001) = 1 state. The lateral strained confinement causes luminescence-peak redshifts and polarization anisotropy, and the anisotropy is more noticeable than that in the unstrained case. The variation of hole energy levels with well widths in the [110] and [001] directions and wave vector along the [110BAR] direction are also obtained.

Impacts of anthropogenic activity on the recent evolution of the Huanghe (Yellow) River delta

Hydrological statistical data, remote sensing images, and bathymetric charts were used to study the recent evolution of the Huanghe (Yellow) River delta under human-induced interventions. It was clear that water and sediment discharge from the Huanghe River had dropped rapidly since 1970, particularly after 1986. The water and sediment discharges for the period of 1986-2000 were found to have been reduced to only 29.2% and 31.2% of those in the period of 1950-69. This was caused by human factors in the upper and middle reaches of the Huanghe River, including water diversion, damming and reservoir construction, and water and soil conservation. Based on the results from visual interpretation of processed Landsat (MSS or TMJETM+) images dated from 1976 to 2001 and two digital elevation models generated from bathymetric charts surveyed in 1976 and 1992, we found that human-induced reduction of water and sediment discharge led to coastline retrogradation, with the maximum mean recession rate of -0.51 km yr-1 over the period of 1976-98, and seabed erosion beyond the -20 m isobath between 1976 and 1992. Other impacts of human activities on the recent evolution of the Huanghe River delta, including tidal flats shrinking, artificial coastline increasing, land surface sinking and so on, were also analyzed. We found that: (i) the whole delta, including subaerial and subaqueous, has turned from a highly constructive period to a destructive phase; (ii) channelization and dredging were two of the main causes of delta destruction; (iii) land loss in the Huanghe River delta caused by submersion will be increased in the near future; (iv) the Huanghe River delta was becoming more fragile and susceptible to natural hazards.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.

Simulating a typhoon storm surge in the East Sea of China using a coupled model

A coupled numerical model with a 2' x 2' resolution grid has been developed and used to simulate five typical typhoon storm surges (5612, 7413, 7910, 8114, and 9711) in the East Sea of China. Three main driving forces have been considered in this coupled model: wave radiation stress, combined wave-current bottom shear stress and wave-state-dependent surface wind stress. This model has then been compared with in situ measurements of the storm set-up. The effect of different driving force components on the total storm surge has also been investigated. This study has found that the coupled model with high resolution is capable of simulating the five typical typhoons better than the uncoupled models, and that the wave-dependent surface wind stress plays an important role in typhoon storm surge-wave coupling in this area and can increase the storm set-up by 1 m. The study of the five typhoon cases has shown that the general coupling effects could increase storm set-up by 20-32%. Thus, it is suggested that to predict typhoon storm surges in the East Sea of China, a storm surge-wave coupled model be adopted. (C) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.