High-precision geoid determination in small areas: A case study in Doñana National Park (Spain)

Doñana National Park is an area of approximately 500 km² located on the SW coast of Spain that shows one of the greatest geoid gradients on the entire Iberian Peninsula, due to its peculiar tectonic characteristics. So, it is necessary to elaborate an accurate geoid model that can be used with GPS for precise surveying, since the existing ones are insufficient, due to their poor resolution and their limited adaptation to a small area with such a strong gradient. The least squares prediction method was tested in order to obtain the undulation from GPS/orthometric points. The results obtained were unsatisfactory because of the strong geoid gradient. In order to improve accuracy the remove-restore technique was used. Global geopotential model EIGEN-CG01C and a Digital Elevation Model (DEM) with a 25 × 25 m resolution and an accuracy better than 3 m were used. Thus, the final geometrical geoid obtained reaches the precision required by other disciplines (3 cm in any point within the Park). Particularly, the geoid model has allowed for the acquisition of a precision DEM that is essential to formulate a hydrodynamic model for the Doñana marsh functions.     

Integrated chronostratigraphy of Proterozoic-Cambrian boundary beds in the western Anabar region, northern Siberia

Carbonate-rich sedimentary rocks of the western Anabar region, northern Siberia, preserve an exceptional record of evolutionary and biogeochemical events near the Proterozoic/Cambrian boundary. Sedimentologically, the boundary succession can be divided into three sequences representing successive episodes of late transgressive to early highstand deposition; four parasequences are recognized in the sequence corresponding lithostratigraphically to the Manykal Formation. Small shelly fossils are abundant and include many taxa that also occur in standard sections of southeastern Siberia. Despite this coincidence of faunal elements, biostratigraphic correlations between the two regions have been controversial because numerous species that first appear at or immediately above the basal Tommotian boundary in southeastern sections have first appearances scattered through more than thirty metres of section in the western Anabar. Carbon- and Sr-isotopic data on petrographically and geochemically screened samples collected at one- to two-metre intervals in a section along the Kotuikan River, favour correlation of the Staraya Reckha Formation and most of the overlying Manykai Formation with sub-Tommotian carbonates in southeastern Siberia. In contrast, isotopic data suggest that the uppermost Manykai Formation and the basal 26 m of the unconformably overlying Medvezhya Formation may have no equivalent in the southeast; they appear to provide a sedimentary and palaeontological record of an evolutionarily significant time interval represented in southeastern Siberia only by the sub-Tommotian unconformity. Correlations with radiometrically dated horizons in the Olenek and Kharaulakh regions of northern Siberia suggest that this interval lasted approximately three to six million years, during which essentially all 'basal Tommotian' small shelly fossils evolved.     

Glyptograptus? persculptus (Salter), its tectonic deformation,and its stratigraphic significance for the Carys Mills Formation of N.E. Maine,U.S.A.

The Carys Mills Formation of NE Maine represents the SW extension of the Matapedia Group and the White Head Formation of southern Gaspé, Quebec. This formation has yielded two main collections, one hitherto considered of Middle Ordovician age (zone 13 of Berry's 1960 Marathon, Texas, zonation) with 7 species referred to 4 genera, and one of the Early Silurian age. A revision of the supposed Middle Ordovician collections reveals that the graptolites are basal Silurian and consist of only one species, Glyptograptus? persculptus. Several hundred specimens, preserved both flattened and in three dimensions, the latter in iron pyrites, allow an analysis of both the rhabdosome changes wrought by diagenetic flattening and those caused by severe tectonic deformation. It is shown that those specimens lying parallel to the bedding plane lineation caused by a high-angle main cleavage have a silhouette appearance close to Glyptograptus tamariscus (Nicholson), whilst those at right angles to this direction superficially resemble robust Amplexograptus spp., or Diplograptus s.s., or Diplograptus s.l. spp., thus explaining the stratigraphic assignment of earlier workers. An analogue video reshaper confirms the above conclusions. This determination assigns the whole Carys Mills Formation to the Lower Silurian, but the possibility that part of the formation is of Late Ordovician age cannot be excluded. This much is suggested by the Matapedia Group and the White Head Formation which are of Late Ordovician to Early Silurian age, and the Honorat Group, which has yielded graptolites of late Middle Ordovician age.     

Evaluating Groundwater Interbasin Flow Using Multiple Models and Multiple Types of Data

The groundwater interbasin flow, Q_y, from the north of Yucca Flat into Yucca Flat simulated using the Death Valley Regional Flow System (DVRFS) model greatly exceeds assessments obtained using other approaches. This study aimed to understand the reasons for the overestimation and to examine whether the Q_y estimate can be reduced. The two problems were tackled from the angle of model uncertainty by considering six models revised from the DVRFS model with different recharge components and hydrogeological frameworks. The two problems were also tackled from the angle of parametric uncertainty for each model by first conducting Morris sensitivity analysis to identify important parameters and then conducting Monte Carlo simulations for the important parameters. The uncertainty analysis is general and suitable for tackling similar problems; the Morris sensitivity analysis has been utilized to date in only a limited number of regional groundwater modeling. The simulated Q_y values were evaluated by using three kinds of calibration data (i.e., hydraulic head observations, discharge estimates, and constant‐head boundary flow estimates). The evaluation results indicate that, within the current DVRFS modeling framework, the Q_y estimate can only be reduced to about half of the original estimate without severely deteriorating the goodness‐of‐fit to the calibration data. The evaluation results also indicate that it is necessary to develop a new hydrogeological framework to produce new flow patterns in the DVRFS model. The issues of hydrogeology and boundary flow are being addressed in a new version of the DVRFS model planned for release by the U.S. Geological Survey.     

Coherent GPS reflections from the sea surface

In this letter, we explore a method to obtain accurate ocean heights using measurements of the global positioning system (GPS) carrier phase after reflection on the sea surface. A carrier tracking algorithm is employed in measuring travel path differences between GPS direct and reflected signals collected from antennas suspended over a marine estuary, when roughness guarantees partially coherent reflections at the GPS L1 frequency. This technique proves to be sensitive to surface roughness and able to follow tide variations with a precision better than 5 cm (1-sigma) for sea states with significant wave heights below 10 cm. It is expected that this technique could be further extended to rougher sea states using GPS frequency combinations with longer synthetic wavelengths.     

Convection in the polar cap after the interplanetary magnetic field becomes northward

As shown by Iwasaki (1971); Maezawa (1976); Kuznetsov and Troshichev (1977) and other investigators, the electric field and the plasma convection in the polar cap change their direction after an appearance of a significant northward component of the interplanetary magnetic field. Two possible mechanisms of this phenomenon may be suggested: (i) the direct penetration of the dusk-to-dawn electric field from the solar wind into the magnetosphere, and (ii) the generation of the observed electric field and convection in a process of the decay of the three-dimensional current system which existed before the appearance of the northward interplanetary field. The latter mechanism implies that the value of the electric field generated in the polar cap will decrease with time after the appearance of the northward interplanetary magnetic field. The results of the experimental investigation show such a decrease which favours the second mechanism.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.