首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 796 毫秒
1.
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.  相似文献   

2.
Summary. A finite element method with uniform and variable resolution meshes is used to model very high Rayleigh number Ra thermal convection in a square box of infinite Prandtl number, Boussinesq fluid with constant viscosity and thermodynamic properties. Heating is either entirely from below or mostly from within and the boundaries are stress free. The variable mesh is coarse in the interior of the convection cell and it is fine in the very thin boundary layers and plumes surrounding the core. The highest resolution variable mesh has a dimensionless grid spacing of 0.027 in the core and 0.0017 in the boundary layers. The boundary layers contain about 10 mesh points even at the highest values of Ra considered and are thus highly resolved. The variable mesh approach is shown to yield reliable simulations of convection as long as the aspect ratio of the most elongated boundary layer elements is not too large; values of about 4 to 6 work well. This aspect ratio also measures the increase in resolution in the boundary layers as compared with the central core. Steady single-cell rolls are computed for bottom heating and Ra up to 5 × 105 times the marginal instability value of the Rayleigh number Racr. One and two-cell roll solutions are calculated for f = 1, 0.8 and 0.6, where f is the fraction of the heat escaping through the top of the box that is generated internally. The values of Racr for f = 1, 0.8 and 0.6 are 1296, 1024 and 864, respectively. The largest of Ra/Racr at which unicellular convection is stable (steady) are approximately 390, 610 and 970, for f = 1, 0.8 and 0.6.  相似文献   

3.
Summary. Numerical convection models are presented in which plates are simulated by imposing piecewise constant horizontal velocities on the upper boundary. A 4 × 1 box of constant viscosity fluid and two-dimensional (2-D) flow is assumed. Four heating modes are compared: the four combinations of internal or bottom heating and prescribed bottom temperature or heat flux. The case with internal heating and an isothermal base is relevant to lower mantle or whole mantle convection, and it yields a lower thermal boundary layer which is laterally variable and can be locally reversed, corresponding to heat flowing back into the core locally. When scaled to the whole mantle, the surface deflections and gravity and geoid perturbations calculated from the models are comparable to those observed at the Earth's surface. For models with migrating ridges and trenches, the flow structure lags well behind the changing surface 'plate'configurations. This may help to explain the poor correlation between the main geoid features and plate boundaries. Trench migration substantially affects the dip of the cool descending fluid because of induced horizontal shear in the vicinity of the trench. Such shear is small for whole mantle convection, but is large for upper mantle convection, and would probably result in the Tonga Benioff zone dipping to the SE, opposite to the observed dip, for the case of upper mantle convection.  相似文献   

4.
b
Spherical harmonics are orthonormalized using the Gram-Schmidt process in a function space. The problem of linear dependence of spherical harmonics over the oceans is studied using the Gram matrices and consequently three sets of orthonormal (ON) functions have been constructed. For the process an efficient formula for computing inner products of spherical harmonics has been developed. Important spectral properties of the ON functions are addressed. The ON functions may be used for representing the sea surface topography (SST) in the analysis of satellite altimeter data. The geoid error can be transformed to a representation by the ON functions and hence the comparison of powers of the geoid error and the SST signal only over the oceans is possible, leading to a better way of determining the cut-off frequency of the SST in the simultaneous solution using satellite altimeter data. As a case study, the modified Levitus SST is expanded into the ON functions. The results show that 99.90 per cent of that signal's energy is contained within degree 24 of the orthonormal functions. Such expansions also render better spectral behaviour of oceanic signals as compared to that from spherical harmonic expansions. The study shows that these generalized Fourier functions are suitable for spectral analyses of oceanic signals and they can be applied to future altimetric mission where the geoid and the SST are to be recovered.  相似文献   

5.
Summary. A total of 3708 1 × 1° free-air gravity anomaly averages have been used to construct a new 1 × 1° gravimetric geoid of the Northwest Pacific Ocean. The 1 × 1° averages are based on a compilation of 147000 surface ship and pendulum gravity measurements. The gravimetric geoid reveals information in the geoid of the Northwest Pacific not present in currently used satellite derived models. The RMS difference between the 1 × 1° geoid and satellite derived models is about ±6 m. Difference geoid undulations range from a maximum of +19 m over the Hawaiian ridge to a minimum of −31 m over the junction of the Kuril and Aleutian trenches. The Hawaiian swell is associated with a geoidal high of up to +15 m with wavelengths of about 2200 km and the topographic rises seaward of deep-sea trenches are associated with geoidal highs of up to 4m with wavelengths of about 220–900 km. The main difference between the gravimetric geoid and the satellite derived models occurs over the Pacific basin where discrepancies reach +10 m with wavelengths of 4000 km. The agreement between the gravi-metric geoid and Skylab-4 and Geos-3 altimeter data is close for wavelengths greater than about 300 km but poor for shorter wavelengths.  相似文献   

6.
Planetary topography can either be modelled as a load supported by the lithosphere, or as a dynamic effect due to lithospheric flexure caused by mantle convection. In both cases the response of the lithosphere to external forces can be calculated with the theory of thin elastic plates or shells. On one-plate planets the spherical geometry of the lithospheric shell plays an important role in the flexure mechanism. So far the equations governing the deformations and stresses of a spherical shell have only been derived under the assumption of a shell of constant thickness. However, local studies of gravity and topography data suggest large variations in the thickness of the lithosphere. In this paper, we obtain the scalar flexure equations governing the deformations of a thin spherical shell with variable thickness or variable Young's modulus. The resulting equations can be solved in succession, except for a system of two simultaneous equations, the solutions of which are the transverse deflection and an associated stress function. In order to include bottom loading generated by mantle convection, we extend the method of stress functions to include loads with a toroidal tangential component. We further show that toroidal tangential displacement always occurs if the shell thickness varies, even in the absence of toroidal loads. We finally prove that the degree-one harmonic components of the transverse deflection and of the toroidal tangential displacement are independent of the elastic properties of the shell and are associated with translational and rotational freedom. While being constrained by the static assumption, degree-one loads can deform the shell and generate stresses. The flexure equations for a shell of variable thickness are useful not only for the prediction of the gravity signal in local admittance studies, but also for the construction of stress maps in tectonic analysis.  相似文献   

7.
Summary. A new set of 1×1° mean free-air anomalies in the Indian Ocean is determined on the basis of previously published free-air anomaly maps (Talwani & Kahle) and the most recent Lamont surface ship gravity measurements. The data are then used to compute a (total) 1×1° gravimetric Indian Ocean geoid. The computation is carried out by combining the Goddard Space Flight Center (GSFC) GEM-6 geoid and a difference geoid that corresponds to the differences between the set of 1×1° surface gravity values and the GEM-6 gravity anomalies. The difference geoid is highest over the Madagascar Ridge (+ 20 m) and lowest over the Timor Trough (-30 m). The total geoid is compared with GEOS-3 radar altimeter derived geoid profiles and geophysical implications are discussed.  相似文献   

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

9.
Applying the infinite Prandtl number approximation, a semi-analytical solution for computing 2-D axisymmetric viscous Stokes flow in a model consisting of two eccentrically nested spheres of different viscosities is derived. Since numerical codes based on spectral or finite techniques for modelling mantle flow in a spherical geometry in the presence of lateral viscosity variation are becoming more and more popular, reliable examples for testing and validating such codes are extremely useful. The eccentrically nested sphere solution was used to test a numerical algorithm based on a mixed spherical-harmonic finite-element formulation of the Stokes problem, and good agreement was obtained.  相似文献   

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

11.
A new spherical mesh generator is described. It represents an efficient, deterministic packing of tetrahedra into a solid sphere, a spherical shell, or both. The mesh can be used for finite-element solutions to a wide variety of global numerical modelling problems in the geosciences. The nodes within the mesh are distributed uniformly, and long, thin tetrahedra are avoided. The method proposed here offers several advantages over 3-D Delaunay algorithms for finite-element mesh generation. For the related problem of trivariate scattered data interpolation, which is not considered here, the 3-D Delaunay algorithms are the method of choice.  相似文献   

12.
In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.  相似文献   

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

14.
The characteristics of a reflected spherical wave at a free surface are investigated by numerical methods; in particular, the polarization angles and amplitude coefficients of a reflected spherical wave are studied. The classical case of the reflection of a plane P wave from a free surface is revisited in order to establish our terminology, and the classical results are recast in a way which is more suited for the study undertaken. The polarization angle of a plane P wave, for a given angle of incidence, is shown to be 90° minus twice the angle of reflection of the reflected S wave. For a Poisson's ratio less than 1/3, there is a non-normal incident angle for which both amplification coefficients are 2 precisely; for this incident angle the direction of the particle motion at the free surface is also the direction of the incident wave. For a wave emanating from a spherical source, the polarization angle, for all angles of incidence, is always less than, or equal to, the polarization angle of a plane P wave. The vector amplification coefficient of a spherical wave, for all angles of incidence, is always greater than the vector amplification coefficient of a plane P wave. As expected, the results for a spherical wave approach the results for a plane P wave in the far field. Furthermore, there was a good agreement between the theoretical modelling and the numerical modelling using the dynamic finite element method (DFEM).  相似文献   

15.
Modelling dynamic rupture for complex geometrical fault structures is performed through a finite volume method. After transformations for building up the partial differential system following explicit conservative law, we design an unstructured bi-dimensional time-domain numerical formulation of the crack problem. As a result, arbitrary non-planar faults can be explicitly represented without extra computational cost. On these complex surfaces, boundary conditions are set on stress fluxes and not on stress values. Prescribed rupture velocity gives accurate solutions with respect to analytical ones depending on the mesh refinement, while solutions for spontaneous propagation are analysed through numerical means. An example of non-planar spontaneous fault growth in heterogeneous media demonstrates the good behaviour of the proposed algorithm as well as specific difficulties of such numerical modelling.  相似文献   

16.
Observations of gravity can be aliased by virtue of the logistics involved in collecting these data in the field. For instance, gravity measurements are often made in more accessible lowland areas where there are roads and tracks, thus omitting areas of higher relief in between. The gravimetric determination of the geoid requires mean terrain-corrected free-air anomalies; however, anomalies based only on the observations in lowland regions are not necessarily representative of the true mean value over the topography. A five-stage approach is taken that uses a digital elevation model, which provides a more accurate representation of the topography than the gravity observation elevations, to reduce the unrepresentative sampling in the gravity observations. When using this approach with the Australian digital elevation model, the terrain-corrected free-air anomalies generated from the Australian gravity data base change by between 77.075 and −84.335 mgal (−0.193 mgal mean and 2.687 mgal standard deviation). Subsequent gravimetric geoid computations are used to illustrate the effect of aliasing in the Australian gravity data upon the geoid. The difference between 'aliased' and 'non-aliased' gravimetric geoid solutions varies by between 0.732 and −1.816 m (−0.058 m mean and 0.122 m standard deviation). Based on these conceptual arguments and numerical results, it is recommended that supplementary digital elevation information be included during the estimation of mean gravity anomalies prior to the computation of a gravimetric geoid model.  相似文献   

17.
Recent studies on the relationship between the Nusselt ( Nu ) and Rayleigh ( Ra ) numbers for base-heated convection in a spherical shell with a constant viscosity show that the power-law index is around 1/4, which is different from the value of 1/3 predicted by a simple boundary layer theory. We show that such a difference may be caused by the flow pattern due to the geometry. The flow pattern of the convection in a spherical shell at relatively low Ra , at least, less than 106, is characterized by narrow upwelling and broad downwelling, which is similar to the opposite flow pattern of internally heated convection. Convection in the internally heated case predicts the power-law index of 1/4. We demonstrate this relationship based on the concept of 'local' Rayleigh ( Ra1 ) and Nusselt ( Nu1 ) numbers  相似文献   

18.
We present a spectral-finite-element approach to the 2-D forward problem for electromagnetic induction in a spherical earth. It represents an alternative to a variety of numerical methods for 2-D global electromagnetic modelling introduced recently (e.g. the perturbation expansion approach, the finite difference scheme). It may be used to estimate the effect of a possible axisymmetric structure of electrical conductivity of the mantle on surface observations, or it may serve as a tool for testing methods and codes for 3-D global electromagnetic modelling. The ultimate goal of these electromagnetic studies is to learn about the Earth's 3-D electrical structure.
Since the spectral-finite-element approach comes from the variational formulation, we formulate the 2-D electromagnetic induction problem in a variational sense. The boundary data used in this formulation consist of the horizontal components of the total magnetic intensity measured on the Earth's surface. In this the variational approach differs from other methods, which usually use spherical harmonic coefficients of external magnetic sources as input data. We verify the assumptions of the Lax-Milgram theorem and show that the variational solution exists and is unique. The spectral-finite-element approach then means that the problem is parametrized by spherical harmonics in the angular direction, whereas finite elements span the radial direction. The solution is searched for by the Galerkin method, which leads to the solving of a system of linear algebraic equations. The method and code have been tested for Everett & Schultz's (1995) model of two eccentrically nested spheres, and good agreement has been obtained.  相似文献   

19.
Traveltime computation methods for strongly heterogeneous 3-D media developed during recent years are well suited for earthquake location. We present here a new method based on the traveltime algorithm of Podvin-Lecomte, related to the inverse problem formulation of Tarantola & Valette. The Podvin-Lecomte method, based on the Huygens principle, is very robust and allows arbitrary surface topography and station placement even for borehole instruments. First arrival traveltimes are computed for each of the recording stations using a fine 3-D velocity mesh (up to 106 cells on a workstation). The traveltime grid allows the use of the Tarantola & Valette formulation, which enables a full non-linear approach. The solution is given as a 3-D probability density function of hypocentre coordinates, which accounts for the arrival time measurements as well as a priori information for the location, the accuracy of both the arrival time readings and the computation of the theoretical traveltimes. This powerful method called 3DGRIDLOC gives the location of the induced seismicity of the gas field of Lacq (France) using 443 520 cells of a 3-D velocity mesh and the observations from nine recording stations, one of which is located at the bottom of a 3880 m deep borehole. Location of synthetic foci as well as more than 500 actual earthquakes shows the real advantages of this new method versus the classical HYPO71. A new insight into the induced seismicity is now possible: induced seismicity may occur as far away as 10 km from the gas reservoir and involve a much greater volume of rock than expected using earlier locations.  相似文献   

20.
We present an adaptive unstructured triangular grid finite element approach for effectively simulating plane-wave diffusive electromagnetic fields in 2-D conductivity structures.
The most striking advantage of irregular grids is their potential to incorporate arbitrary geometries including surface and seafloor topography. Adaptive mesh refinement strategies using an a posteriori error estimator yield most efficient numerical solutions since meshes are only refined where required.
We demonstrate the robustness of this approach by comparison with analytical solutions and previously published numerical simulations. Maximum errors may systematically be reduced to, for example, 0.8 per cent for the apparent resistivity and 0.2° in the phase.
An additional accuracy study of the thickness of the air layer in E-polarization suggests to keep a minimum thickness depending on lateral conductivity contrasts within the earth.
Furthermore, we point out the new quality and flexibility of our simulation technique by addressing two marine magnetotelluric applications. In the first case, we discuss topographic effects associated with a synthetic sinusoidal sea bottom model and in the second case, we show a close-to-reality scenario using real bathymetry data from the East Pacific Rise at 17°S.  相似文献   

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

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