Application of Gravity Gradients in the Process of GOCE Orbit Determination

The possibility of improving the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) mission satellite orbit using gravity gradient observations was investigated. The orbit improvement is performed by a dedicated software package, called the Orbital Computation System (OCS), which is based on the classical least squares method. The corrections to the initial satellite state vector components are estimated in an iterative process, using dynamic models describing gravitational perturbations. An important component implemented in the OCS package is the 8th order Cowell numerical integration procedure, which directly generates the satellite orbit. Taking into account the real and simulated GOCE gravity gradients, different variants of the solution of the orbit improvement process were obtained. The improved orbits were compared to the GOCE reference orbits (Precise Science Orbits for the GOCE satellite provided by the European Space Agency) using the root mean squares (RMS) of the differences between the satellite positions in these orbits. The comparison between the improved orbits and the reference orbits was performed with respect to the inertial reference frame (IRF) at J2000.0 epoch. The RMS values for the solutions based on the real gravity gradient measurements are at a level of hundreds of kilometers and more. This means that orbit improvement using the real gravity gradients is ineffective. However, all solutions using simulated gravity gradients have RMS values below the threshold determined by the RMS values for the computed orbits (without the improvement). The most promising results were achieved when short orbital arcs with lengths up to tens of minutes were improved. For these short arcs, the RMS values reach the level of centimeters, which is close to the accuracy of the Precise Science Orbit for the GOCE satellite. Additional research has provided requirements for efficient orbit improvement in terms of the accuracy and spectral content of the measured gravity gradients.     

Determination of the Croatian geomagnetic observatory location

Ground survey within the Nature Park Lonjsko Polje, placed in the middle-northern Croatia was performed during the time interval 2007–2010 in order to find the best location for installing the geomagnetic observatory. The total magnetic field has been measured a few times using the Overhauser proton magnetometers. The horizontal and vertical gradients of the total field, and its temporal behaviour were investigated over the restricted region that we estimated as suitable for the observatory. The results obtained from thoroughly conducted measurements allowed us to find definitive positions for the instrument pillars. These results are in agreement with previously suggested location found based on combination of Comprehensive CM4 model prediction and measurements conducted from 2003 to 2005. This study contributes to the development of geomagnetism in Croatia and paves a way to install the first geomagnetic observatory in Croatia.     

Dynamic response analysis of solitary flexible rocking bodies: modeling and behavior under pulse‐like ground excitation

Results obtained for rigid structures suggest that rocking can be used as seismic response modification strategy. However, actual structures are not rigid: structural elements where rocking is expected to occur are often slender and flexible. Modeling of the rocking motion and impact of flexible bodies is a challenging task. A non‐linear elastic viscously damped zero‐length spring rocking model, directly usable in conventional finite element software, is presented in this paper. The flexible rocking body is modeled using a conventional beam‐column element with distributed masses. This model is verified by comparing its pulse excitation response to the corresponding analytical solution and validated by overturning analysis of rocking blocks subjected to a recorded ground motion excitation. The rigid rocking block model provides a good approximation of the seismic response of solitary flexible columns designed to uplift when excited by pulse‐like ground motions. Guidance for development of rocking column models in ordinary finite element software is provided. Copyright © 2014 John Wiley & Sons, Ltd.     

Deficit and surplus of precipitation as a continuous function of time

The quantitative measure of drought is a complex problem. A common factor in all types of drought is that they originate from a deficiency of precipitation. The difference between the real and potentially expected daily precipitation amount is what we have used to define the continuous function of daily deficit and surplus of precipitation (DSP). We have also defined the monthly and yearly DSP. We tested the precipitation data series for the normal climatic period 1961–1990 observed at the Meteorological Observatory of Belgrade. The DSP considers that the intensity of drought during the current period is dependent upon the current precipitation plus the cumulative precipitation of the previous period. When applied to longer periods, this function denotes dry and wet periods.     

Absolute Sea Level Surface Modeling for the Mediterranean from Satellite Altimeter and Tide Gauge Measurements

Changes in the height of the ocean can be described through the relative and absolute sea level changes depending on the geodetic reference the sea level records are related to. Satellite altimetry provides absolute sea level (ASL) measurements related to the global geodetic reference, whereas tide gauges provide relative sea level (RSL) measurements related to the adjacent land. This study aims at computing the ASL surfaces for different time epochs from combined satellite altimeter and tide gauge records. A method of sea level data fusion is proposed to enable modeling of the impact of present and future sea level changes on the coast. Sea surface modeling was investigated for ten different gridding methods commonly used for the interpolation of altimeter data over the open ocean and extrapolation over the coastal zones. The performance of gridding methods was assessed based on the comparison of the gridded altimeter data and corrected tide gauge measurements. Finally, the sea level surfaces related to the GRS80 global reference ellipsoid were computed for the Mediterranean Sea over the altimeter period. In addition, the current sea level trends were estimated from both sea level measurements.     

Ophiolitic detritus in Cretaceous clastic formations of the Dinarides (NW Croatia): evidence from Cr-spinel chemistry

Cr-spinel is a common heavy mineral in the sandstones of Cretaceous synorogenic sedimentary formations of the NW Dinarides, Croatia. The rocks occur in isolated exposures in the uplifted basement units of Medvednica, Ivanščica, Žumberak and Samobor Mountains near Zagreb. In this area, evidence of the early Alpine evolution of the Dinarides is obscured due to strong dismemberment of pre-Tertiary tectonostratigraphic units resulting from an intense tectonic history, as well as due to the widespread sedimentary cover of the Pannonian Basin. Electron microprobe analyses of detrital Cr-spinels from the Oštrc Formation reveal that in the Early Cretaceous the ophiolitic source area was predominantly composed of harzburgite peridotites and associated cumulate rocks, which developed in a supra-subduction zone setting. The supply of Cr-spinels with the same chemical signature remained dominant until the end of the Cretaceous, suggesting that exposed remnants of the same ophiolite belt persisted through the Cretaceous and/or that recycling was significant. Similarities with data reported from the Northern Calcareous Alps and the Transdanubian Central Range imply that a rather extensive harzburgitic ophiolite belt probably extended along the Adriatic margin during the Early Cretaceous. A slight trend of increasing variation in the Cr# is observed from the Early to the latest Cretaceous, suggesting that the source areas became more heterogeneous with the ongoing Cretaceous tectonic evolution. Differences in Cr-spinel compositions in two contemporaneous latest Cretaceous formations are well in line with existing data on heavy mineral proportions, which together identify contrasting hinterland geology for these formations and strongly suggest the coeval existence of two separate basins. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development of Productivity Models for the Northern Gulf of Mexico Based on Oxygen Concentrations and Stable Isotopes

Oxygen concentrations have been used for decades to estimate primary production (P) and respiration (R) in aquatic ecosystems. Yet, this approach cannot separate the effects of biological and physical processes affecting oxygen dynamics; therefore, it is now often complemented with the analysis of stable oxygen isotopes. Existing algorithms for calculating primary production and P/R have been developed for closed systems and steady-state open systems. None of these formulations are applicable to productive aquatic ecosystems where diurnal changes in oxygen concentrations and isotope values are usually large. Here, we describe a novel P/R model that includes algorithms for air–sea gas exchange and is not constrained by steady-state conditions. Our objective was to test model sensitivity to variations in input parameters for values commonly observed in coastal waters of the northern Gulf of Mexico. The model was highly sensitive to variations in fractionation factor for respiration (ε) but less sensitive to variations in wind speed, oxygen isotope values of source-water, or oxygen flux through the picnocline. This model is easily transferable to other coastal ecosystems, with a caveat that system-specific values for ε are needed to obtain realistic estimates of P/R.