首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 34 毫秒
1.
This article is a continuation of the previous one considering the evolution of a long-period comet. The changes of the local radius and the moment of inertia as well as the orientation of the rotation axis are calculated. With reference to the previous model, the time-dependent orbital parameters are introduced. The procedure of sublimation-driven evolution is implemented to the comets 67P/Churyumov-Gerasimenko, 9P/Tempel 1, and 81P/Wild. The inclination of the rotation axis of Comet Ch-G is calculated. The position of the rotation axis is a result of the best fit of the water production rate curve received from modelling to the observational water production rate curve. The method is verified through application to Comets 9P/Tempel 1 and 81P/Wild with well-known positions of the rotation axes. The best fit is for inclination of the rotation axis which is close to the position of the rotation axis with minimum energy (maximum of the moment of inertia). For the best fitted position of the rotation axis I = 90° and Φ = 60°, the largest decrease of radius was about 5.6 m in the northern polar region. The smallest decrease of radius, of about 0.3 m, was noticed on cometographic latitudes between 39° and 46°.  相似文献   

2.
3.
Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 104  s−1m−2), followed by the hot central core (0.1–0.4 × 104 J s−1m−2) and cold annulus (0.04–0.1 × 104 J s−1m−2). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s−1), but the cold annulus contributed equal power (2.41–7.07 MJ s−1) as the hot central core (2.68–6.92 MJ s−1) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.  相似文献   

4.
There is a general agreement that forest ecosystems in the Northern Hemisphere function as signifi-cant sinks for atmospheric CO2; however, their magnitude and distribution remain large uncertainties. In this paper, we report the carbon (C) stock and its change of vegetation, forest floor detritus, and mineral soil, annual net biomass increment and litterfall production, and respiration of vegetation and soils between 1992 to 1994, for three temperate forest ecosystems, birch (Betula platyphylla) forest, oak (Quercus liaotungensis) forest and pine (Pinus tabulaeformis) plantation in Mt. Dongling, Beijing, China. We then evaluate the C budgets of these forest ecosystems. Our results indicated that total C density (organic C per hectare) of these forests ranged from 250 to 300 t C ha-1, of which 35―54 t C ha-1 from vegetation biomass C and 209―244 t C ha-1 from soil organic C (1 m depth, including forest floor detritus). Biomass C of all three forests showed a net increase, with 1.33―3.55 t C ha-1 a-1 during the study period. Litterfall production, vegetation autotrophic respiration, and soil heterotrophic respira-tion were estimated at 1.63―2.34, 2.19―6.93, and 1.81―3.49 t C ha-1 a-1, respectively. Ecosystem gross primary production fluctuated between 5.39 and 12.82 t C ha-1 a-1, about half of which (46%―59%, 3.20―5.89 t C ha-1 a-1) was converted to net primary production. Our results suggested that pine forest fixed C of 4.08 t ha-1 a-1, whereas secondary forests (birch and oak forest) were nearly in balance in CO2 exchange between the atmosphere and ecosystems.  相似文献   

5.
Receiver function study in northern Sumatra and the Malaysian peninsula   总被引:1,自引:0,他引:1  
In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V S) near the surface in the range of 3.4–3.6 km s − 1 attributed to crystalline rocks and 3.6–4.0 km s − 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V S of 2.6–2.9 km s − 1 indicate sediment layers. High V S of 4.2 km s − 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V S between 10 and 20 km depth. Within the subducting slab V S ≈ 4.7 km s − 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.  相似文献   

6.
Besides generating seismic waves, which eventually dissipate, an earthquake also generates a static displacement field everywhere within the Earth. This global displacement field rearranges the Earth’s mass thereby causing the Earth’s rotation and gravitational field to change. The size of this change depends upon the magnitude, focal mechanism, and location of the earthquake. The Sumatran earthquake of December 26, 2004 is the largest earthquake to have occurred since the 1960 Chilean earthquake. Using a spherical, layered Earth model, the coseismic effect of the Sumatran earthquake upon the Earth’s length-of-day, polar motion, and low-degree harmonic coefficients of the gravitational field are computed. Using a model of the earthquake source that is composed of five subevents having a total moment-magnitude M w of 9.3, it is found that this earthquake should have caused the length-of-day to decrease by 6.8 microseconds, the position of the Earth’s generalized figure axis to shift 2.32 milliarcseconds towards 127° E longitude, the Earth’s oblateness J 2 to decrease by 2.37 × 10−11 and the Earth’s pear-shapedness J 3 to decrease by 0.63 × 10−11. The predicted change in the length-of-day, position of the generalized figure axis, and J 3 are probably not detectable by current measurement systems. But the predicted change in oblateness is perhaps detectable if other effects, such as those of the atmosphere, oceans, and continental water storage, can be adequately removed from the observations.  相似文献   

7.
The characteristic of seismic coda wave attenuation in Yunnan area in 7 frequency-bands range from 1 Hz to 20 Hz was estimated by using the local earthquake's waveform data recorded from 22 Yunnan digital seismic stations.Coda attenuation Q-c1 of each station was firstly calculated by single scattering method. Then, mean free path Le and seismic albedo Bo of each station were calculated, and scattering attenuation Q-1s and intrinsic attenuation Q-1i were separated from total attenuation Q-1t by multiple lapse time window analysis based on the multiple scattering model in uniform random isotropic scattering medium. The attenuating characteristics in Yunnan show that most value of Le are in 10~30 km, with maximal within 2~6 Hz;Bo are about 0.5 at 1~2 Hz, but less than 0.5at other frequency-bands, which means Q-1i is comparable with Q-1s at 1~2 Hz, and after 1~2 Hz, Q-1i is greater than Q-1s and dominates the attenuation process. Q-1c is close to Q-1i at other frequency bands except 1~2 Hz.Results show that Q-1 especially Qs-1 varies spatially, Q-1 in eastern Yunnan zone is a bit higher than in northwestern Yunnan zone;northwestern Yunnan zone higher than southwestern Yunnan zone. Comparing with other results in global, Qs-1 in Yunnan is lower than the global average value among these results, Q-1i is higher than the global average value, and Q-1t lies the middle among these results.  相似文献   

8.
The results from Raman spectroscopy analysis of salt aqueous solutions at −170°C demonstrate that for those clearly sharp iron peaks whose Raman wavenumber is close to each other such as NO 3 and CO 3 2− , their original shape could be restorable by the stripping technique, and that ice’s sharp characteristic peak 3090–3109 cm−1) is steady, while the spectrum band of the complex compound (nCl-[H+-OH]n) chlorine ion combined chemically with water molecule is 3401–3413 cm−1. On the other hand, the research shows that the higher the negative iron concentration, the stronger its Raman characteristic peak intensity and the smaller the ice’s. Based on the number of data and theoretical work, the strong correlation of the molar concentration of negative ion with the S i/S H 2O band area ratio is built up. Moreover, the developed Raman method is successfully used in the component analysis of the field fluid inclusions from Silurian sandstone in Tarim basin.  相似文献   

9.
We used in situ measurements and remote-sensing data sets to evaluate the mass budgets of the Lambert, Mellor and Fisher Glaciers and the basal melting and freezing rates beneath their flowbands on the Amery Ice Shelf. Our findings show the Lambert and Mellor Glaciers upstream of the ANARE Lambert Glacier Basin (LGB) traverse may have positive imbalances of 3.9±2.1 Gt a-1 and 2.1±2.4 Gt a-1, respectively, while the Fisher Glacier is approximately in balance. The upstream region as a whole has a positive imbalance of 5.9±4.9 Gt a-1. The three same glaciers downstream of the ANARE LGB traverse line are in negative imbalance, where the whole downstream region has a negative imbalance of -8.5±5.8 Gt a-1. Overall the mass budgets of the Lambert, Mellor, and Fisher Glaciers are close to bal-ance, and the collective three-glacier system is also nearly in balance with a mass budget of -2.6±6.5 Gt a-1. The significant positive imbalances for the interior basin upstream of the ice-movement stations established in the early 1970s (GL line) reported previously are possibly due to an overestimate of the total accumulation and an underestimate of the ice flux through the GL line. The mean melting rate is -23.0±3.5 m ice a-1 near the southern grounding line, which decreases rapidly downstream, and transitions to refreezing at around 300 km from the southern extremity of the Amery Ice Shelf. Freezing rates along the flowbands are around 0.5±0.1 to 1.5±0.2 m ice a-1. The per-centage of ice lost from the interior by basal melting beneath the flowbands is about 80%±5%. The total basal melting and refreezing beneath the three flowbands is 50.3±7.5 Gt ice a-1 and 7.0±1.1 Gt ice a-1, respectively. We find a much larger total basal melting and net melting than the results for the whole Amery Ice Shelf derived from previous modeling and oceanographic measurements.  相似文献   

10.
In this paper we suggest that conditional estimator/predictor of rockburst probability (and rockburst hazard, P T (t)) can be approximated with the formula P T (t) = P 1(θ 1)…P N (θ N P dyn T (t), where P dyn T (t) is a time-dependent probability of rockburst given only the predicted seismic energy parameters, while P i (θ i ) are amplifying coefficients due to local geologic and mining conditions, as defined by the Expert Method of (rockburst) Hazard Evaluation (MRG) known in the Polish mining industry. All the elements of the formula are (approximately) calculable (on-line) and the resulting P T value satisfies inequalities 0 ≤ P T (t) ≤ 1. As a result, the hazard space (0–1) can be always divided into smaller subspaces (e.g., 0–10−5, 10−5–10−4, 10−4–10−3, 10−3–1), possibly named with symbols (e.g., A, B, C, D, …) called “hazard states” — which saves the prediction users from worrying of probabilities. The estimator P T can be interpreted as a formal statement of (reformulated) Comprehensive Method of Rockburst State of Hazard Evaluation, well known in Polish mining industry. The estimator P T is natural, logically consistent and physically interpretable. Due to full formalization, it can be easily generalized, incorporating relevant information from other sources/methods.  相似文献   

11.
Observations of thermohaline properties and currents were undertaken in the Curimataú River estuary (6°18′S), Rio Grande do Norte state (RN), Brazil, during consecutive neap–spring tidal cycles in the austral autumn rainy season. Highly asymmetric neap tide along channel velocities (−0.4 to 0.9 m s−1) and highly stratified conditions were generated by an increase of the buoyancy energy from the freshwater input (R iE≈5.6). During the spring-tidal cycle the river discharge decreased and the longitudinal velocity components were higher, less asymmetrical (−0.8 to 1.1 m s−1) and semidiurnal, associated with moderately stratified conditions (R iE≈0.1) due to the increase of the kinetic tidal energy forcing mechanism. The overall salinity variation from surface to bottom during two tidal cycles was from 20.5 to 36.3 and 29 to 36.7 in the neap and spring tide experiments, respectively; in the last experiment, the tropical water (TW) mass intrusion was enhanced. The net salt transport reversed from down to up estuary during the neap and spring tide experiments, respectively, varied from 6.0 to –2.0 kg m−1 s−1, an indication of changes in the main forcing of the estuary dynamics. Evaluation of a classical steady analytical model, in comparison with nearly steady experimental vertical profiles of velocity, shows an agreement classifiable as reasonably fair.  相似文献   

12.
Based on the stem analysis of 59 individuals of Pinus elliottii in combination with tree biomass models, we calculated annual biomass increment of forest plots at Qianyanzhou Ecological Station, Chinese Academy of Sciences in subtropical China. In addition, canopy layer and community NPP were calculated based on 12 years’ litter fall data. NPP of the 21-year-old forest was estimated by using the BIOME BGC model; and both measured NPP and estimated NPP were compared with flux data. Community biomass was 10574 g · m−2; its distribution patterns in tree layer, shrub layer, herbaceous layer, tree root, herbaceous and shrub roots and fine roots were 7542, 480, 239, 1810, 230, 274 and 239 g · m−2, respectively. From 1999 to 2004, the average annual growth rate and litter fall were 741 g · m−2 · a−1 (381.31 gC · m−2 · a−1) and 849 g · m−2 · a−1 (463 gC · m−2 · a−1), respectively. There was a significant correlation between annual litter fall and annual biomass increment; and the litter fall was 1.19 times the biomass increment of living trees. From 1985 to 2005, average NPP and GPP values based on BGC modeling were 630.88 (343.31–906.42 gC · m−2 · a−1) and 1 800 gC · m−2 · a−1 (1351.62–2318.26 gC · m−2 · a−1). Regression analysis showed a linear relationship (R 2=0.48) between the measured and simulated tree layer NPP values. NPP accounted for 30.2% (25.6%–32.9%) of GPP, while NEP accounted for 57.5% (48.1%–66.5%) of tree-layer NPP and 41.74% (37%–52%) of stand NPP. Soil respiration accounted for 77.0% of measured tree NPP and 55.9% of the measured stand NPP. NEE based on eddy covariance method was 12.97% higher than the observed NEP. Supported by the National Key Basic Research Special Foundation of China (Grant No. 2002CB4125), International Joint Research Project under Ministry of Science and Technology of China (Grant No. 2006DFB91920)  相似文献   

13.
The Luliang and Baoshan basins are two small ba- sins in Yunnan Province. In the recent ten years or so, there have been found a number of natural gas pools of commercial importance in the two basins. Although the gas pools are small in size, the natural …  相似文献   

14.
Deep water originating in the North Atlantic is transported across the Antarctic Circumpolar Current by eddies and, after circumnavigating of the Antarctic, enters the Weddell Gyre south of Africa. As it does so, it rises up from mid-depth towards the surface. The separate temperature and salinity maxima, the Upper and Lower Circumpolar Deep Waters, converge to form the Warm Deep Water. Cores of this water mass on the southern flank of the eastern Weddell Gyre show a change in characteristic as they flow westward in the Lazarev Sea. Observations have been made along four meridional sections at 3° E, 0°, 3° W and 6° W between 60 and 70° S during the Polarstern Cruise ANTXXIII/2 in 2005/2006. These show that a heterogeneous series of warm and salty cores entering the region from the east both north and south of Maud Rise (65° S, 3° W) gradually merge and become more homogeneous towards the west. The gradual reduction in the variance of potential temperature on isopycnals is indicative of isopycnic mixing processes. A multiple regression technique allows diagnosis of the eddy diffusivities and, thus, the relative importance of isopycnic and diapycnic mixing. The method shows that the isopycnic diffusivity lies in the range 70–140 m2 s−1 and the diapycnic diffusivity reaches about 3 × 10−6 m2 s−1. Scale analysis suggests that isopycnic diffusion dominates over diapycnic diffusion in the erosion of the Warm Deep Water cores.  相似文献   

15.
Liquid conductivity (EC) measurement was conducted for the samples collected from several snow pits and ice cores over the Qinghai-Xizang (Tibet) Plateau, with their time range covering seasonal, decadal and centennial scales. Unlike the previous attention mostly focused on the acidity (H+) responding to the solid conductance (ECM) of glacial ice, we introduce the alkalinity (OH) of snow and ice to show how it responds to EC. Strong linear relationship was established between EC and OH for these snow pits and ice cores. Positive correlation is also established between EC and major cations (Ca2+, Mg2+, Na+ and K+). Since the cations are known as the proxies for the intensity of mineral dust influx onto glaciers of the northern Qinghai-Xizang Plateau, we believe that EC could be used as an indicator for the history of dust input in deep ice core study. In fact, records in Guliya ice core since the Little Ice Age (LIA) indicate that dust load in glacier may depend on the combination of temperature and humidity. “Cold-dry” combination favors the dust arising, and results in higher EC and OH values, while “warm-wet” combination prevents dust form and EC and OH values are lower. In the past century, with the atmospheric warming and precipitation increasing over the northern plateau, which means an atmospheric condition of dust decreasing, both EC and OH displayed rapid decline.  相似文献   

16.
Data collected at Somma-Vesuvius during the 1998–1999 radon surveys have been revisited and reinterpreted in light of recent geophysical and geochemical information. The duration of selected radon anomalies, together with the decay properties of radon, have been used to estimate the permeability and porosity of rocks of the deep hydrothermal system. The current local cyclic seismicity is explained by means of a double convective-cell model. Convective cells are separated by a low-permeability horizon located at about 2–2.5 km below sea level. Fluids convecting within the upper cells show temperatures ranging 300–350°C. Rock permeabilities in this sector are estimated on the order of 10−12 m2, for porosities (ϕ) of about 10−5 typical of a brittle environment where fluid velocities may reach ∼800 m/day. Fluid temperatures within the lower cells may be as high as 400–450°C, consistent with supercritical regimes. The hydrodynamic parameters for these cells are lower, with permeability k ∼ 10−15 m2, and porosity ranging from 10−6 to 10−7. Here, fluid motion toward the surface is controlled by the fracture network within a porous medium approaching brittle–ductile behaviour, and fluid velocities may reach ∼1,800 m/day. The low-permeability horizon is a layer where upper and lower convecting cells converge. In this region, fluids (convecting both at upper and lower levels) percolate through the wallrock and release their brines. Due to self-sealing processes, permeability within this horizon reaches critical values to keep the fluid pressure near lithostatic pressure (for k ∼ 10−18 m2). Deep fluid pressure buildups precede the onset of hydrothermally induced earthquakes. Permeability distribution and rock strength do not exclude that the next eruption at Somma-Vesuvius could be preceded by a seismic crisis, eventually leading to a precursory phreatic explosion. The coupling of these mechanisms has the potential of inducing pervasive failure within rocks of the hydrothermal shell, and may be a prelude to a magmatic eruption. It is finally emphasised that the integrated analysis of seismic and geochemical data, including radon emissions, could be successfully used in testing temperature distributions and variations of porosity and permeability in active geothermal reservoirs.  相似文献   

17.
The long-term variations in the second degree sectorial Stokes parameters of the geopotential have been determined from TOPEX-POSEIDON (T/P) satellite altimeter data, covering the period of January 1, 1993 to January 3, 2001 (T/P cycles 11-305). It is the first attempt to determine the variations in the second sectorial harmonics in the Earth’s inertia tensor due to the ocean dynamics. The variations amount to about 1 × 10−10 (J 2 (2) ≈ 1.6 × 10−6 and S 2 (2) ≈ −0.9 × 10−6). The variations are about 5% of the tidal effect. This corresponds to variations in the directions of the equatorial axes of the Earth’s inertia ellipsoid of up to 10 arc-seconds. Consequently, the annual and semi-annual variations of the Earth’s equatorial flattening is about 10−9; i.e. it corresponds to a change of 8 units of its denominator of 91 030. (The equatorial flattening ≈ 1/91 030). Since the coverage of the Earth’s ocean surface is not worldwide, and the inclination of T/P is i = 66°, it is only 58.2% (min. depth of the ocean 2 000 m) of the Earth’s surface which is processed, the torque, resulting from the seasonal transfer of masses within a sea surface layer, is not zero. It amounts up to 1016 kg m2s−2, which is comparable to the total indirect tidal torque due to the Moon and the Sun, ∼ 3.9 × 1016 kg m2s−2. However, the above estimate strongly depends on the adopted thickness of the sea surface layer, ΔR = 50 m. For a larger thickness of ΔR = 100 m, the seasonal torque amounts to about ∼ 2.3 × 1016 kg m2s−2.  相似文献   

18.
During the Pamir Himalayan project in the year 1975 seismic refraction and wide-angle reflection data were recorded along a 270 km long Lawrencepur-Astor (Sango Sar) profile in the northwest Himalayas. The profile starts in the Indus plains and crosses the Main Central Thrust (MCT), the Hazara Syntaxis, the Main Mantle Thrust (MMT) and ends to the east of Nanga Parbat. The seismic data, as published by Guerra et al. (1983), are reinterpreted using the travel-time ray inversion method of Zelt and Smith (1992) and the results of inversion are constrained in terms of parameter resolution and uncertainty estimation. The present model shows that the High Himalayan Crystallines (HHC, velocity 5.4 km s−1) overlie the Indian basement (velocity 5.8–6.0 km s−1). The crust consists of four layers of velocity 5.8–6.0, 6.2, 6.4 and 6.8 km s−1 followed by the upper mantle velocity of 8.2 km s−1 at a depth of about 60 km.  相似文献   

19.
Lower crustal xenoliths brought up rapidly by basaltic magma onto the earth surface may provide di-rect information on the lower crust. The main purpose of this research is to gain an insight into the rheology of the lower crust through the detailed study of lower crustal xenoliths collected from the Hannuoba basalt, North China. The lower crustal xenoliths in this area consist mainly of two pyroxene granulite, garnet granulite, and light-colored granulite, with a few exception of felsic granulite. The equilibration temperature and pressure of these xenoliths are estimated by using geothermometers and geobarometers suitable for lower crustal xenoliths. The obtained results show that the equilibration temperature of these xenoliths is within the range of 785―900℃, and the equilibrium pressure is within the range of 0.8―1.2 GPa, corresponding to a depth range of 28―42 km. These results have been used to modify the previously constructed lower crust-upper mantle geotherm for the studied area. The dif-ferential stress during the deformation process of the lower crustal xenoliths is estimated by using recrystallized grain-size paleo-piezometer to be in the range of 14―20 MPa. Comparing the available steady state flow laws for lower crustal rocks, it is confirmed that the flow law proposed by Wilks et al. in 1990 is applicable to the lower crustal xenoliths studied in this paper. The strain rate of the lower crust estimated by using this flow law is within the range of 10-13―10-11 s-1, higher than the strain rate of the upper mantle estimated previously for the studied area (10-17―10-13 s-1); the equivalent viscosity is estimated to be within the range of 1017―1019Pa·s, lower than that of the upper mantle (1019―1021 Pa·s). The constructed rheological profiles of the lower crust indicate that the differential stress shows no significant linear relation with depth, while the strain rate increases with depth and equivalent vis-cosity decrease with depth. The results support the viewpoint of weak lower continental crust.  相似文献   

20.
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.  相似文献   

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

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