Lidar study of the characteristics of gravity waves in the mesopause region at a southern low-latitude location

Using the technique developed by Gardner and Voelz [1987. Lidar studies of the nighttime sodium layer over Urbana, Illinois, 2. Gravity waves. Journal of Geophysical Research 92, 4673–4693] and a new method proposed by us, two groups of gravity wave parameters are extracted from 11 years sodium lidar measurements made at a southern low-latitude location. The wave occurrence frequencies, wave parameter distributions, and wave parameter relationships are given and compared with other lidar observations. The different lidar derived results can be attributed to different atmospheric parameters. The characteristics of these two groups of gravity waves are also compared, as well as that predicted by diffusive filtering theory, and we find the gravity wave relationships derived from our method are in better agreement with diffusive filtering theory predictions.     

Parameterization of gravity wave momentum deposition based on nonlinear wave interactions: basic formulation and sensitivity tests

We present a practical method for parameterization of gravity wave drag based on the Medvedev and Klaassen (1995. Journal of Geophysical Research 100, 25,841–25,853) theory of gravity wave spectral evolution and saturation. The only tuning necessary for the scheme involves assumptions about the nature of the source spectrum of subgrid-scale gravity waves, i.e. the wavenumbers, launch heights and amplitudes of the spectral components. In this paper we employ a column model with representative distributions of mean wind and temperature to examine the sensitivity of the parameterized wave drag to the source spectra. For the range of anticipated variability of source spectra in the troposphere the scheme produces plausible results consistent with observations and with theoretical estimates. Computationally, the scheme is as efficient as a multiple-wave Lindzen scheme, and suitable for use in general circulation models (GCMs) of the atmosphere.     

Anisotropy of the horizontal velocity fluctuation field in the large wavenumber region

Mesoscale atmospheric wind, temperature, and density fluctuations are described in terms of their power spectra. In the early 1980s, VanZandt[1] showed that the observed spectrum of the atmospheric fluctuations can result, as in the ocean, from an interna…     

Wavenumber Sampling Issues in 2.5D Frequency Domain Seismic Modelling

There are several important wavenumber sampling issues associated with 2.5D seismic modelling in the frequency domain, which need careful attention if accurate results are to be obtained. At certain critical wavenumbers there exist rapid disruptions in the mainly smooth oscillatory spectra. The amplitudes of these disruptions can be very large, and this affects the accuracy of the inverse Fourier transformed frequency-space domain solution. In anisotropic elastic media there are critical wavenumbers associated with each wave mode—the quasi-P (qP) wave, and the two quasi-shear (qS1 and qS2) waves. A small wavenumber sampling interval is desirable in order to capture the highly oscillatory nature of the wavenumber spectrum, especially at increasing distance from the source. Obviously a small wavenumber sampling interval adds greatly to the computational effort because a 2D problem must be solved for every wavenumber and every frequency. The discretisation should be carried out up to some maximum wavenumber, beyond which the field becomes evanescent (exponentially decaying or diffusive). For receivers close to the source, activity persists beyond the critical wavenumber associated with the minimum shear wave velocity in the model. Fortunately, for receivers well removed from the source, the contribution from the evanescent energy is negligible and so there is no need to sample beyond this critical wavenumber. Sampling at Gauss–Legendre spacings is a satisfactory approach for acoustic media, but it is not practical in elastic media due to the difficulty of partitioning the integration around the different critical wavenumbers. We found to our surprise that in transversely isotropic media, the critical wavenumbers are independent of wave direction, but always occur at those wavenumbers corresponding to the maximum phase velocities of the three wave modes (qP, qS1 and qS2), which depend only on the elastic constants and the density. Additionally, we have observed that intermediate layers between source and receiver can filter out to a large degree, the sharp irregularities around the critical wavenumbers in the ω–k _y spectra. We have found that, using the spectral element method, the singularities (poles) at the critical wavenumbers which exist with analytic solutions, do not arise. However, the troublesome spike-like behaviour still occurs and can be damped out without distorting the spectrum elsewhere, through the introduction of slight attenuation.     

Inertia gravity wave activity in the troposphere and lower stratosphere observed by Wuhan MST radar

The troposphere and lower stratosphere (TLS) is a region with active atmospheric fluctuations. The Wuhan Mesosphere-Stratosphere-Troposphere (MST) radar is the first MST radar to have become operational in Mainland China. It is dedicated to real-time atmospheric observations. In this paper, two case studies about inertia gravity waves (IGWs) derived from three-dimensional wind field data collected with the Wuhan MST radar are presented. The intrinsic frequencies, vertical wavelengths, horizontal wavelengths, vertical wavenumber spectra, and energy density are calculated and analyzed. In this paper, we also report on multiple waves existing in the lower stratosphere observed by the Wuhan MST radar. Lomb-Scargle spectral analysis and the hodograph method were used to derive the vertical wavenumber and propagation direction. Meanwhile, an identical IGW is observed by Wuhan MST radar both in troposphere and lower stratosphere regions. Combining the observations, the source of the latter IGW detected in the TLS would be the jet streams located in the tropopause region, which also produced wind shear above and below the tropopause.     

火箭观测到的中层垂直速度扰动细微结构

使用一枚携带ｃｈａｆｆ的火箭观测到的垂直速度数据，研究中层大气垂直速度扰动细微结构。结果表明，垂直速度扰动有明显的层状结构，各层厚度约１５０－５００ｍ，层间垂直分离距离在几百米到１－３ｋｍ间。在垂直速度方差和低Ｒｉ数峰值间以及动力不稳定和ＭＳＴ雷达回波区域间有好的对应关系。水平速度垂直波数谱与饱和谱在谱斜率和谱振幅上存在好的一致．这些观测结果表明，火箭测量到的细微结构可以用波场饱和解释。     

Latitudinal and seasonal variability of the micrometeor input function: A study using model predictions and observations from Arecibo and PFISR

In this work, we use a semi-empirical model of the micrometeor input function (MIF) together with meteor head-echo observations obtained with two high power and large aperture (HPLA) radars, the 430 MHz Arecibo Observatory (AO) radar in Puerto Rico (18°N, 67°W) and the 450 MHz Poker flat incoherent scatter radar (PFISR) in Alaska (65°N, 147°W), to study the seasonal and geographical dependence of the meteoric flux in the upper atmosphere. The model, recently developed by Janches et al. [2006a. Modeling the global micrometeor input function in the upper atmosphere observed by high power and large aperture radars. Journal of Geophysical Research 111] and Fentzke and Janches [2008. A semi-empirical model of the contribution from sporadic meteoroid sources on the meteor input function observed at arecibo. Journal of Geophysical Research (Space Physics) 113 (A03304)], includes an initial mass flux that is provided by the six known meteor sources (i.e. orbital families of dust) as well as detailed modeling of meteoroid atmospheric entry and ablation physics. In addition, we use a simple ionization model to treat radar sensitivity issues by defining minimum electron volume density production thresholds required in the meteor head-echo plasma for detection. This simplified approach works well because we use observations from two radars with similar frequencies, but different sensitivities and locations. This methodology allows us to explore the initial input of particles and how it manifests in different parts of the MLT as observed by these instruments without the need to invoke more sophisticated plasma models, which are under current development. The comparisons between model predictions and radar observations show excellent agreement between diurnal, seasonal, and latitudinal variability of the detected meteor rate and radial velocity distributions, allowing us to understand how individual meteoroid populations contribute to the overall flux at a particular location and season.     

A review of: “Problems of Solar and Stellar Oscillations. Iau Colloquium No. 66. Edited by D. O. Gough. D. Reidel Publishing Company, 493 pp. Df1195.00/U.S. $84.50 (ISBN 90-277-1554-8). 1983.”

Abstract

The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed.     

广义散射层析成像反演

本文详细地给出了基于非均匀介质的体散射广义散射层析成像反演的基本理论.广义散射层析成像反演可以描述为波场的反向传播和对成像场进行局部波数域滤波的过程.在数值算例中,利用背景速度沿深度方向均匀变化的v（z）介质中的简单的方块作为速度异常体的模型,通过对该模型产生的低频的Born数据和声波的正演数据的测试,在对采集系统进行有限频率带宽和空间孔径的校正来进行局部成像矩阵谱的恢复中,可以看出模型中各点的谱在恢复后的质量无论从覆盖的面积范围还是幅值的均一性上都有着明显的提高;在对速度模型的重建中,广义散射层析成像反演能够很好地恢复速度模型的低频分量,即便是方块速度异常体相对于背景速度的平均速度扰动是23%也能很好地重建模型中的速度,且对于不同的背景速度模型基本上都能很好地恢复Marmousi速度模型的低频分量.所以该方法将基于Born模型的层析成像反演适应范围进行了一定程度的扩展.     

An alternative construction of normalized seismic design spectra for near-fault regions

This paper presents a methodology for constructing seismic design spectra in near-fault regions.By analyzing the characteristics of near-fault pulse-type ground motions,an equivalent pulse model is proposed,which can well represent the characteristics of the near-fault forward-directivity and fling-step pulse-type ground motions.The normalized horizontal seismic design spectra for near-fault regions are presented using recorded near-fault pulse-type ground motions and equivalent pulse-type ground motions,which are derived based on the equivalent pulse model coupled with ground motion parameter attenuation relations.The normalized vertical seismic design spectra for near-fault regions are obtained by scaling the corresponding horizontal spectra with the vertical-to-horizontal acceleration spectral ratios of near-fault pulse-type ground motions.The proposed seismic design spectra appear to have relatively small dispersion in a statistical sense.The seismic design spectra for both horizontal and vertical directions can provide alternative spectral shapes for seismic design codes.     

板内应力场研究新成果的总结：JGR世界应力图项目专辑简介

ＪＧＲ杂志９７卷Ｂ８期是世界应力图项目的专辑。该专辑总结了１９８６年至１９９０年期间国际合作研究世界板内应力场的新成果，本文对这些成果作了简要介绍。     

THE RELATIVE CONTRIBUTION OF SOIL TILLAGE AND OVERLAND FLOW EROSION TO SOIL REDISTRIBUTION ON AGRICULTURAL LAND

This study uses evidence for the long-term (35 years) pattern of soil redistribution within two agricultural fields in the UK to identify the relative importance of tillage and overland flow erosion. Spatially distributed long-term total soil redistribution data for the fields (Dalicott Farm and Rufford Forest Farm) were obtained using the caesium-137 (¹³⁷Cs) technique. These data were compared with predicted patterns of soil redistribution. Recent studies have demonstrated that the redistribution of soil by tillage may be described as a diffusive process. A two-component model was, therefore, developed which accounts for soil redistribution by both overland flow and diffusive processes. Comparison of the predicted patterns of overland flow erosion alone with the observed (¹³⁷Cs-derived) data indicated a poor agreement (r² = 0.17 and 0.11). In contrast, a good agreement exists between the predicted pattern of diffusive redistribution and the observed data (r² = 0.43 and 0.41). These results give a clear indication that diffusive processes are dominant in soil redistribution within these fields. Possible diffusive processes include splash erosion, soil creep and tillage. However, the magnitude of the diffusion coefficients for the optimum predicted pattern (c. 350–400 kg m⁻¹ a⁻¹) demonstrates that tillage is the only process capable of explaining the very significant soil redistribution which is indicated by the ¹³⁷Cs data. Consideration is given to the implications of these results for both soil erosion prediction and landscape interpretation.     

Statistics of gravity wave spectra in the troposphere and lower stratosphere over Beijing

We utilize the temperature profiles with a height resolution of 50-m obtained over the Beijing Observatory in the period between January of 2002 and December of 2002 to study vertical wavenumber spectra of normalized temperature fluctuations in the 1.67-8.02 km and 13.57-19.92 km altitude ranges and compare them with linear saturation model.Results indicate that individual vertical wavenumber spectra reveal a considerable variability in both slope and amplitude.The observed variability is not consistent wit...     

集合资料同化中方差滤波技术研究及试验

本文基于YH4DVAR业务系统构建了集合资料同化试验平台,利用10个集合样本统计得到的流依赖背景误差能显著改进业务应用中背景误差方差的结构和大小.但是受样本数的限制,背景误差方差的集合估计值中引入了大量的随机取样噪声.为了降低噪声对估计值的影响,本文采用谱滤波方法,根据信号和噪声尺度的统计特征构造一个低通滤波器来滤除背景误差方差估计值中的大部分随机取样噪声.在2013年第九号台风"飞燕"的集合方差滤波试验中,10个样本的滤波结果优于30个样本的集合估计值.谱滤波方法的成功应用有效降低了集合资料同化系统对集合样本数的要求,将是集合资料同化系统未来业务化运行的一项不可或缺的关键技术.     

Distance scaling of vertical and horizontal response spectra of the Loma Prieta earthquake

Attenuations of the vertical and horizontal response spectra of the 17 October 1989 Loma Prieta, California, earthquake are developed through analyses of the ground motion at 53 sites within a 100 km radius of the source. The analyses are performed on the spectral ordinates for 16 incremental periods ranging from 0.05 to 2.0 sec. The response spectra are modelled empirically for two different site conditions characterized by rock and stiff-soil geologies. Data analysis is performed by the application of a non-linear multivariate regression procedure allowing for distance and site factor as independent variables. Variation of the vertical-to-horizontal (V/H) spectral ratios with wave frequency and distance shows the same behaviour as observed previously in the widely separated geographic regions of northeastern Taiwan and east-central Iran. The predicted ratios at sites underlain by stiff soil are generally higher than the commonly used value of 2/3 at high frequencies ( > 5 Hz) in the near-source region (R < 30 km), but reduce to 1/2 or less at longer periods and farther distances. This behaviour is also observed at rock sites; however, it is somewhat less pronounced. With a faster attenuation of spectral ordinates at higher frequencies, the shape of the response spectrum is found to change with distance. As expected, the spectral attenuation with distance is generally higher for the vertical spectrum than for the horizontal spectrum. The difference is particularly significant at the higher-frequency end of spectrum. Site amplification factors for stiff soil with respect to rock geology varies between 1.17 and 1.72 for horizontal spectrum and 1.01 and 1.81 for vertical spectrum. Spectral amplifications at four sites underlain by soft soil and artificial fill, are also evaluated. This is done by a comparison of the observed spectra with those predicted for rock geology at corresponding distances. As expected, the resulting amplification factors at soft-soil sites show significant increase relative to those at sites underlain by rock.     

Solar–geomagnetic activity and Aa indices toward a standard classification

Legrand and Simon [1989. Solar cycle and geomagnetic activity: a review for geophysicists. Part I. The contributions to geomagnetic activity of shock waves and of the solar wind. Annales Geophysicae 7(6), 565–578] classified one century (1868–1978) of geomagnetic activity, using the Mayaud's Aa index, in four classes related to solar activity: (1) the magnetic quiet activity due to slow solar wind flowing around the magnetosphere, (2) the recurrent activity related to high wind speed solar wind, (3) the fluctuating activity related to fluctuating solar wind and (4) the shock activity due to shock events (CME). In this paper, we use this classification to analyse the solar–geomagnetic activity from 1978 to 2005. We found that during the last three decades the level of geomagnetic quiet activity estimated by Aa indices is decreasing: 2003 is the year of the smallest level of quiet geomagnetic activity since 1868. We compare Legrand and Simon's classification with new in situ solar wind data [Richardson, I.G., Cliver, E.W., Cane, H.V., 2000. Sources of geomagnetic activity over the solar cycle: relative importance of coronal mass ejections, high-speed streams, and slow solar wind. Journal of Geophysical research 105(A8), 18,200–18,213; Richardson, I.G., Cane, H.V., 2002. Sources of geomagnetic activity during nearly three solar cycles (1972–2000). Journal of Geophysical Research 107(A8), 1187] and find a rather good agreement. The differences are only due to minor definitions of the extent of the classes. An attempt is made at defining a more precise standard classification of solar phenomena and at defining time scales of these to understand more precisely the geomagnetic signatures of solar activity.     

Extended reflectivity method for modelling the propagation of diffusive–viscous wave in dip‐layered media

The reflectivity method plays an important role in seismic modelling. It has been used to model different types of waves propagating in elastic and anelastic media. The diffusive–viscous wave equation was proposed to investigate the relationship between frequency dependence of reflections and fluid saturation. It is also used to describe the attenuation property of seismic wave in a fluid‐saturated medium. The attenuation of diffusive–viscous wave is mainly characterised by the effective attenuation parameters in the equation. Thus, it is essential to obtain those parameters and further characterise the features of the diffusive–viscous wave. In this work, we use inversion method to obtain the effective attenuation parameters through quality factor to investigate the characteristics of diffusive–viscous wave by comparing with those of the viscoacoustic wave. Then, the reflection/transmission coefficients in a dip plane‐layered medium are studied through coordinate transform and plane‐wave theory. Consequently, the reflectivity method is extended to compute seismograms of diffusive–viscous wave in a dip plane multi‐layered medium. Finally, we present two models to simulate the propagation of diffusive–viscous wave in a dip plane multi‐layered medium by comparing the results with those in a viscoacoustic medium. The numerical results demonstrate the validity of our extension of reflectivity method to the diffusive–viscous medium. The numerical examples in both time domain and time–frequency domain show that the reflections from a dip plane interface have significant phase shift and amplitude change compared with the results of horizontal plane interface due to the differences in reflection/transmission coefficients. Moreover, the modelling results show strong attenuation and phase shift in the diffusive–viscous wave compared to those of the viscoacoustic wave.     

Array performances for ambient vibrations on a shallow structure and consequences over <Emphasis Type="Italic">V</Emphasis><Subscript><Emphasis Type="Italic">s</Emphasis></Subscript> inversion

Valuable information about one-dimensional soil structures can be obtained by recording ambient vibrations at the surface, in which the energy contribution of surface waves predominates over the one of other types of waves. The dispersion characteristics of surface waves allow the retrieval of the shear-wave velocity as a function of depth. Microtremor studies are usually divided in two stages: deriving the dispersion (or auto-correlation) curve from the recorded signals and inverting it to obtain the site velocity profile. The possibility to determine the dispersion curve over the adequate frequency range at one site depends on the array aperture and on the wavefield spectra amplitude that can be altered by filtering effects due to the ground structure. Microtremors are usually recorded with several arrays of various apertures to get the spectral curves over a wide frequency band, and different methods also exist for processing the raw signals. With the objective of defining a strategy to achieve reliable results for microtremor on a shallow structure, we analyse synthetic ambient vibrations (vertical component) simulated with 333 broadband sources for a 25-m deep soil layer overlying a bedrock. The first part of our study is focused on the determination of the reliable frequency range of the spectral curves (dispersion or auto-correlation) for a given array geometry. We find that the wavenumber limits deduced from the theoretical array re sponse are good estimates of the valid spectral curve range. In the second part, the spectral curves are calculated with the three most popular noise-processing techniques (frequency–wavenumber, high-resolution frequency–wavenumber and spa tial auto-correlation methods) and inverted indi vidually in each case. The inversions are performed with a tool based on the neighbour hood algorithm that offers a better estimation of the global uncertainties than classical linearised methods, especially if the solution is not unique. Several array apertures are necessary to construct the dispersion (auto-correlation) curves in the appropriate frequency range. Considering the final velocity profiles, the three tested methods are almost equivalent, and no significant advantage was found for one particular method. With the chosen model, all methods exhibit a penetration limited to the bedrock depth, as a consequence of the filtering effect of the ground structure on the vertical component, which was observed in numerous shallow sites.     

Spectra of atmospheric waves in the dynamo region of the ionosphere at Kamchatka

A spectral analysis of the diurnal variations in the geomagnetic field horizontal component, observed at Kamchatka and Barrow polar observatory in September–October 1999, has been performed. The complete set of oscillations of thermal tidal atmospheric waves with T = 24, 12, 8, and 4 h has been detected in the variation spectral power (Sq) at Kamchatka, and only the fundamental harmonic with T = 24 h has been distinguished at Barrow. The above periods vary in both directions relative to stable maximums during strong geomagnetic disturbances. The relative spectral intensity at subharmonics also vary toward the fundamental harmonic with a period of 24 h. In the frequency band 0.5–3 h (IGW periods), the maximal intensity in the background spectra is observed at T ～ 2 h and increases by an order of magnitude with increasing geomagnetic activity at both Kamchatka and Barrow. A day before earthquakes, the intensity of this maximum is below the rms background values, and the spectra widen toward the region of periods shorter than 2 h. A similar effect was previously observed in the power spectra of the diurnal variations in the quasistatic electric field and VLF noise, simultaneously measured in September–October 1999.     

Rotational and magnetic instability in the diffusive tachocline

In this article we study the linear instability of the two-dimensional strongly stratified model for global MHD in the diffusive solar tachocline. Gilman and Fox [Gilman, P.A. and Fox, P., Joint instability of the latitudinal differential rotation and toroidal magnetic fields below the solar convection zone. Astrophys. J., 1997, 484, 439–454] showed that for ideal MHD, the observed surface differential rotation becomes more unstable than is predicted by Watson's [Watson, M., Shear instability of differential rotation in stars. Geophys. Astrophys. Fluid Dyn., 1981, 16, 285–298] nonmagnetic analysis. They showed that the solar differential rotation is unstable for essentially all reasonable values of the differential rotation in the presence of an antisymmetric toroidal field. They found that for the broad field case B _φ～sinθcosθ, θ being the co-latitude, instability occurs only for the azimuthal m?=?1 mode, and concluded that modes which are symmetric (meridional flow in the same direction) about the equator onset at lower field strengths than the antisymmetric modes. We study the effect of viscosity and magnetic diffusivity in the strongly stably stratified case where diffusion is primarily along the level surfaces. We show that antisymmetric modes are now strongly preferred over symmetric modes, and that diffusion can sometimes be destabilising. Even solid body rotation can be destabilised through the action of magnetic field. In addition, we show that when diffusion is present, instability can occur when the longitudinal wavenumber m?=?2.