Nonlinear ray perturbation theory with its applications to ray tracing and inversion in anisotropic media

A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters.     

Isotropic and anisotropic correlations in turbulent wind speed data

The correlations in turbulent wind speed data from the open ocean are studied. They are calculated as Fourier transforms of the spectraldensities, which are computed from Fourier transforms of the measured windspeed components. All correlations are found to decay within a few seconds.It is shown that the spectra are most efficiently calculated using data sections whose length does not exceed the observed correlation decay time.Longer runs include more random noise without containing more information.The ratio of the along-wind to transverse spectral levels is investigatedand found to be in excellent agreement with the assumption of isotropy inthe inertial wave number range.     

圆内解析函数在闭圆上连续的条件

研究解析函数的边界性质。通过将Ｐｌｅｓｓｎｅｒ定理的可积条件加强为连续，得到了国内解析函数在闭圆上实部连续导致虚部连续的充分必要条件。     

Modern saltmarsh diatom distributions of the Outer Banks, North Carolina, and the development of a transfer function for high resolution reconstructions of sea level

We collected modern diatom samples from Currituck Barrier Island, Oregon Inlet and Pea Island marshes, Outer Banks, North Carolina, USA, which have different salinity regimes due to their varying distances from a major barrier island inlet. Multivariate analyses separate the saltmarsh diatom assemblages into distinct elevational zones, dominated by differing abundances of polyhalobous, mesohalobous and oligohalobous taxa, suggesting that the distribution of saltmarsh diatoms is a direct function of elevation, with the most important controlling factors being the duration and frequency of subaerial exposure.We developed the first diatom-based transfer function for the east coast of North America to reconstruct former sea levels based upon the relationship between diatom assemblage and elevation. Results imply that this is possible to a precision of ±0.08 m, superior to most similar studies from temperate, mid-latitude environments. The transfer function is used to construct a relative sea-level curve from fossil assemblages from Salvo, North Carolina. These results suggest a sea-level rise of 0.7 m over the last c. 150 years, at an average of c. 3.7 mm year⁻¹. This is consistent with existing sea-level data, and illustrates the utility of the transfer function approach.     

Wave parameters and functions in wavelet analysis

A preliminary study of wave parameters and functions in wavelet analysis is conducted in this paper. The Morlet wavelet transform is used to calculate the time–frequency wavelet energy density function, its volume (i.e. the total energy), its frequency-integral (i.e. the wavelet smoothed instantaneous wave energy history), its time-integral (i.e. the wavelet spectral density function), and two non-dimensional wave indices (NIF, NIT). The processing of the measured wave data obtained from the Chi-Gu coastal observation tower during the period August 2000 to July 2001 indicates that the inter-comparison of wavelet smoothed instantaneous wave energy history and smoothed instantaneous wave energy history (SIWEH) as defined by Funke and Mansard (Proc 17th Int Conf on Ocean Eng, 1980) can reveal the noise structure of the wave signal. The wave data with index NIF greater than 2 is always accompanied with noise, therefore NIF can be used as one of the data quality criteria. The index NIT is linearly correlated with the significant wave period and with the significant wave height, therefore NIT can be used to study the wave growth and decaying phenomena.     

Thermodynamics of gas and steam-blast eruptions

Eruptions of gas or steam and non-juvenile debris are common in volcanic and hydrothermal areas. From reports of non-juvenile eruptions or eruptive sequences world-wide, at least three types (or end-members) can be identified: (1) those involving rock and liquid water initially at boiling-point temperatures (boiling-point eruptions); (2) those powered by gas (primarily water vapor) at initial temperatures approaching magmatic (gas eruptions); and (3) those caused by rapid mixing of hot rock and ground- or surface water (mixing eruptions). For these eruption types, the mechanical energy released, final temperatures, liquid water contents and maximum theoretical velocities are compared by assuming that the erupting mixtures of rock and fluid thermally equilibrate, then decompress isentropically from initial, near-surface pressure (10 MPa) to atmospheric pressure. Maximum mechanical energy release is by far greatest for gas eruptions (1.3 MJ/kg of fluid-rock mixture)-about one-half that of an equivalent mass of gunpowder and one-fourth that of TNT. It is somewhat less for mixing eruptions (0.4 MJ/kg), and least for boiling-point eruptions (0.25 MJ/kg). The final water contents of crupted boiling-point mixtures are usually high, producing wet, sloppy deposits. Final erupted mixtures from gas eruptions are nearly always dry, whereas those from mixing eruptions vary from wet to dry. If all the enthalpy released in the eruptions were converted to kinetic energy, the final velocity (v _max) of these mixtures could range up to 670 m/s for boiling-point eruptions and 1820 m/s for gas eruptions (highest for high initial pressure and mass fractions of rock (m _r) near zero). For mixing eruptions, v _max ranges up to 1150 m/s. All observed eruption velocities are less than 400 m/s, largely because (1) most solid material is expelled when m _r is high, hence v _max is low; (2) observations are made of large blocks the velocities of which may be less than the average for the mixture; (3) heat from solid particles is not efficiently transferred to the fluid during the eruptions; and (4) maximum velocities are reduced by choked flow or friction in the conduit.     

运用正交函数拟合预测新疆阿舍勒地区块状硫化物矿床

新疆阿舍勒块状硫化物铜锌矿床,是一个正在勘探的大型富铜矿床,矿体呈厚大透镜状,具层控特征,矿体和围岩之间密度差较大。地面重力方法,理应是最有效的勘探方法,但是,由于阿舍勒矿区地质构造过分复杂,地下不同深度的多种密度体叠合在一起,往往掩盖了矿体所产生的重力异常。为了解决这个疑难问题,通过研究,选择适当的正交函数去进行拟合,实现目标信息分离,将成矿层位的重力异常从实测重力异常值中分离出来,勾绘出了矿区线状构造成矿带的轮廓,揭示了矿区构造;指出找矿方向,圈出１１个重力异常,根据矿区矿石的最低密度大于围岩最高密度的实事,可以说：“凡有矿的地方必然显示出重力异常”。故使用本方法处理重力资料,在阿舍勒这类矿区,可具直接找矿效果。     

Source history of the 1905 great Mongolian earthquakes (Tsetserleg, Bolnay)

Two great Mongolian earthquakes, Tsetserleg and Bolnay, occurred on 1905 July 9 and 23. We determined the source history of these events using body waveform inversion. The Tsetserleg rupture (azimuth N60°) correspond to a N60° oriented branch of the long EW oriented Bolnay fault.
Historical seismograms recorded by Wiechert instruments are digitized and corrected for the geometrical deformation due to the recording system. We use predictive filters to recover the signals lost at the minute marks.
The total rupture length for the Tsetserleg earthquake may reach up to 190 km, in order to explain the width of the recorded body waves. This implies adding 60 km to the previously mapped fault. The rupture propagation is mainly eastward. It starts at the southwest of the central subsegment, showing a left lateral strike-slip with a reverse component. The total duration of the modelled source function is 65 s. The seismic moment deduced from the inversion is 10²¹ N m, giving a magnitude   M _w = 8  .
The nucleation of the Bolnay earthquake was at the intersection between the main fault (375 km left lateral strike-slip) and the Teregtiin fault (N160°, 80 km long right lateral strike-slip with a vertical component near the main fault). The rupture was bilateral along the main fault: 100 km to the west and 275 km to east. It also propagated 80 km to the southeast along the Teregtiin fault. The source duration was 115 s. The moment magnitude Mw varies between 8.3 and 8.5.
The nucleation and rupture depths remain uncertain. We tested three cases: (1) nucleation and rupture depth limited to the seismogenic zone; (2) nucleation in the seismogenic zone and rupture propagation going to the base of the crust and (3) nucleation within the crust–upper mantle interface and rupture propagation within the upper mantle.     

Crustal thickness variations in the margins of the Gulf of California from receiver functions

Receiver functions (RFs) from teleseismic events recorded by the NARS-Baja array were used to map crustal thickness in the continental margins of the Gulf of California, a newly forming ocean basin. Although the upper crust is known to have split apart simultaneously along the entire length of the Gulf, little is known about the behaviour of the lower crust in this region. The RFs show clear P -to- S wave conversions from the Moho beneath the stations. The delay times between the direct P and P -to- S waves indicate thinner crust closer to the Gulf along the entire Baja California peninsula. The thinner crust is associated with the eastern Peninsular Ranges batholith (PRB). Crustal thickness is uncorrelated with topography in the PRB and the Moho is not flat, suggesting mantle compensation by a weaker than normal mantle based on seismological evidence. The approximately W–E shallowing in Moho depths is significant with extremes in crustal thickness of ∼21 and 37 km. Similar results have been obtained at the northern end of the Gulf by Lewis et al., who proposed a mechanism of lower crustal flow associated with rifting in the Gulf Extensional Province for thinning of the crust. Based on the amount of pre-Pliocene extension possible in the continental margins, if the lower crust did thin in concert with the upper crust, it is possible that the crust was thinned during the early stages of rifting before the opening of the ocean basin. In this case, we suggest that when breakup occurred, the lower crust in the margins of the Gulf was still behaving ductilely. Alternatively, the lower crust may have thinned after the Gulf opened. The implications of these mechanisms are discussed.