首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Summary. Most of the Earth's magnetic field and its secular change originate in the core. Provided the mantle can be treated as an electrical insulator, stochastic inversion enables surface observations to be analysed for the core field. A priori information about the variation of the field at the core boundary leads to very stringent conditions at the Earth's surface. The field models are identical with those derived from the method of harmonic splines (Shure, Parker & Backus) provided the a priori information is specified appropriately.
The method is applied to secular variation data from 106 magnetic observatories. Model predictions for fields at the Earth's surface have error estimates associated with them that appear realistic. For plausible choices of a priori information the error of the field at the core is unbounded, but integrals over patches of the core surface can have finite errors. The hypothesis that magnetic fields are frozen to the core fluid implies that certain integrals of the secular variation vanish. This idea is tested by computing the integrals and their standard and maximum errors. Most of the integrals are within one standard deviation of zero, but those over the large patches to the north and south of the magnetic equator are many times their standard error, because of the dominating influence of the decaying dipole. All integrals are well within their maximum error, indicating that it will be possible to construct core fields, consistent with frozen flux, that satisfy the observations.  相似文献   

2.
Summary. Laplace and Bessel Transforms are used to solve for the transient behaviour of the electromagnetic fields after switching off a steady current in a grounded infinitesimal horizontal dipole on the surface of a uniformly conducting half-space. Simple analytic expressions, which are valid for times sufficiently long after the switch that displacement terms can be ignored, are obtained on the surface of the half-space for the electric field and the time derivative of the magnetic field. At the instant of switching an infinitesimally long image becomes established directly under the source dipole. It is the diffusion of this image which gives the vertical magnetic field and horizontal electric fields their transient behaviour. During the transient, there is also a decaying charge distribution on the surface.  相似文献   

3.
To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory (MATLAB) software, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by comparing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatially random fields.  相似文献   

4.
Summary. A conducting slab of finite thickness divided into three segments of different conductivities and overlying a perfect conductor is proposed as a suitable two-dimensional 'control' model for testing the accuracy of the various numerical modelling programs that are available for calculating the fields induced in the Earth by an external, time-varying magnetic source. An analytic solution is obtained for this control model for the case of the magnetic field everywhere parallel to the conductivity boundaries ( B -polarization). Values of the field given by this solution for a particular set of model parameters are calculated at selected points on the surface and on a horizontal plane inside the conductor, and are tabulated to three figure accuracy for reference. They are used to check the accuracy of the results given by the finite difference program of Brewitt-Taylor & Weaver and the finite element program of Kisak & Silvester for the same model. Improved formulae for calculating the derived electric field components in B -polarization are first developed for incorporation in the finite difference program, and these give surface electric fields within 1 per cent of the analytic values, while all three field components inside the conductor are calculated to better than 96 per cent accuracy by the finite difference program. The results given by the finite element program are not quite so satisfactory. Errors somewhat greater than 10 per cent are present and although the program requires much less disk space it takes rather more CPU time to complete the calculations.  相似文献   

5.
Surface integral formulae are derived expressing any one of certain field quantities, namely current functions, magnetic potentials and normal components of magnetic fields, in terms of any one other, for current systems flowing in concentric spherical surfaces. In all, 36 such formulae are obtained, which should prove useful in many geomagnetic studies, especially in geomagnetic induction problems.  相似文献   

6.
Summary. A problem in modelling electromagnetic fields used in exploration geophysics is treated mathematically. Analytical expressions are obtained for the electric field due to a harmonic current in a horizontal loop on or above a conducting ground in which is buried a conductive and permeable sphere (ore body). The loop is coaxial with the sphere. For a general time-varying current in the loop, the analysis is carried to the stage where a Fourier inversion can be used to obtain readily the electric field in the time-domain. A new relationship between spherical and cylindrical wave functions is obtained as a transformation of local elements.
Solution of this problem has not been presented before in this form. Lee's solution of 1975 which uses an integral-equation formulation treats a similar problem without taking account of differences in magnetic permeability. The effects of magnetic permeability may have important and useful implications for geophysical explorations.  相似文献   

7.
We consider the kinematic production of magnetic fields in a sphere by velocity fields dominated by differential rotation and spiralling convective cells. The high magnetic Reynolds number limit of Braginsky (1964) is considered and formulae are derived allowing an α-effect parametrization of such flows to be easily calculated. This permits an axisymmetric system to be investigated in parallel with the direct 3-D numerical computations. Good agreement between the asymptotic and 3-D calculations is found. The 'spiralling' property typical of convective motion in rotating spheres is important in terms of dynamo action; the differential rotation coexisting with this feature is also vital. Indeed, it is the presence of both features which allows the analysis of Braginsky to be employed. With flows approximating the columnar form anticipated for rapidly rotating convection, dynamo action is relatively easily achieved for all azimuthal wavenumbers; modes of differing wavenumbers interact almost by a simple superposition. With flows of more complex latitudinal form, the mutual interactions between modes become more complicated. For columnar-type flows, dipole magnetic fields are favoured when the sense of outward spiralling is prograde and the zonal flow is eastwards, as is physically preferred.  相似文献   

8.
Many geophysical inverse problems derive from governing partial differential equations with unknown coefficients. Alternatively, inverse problems often arise from integral equations associated with a Green's function solution to a governing differential equation. In their discrete form such equations reduce to systems of polynomial equations, known as algebraic equations. Using techniques from computational algebra one can address questions of the existence of solutions to such equations as well as the uniqueness of the solutions. The techniques are enumerative and exhaustive, requiring a finite number of computer operations. For example, calculating a bound to the total number of solutions reduces to computing the dimension of a linear vector space. The solution set itself may be constructed through the solution of an eigenvalue problem. The techniques are applied to a set of synthetic magnetotelluric values generated by conductivity variations within a layer. We find that the estimation of the conductivity and the electric field in the subsurface, based upon single-frequency magnetotelluric field values, is equivalent to a linear inverse problem. The techniques are also illustrated by an application to a magnetotelluric data set gathered at Battle Mountain, Nevada. Surface observations of the electric ( E y ) and magnetic ( H x ) fields are used to construct a model of subsurface electrical structure. Using techniques for algebraic equations it is shown that solutions exist, and that the set of solutions is finite. The total number of solutions is bounded above at 134 217 728. A numerical solution of the algebraic equations generates a conductivity structure in accordance with the current geological model for the area.  相似文献   

9.
i
Oriented cores have been secured from fourteen sites in the Great Rhodesian Dyke, by means of a portable sampling drill. The natural remanent magnetizations showed high dispersion at all sites except one. After demagnetization in alternating magnetic fields, nine sites gave well grouped directions of primary magnetization. These sites include five rock types distributed among three Complexes of the Great Dyke and two satellite dykes, over 200 miles of the length of the Dyke and through several thousand feet in depth as the rocks were originally intruded. The nine site mean directions of primary magnetization are closely grouped and are believed to represent directions of thermo-remanent magnetization at the date of intrusion of the Great Dyke. It is suggested that the dates of magnetization at the sites must cover a sufficient time interval to give a mean pole position close to the axial geocentric dipole freed from secular variation. On the assumption of a geocentric dipole field, the position of the mean South magnetic pole is 211/2 °N, 611/2 °E, with radius of 95 per cent confidence 9°. This pole position is close to positions of North magnetic poles given by studies of the palaeo-magnetism of the Pilansberg Dykes and Bushveld gabbro.  相似文献   

10.
Summary. The method of stochastic inversion, previously applied to secular variation data, is applied to main field data. Adaptations to the method are required: non-linear, as well as linear, data are used; allowance is made for crustal components in the observatory data; and the prior information is specified differently. The requirement that the models should satisfy a finite lower bound on the Ohmic heating in the core provides strong prior information and gives finite error estimates at the core—mantle boundary.
The new method is applied to data from the epochs 1969.5 and 1980.0. The resulting field models are very much more complex than other models, such as the IGRF models extrapolated to the core, and show considerable small-scale detail which, on the basis of the error analysis, can be believed.
The flux integral over the northern hemisphere is computed at each epoch; the difference between the two epochs is approximately one standard deviation, suggesting that the question as to whether the decay of the dipole is consistent with the frozen-flux hypothesis has been resolved in favour of the hypothesis.  相似文献   

11.
A physical model for palaeosecular variation   总被引：1，自引：0，他引：1  
Summary. A new model to describe the latitude dependence of the angular dispersion of the palaeomagnetic field (palaeosecular variation) is developed following previous models, but with crucial differences. It is shown that if the probability distribution of virtual geomagnetic poles (VGPs) is circularly symmétric about the rotation axis then the geométry of the distribution of field directions is latitude dependent. This has a significant effect on the latitude dependence of dispersion and is accounted for in the model. The dipole and non-dipole parts of the field are not artificially separated but are intimately linked through an observationally determined relation that the time averaged intensity of the non-dipole field is dependent upon the intensity of the dipole field. It is shown that a consequence of this relation is that no knowledge of the probability distribution of the geomagnetic dipole moment is required. This is a fundamental improvement over previous models.
The model provides excellent fits to the palaeodata and, unlike previous models, is not inconsistent with the latitude variation of the non-dipole field dispersion determined from the present field. For the past 5 Ma the point estimate of the VGP dispersion due to dipole wobble is 7.2° and of the VGP dispersion at the equator due to variation in the non-dipole field is 10.6°. This estimate of the dispersion due to variation in the non-dipole field is in excellent agreement with the value predicted from an analysis of the variation in field intensities over the same period. Fits of the model to data from earlier periods indicate that dispersion due to variation in the non-dipole field is essentially independent of the geomagnetic reversal rate while dipole wobble is positively correlated with reversal rate.  相似文献   

12.
13 lava flows of known age (ages from 14C dating), which have been erupted in the last 30 000 years, have been studied to determine the palaeosecular variation of the geomagnetic field in Central Mexico. Samples were taken from two different monogenetic volcanic fields: the Michoacan-Guanajuato volcanic field (six sites) and the Chichinautzin Formation (seven sites), both part of the Transmexican Volcanic Belt. The lavas were studied in detail using rock magnetic methods (magnetic susceptibility at room temperature, low-temperature susceptibility behaviour, hysteresis loops, Curie temperatures), combined with reflected light microscopy, in order to deduce their magnetic mineralogy and the domain states of the magnetic minerals. The magnetic carriers are titanomagnetites, which show differing degrees of high-temperature deuteric oxidation, and seem to be predominantly pseudo-single domain (PSD), though in many cases are probably a mixture of domain states. Mean palaeomagnetic directions and palaeointensity values using Shaw and Thellier techniques were obtained using several specimens from each flow. Our data seem to indicate a sharp easterly swing in declination about 5000 years ago, which is also observed in lake sediments from Central Mexico. The calculated values of the virtual dipole moment (VDM) range from 3.1 to 14.9 × 1022 A m2. Our data indicate that the virtual dipole moment seems to have increased gradually in magnitude over the last 30 kyr, with a peak at about 9000 years BP. These are features that have been observed in other parts of the globe and are probably caused by variations in the dipole part of the geomagnetic field.  相似文献   

13.
Summary New formulae are obtained for the displacement potentials and displacements due to a point source with moment tensor, and for a fault with prescribed slip. These formulae, unlike previous formulae, are invariant, i.e. they are valid in any coordinate system, not just Cartesian coordinates or orthogonal curvilinear coordinates.
Apart from their invariance, these formulae have other advantages: they are exceedingly simple; the expressions for P -motion are nearly the same as those for S -motion; and the separation of far field motion from final static displacement is automatic.  相似文献   

14.
Geographic representation has become more complex through time as researchers have added new concepts, leading to apparently endless proliferation and creating a need for simplification. We show that many of these concepts can be derived from a single foundation that we term the atomic form of geographic information. The familiar concepts of continuous fields and discrete objects can be derived under suitable rules applied to the properties and values of the atomic form. Fields and objects are further integrated through the concept of phase space, and in the form of field objects. A second atomic concept is introduced, termed the geo‐dipole, and shown to provide a foundation for object fields, metamaps, and the association classes of object‐oriented data modelling. Geographic dynamics are synthesized in a three‐dimensional space defined by static or dynamic object shape, the possibility of movement, and the possibility of dynamic internal structure. The atomic form also provides a tentative argument that discrete objects and continuous fields are the only possible bases for geographic representation.  相似文献   

15.
Summary. Simple and convenient formulae are derived which describe the transformation of a multipolar expansion under an arbitrary proper rotation of the reference frame. When combined with the corresponding formulae for a translation, these results show how multipolar representations of source fields transform under any proper displacement of the reference frame. Particular emphasis is placed on the seismic source problem; however, these results find applications in many other physical problems.  相似文献   

16.
We developed an inversion method to estimate the stress fields related to earthquake generation (seismogenic stress fields) from the centroid moment tensors (CMT) of seismic events by using Akaike's Bayesian information criterion (ABIC). On the idea that the occurrence of an earthquake releases some part of the seismogenic stress field around its hypocentre, we define the CMT of a seismic event by a weighted volume integral of the true but unknown seismogenic stress field. Representing each component of the seismogenic stress field by the superposition of a finite number of 3-D basis functions (tri-cubic B-splines), we obtain a set of linear observation equations to be solved for the expansion coefficients (model parameters). We introduce prior constraint on the roughness of the seismogenic stress field and combine it with observed data to construct a Bayesian model with hierarchic, highly flexible structure controlled by hyper-parameters. The optimum values of the hyper-parameters are objectively determined form observed data by using ABIC. Given the optimum values of the hyper-parameters, we can obtain the best estimates of model parameters by using a maximum likelihood algorithm. We tested the validity of the inversion method through numerical experiments on two synthetic CMT data sets, assuming the distribution of fault orientations to be aligned with the maximum shear stress plane in one case and to be random in the other case. Then we applied the inversion method to actual CMT data in northeast Japan, and obtained the pattern of the seismogenic stress field consistent with geophysical and geological observations.  相似文献   

17.
Palaeomagnetic results are presented from a volcanic sequence in Zhijin County, Guizhou Province, Southern China. The lavas and associated volcanic breccias comprising the sequence represent a southern extension of the Emeishan volcanic province. Biostratigraphic dating of interbedded limestone units and stratigraphic constraints indicate that the section formed during the late Permian (∼263–255  Ma), and is thus somewhat older than the Emei stratotype section in Sichuan Province, and close in age to reported estimates of the termination of the Permo-Carboniferous (Kiaman) reverse superchron. Rock magnetic analyses and reflected light microscopy indicate that the magnetic mineralogy of the lava units is dominated by fresh, primary magnetites containing a significant fraction of single-domain grains. Thermal demagnetization behaviour of the breccia units is poor, but most lava samples have one or two components of remanence above 250 °C. The normal polarity characteristic remanence held by the lavas implies a post-Kiaman age for this succession and suggests that the termination of the Kiaman occurred prior to 263  Ma, supporting recently published estimates. According to standard criteria, Thellier palaeointensity results from the lavas are of good quality and reveal that the dipole field strength was comparatively low shortly after the termination of the superchron. 80 per cent of samples record relative VDM values in the range 42–52 per cent of the present-day value, supporting recent studies of mid-Kiaman field intensity. This suggests that a low-energy dipole existed at least between 300 and 255  Ma and does not appear to have been confined to the stable reverse polarity interval.  相似文献   

18.
A sandbox model consisting of two adjacent mechanically different decollements (frictional and viscous) loosely simulated the southeastern part of the Zagros fold-thrust-belt. Digital images of the model surface are used to coordinate passive markers on the surface and quantify displacement fields and estimate 2-D finite strains. These analyses show that, mapped in a fixed coordinate system, the deformation front propagates at different rates above the two decollements. Strain analysis of the model surface at different stages of deformation also shows that cumulative strain is more heterogeneous above the viscous decollement where strain domains are separated by fault zones. Maps of displacement fields, finite strain ellipses and dilatation also differ in character above the two decollements. Displacements above a viscous decollement decrease gradually towards the foreland, whereas they decrease sharply in front of the frontal thrust above the frictional decollement. Our analyses also show that the estimated finite strain depends not only on the density of the marker points chosen for the analysis, but also their initial distribution relative to the structures. This comparison shows that marker density limits measuring the actual strains in a heterogeneously deforming fold-thrust-belt and marker density and distribution have a strong impact on the strain analysed in nature. The similarity of our model with nature is examined with recent GPS study in the Zagros fold-thrust-belt (SW Iran) and shows, similar to the model results, that a weak salt decollement causes divergent movement in the sedimentary cover in SE Zagros.  相似文献   

19.
A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g 12  and   h 12  , suggestive of long-term non-zonal structure in the secular variations.  相似文献   

20.
  相似文献   

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

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