首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
The Bering Sea circulation is derived as a variational inverse of hydrographic profiles(temperature and salinity),atmospheric climatologies and historical observation of ocean curents.The important result of this study is estimate of the mean climatological sea surface height(SSH) that can be used as a reference for satellite altimetry sea level anomaly data in the Bering Sea region.Numerical experiments reveal that,when combined with satellite altimetry,the obtained reference SSH effectively constrains a realistic reconstruction of the Amukta Pass circulation.  相似文献   

2.
The inversion of high-resolution geoid anomaly maps derived from satellite altimetry should allow one to retrieve the lithospheric elastic thickness, T e , and crustal density, c . Indeed, the bending of a lithospheric plate under the load of a seamount depends on both parameters, and the associated geoid anomaly is correspondingly dependent on the two parameters. The difference between the observed and modelled geoid signatures is estimated by a cost function, J , of the two variables, T e and c . We show that this cost function forms a valley structure along which many local minima appear, the global minimum of J corresponding to the true values of the lithospheric parameters. Classical gradient methods fail to find this global minimum because they converge to the first local minimum of J encountered, so that the final parameter estimate strongly depends on the starting pair of values ( T e ,   c ). We here implement a non-linear optimization algorithm to recover these two parameters from altimetry data. We demonstrate from the inversion of synthetic data that this approach ensures robust estimates of T e and c by activating two search phases alternately: a gradient phase to find a local minimum of J , and a tunnelling phase through high values of the cost function. The accuracy of the solution can be improved by a search in an iteratively restricted parameter subspace. Applying our non-linear inversion to the Great Meteor Seamount geoid data, we further show that the inverse problem is intrinsically ill-posed. As a consequence, minute geoid (or gravity) data errors can induce large changes in any recovery of lithospheric elastic thickness and crustal density.  相似文献   

3.
联合基于GRACE重力卫星观测资料恢复的重力场模型(EIGEN-GL04S1)和卫星测高推求的平均海面高模型(KMSS04)来构造南极绕极流区域的平均海面动力地形,并利用小波滤波方法去掉短波及噪声信号,进而推算大、中尺度的绕极流。与非卫星重力场模型、同化资料及海洋水文资料确定相应结果的验证分析表明:基于新的卫星重力场模型推算的南极绕极流区域的海面动力地形、PF、SAF和表层流场等都与海洋学结果相吻合,且局部特征更加清晰。表明卫-卫跟踪重力卫星计划确定的地球重力场模型较之以前存在的重力场模型在中长波部分精度有较大提高,从大地测量(从空间)角度来探测南极绕极流已达到较高的精度。  相似文献   

4.
The relationship between gravity and bathymetry in the Pacific Ocean   总被引：2，自引：0，他引：2  
Summary. Surface-ship and satellite derived data have been compiled in new free-air gravity anomaly, bathymetry and geoid anomaly maps of the Pacific Ocean basin and its margin. The maps are based on smoothed values of the gravity anomaly, bathymetry and geoid interpolated on to a 90 × 90 km grid. Each smoothed value was obtained by Gaussian filtering measurements along individual ship and subsatellite tracks. The resulting maps resolve features in the gravity, bathymetry and geoid with wavelengths that range from a few hundred to a few thousand kilometres. The smoothed values of bathymetry and geoid anomaly have been corrected for age. The resulting maps show the Pacific ocean basin is associated with a number of ENE–WSW-trending geoid anomaly highs with amplitudes of about ± 5 m and wavelengths of about 3000 km. The most prominent of these highs correlate with the Magellan seamounts–Marshall Gilbert Islands–Magellan rise and the Hess rise–Hawaiian ridge regions. The correlation between geoid anomaly and bathymetry cannot be explained by models of static compensation, but is consistent with a model in which the geoid anomaly and bathymetry are supported by some form of dynamic compensation. We suggest that the dynamic compensation, which characterizes oceanic lithosphere older than 80 Myr, is the result of mantle convection on scales that are smaller than the lithospheric plates themselves.  相似文献   

5.
This paper describes a method for determining Moho depth, lithosphere thinning factor (γ= 1 − 1/β) and the location of the ocean–continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere. The lateral density changes caused by the elevated geotherm in thinned continental margin and adjacent ocean basin lithosphere produce a significant lithosphere thermal gravity anomaly which may be in excess of −100 mGal, and for which a correction must be made in order to determine Moho depth accurately from gravity inversion. We describe a method of iteratively calculating the lithosphere thermal gravity anomaly using a lithosphere thermal model to give the present-day temperature field from which we calculate the lithosphere thermal density and gravity anomalies. For continental margin lithosphere, the lithosphere thermal perturbation is calculated from the lithosphere thinning factor (γ= 1 − 1/β) obtained from crustal thinning determined by gravity inversion and breakup age for thermal re-equilibration time. For oceanic lithosphere, the lithosphere thermal model used to predict the lithosphere thermal gravity anomaly may be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere. A correction is made for crustal melt addition due to decompression melting during continental breakup and seafloor spreading. We investigate the sensitivity of the lithosphere thermal gravity anomaly and the predicted Moho depth from gravity inversion at continental rifted margins to the methods used to calculate and condition the lithosphere thermal model using both synthetic models and examples from the North Atlantic.  相似文献   

6.
Summary. Three complementary methods for the extraction of the M2 ocean tide using SEASAT altimetry are presented and compared. The first method (that developed by Cartwright & Alcock), which provides 'point measurements'of the tide at the crossovers of the SEASAT repeat orbit ground track, has been applied to a study of the tide in tropical ocean areas. The other two methods involve spatial expansions of M2 in terms of either surface spherical harmonics (in the case of the method developed by Mazzega) or Platzman normal modes of the world ocean. The results obtained by each method from only one month of SEASAT data reproduce many features of the tide represented in recent tidal models, and promise well for satellite altimetry as a future source of tidal knowledge.  相似文献   

7.
Estimating the amount of erosion experienced by a sedimentary basin during its geological history plays a key role in basin modelling. In this paper, we present a novel probabilistic approach to estimate net erosion from porosity–depth data from a single well. Our approach uses a Markov chain Monte Carlo algorithm which readily allows us to deal with imprecise knowledge of the lithology-dependent compaction parameters in a joint inversion scheme using multiple lithologies. The results using synthetic data highlight the advantages of our approach over conventional techniques for net erosion estimation: (a) uncertainties on compaction parameters can be effectively mapped into a probabilistic solution for net erosion; (b) posterior uncertainties are easy to quantify; (c) the joint inversion scheme can automatically reconcile porosity data from different lithologies. Our results also underscore the critical role of prior assumptions on controlling the retrieved estimates for net erosion. Using real data from a well in the Barents Sea, we simulate three possible scenarios of variable prior assumptions on compaction parameters to demonstrate the general applicability of our approach. Strong prior assumptions on the compaction parameters led to unrealistic estimates of net erosion for the target well, indicating the assumptions are probably inappropriate. Our preferred strategy for this dataset is to include additional data to constrain the normal compaction trend of the sediments. This provides a net erosion estimate for the target well of about 2300 m with a standard deviation of 140 m which is in line with previous studies. Finally, we discuss potential guidelines to deal with real applications in which data from normally compacted sediments are not available. One is to use our algorithm as a hypothesis-testing tool to evaluate the results under a large set of assumed compaction parameters. A second is to infer compaction parameters and net erosion simultaneously from the target well porosity data. Although appealing and successful with synthetic data, this strategy provides results which are strongly dependent on the calibration data and the geological history of the sediments sampled by the target well.  相似文献   

8.
The VLF-R (very low frequency-resistivity) data, i.e. the apparent resistivity ( ρ a ) and phase ( φ ) data, were inverted individually and jointly using the VFSA (very fast simulated annealing) global inversion approach. Global inversion results for synthetic data without and with various amounts of random and normally distributed Gaussian noise reveal that the inversion of neither the ρ a nor φ data alone yields the true parameters of the structures. However, the joint inversion of the ρ a and φ data yields very good estimates of the model parameters. Five models, representing typical subsurface structures in the shield areas, are studied here. Various models achieved after 10 VFSA runs were used to compute the mean model and the corresponding covariance and correlation matrices, which were used to estimate the uncertainties in the mean model parameters and correlations between the model parameters. We observe that these correlations follow the physics associated with the problem. VLF-R field data due to a nearly vertical contact structure and a very thick dyke-like structure were also inverted to demonstrate the efficacy of the approach in the delineation of the parameters of 2-D structures.  相似文献   

9.
Among the mantle hotspots present under oceanic areas, a large number are located on—or close to—active oceanic ridges. This is especially true in the slow-spreading Atlantic and Indian oceans. The recent availability of worldwide gravity grids and the increasing coverage of geochemical data sets along active spreading centres allow a fruitful comparison of these data with global geoid and seismic tomography models, and allow one to study interactions between mantle plumes and active slow-spreading ridges. The observed correlations allow us to draw preliminary conclusions on the general links between surficial processes, which shape the detailed morphology of the ridge axes, and deeper processes, active in the upper mantle below the ridge axial domains as a whole. The interactions are first studied at the scale of the Atlantic (the Mid-Atlantic Ridge from Iceland to Bouvet Island) from the correlation between the zero-age free-air gravity anomaly, which reflects the zero-age depth of the ridge axis, and Sr isotopic ratios of ridge axis basalts. The study is then extended to a more global scale (the slow ridges from Iceland to the Gulf of Aden) by including geoid and upper-mantle tomography models. The interactions appear complex, ranging from the effect of large and very productive plumes, almost totally overprinting the long-wavelength segmentation pattern of the ridge, to that of weaker hotspots, barely marking some of the observables in the ridge axial domain. Intermediate cases are observed, in which hotspots of medium activity (or whose activity has gradually decreased) located at some distance from the ridge axis produce geophysical or geochemical signals whose variation along the axis can be correlated with the geometry of the plume head in the upper mantle. Such observations tend to preclude the use of a single hotspot/ridge interaction model and stress the need for additional observations in various plume/ridge configurations.  相似文献   

10.
Velocity and mass fields from a constant-density, near-global ocean model, driven with observed twice-daily surface wind stresses and atmospheric pressures for the period October 1992-September 1993, are used to calculate oceanic excitation functions for the length of day (LOD) and for polar motion (PM), and results are analysed as a function of the frequency band. Variable currents and mass redistributions are both important in determining oceanic excitation functions. For bands with periods longer than one month, wind-driven variability is the primary cause of oceanic excitation signals. At higher frequency bands, larger deviations from the inverted barometer response occur, and pressure-driven signals contribute more significantly to the variance in the excitation functions. Oceanic LOD excitation is generally small compared to that of the atmosphere, except for the 2-10 day band. At these scales, adding oceanic to atmospheric excitation series does not lead to better agreement with the observed LOD, although this result may be related to data quality issues. With regard to the excitation of PM, the ocean is in general as important as the atmosphere at most time scales. Combined oceanic and atmospheric excitation series compare visibly better with geodetic series than do atmospheric series alone, pointing to the ocean as a source of measurable signals in PM.  相似文献   

11.
The ocean geoid can be inferred from the topography of the mean sea surface. Satellite altimeters transmit radar pulses and determine the return traveltime to measure sea-surface height. The ERS-1 altimeter stacks 51 consecutive radar reflections on board the satellite to a single waveform. Tracking the time shift of the waveform gives an estimate of the distance to the sea surface. We retrack the ERS-1 radar traveltimes using a model that is focused on the leading edge of the waveforms. While earlier methods regarded adjacent waveforms as independent statistical events, we invert a whole sequence of waveforms simultaneously for a spline geoid solution. Smoothness is controlled by spectral constraints on the spline coefficients. Our geoid solutions have an average spectral density equal to the expected power spectrum of the true geoid. The coherence of repeat track solutions indicates a spatial resolution of 31  km, as compared to 41  km resolution for the ERS-1 Ocean Product. While the resolution of the latter deteriorates to 47  km for wave heights above 2  m, our geoid solution maintains its resolution of 31  km for rough sea. Retracking altimeter waveform data and constraining the solution by a spectral model leads to a realistic geoid solution with significantly improved along-track resolution.  相似文献   

12.
Summary. Rayleigh-wave phase velocities at very long periods (185–290 s) are investigated and regionalized, taking into account the lateral heterogeneities within ocean plates revealed by earlier studies at shorter periods. The two-station method is applied to a few 'pure-age' oceanic paths, and is shown to be compatible with the average Earth model C2 (Anderson & Hart 1976) below depths of 180 km. Under this assumed oceanic model, regionalized for age above 180 km, continental velocities are then derived from a set of experimental great-circle values, both new or taken from previously published studies. The results basically agree with earlier studies (Dziewonski 1970; Kanamori 1970), although they exhibit less scatter than Kanamori's model. Results are successfully checked against a set of values derived by the two-station method from a pure continental path.
Although the shield velocities are substantially different from the mean oceanic ones, they still fall within the range of variation of oceanic velocities with the age of the plate. This makes velocities derived theoretically from Jordan's (1975a, b) models of deep continent—ocean lateral heterogeneities, inconsistent with the present set of experimental data. Finally, we show that Dziewonski's (1971) model S2 reconciles all experimental seismic data relative to shields, without being significantly different from oceanic models below 240 km.  相似文献   

13.
Viscous and viscoelastic models for a subduction zone with a faulted lithosphere and internal buoyancy can self-consistently and simultaneously predict long-wavelength geoid highs over slabs, short-wavelength gravity lows over trenches, trench-forebulge morphology, and explain the high apparent strength of oceanic lithosphere in trench environments. The models use two different free-surface formulations of buoyancy-driven flows (see, for example, Part I): Lagrangian viscoelastic and pseudo-free-surface viscous formulations. The lower mantle must be stronger than the upper in order to obtain geoid highs at long wavelengths. Trenches are a simple consequence of the negative buoyancy of slabs and a large thrust fault, decoupling the overriding from underthrusting plates. The lower oceanic lithosphere must have a viscosity of less than to24 Pa s in order to be consistent with the flexural wavelength of forebulges. Forebulges are dynamically maintained by viscous flow in the lower lithosphere and mantle, and give rise to apparently stiffer oceanic lithosphere at trenches. With purely viscous models using a pseudo-free-surface formulation, we find that viscous relaxation of oceanic lithosphere, in the presence of rapid trench rollback, leads to wider and shallower back-arc basins when compared to cases without viscous relaxation. Moreover, in agreement with earlier studies, the stresses necessary to generate forebulges are small (∼ 100 bars) compared to the unrealistically high stresses needed in classic thin elastic plate models.  相似文献   

14.
An inverse model is used to infer the circulation in Prydz Bay and its adjacent open ocean using hydrographic data obtained by the cruise of the 7th Chinese National Antarctic Research Expedition (CHINARE-7), 1990/91. Barotropie components are found to be strong in the study area, esp. at the Antarctic Divergence, and from a whole view, the velocity is rather small. In the open ocean, the flow is quasizonal, but outside the bay it shows a tendency of pressing onto the shelf from surface to bottom, and a feature of intensification just east of Fram Bank. We suggest here be the most important place to detect the possibility of the Antarctic Bottom Water formation. The meridional profiles of the distribution indicate a strong (relative to the ambient) core and a slope-trapped part into the bargain. In the southeastern part of the bay, there seems to exist a strong coastal current flowing westward. The computed upwelling centers are mainly situated in the west of the study region, as agrees quite well with  相似文献   

15.
海冰是全球气候变化的敏感指示器,对全球热平衡、大气环流、海洋水循环和温盐平衡起到至关重要的作用。海冰厚度作为海冰中最重要的参数之一,也是最难探测的地球物理参数。在回顾和分析基于仰视声纳、走航观测、电磁感应、微波遥感等方法进行极地海冰厚度探测和应用研究的基础上,着重阐述近年来利用卫星测高技术估算极地海冰厚度的研究现状与趋势,为我国开展相关科研工作提供参考。  相似文献   

16.
Summary A technique is presented for calculating geoid height anomalies over two-dimensional models of Earth structure. The method consists of convolving gravity anomalies over the structure with filters which take into account the finite size of the structure in the third dimension and the curvature of the Earth. Similar filters are also developed for a flat earth case. The method is applied to a sea-surface gravity profile crossing the Tonga-Kermadec trench and is found to give good agreement with a Geos-3 radar altimetry profile in the same region. The example demonstrates that introducing arbitrary offsets in computing gravity anomalies can result in spurious long-wavelength effects in the computed geoid. Comparison of the results obtained using flat earth and spherical earth filters suggests that the effects of the curvature of the Earth only become significant for wavelengths in the gravity field greater than about 1000 km.  相似文献   

17.
本文讨论了大洋风场,风应力和某些大气动力学特征,并介绍了需要计算的近海面大气运动学和动力学的量值,特别是风应力的计算方法,最后给出了利用迄今全球大洋最完善的COADS资料求得的大洋风应力和合成风的计算与分析个例.  相似文献   

18.
Inverse methods are useful tools not only for deriving estimates of unknown parameters of the subsurface, but also for appraisal of the thus obtained models. While not being neither the most general nor the most efficient methods, Bayesian inversion based on the calculation of the Jacobian of a given forward model can be used to evaluate many quantities useful in this process. The calculation of the Jacobian, however, is computationally expensive and, if done by divided differences, prone to truncation error. Here, automatic differentiation can be used to produce derivative code by source transformation of an existing forward model. We describe this process for a coupled fluid flow and heat transport finite difference code, which is used in a Bayesian inverse scheme to estimate thermal and hydraulic properties and boundary conditions form measured hydraulic potentials and temperatures. The resulting derivative code was validated by comparison to simple analytical solutions and divided differences. Synthetic examples from different flow regimes demonstrate the use of the inverse scheme, and its behaviour in different configurations.  相似文献   

19.
Collocated magnetotelluric (MT) and seismic profiling is emerging as a necessary combined approach for deep and near-surface imaging but the resulting experimental data are typically interpreted separately since no production programs exist for multidimensional joint inversion of MT and seismic data. We present a joint 2-D inversion approach for imaging collocated MT and seismic refraction data with cross-gradient structural constraints. We describe the main features of the algorithm and first apply it to synthetic data generated for a hypothetical complex geological model. For the synthetic data, we find that the scheme leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and seismic velocity. We apply the scheme to near-surface field data to test the consistency of a previously suggested resistivity–velocity interrelationship and its potential use for subsurface lithofacies discrimination or structural classification. The MT-seismic relationship is found to be in excellent accord with that derived previously for DC resistivity and seismic data set at the test site. Our results suggest that joint MT-seismic cross-gradient imaging leads to improved characterization of heterogeneous geological targets at near-surface to mantle depths.  相似文献   

20.
The deep seismic reflection profile Western Approaches Margin (WAM) cuts across the Goban Spur continental margin, located southwest of Ireland- This non-volcanic margin is characterized by a few tilted blocks parallel to the margin. A volcanic sill has been emplaced on the westernmost tilted block. The shape of the eastern part of this sill is known from seismic data, but neither seismic nor gravity data allow a precise determination of the extent and shape of the volcanic body at depth. Forward modelling and inversion of magnetic data constrain the shape of this volcanic sill and the location of the ocean-continent transition. The volcanic body thickens towards the ocean, and seems to be in direct contact with the oceanic crust. In the contact zone, the volcanic body and the oceanic magnetic layer display approximately the same thickness. The oceanic magnetic layer is anomalously thick immediately west of the volcanic body, and gradually thins to reach more typical values 40 km further to the west. The volcanic sill would therefore represent the very first formation of oceanic crust, just before or at the continental break-up. The ocean-continent transition is limited to a zone 15 km wide. The continental magnetic layer seems to thin gradually oceanwards, as does the continental crust, but no simple relation is observed between their respective thinnings.  相似文献   

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

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