Global methane flux into the atmosphere and its seasonal variations

We have analyzed the regularities in the spatiotemporal behavior of the methane concentration and mass, recorded in different regions of the globe. In the southern hemisphere (SH), the methane concentration does not depend on latitude. In the northern hemisphere (NH), the methane concentration increases towards the Arctic zone. The maximum CH₄ concentration in the fall is due to the action of the additional methane sources in the NH Arctic region. The methane flux calculations showed that in the SH the methane flux into the Earth’s atmosphere is barely season dependent. In the NH during the fall season, the methane flux considerably increases, and the difference between the maximum flux in August-September and the flux in December-January reaches three. The additional, still disregarded sources, sustaining high methane emission in the spring season at high SH latitudes are estimated to be as strong as 100–120 Tg per year. To equilibrate the methane supply to the atmosphere and methane sink in SH, additional, yet unidentified sources of over 300 Tg/yr are required. The global methane flux into the Earth’s atmosphere should be no less than 1050–1100 Tg yearly.     

Shear waves at a corrugated interface between anisotropic elastic and visco-elastic solid half-spaces

Rayleigh’s method of approximation is employed to find out the reflection and transmission coefficients due to an incident plane SH wave at a corrugated interface between a laterally and vertically inhomogeneous anisotropic elastic solid half-space and a laterally and vertically inhomogeneous isotropic visco-elastic solid half-space. The lateral and vertical inhomogeneities are described by the exponential variations of elastic parameters. The formulae of reflection and transmission coefficients are derived in closed form for the first-order approximation of the corrugation. The effects of the corrugation of the interface, the inhomogeneity, the anisotropy, the visco-elasticity and the frequency of the incident wave on these coefficients are studied analytically and numerically for a specific model containing a periodic interface. The results of earlier workers have been reduced as particular cases from the present formulation.     

Macroscopic equations for vapor transport in a multi-layered unsaturated zone

Using the method of homogenization, we present a systematic derivation of the macroscopic equations for air flow and chemical vapor transport in an unsaturated zone with a periodic structure of heterogeneity. The effective specific discharge and hydrodynamic dispersion coefficient are expressible in terms of some cell functions, whose analytical solutions are sought for the simple case of alternate stacking of two strictly plane layers of different properties. For this kind of bi-layered composite, the effects of convection velocity and layer property contrasts on the longitudinal and transverse components of the hydrodynamic dispersion coefficient are investigated.     

Analytical solutions of one-dimensional macrodispersion in stratified porous media by the quadrupole method: convergence to an equivalent homogeneous porous medium

In this article, the quadrupole method is implemented in order to simulate the effects of heterogeneities on one dimensional advective and diffusive transport of a passive solute in porous media. Theoretical studies of dispersion in heterogeneous stratified media can bring insight into transport artefacts linked to scale effects and apparent dispersion coefficients. The quadrupole method is an efficient method for the calculation of transient response of linear systems. It is based here on the Laplace transform technique. The analytical solutions that can be derived by this method assists understanding of upscaled parameters relevant to heterogeneous porous media.First, the method is developed for an infinite homogeneous porous medium. Then, it is adapted to a stratified medium where the fluid flow is perpendicular to the interfaces. The first heterogeneous medium studied is composed of two semi-infinite layers perpendicular to the flow direction each having different transport properties. The concentration response of the medium to a Dirac injection is evaluated. The case studied emphasises the importance in the choice of the boundary conditions.In the case of a periodic heterogeneous porous medium, the concentration response of the medium is evaluated for different numbers of unit-cells. When the number of unit cells is great enough, depending on the transport properties of each layer in the unit cell, an equivalent homogeneous behaviour is reached. An exact determination of the transport properties (equivalent dispersion coefficient) of the equivalent homogeneous porous medium is given.     

Effect of point source and heterogeneity on the propagation of SH-Waves in a viscoelastic layer over a viscoelastic half space

The present paper is concerned with the propagation of shear waves in a homogeneous viscoelastic isotropic layer lying over a semi-infinite heterogeneous viscoelastic isotropic half-space due to point source. The inhomogeneity parameters associated to rigidity, internal friction and density are assumed to be functions of depth. The dispersion equation of shear waves has been obtained using Green’s function technique. The dimensionless angular frequency has been plotted against dimensionless wave number for different values of inhomogeneity parameters. The effects of inhomogeneity have been shown in the dispersion curves. graphical user interface (GUI) software in MATLAB has been developed to show the effect of various inhomogeneity parameters on angular frequency. The topic can be of interest for geophysical applications in propagation of shear waves on the Earth’s crust.     

声波测井仪器的等效理论及其应用(英文)

在处理实际声波测井模式波数据或用模式波的反演地层参数时,必须考虑测井仪器对频散模式波传播的影响。本文介绍了一种用具有等效弹性模量的圆柱体来模拟声波测井仪器的等效理论,并将该理论推广应用到电缆、随钻多极子声波测井。在波长大于仪器半径的前提下,将井中的多极子声场与该等效仪器模型进行声导纳/阻抗匹配,便可导出这一理论。通过数值模拟和现场资料对该理论的有效性和实用性进行了验证,验证结果证明该理论既简单又足够精确。利用该理论模型处理实际资料时,通过调整仪器等效弹性模量和半径在给定频段标定仪器的影响,而不必考虑仪器的具体材料和结构。本文以电缆偶极子声波测井现场资料的频散校正为例,进一步证明了该等效理论处理现场资料的准确性和高效性。     

Structure and physical properties of the garth＇s interior

The studies on the structure and physical properties of the Earth‘s interior done by Chinese geophysicists from 1999 to 2002 were reviewed in this paper. It includes several research areas: the structure of the Earth‘s interiors using seismic tomography, anisotropy of the upper mantle in China and its adjacent areas, quality factor Qβ for S waves, subduction zone, mantle discontinuities, physical properties of Earth‘s materials and others. The review concerns mainly the contents, the methods and the results of the studies. It can be seen that new progress in the study on the structure and physical properties of the Earth‘s interior has been made in the last 4 years in China. It is shown on three aspects: advancement made on some preexistent areas; pioneering on some new fields and new methods adopted.     

Models of sub-grid variability in numerical simulations of solute transport in heterogeneous porous formations: three-dimensional flow and effect of pore-scale dispersion

Modeling transport of contaminants in the earths subsurface relies on numerical solutions over grids with blocks larger than Darcys scale. The hydraulic conductivity is homogenized over the grid blocks and the plumes spreading is reduced as a consequence of the wiped-out variability. To compensate for this loss Rubin et al. (1999) proposed to augment mixing by block-effective dispersion coefficients, and Rubin et al. (2003) showed, by means of two dimensional simulations, how this concept can be applied in practice. In this paper, we present new solutions of the block-effective dispersion tensor for an axisymmetric exponential covariance model. In addition, we discuss the influence of pore-scale dispersion in both two- and three-dimensional applications.     

Study on evaluation of cities’ ability reducing earthquake disasters

Cities’ ability reducing earthquake disasters is a complex system involving numerous factors, moreover the research on evaluating cities’ ability reducing earthquake disasters relates to multi-subject, such as earthquake science, social science, economical science and so on. In this paper, firstly, the conception of cities’ ability reducing earthquake disasters is presented, and the ability could be evaluated with three basic elements — the possible seismic casualty and economic loss during the future earthquakes that are likely to occur in the city and its surroundings and time required for recovery after earthquake; based upon these three basic elements, a framework, which consists of six main components, for evaluating city’s ability reducing earthquake disasters is proposed; then the statistical relations between the index system and the ratio of seismic casualty, the ratio of economic loss and recovery time are gained utilizing the cities’ prediction results of earthquake disasters which were made during the ninth five-year plan; at last, the method defining the comprehensive index of cities’ ability reducing earthquake disasters is presented. Thus the relatively comprehensive theory frame is set up. The frame can evaluate cities’ ability reducing earthquake disasters absolutely and quantitatively and consequently instruct the decision-making on reducing cities’ earthquake disasters loss. Foundation item: State Important Research Project of China (95130603).     

Solar EUV/FUV irradiance variations: analysis and observational strategy

The knowledge of solar extreme and far ultraviolet (EUV) irradiance variations is essential for the characterization of the Earth’s upper atmosphere. For a long time, this knowledge has been based on empirical models, which are themselves based on proxies of the solar activity. However, the accurate modeling and prediction of the Earth’s upper atmosphere necessitate to improve the precision on the irradiance and its variations below about 200 nm. Here, we present a review of recent works made by the authors that aim at quantifying the irradiance variability at these wavelengths, and that lead to new way of monitoring the solar EUV/FUV irradiance spectrum. In more details, it is shown that the quantification of the high level of redundancy in the solar spectrum variability allows to envisage measuring only a small portion of the spectrum without losing essential knowledge. Finally, we discuss what should and could be measured in order to retrieve the solar extreme and far ultraviolet spectrum.     

Solar limb darkening in the wings of the CaII H and K lines

The method for measuring solar limb darkening has been proposed, and the formulas, describing the law of solar limb darkening in the wings of the CaII H (396.849 nm) and K (393.369 nm) lines, have been derived. To operate at arbitrary points on the Sun’s surface, it is necessary to know the law of solar limb darkening, which is specific in different spectral regions. The procedure of spectrum correction for a flat field, proposed in our previous works, is based on comparing line-free spectral regions with the solar center reference spectrum from the atlas by Brault and Neckel [1994].     

Seismic response of L��Aquila downtown from comparison between 2D synthetics spectral ratios of SH, P-SV and Rayleigh waves and observations of the 2009 earthquake sequence

This is the first part of a study on the seismic response of the L’Aquila city using 2D simulation and experimental data. We have studied two velocity-depth models with the aim of outlining the behavior of a velocity reversal in the top layer, which is associated with the stiff Brecce de L’Aquila unit (BrA). In this setting, the SMTH model is topped by a layer with about 2:1 impedance contrast with the underlying layer while the NORV model has no velocity reversal. We have simulated the propagation of SH and P-SV wavefields in the range 0–10 Hz for incidence 0°–90°. Earthquake spectral ratios of the horizontal and vertical components at six sites in L’Aquila downtown are compared to corresponding synthetics spectral ratios. The vertical component of P-SV synthetics enables us to investigate a remarkable amplification effect seen in the vertical component of the recorded strong motion. Sites AQ04 and AQ05 are best matched by synthetics from the NORV model while FAQ5 and AQ06 have a better match with synthetics spectral ratios from the SMTH model. All simulations show this behavior systematically, with horizontal and near-horizontal incident waves predicting the overall pattern of matches more clearly than vertical and near-vertical incidence. The model inferences are in agreement with new geological data reporting lateral passages in the top layer from the stiff BrA to softer sediments. Matches are good in terms of frequency of the first amplification peak and of spectral amplitude: the horizontal components have spectral ratio peaks predominantly at 0.5 Hz in the simulations and at 0.7 Hz in the data, both with amplitudes of 4, while the vertical component spectral ratios reach values of 6 at frequencies of about 1 Hz in both data and simulations. The vertical component spectral ratios are very well matched using Rayleigh waves with incidence at 90°. The NORV model without the velocity reversal predicts spectral ratio peaks for the horizontal components at frequencies up to 6 Hz. The reversal of velocity acts as a low-pass frequency filter on the horizontal components reducing the amplification effect of the sediment filled valley.     

用面波方法研究上扬子克拉通壳幔速度结构

本文研究采用单台法和双台法提取了穿越上扬子的基阶面波相速度和群速度频散;通过对提取的面波群速度和相速度频散进行联合反演,得到的1-D SV速度模型显示上扬子块体下地壳S波速度与典型克拉通区域相当,其上地幔顶部80~170 km深处存在高速的岩石圈盖层,较AK135模型要快2%~3%,其岩石圈厚度约为180 km.在上扬子地区,径向各向异性集中分布在300 km以浅的岩石圈与软流圈部分,其中岩石圈部分SH波比SV波波速要快2%~4%,软流圈部分SH波比SV波波速要快3%~5%;Rayleigh波相速度方位各向异性分析结果显示,上扬子块体周期为25~45 s（大致相当于30~70 km深度范围内）的Rayleigh波相速度存在1.8%~2.7%不等的方位各向异性,其快波方向介于147°~174°.我们认为上扬子块体径向各向异性集中分布在岩石圈、软流圈部分,且各向异性随深度变化, 其岩石圈部分各向异性为大陆克拉通化的遗迹,软流圈部分各向异性与现今板块运动相关.     

Suspended sediment transport along an idealised tidal embayment: settling lag, residual transport and the interpretation of tidal signals

We present semi-analytical solutions for suspended sediment concentration (SSC) and residual sediment transport in a simple mathematical model of a short tidal embayment. These solutions allow us to investigate in some detail the characteristic tidal and semi-tidal variation of SSC and the processes leading to residual sediment transport, including settling and scour lags, the roles of ‘local’ and ‘advective’ contributions, and the presence of internally or externally generated overtides. By interpreting the transport mechanisms in terms of the classic conceptual models of settling lag we clarify how these models may be expressed in mathematical terms. Our results suggest that settling lag is usually a more important process than scour lag, and that a local model which neglects advection may predict the direction of net sediment transport incorrectly. Finally, we discuss our results in the context of other transport processes and morphodynamic feedback.     

Far-zone contributions to the gravity field quantities by means of Molodensky’s truncation coefficients

To reduce the numerical complexity of inverse solutions to large systems of discretised integral equations in gravimetric geoid/quasigeoid modelling, the surface domain of Green’s integrals is subdivided into the near-zone and far-zone integration sub-domains. The inversion is performed for the near zone using regional detailed gravity data. The farzone contributions to the gravity field quantities are estimated from an available global geopotential model using techniques for a spherical harmonic analysis of the gravity field. For computing the far-zone contributions by means of Green’s integrals, truncation coefficients are applied. Different forms of truncation coefficients have been derived depending on a type of integrals in solving various geodetic boundary-value problems. In this study, we utilise Molodensky’s truncation coefficients to Green’s integrals for computing the far-zone contributions to the disturbing potential, the gravity disturbance, and the gravity anomaly. We also demonstrate that Molodensky’s truncation coefficients can be uniformly applied to all types of Green’s integrals used in solving the boundaryvalue problems. The numerical example of the far-zone contributions to the gravity field quantities is given over the area of study which comprises the Canadian Rocky Mountains. The coefficients of a global geopotential model and a detailed digital terrain model are used as input data.     

27.3-day and 13.6-day atmospheric tide

An analysis of time variations in the earth's length of day (LOD) for 25 years (1973-1998) versus at- mospheric circulation changes and lunar phase is presented. It is found that, on the average, there is a 27.3-day and 13.6-day period oscillation in global zonal wind speed, atmospheric geopotential height, and LOD following alternating changes in lunar phase. Every 5-9 days (6.8 days on average), the fields of global atmospheric zonal wind and geopotential height and LOD undergo a sudden change in rela- tion to a change in lunar declination. The observed atmospheric oscillation with this time period may be viewed as a type of atmospheric tide. Ten atmospheric tidal cases have been analyzed by comparing changes in LOD, global zonal wind speed and atmospheric geopotential height versus change in lunar declination. Taken together these cases reveal prominent 27.3-day and 13.6-day tides. The lunar forcing on the earth's atmosphere is great and obvious changes occur in global fields of zonal wind speed and atmospheric geopotential height over the equatorial and low latitude areas. The driving force for the 27.3-day and 13.6-day atmospheric tides is the periodic change in lunar forcing during the moon's revolution around the earth. When the moon is located on the celestial equator the lunar declination equals zero and the lunar tidal forcing on the atmosphere reaches its maximum, at this time the global zonal wind speed increases and the earth's rotation rate decreases and LOD increases. Conversely, when the moon reaches its most northern or southern positions the lunar declination is maximized, lunar tidal forcing decreases, global zonal wind speed decreases, earth's rotation rate increases and LOD decreases. 27.3-day and 13.6-day period atmospheric tides deserve deeper study. Lunar tidal forcing should be considered in models of atmospheric circulation and in short and medium range weather forecasting.     

The structure of the Vladikavkaz Fault Zone based on the study utilizing a complex of geological-geophysical methods

The Vladikavkaz Fault Zone is a regional boundary between the Alpine folded-block system of the Greater Caucasus and Pre-Caucasus foredeep (Ossetian Depression). A complex of geological and geomorphological techniques was used for its study. The data of recent activity of discontinuous and folded dislocations were obtained as a result. Additional application of the microseismic sounding method allowed us to trace the behavior at depth and along the direction of the Vladikavkaz Fault Zone’s branches known on the Earth’s surface from geological and geophysical surveys. An idea of the three-dimensional structure of this system of tectonic faults was made up. It was shown that the structure of the upper section of the fault zone has a form of listric foldthrust, and in the deep layers of the Earth’s crust it is represented by a well-formed contrast zone with subvertical plunging. Expressive deformations of young sediments and a clear manifestation of the fault zone at a depth make it possible to consider the Vladikavkaz Fault as a major seismogenerating structure.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.     

Non-stationary spatial covariance structure estimation in oversampled domains by cluster differences scaling with spatial constraints

In the analysis of spatiotemporal processes underlying environmental studies, the estimation of the non-stationary spatial covariance structure is a well known issue in which multidimensional scaling (MDS) provides an important methodological approach (Sampson and Guttorp in J Am Stat Assoc 87:108–119, 1992). It is also well known that approximating dispersion by a non-metric MDS procedure offers, in general, low precision when accurate differences in spatial dispersion are needed for interpolation purposes, specially if a low dimensional configuration is employed besides a high number of stations in oversampled domains. This paper presents a modification, consisting of including geographical spatial constraints, of Heiser and Groenen’s (Psychometrika 62:63–83, 1997) cluster differences scaling algorithm by which not the original stations but the cluster centres can be represented, while the stations and clusters retain their spatial relationships. A decomposition of the sum of squared dissimilarities into contributions from several sources of variation can be employed for an exploratory diagnosis of the model. Real data are analyzed and differences between several cluster-MDS strategies are discussed.     

The relation between the response spectrum and the phase diference spectrum

ＴｈｅｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｈｅｒｅｓｐｏｎｓｅｓｐｅｃｔｒｕｍａｎｄｔｈｅｐｈａｓｅｄｉｆｅｒｅｎｃｅｓｐｅｃｔｒｕｍＦＥＮＧ－ＸＩＮＺＨＡＯ（赵凤新）ａｎｄＹＵ－ＸＩＡＮＨＵ（胡聿贤）ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，Ｓｔａ...