首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 3 毫秒
1.
  相似文献   

2.
A simple method is described, based on standard VHF wind-profiler data, where imbalances of echo power between four off-vertical radar beams, caused by mountain waves, can be used to calculate the orientation of the wave pattern. It is shown that the mountain wave azimuth (direction of the horizontal component of the wavevector), is given by the vector are radar echo powers, measured in dB, in beams pointed away from vertical by the same angle towards north, south, east and west respectively, and W is the vertical wind velocity. The method is applied to Aberystwyth MST radar data, and the calculated wave vector usually, but not always, points into the low-level wind direction. The mean vertical wind at Aberystwyth, which may also be affected by tilted aspect-sensitive layers, is investigated briefly using the entire radar output 1990–1997. The mean vertical-wind profile is inconsistent with existing theories, but a new mountain-wave interpretation is proposed.  相似文献   

3.
To determine maximum-likelihood magnitudes for a given seismological network requires estimates of the detection thresholds of short-period P for the individual stations. Ideally the detection threshold of a station with fixed instrumentation and operating procedures would be constant. However, some stations reporting to the International Seismological Centre (ISC) show systematic variations with time. Other stations are unreliable or report amplitudes (A) and periods (T) only intermittently. Ways are therefore needed for choosing stations that have a consistent performance over many years. Two rough measures are described here that can be used to look at variation in the apparent detection thresholds with time: one is the number of A/T observations (at distances greater than 20∘) and the other , for the P detections associated each month by the ISC; being the average . Some of the factors that control detection thresholds and clearly revealed by these rough measures are: array stations have lower thresholds than conventional observatories; installation of automatic detectors can lower the threshold; and annual variations in noise level are reflected in annual variations in detection thresholds. For some stations an apparent lowering of the threshold appears not to be due to changes at the station, but is a consequence of the lowering of the ISC reporting thresholds in the region around the station. For stations of the Swiss seismological network there is a slow rise in the threshold from the mid-1970s to 1992 which is difficult to explain but seems to be due to a gradual changes in the detail with which the seismograms are analysed.  相似文献   

4.
The 1999 basaltic eruption of Shishaldin volcano (Alaska, USA) included both Strombolian and Subplinian activity, as well as a “pre-Subplinian” phase interpreted as the local coalescence within a long foam in the conduit. Although few visual observations were made of the eruption, a great deal of information regarding gas velocity, gas flux at the vent and plume height may be inferred by using acoustic recordings of the eruption. By relating acoustic power to gas velocity, a time series of gas velocity is calculated for the Subplinian and pre-Subplinian phases. These time series show trends in gas velocity that are interpreted as plumes or, for those signals lasting only a short time, thermals. The Subplinian phase is shown to be composed of a thermal followed by five plumes with a total expelled gas volume of .The initiation of the Subplinian activity is probably related to the arrival of a large overpressurised bubble close to the top of the magma column. A gradual increase in low-frequency (0.01–0.5 Hz) signal prior to this “trigger bubble” may be due to the rise of the bubble in the conduit. This delay corresponds to a reservoir located at ≈3.9 km below the surface, in good agreement with studies on other volcanoes.The presence of two thermal phases is also identified in the middle of the pre-Subplinian phase with a total gas release of and . Gas velocity at the vent is found to be and for the Subplinian plumes and the pre-Subplinian thermals respectively.The agreement is very good between estimates of the gas flux from modelling the plume height and those obtained from acoustic measurements, leading to a new method by which eruption physical parameters may be quantified. Furthermore, direct measurements of gas velocity can be used for better estimates of the flux released during the eruption.  相似文献   

5.
Krainov  S. R.  Ryzhenko  B. N.  Pavlov  Yu. S. 《Water Resources》2003,30(2):154-176
Synthesis of empirical natural materials and thermodynamic computer modeling of geochemical processes in water–rock systems at different boundary conditions (solid-to-liquid ratio, , T) were used to determine the genetic causes of the inverse geochemical zonality that forms in deep horizons of oil-and-gas bearing structures. The geochemical pattern of inversion water was found to form chiefly because of changes in the Eh–pH-conditions of the original groundwater under the effect of organic components of rocks and because of an increase in temperature to 100°C at low values of solid-to-liquid ratios and at no higher than 10–2 bar.  相似文献   

6.
We present a comparative study of soil CO2 flux () measured by five groups (Groups 1–5) at the IAVCEI-CCVG Eighth Workshop on Volcanic Gases on Masaya volcano, Nicaragua. Groups 1–5 measured using the accumulation chamber method at 5-m spacing within a 900 m2 grid during a morning (AM) period. These measurements were repeated by Groups 1–3 during an afternoon (PM) period. Measured ranged from 218 to 14,719 g m−2 day−1. The variability of the five measurements made at each grid point ranged from ±5 to 167%. However, the arithmetic means of fluxes measured over the entire grid and associated total CO2 emission rate estimates varied between groups by only ±22%. All three groups that made PM measurements reported an 8–19% increase in total emissions over the AM results. Based on a comparison of measurements made during AM and PM times, we argue that this change is due in large part to natural temporal variability of gas flow, rather than to measurement error. In order to estimate the mean and associated CO2 emission rate of one data set and to map the spatial distribution, we compared six geostatistical methods: arithmetic and minimum variance unbiased estimator means of uninterpolated data, and arithmetic means of data interpolated by the multiquadric radial basis function, ordinary kriging, multi-Gaussian kriging, and sequential Gaussian simulation methods. While the total CO2 emission rates estimated using the different techniques only varied by ±4.4%, the maps showed important differences. We suggest that the sequential Gaussian simulation method yields the most realistic representation of the spatial distribution of , but a variety of geostatistical methods are appropriate to estimate the total CO2 emission rate from a study area, which is a primary goal in volcano monitoring research.Editorial responsibility: H Shinohara  相似文献   

7.
Based on P- and S-wave amplitudes and some clear initial P-wave motion data, we calculated focal mechanism solutions of 928 M≥2.5 earthquakes （1994-2005） in four sub-blocks of Sichuan and Yunnan Provinces, namely Sichuan-Qinghai, Yajiang, Central Sichuan and Central Yunnan blocks. Combining these calculation results with those of the focal mechanism solutions of moderately strong earthquakes, we analyzed the stress field characteristics and dislocation types of seismogenic faults that are distributed in the four sub-blocks. The orientation of principal compressive stress for each block is： EW in Sichuan-Qinghai, ESE or SE in Yajiang, Central Sichuan and Central Yunnan blocks. Based on a great deal of focal mechanism data, we designed a program and calculated the directions of the principal stress tensors, σ1, σ2 and σ3, for the four blocks. Meanwhile, we estimated the difference （also referred to as consistency parameter θ^- ） between the force axis direction of focal mechanism solution and the direction of the mean stress tensor of each block. Then we further analyzed the variation of θ^- versus time and the dislocation types of seismogenic faults. Through determination of focal mechanism solutions for each block, we present information on the variation in θ^- value and dislocation types of seismogenic faults.  相似文献   

8.
The magnetoconvection problem under the magnetostrophic approximation is investigated as the nonlinear regime is entered. The model consists of a fluid filled sphere, internally heated, and rapidly rotating in the presence of a prescribed, axisymmetric, toroidal magnetic field. For simplicity only a dipole parity and a single azimuthal wavenumber (m = 2) is considered here. The leading order nonlinearity at small amplitude is the geostrophic flow U g which is introduced to the previously linear model (Walker and Barenghi, 1997a, b). Walker and Barenghi (1997c) considered parameter space above critical and found that U g acts as an equilibration mechanism for moderately supercritical solutions. However, for solutions well above critical a Taylor state is approached and the system can no longer equilibrate. More importantly though, in the context of this paper, is that subcritical solutions were found. Here subcritical solutions are considered in more detail. It was found that, at is strongly dependent on . ( is the critical value of the modified Rayleigh number is a measure of the maximum amplitude of the generated geostrophic flow while , the Elsasser number, defines the strength of the prescribed toroidal field.) Rm at proves to be the key measure in determining how far into the subcritical regime the system can advance.  相似文献   

9.
A generalized turbulent diffusion model has been developed which evaluates the time rate of growth of a simulated cloud of particles released into a turbulent (i.e. diffusive) atmosphere. The general model, in the form of second-order differential equations, computes the three-dimensional size of the cloud as a function of time. Parameters which influence the cloud growth, and which are accounted for in the model equations, are: (1) length scales and velocity magnitudes of the diffusive field, (2) rate of viscous dissipation , (3) vertical stability as characterized by the relative adiabatic lapse rate (1/T)(g/C p +T/z), and (4) vertical shear in the mean horizontal winds , and , for a given height and of spatial extent equal to that of the diffusing cloud. Sample results for near ground level and for upper stratospheric heights are given. For the atmospheric boundary layer case, the diffusive field is microscale turbulence. In the upper stratospheric case it is considered to be a field of highly interactive and dispersive gravity waves.  相似文献   

10.
Riassunto Data una variabile casuale X che segue la legge normale di probabilitl con valor medio a ed error medio y 1'A. considera un'altra variabile casuale che prende il valore intero r quando r–1/2
Summary Given a random variable X following the normal probability law, with expectation a and standard error p, the author considers another random variable , that takes the entire value r when r–1/2  相似文献   

11.
ＬｙａｐｕｎｏｖｅｘｐｏｎｅｎｔａｎｄｄｉｍｅｎｓｉｏｎｏｆｔｈｅｓｔｒａｎｅａｔｔｒａｃｔｏｒｏｆｅｌａｓｔｉｃｆｒｉｃｔｉｏｎａｌｓｙｓｔｅｍＺｈｉ－ＲｅｎＮＩＵ（牛志仁）ａｎｄＤａｎｇ－ＭｉｎＣＨＥＮ（陈党民）（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒ...  相似文献   

12.
This study examined ${\rm NH}_{{\rm 4}}^{{\rm + }} $ , ${\rm PO}_{{\rm 4}}^{{\rm 3}- } $ recovery and the concentration of residual ions from anaerobic effluent of the potato processing industry through magnesium ammonium phosphate (MAP) precipitation using a Box–Behnken design. The regression model was statistically significant in terms of ${\rm NH}_{{\rm 4}}^{{\rm + }} $ and ${\rm PO}_{{\rm 4}}^{{\rm 3}- } $ removal efficiency and residual ion concentrations. Optimum ${\rm NH}_{{\rm 4}}^{{\rm + }} $ and ${\rm PO}_{{\rm 4}}^{{\rm 3}- } $ removal was obtained at pH 9.50 and at Mg2+/${\rm NH}_{{\rm 4}}^{{\rm + }} $ /${\rm PO}_{{\rm 4}}^{{\rm 3}- } $ molar ratio of 1.8:1:1.8. Under these conditions, Mg, Ca, K, Fe, and Cl concentrations required for plant growth significantly decreased with MAP precipitation, which was supported by EDX analysis of dry MAP precipitate. The fertilizer effect of MAP on the growth of corn and tomato plants was compared with chemical fertilizers through pot trials. Nutrient element uptake levels of plants were examined in different fertilizer sources. While Mg, Fe, Cu, Mn, and Zn nutrient element uptake levels were sufficient in MAP pots, Ca uptake exceeded sufficient level. Average levels of N, P, K, Mg, Cu, and Mn of corn plant were higher in MAP than other pots. The average N, P, and Mg levels of tomato plant in MAP pots were higher than other pots. N/K ratio, which is important in tomato plants, was better optimized in MAP pots. Only Ni, Cr, and Pb heavy metals were found in plants.  相似文献   

13.
The intrinsic dissipation and scattering attenuation in southwestern (SW) Anatolia, which is a tectonically active region, is studied using the coda waves. First the coda quality factor (Qc) assuming single scattering is estimated from the slope of the coda-wave amplitude decay. Then the Multiple Lapse Time Window (MLTW) analysis is performed with a uniform earth model. Three non-overlapping temporal data windows are used to calculate the scattered seismic energy densities against the source-receiver distances, which, in turn, are used to calculate separate estimates of the intrinsic and scattering factors. In order to explore the frequency dependency, the observed seismograms are band pass-filtered at the center frequencies of 0.75, 1.5, 3.0, 6.0 and 12.0. The scattering attenuation (Qs−1) is found lower than the intrinsic attenuation (Qi−1) at all frequencies except at 0.75 Hz where the opposite is observed. Overall the intrinsic attenuation dominates over the scattering attenuation in the SW Anatolia region. The integrated energy curves obtained for the first energy window (i.e., 0–15 s) are somewhat irregular with distance while the second (i.e., 15–30 s) and third (i.e., 30–45 s) data windows exhibit more regular change with distance at most frequencies. The seismic albedo B0 is determined as 0.61 at 0.75 Hz and 0.34 at 12.0 Hz while the total attenuation factor denoted by Le−1 changes in the range 0.034–0.017. For the source-station range 20–180 km considered the scattering attenuation is found strongly frequency dependent given by the power law Qs−1 = 0.010*f−1.508. The same relations for Qi−1, Qt−1 (total), Qc−1 and (expected) hold as Qi−1 = 0.0090*f−1.17, Qt−1 = 0.019*f−1.31, Qc−1 = 0.008*f−0.84 and respectively. Compared to the other attenuation factors Qc−1 and are less dependent on the frequency.  相似文献   

14.
Fermat's variational principle states that the signal propagates from point S to R along a curve which renders Fermat's functional (l) stationary. Fermat's functional (l) depends on curves l which connect points S and R, and represents the travel times from S to R along l. In seismology, it is mostly expressed by the integral (l) = (x k,x k ')du, taken along curve l, where (x k,x k ') is the relevant Lagrangian, x k are coordinates, u is a parameter used to specify the position of points along l, and x k ' = dx k÷du. If Lagrangian (x k,x k ') is a homogeneous function of the first degree in x k ', Fermat's principle is valid for arbitrary monotonic parameter u. We than speak of the first-degree Lagrangian (1)(x k,x k '). It is shown that the conventional Legendre transform cannot be applied to the first-degree Lagrangian (1)(x k,x k ') to derive the relevant Hamiltonian (1)(x k,p k), and Hamiltonian ray equations. The reason is that the Hessian determinant of the transform vanishes identically for first-degree Lagrangians (1)(x k,x k '). The Lagrangians must be modified so that the Hessian determinant is different from zero. A modification to overcome this difficulty is proposed in this article, and is based on second-degree Lagrangians (2). Parameter u along the curves is taken to correspond to travel time , and the second-degree Lagrangian (2)(x k, k ) is then introduced by the relation (2)(x k, k ) = [(1)(x k, k )]2, with k = dx k÷d. The second-degree Lagrangian (2)(x k, k ) yields the same Euler/Lagrange equations for rays as the first-degree Lagrangian (1)(x k, k ). The relevant Hessian determinant, however, does not vanish identically. Consequently, the Legendre transform can then be used to compute Hamiltonian (2)(x k,p k) from Lagrangian (2)(x k, k ), and vice versa, and the Hamiltonian canonical equations can be derived from the Euler-Lagrange equations. Both (2)(x k, k ) and (2)(x k,p k) can be expressed in terms of the wave propagation metric tensor g ij(x k, k ), which depends not only on position x k, but also on the direction of vector k . It is defined in a Finsler space, in which the distance is measured by the travel time. It is shown that the standard form of the Hamiltonian, derived from the elastodynamic equation and representing the eikonal equation, which has been broadly used in the seismic ray method, corresponds to the second-degree Lagrangian (2)(x k, k ), not to the first-degree Lagrangian (1)(x k, k ). It is also shown that relations (2)(x k, k ) = ; and (2)(x k,p k) = are valid at any point of the ray and that they represent the group velocity surface and the slowness surface, respectively. All procedures and derived equations are valid for general anisotropic inhomogeneous media, and for general curvilinear coordinates x i. To make certain procedures and equations more transparent and objective, the simpler cases of isotropic and ellipsoidally anisotropic media are briefly discussed as special cases.  相似文献   

15.
Summary The mechanism of beating of Pc3 type pulsations is studied. Using the method of numerical computation of a sonagram (the method of frequency-time analysis) a set of samples of pulsations from the Budkov Observatory is treated (1968–1969) mostly at K-indices equal to 2–3. By comparing f–t diagrams with the spectra of the samples an attempt has been made at interpreting the beating as a superposition of the frequency components, contained in the pulsation signal. In most observed cases it is possible to determine two close frequencies, the difference of which is on the average =5.4 mHz. The average carrier frequency of the samples was =37.6 mHz, and the average frequency of the beating =2.7 mHz. The interval of observed values of fB amounted to 1–5 mHz. A tendency was observed for fB to increase with increasing degree of disturbance of the geomagnetic field.  相似文献   

16.
Summary Seven optimal networks consisting of 4 to 10 stations are compared for a given region, where velocity-depth profiles and the distribution of seismic intensity are known. Assuming that the standard error of arrival time is t =0.05 s and the standard errors of the parameters of velocity-depth profiles are equal to 5% of their values, the average standard errors of the origin time and focus coordinates are estimated. The application of optimum methods to the planning of seismic networks in the Lublin Coal Basin is presented, and maps of standard errors of origin time , depth and epicenter ( xy ) for the case of an optimum network of 6 seismic stations are given.  相似文献   

17.
Predictive relations are developed for peak ground acceleration (PGA) from the engineering seismoscope (SRR) records of the 2001 Mw 7.7 Bhuj earthquake and 239 strong-motion records of 32 significant aftershocks of 3.1 ≤ Mw ≤ 5.6 at epicentral distances of 1 ≤ R ≤ 288 km. We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relation for peak horizontal acceleration in the Kachchh seismic zone, Gujarat. This new analysis uses the Joyner-Boore’s method for a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. The resulting attenuation equation is,
where, Y is peak horizontal acceleration in g, Mw is moment magnitude, rjb is the closest distance to the surface projection of the fault rupture in kilometers, and S is a variable taking the values of 0 and 1 according to the local site geology. S is 0 for a rock site, and, S is 1 for a soil site. The relation differs from previous work in the improved reliability of input parameters and large numbers of strong-motion PGA data recorded at short distances (0–50 km) from the source. The relation is in demonstrable agreement with the recorded strong-ground motion data from earthquakes of Mw 3.5, 4.1, 4.5, 5.6, and 7.7. There are insufficient data from the Kachchh region to adequately judge the relation for the magnitude range 5.7 ≤ Mw ≤ 7.7. But, our ground-motion prediction model shows a reasonable correlation with the PGA data of the 29 March, 1999 Chamoli main shock (Mw 6.5), validating our ground-motion attenuation model for an Mw6.5 event. However, our ground-motion prediction shows no correlation with the PGA data of the 10 December, 1967 Koyna main shock (Mw 6.3). Our ground-motion predictions show more scatter in estimated residual for the distance range (0–30 km), which could be due to the amplification/noise at near stations situated in the Kachchh sedimentary basin. We also noticed smaller residuals for the distance range (30–300 km), which could be due to less amplification/noise at sites distant from the Kachchh basin. However, the observed less residuals for the longer distance range (100–300 km) are less reliable due to the lack of available PGA values in the same distance range.  相似文献   

18.
This article presents a new approach for the implementation of a planar-free surface boundary condition. It is based on a vertical grid-size reduction above the free surface during the explicit computation of a free surface boundary condition. This technique is very much similar to the well-known stress imaging technique. VGR-stress imaging technique name is proposed for this new free surface boundary condition (VGR stands for ‘vertical grid-size reduction’). To study the performance of the proposed VGR-stress imaging technique, it was implemented in a newly developed second order accurate in time and fourth-order accurate in space (2, 4) staggered grid SH-wave finite-difference (FD) algorithm with variable grid size. It was confirmed that the effective thickness (ETH) of first soil layer becomes less by one-half of vertical grid size than the assigned thickness (ATH), if stress imaging technique is used as a free surface boundary condition. The qualitative and quantitative results of various numerical experiments revealed that the proposed VGR-stress imaging technique is better than the stress imaging technique since it is free from the thickness discrepancy arising due to the use of images of stress components across the free surface. On the basis of iterative numerical experiments, it was confirmed that the stability condition for this FD scheme with variable grid size is It was also inferred that at least five to six grid points per shortest wavelength are required to avoid the grid dispersion. The maximum grid-spacing ratio up to 12.5 or even more did not affect the accuracy of (2,4) SH-wave algorithm. The obtained reduction of 10.46 and 5.38 folds in the requirement of computational memory and time for a particular basin-edge model, as compared with the homogeneous grid size, reflects the efficacy of the new FD algorithm.  相似文献   

19.
Summary The application of Euler's postulate , which determines rotation (here denotes the derivative of the angular momentum with respect to time, L the resultant torque of exterior forces relative to the centre of inertia), to a celestial body, which is not spherically symmetric, had aroused doubts in Euler himself because of the displacement of the resultant attractive force of other celestial bodies from the said centre. As Minding had noted, it is more logical to compute the resultant torque relative to a point at which the torque has a minimum. This point should be named the attrahentis centrum after Euler. Newton was aware of this displacement of the resultant force and had mentioned it in the Principia. A review of Euler's works connected with this subject is given.  相似文献   

20.
A new technique for the treatment of the kinematic dynamo problem is presented. The method is applicable when the dynamo is surrounded by a medium of finite conductivity and is based on a reformulation of the induction equation and boundary conditions at infinity into an integral equation. We show that the integral operator involved here is compact in the case of homogeneous conductivity, which is important for both mathematical and numerical treatment. A lower bound for the norm of then yields a necessary condition for the generation of magnetic fields by kinematic dynamos. Numerical results are presented for some simple 2-dynamo models. The far-field asymptotics for stationary and time-dependent field modes are discussed.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号