The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model

We apply a coupled thermomechanical ice sheet—self-gravitating viscoelastic solid Earth model (SGVEM), allowing for the dynamic exchange of ice thickness and bedrock deformation, in order to investigate the effect of viscoelastic deformation on ice dynamics and vice versa. In a synthetic glaciation scenario, we investigate the interaction between the ice sheet and the solid Earth deformation, the glacial-isostatic adjustment (GIA), accounting for an atmospheric forcing depending on the ice sheet surface altitude. We compare the results from the coupled model to runs with the common elastic lithosphere/relaxing asthenosphere (ELRA) model, where the lithosphere is represented by a thin plate and the mantle relaxes with one characteristic relaxation time, as well as to a rigid Earth without any deformation. We find that the deformational behaviour of the SGVEM on ice dynamics (i.e. stored ice volume, ice thickness and velocity field) is comparable to the ELRA for an optimal choice of the parameters in steady state, but exhibits differences in the transient behaviour. Beyond the ice sheet, in the region of peripheral forebulge, the differences in the transient surface deformation between ELRA and SGVEM are substantial, demonstrating the inadequacy of the ELRA model for interpreting constraints on GIA in the periphery of the ice sheet, such as sea-level indicators and GPS uplift rates.     

Manifestations of the North – South Asymmetry in the Photosphere and in the Green Line Corona

The north – south asymmetries (NSA) of three solar activity indices are derived and mutually compared over a period of more than five solar cycles (1945 – 2001). A catalogue of the hemispheric sunspot numbers, the data set of the coronal green line brightness developed by us, and the magnetic flux derived from the NSO/KP data (1975 – 2001) are treated separately within the discrete low- and mid-latitude zones (5° – 30°, 35° – 60°). The calculated autocorrelations, cross-correlations, and regressions between the long-term NSA data sets reveal regularities in the solar activity phenomenon. Namely, the appearance of a distinct quasi-biennial oscillation (QBO) is evident in all selected activity indices. Nevertheless, a smooth behavior of QBO is derived only when sufficient temporal averaging is performed over solar cycles. The variation in the significance and periodicity of QBO allows us to conclude that the QBO is not persistent over the whole solar cycle. A similarity in the photospheric and coronal manifestations of the NSA implies that their mutual relation will also show the QBO. A roughly two-year periodicity is actually obtained, but again only after significant averaging over solar cycles. The derived cross-correlations are in fact variable in degree of correlation as well as in changing periodicity. A clear and significant temporal shift of 1 – 2 months in the coronal manifestation of the magnetic flux asymmetry relative to the photospheric manifestation is revealed as a main property of their mutual correlation. This shift can be explained by the delayed large-scale coronal manifestation in responding to the emergence of the magnetic flux in the photosphere. The reliability of the derived results was confirmed by numerical tests performed by selecting different numerical values of the used parameters.     

Tumbling asteroids

We present both a review of earlier data and new results on non-principal axis rotators (tumblers) among asteroids. Among new tumblers found, the best data we have are for 2002 TD₆₀, 2000 WL₁₀₇, and (54789) 2001 MZ₇—each of them shows a lightcurve with two frequencies (full terms with linear combinations of the two frequencies are present in the lightcurve). For 2002 TD₆₀, we have constructed a physical model of the NPA rotation. Other recent objects which have been found to be likely tumblers based on their lightcurves that do not fit with a single periodicity are 2002 NY₄₀, (16067) 1999 RH₂₇, and (5645) 1990 SP. We have done a statistical analysis of the present sample of the population of NPA rotators. It appears that most asteroids larger than ∼0.4 km with estimated damping timescales (Harris, 1994, Icarus 107, 209) of 4.5 byr and longer are NPA rotators. The statistic of two short-period tumblers (D=0.04 and 0.4 km) with non-zero tensile strength suggests that for them the quantity μQ/T, where μ is the mechanical rigidity, Q is the elastic dissipation factor, and T is a spin excitation age (i.e., a time elapsed since the last significant spin excitation event), is greater by two to four orders of magnitude than the larger, likely rubble-pile tumblers. Among observational conditions and selection effects affecting detections of NPA rotations, there is a bias against detection of low-amplitude (small elongation) tumblers.     

Review: The Yucat��n Peninsula karst aquifer, Mexico

The Yucatán Peninsula karst aquifer is one of the most extensive and spectacular karst aquifer systems on the planet. This transboundary aquifer system extends over an area of approximately 165,000? km² in México, Guatemala and Belize. The Triassic to Holocene Yucatán limestone platform is located in the vicinity of the North American/Caribbean plate boundary and has been reshaped by a series of tectonic events over its long geologic history. At the end of the Cretaceous period, the Yucatán Peninsula was hit by a large asteroid, which formed the Chicxulub impact crater. The Yucatán Peninsula karst aquifer hosts large amounts of groundwater resources which maintain highly diverse groundwater-dependent ecosystems. Large parts of the aquifer are affected by seawater intrusion. Anthropogenic pollution of the aquifer has been increasing over the past few decades, owing to relentless economic development and population growth on the Peninsula. This review summarizes the state of knowledge on the Yucatán Peninsula karst aquifer and outlines the main challenges for hydrologic research and practical groundwater-resources management on the Peninsula.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

������С��������Ԥ��ȫ�������ܵ��Ӻ���

??????????????????(VTEC)??????????????????????????????????4?????????????????????Э????????????????С???????÷????????????????????????г????????????????????????????????????????????????????????(CODE)????????????????????????0.812???????????ο???в???4 TECu????????????????????????????????????????????γ20°?????????????????????????????????????1??7 TECu?????????????????2.4 TECu?????CODE????????20.14%????????????????3.5 TECu?????CODE???????????????????3.0 TECu?????CODE????????8.25%??     

Comparison of two different solutions to Molodensky’s G1 term

The main purpose of this paper is comparison of two different approaches of solution to the Simple Molodensky’s Problem, the Molodensky’s Approach and the Analytical Continuation Approach, based on numerical computation. Although these approaches have been described theoretically by several authors, e.g. Molodensky et al. (1960), Heiskanen and Moritz (1967), Vaníček (1974), Moritz (1980) and Holota (1991, unpublished results) and theoretical proof of equivalence was given by Heiskanen and Moritz (1967), Moritz (1971), Ecker (1971) and Pellinen (1972, unpublished results), only very few practical experiences about the differences between particular solutions and computational efficiency exist. In this paper we compare the above two mentioned approaches in terms of the G₁-effect on quasigeoid. Both quasigeoid solutions are tested by the independent GPS/levelling approach and are also compared with the previous quasigeoid model of Slovakia where the G₁-term has been approximated using the classical terrain correction. The effect of the G₂-term is also numerically estimated, revealing that it might be significant for the precise quasigeoid determination. Some practical comments and recommendations are given at the end of the paper.     

Stability of the solar system

If the solar system is considered as a mechanical clockwork consisting of its present members which attract each other as mass-points, the extent of its present approach to secular stability (i.e., the state of minimum potential energy) — manifested by the existence of a number of nearcommensurabilities of the present orbital periods, not only of the planets, but also of their satellites —could not have been attained in a time-span of 4.6×10⁹ yr of its age by gravitational perturbations alone.The existence of such commensurabilities — striking in many instances— could then be understood only on the assumption that either (a) the solar system was actually born with the present 2-, 3- and 4-term couplings between the orbital period of the planets already built-in from the outset (which is improbable on any known grounds); or (b) that these couplings — in particular, the 25 Jupiter-Saturn commensurability — have arisen as a result of tidal interaction between proto-planetary globes of much larger dimensions than these planets possess today. For the present dimensions and mutual distances of these planets, their tidal interaction in 10⁹ yr would exert but negligible effects; and during that time neither their masses, nor the scale of the solar system underwent any essential change.Therefore, a hypothesis is proposed that the situation now obtaining had its origin in the early days of the formation of the solar system, when the planetary globes — in particular, those of Jupiter and Saturn (now in the terminal stage of Kelvin contraction) — were very much larger than they are today; and when, as a result, the tidal coupling between them operated at a much higher rate than at the present time.Paper presented at the European Workshop on Planetary Sciences, organised by the Laboratorio di Astrofisica Spaziale di Frascati, and held between April 23–27, 1979, at the Accademia Nazionale del Lincei in Rome, Italy.     

Morphology of the Wadati-Benioff zone,andesitic volcanism,and active fracture zones in Central America

Summary The morphology of the Wadati-Benioff zone in the region of Central America, based on the distribution of 1377 earthquake foci, verified the existence of an intermediate aseismic gap and its relation to active andesitic volcanism, and the non-uniformity of subduction due to the hampering effect of the main structural features of the subducting Cocos plate. Four deep seismically active fracture zones, genetically connected with the process of subduction, and three fracture zones manifesting the possible boundary between the Americas and Caribbean plates were identified in the continental wedge.     

Mercury monitoring: Mercury stability in bird feathers

The influence of ultraviolet light, heating, freezing and weathering on the mercury concentration in the primary feathers from Guillemot and Black Guillemot has been examined. Even within 8 months of exposure variation in mercury concentration due to either loss of mercury or weight loss of the feathers has been found to be less than 10% relative.