Differentiation for High-Precision GPS Velocity and Acceleration Determination

Accurate estimates of the velocity and acceleration of a platform are often needed in high dynamic positioning, airborne gravimetry, and geophysics. In turn, differentiation of GPS signals is a crucial process for obtaining these estimates. It is important in the measurement domain where, for example, the phase measurements are used along with their instantaneous derivative (Doppler) to estimate position and velocity. It is also important in postprocessing, where acceleration is usually estimated by differentiating estimates of position and velocity. Various methods of differentiating a signal can have very different effects on the resulting derivative, and their suitability varies from situation to situation. These comments set the stage for the investigations in this article. The objective is twofold: (1) to carry out a comprehensive study of possible differentiation methods, characterizing each in the frequency domain; and (2) to use real data to demonstrate each of these methods in both of the measurement and position domains, in conditions of variable, high, or unknown dynamics. Examples are given using real GPS data in both the measurement domain and in the position and velocity domain. The appropriate differentiator is used in several cases of varying dynamics to derive a Doppler signal from carrier phase measurements (rather than using the raw Doppler generated by the receiver). In the statistic case, it is seen that the accuracy of velocity estimates can be improved from 4.0 mm/s to 0.7 mm/s by using the correct filter. In conditions of medium dynamics experienced in an airborne gravity survey, it is demonstrated that accelerations as the 2–4 mGal level (1 mGal = 0.00001 m/s²) can be obtained at the required filtering periods. Finally, a precision motion table is used to show that when using the correct filter, velocity estimates under high dynamics can be improved by an order of magnitude to 27.0 mm/s. ? 1999 John Wiley & Sons, Inc.     

Ocean–atmosphere coupling and the boreal winter MJO

The influence of ocean–atmosphere coupling on the simulation and prediction of the boreal winter Madden–Julian Oscillation (MJO) is examined using the Seoul National University coupled general circulation model (CGCM) and atmospheric—only model (AGCM). The AGCM is forced with daily SSTs interpolated from pentad mean CGCM SSTs. Forecast skill is examined using serial extended simulations spanning 26 different winter seasons with 30-day forecasts commencing every 5 days providing a total of 598 30-day simulations. By comparing both sets of experiments, which share the same atmospheric components, the influence of coupled ocean–atmosphere processes on the simulation and prediction of MJO can be studied. The mean MJO intensity possesses more realistic amplitude in the CGCM than in AGCM. In general, the ocean–atmosphere coupling acts to improve the simulation of the spatio-temporal evolution of the eastward propagating MJO and the phase relationship between convection (OLR) and SST over the equatorial Indian Ocean and the western Pacific. Both the CGCM and observations exhibit a near-quadrature relationship between OLR and SST, with the former lagging by about two pentads. However, the AGCM shows a less realistic phase relationship. As the initial conditions are the same in both models, the additional forcing by SST anomalies in the CGCM extends the prediction skill beyond that of the AGCM. To test the applicability of the CGCM to real-time prediction, we compute the Real-time Multivariate MJO (RMM) index and compared it with the index computed from observations. RMM1 (RMM2) falls away rapidly to 0.5 after 17–18 (15–16) days in the AGCM and 18–19 (16–17) days in the CGCM. The prediction skill is phase dependent in both the CGCM and AGCM.     

A combinatorial data model for representing topological relations among 3D geographical features in micro‐spatial environments

This research is motivated by the need for 3D GIS data models that allow for 3D spatial query, analysis and visualization of the subunits and internal network structure of ‘micro‐spatial environments’ (the 3D spatial structure within buildings). It explores a new way of representing the topological relationships among 3D geographical features such as buildings and their internal partitions or subunits. The 3D topological data model is called the combinatorial data model (CDM). It is a logical data model that simplifies and abstracts the complex topological relationships among 3D features through a hierarchical network structure called the node‐relation structure (NRS). This logical network structure is abstracted by using the property of Poincaré duality. It is modelled and presented in the paper using graph‐theoretic formalisms. The model was implemented with real data for evaluating its effectiveness for performing 3D spatial queries and visualization.     

Occurrence and characteristics of sporadic sodium layer observed by lidar at a mid-latitude location

Na lidar observations of SSL during the past 5 years at a mid-latitude location (Wuhan, China, 31°N, 114°E) are reported in this paper. From 26 SSL events detected in about 230 h of observation, an SSL occurrence rate of 1 event/9 h at our location was obtained. This result, combined with that reported by Nagasawa and Abo (Geophys. Res. Lett. 22 (1995) 263) at Tokyo, Japan, reveals that the SSL occurrence can be relatively frequent at some mid-latitude locations. The statistical analyses of main parameters for the 26 SSL events were made, and the results were compared with those previously reported. By examining the corresponding data from a nearby ionosonde, it was found that there was a fairly good correlation between SSL and Es. Of the 18 pairs of SSL and Es events checked, 15 of SSL were accompanied by Es, and 8 pairs of them were correlated within 1 h in time and within 5 km in altitude. From the analyses of observed perturbations during SSL development, the role of dynamic processes of atmosphere in the SSL formations were emphasized.     

Petrogenesis of Peraluminous Granites, Monashee Mountains, Southeastern Canadian Cordillera

Two-mica granites that locally contain garnet and sillimaniteoccur as dikes, sills, and sheets up to 50 m thick within thesillimanite zonc of the Monashee Mountains in the southeasternCanadian Cordillera of British Columbia. Syn-kinematic and post-kinematicgranites are recognized. U-Pb dating of zircon demonstrates that the syn-kinematic granitesare 100.4?0.3 Ma old, based on duplicate concordant single abradedzircon analyses. Other zircons have slightly older Pb/Pb dates,indicating small amounts of inherited zircons. Monazites are99?10 Ma old. Post-kinematic granites have 62.5?0.2 Ma zirconages and 634+0.1 Ma monazite ages. High initial 87 ratios (0.71492–0.74181)and evidence of Precambrian Pb inheritance indicates that bothsyn- and post-kinematic granites were derived from a crustalsource. Geobarometric estimates suggest that both generationsof granites equilibrated at 6–8 kb (22–30 km). Zirconand monazite saturation temperatures range from 660–824?Cand indicate that these minerals were liquidus phases earlyin the crystallization history of the granites. Because monazitesaturation temperatures generally exceed those of zircon, itis possible that some monazites may be inherited. Apatite saturationtemperatures in excess of 900?C suggest that both generationsof granites contain source inherited apatite. Syn- and post-kinematic granites have essentially identicalmajor and trace element chemistries. Syn-kinematic graniteshave steep light rare earth element (LREE) enriched patternswith pronounced negative Eu anomalies. The REE patterns of post-kinematicgranites range from steep LREE enriched patterns with negativeNd and Eu anomalies to flat patterns with low LREE contents,negative Nd anomalies, and both positive and negative Eu anomalies.Modelling of REE, Rb, Sr, and Ba contents demonstrates thatsyn-kinematic gramtes could have been generated by a low degreeof partial melting (with 10–25% feldspar fractionationof the melt) of Late Proterozoic Horsethief Creek Group metapelitesleaving a monazite-bearing upper amphibolite facies residue.Post-kinematic granites were produced by partial melting ofa geochemically and isotopically similar metapelitic source.The suite of post-kinematic granites can be related by a smallamount (up to 0.1%) of monazite crystal fractionation.     

Temperature dependence of the cation distribution in nickel aluminate (NiAl2O4) spinel: a powder XRD study

NiAl₂O₄ is a largely inverse spinel, which in detail shows increasing randomisation with temperature of Ni and Al between the octahedral and tetrahedral cation sites of the spinel structure. We have used powder XRD to determine this cation distribution in various samples of NiAl₂O₄ quenched after annealing between 700 and 1400° C. The inversion parameter (x) can be measured with a precision of ± 0.004 (one standard deviation), and a comparison of different methods of synthesis, X-ray diffraction and refinement techniques, suggests a probable accuracy of better than 0.01. The results are supported by some preliminary single crystal refinements on flux-grown samples.Below 800° C the rate of cation ordering becomes very slow, and, despite reaching an apparently steady state, it is doubtful if our samples attained complete internal equilibrium. Above 1250° C the cation redistribution becomes so fast that the quenching method becomes unreliable. Between 800 and 1250° C inclusive, the degree of inversion changes smoothly from 0.87 at 800° C to 0.79 at 1250° C, and is accompanied by linear changes in u, the oxygen parameter, from 0.2555 to 0.2563 (±0.0002), and a₀, the lattice parameter, from 8.0462 to 8.0522 Å (±0.0002 Å).     

Chromian spinel as a petrogenetic indicator in abyssal and alpine-type peridotites and spatially associated lavas

The composition of chromian spinel in alpine-type peridotites has a large reciprocal range of Cr and Al, with increasing Cr# (Cr/(Cr+Al)) reflecting increasing degrees of partial melting in the mantle. Using spinel compositions, alpine-type peridotites can be divided into three groups. Type I peridotites and associated volcanic rocks contain spinels with Cr#<0.60; Type III peridotites and associated volcanics contain spinels with Cr#>0.60, and Type II peridotites and volcanics are a transitional group and contain spinels spanning the full range of spinel compositions in Type I and Type II peridotites. Spinels in abyssal peridotites lie entirely within the Type I spinel field, making ophiolites with Type I alpine-type peridotites the most likely candidates for sections of ocean lithosphere formed at a midocean ridge. The only modern analogs for Type III peridotites and associated volcanic rocks are found in arc-related volcanic and intrusive rocks, continental intrusive assemblages, and oceanic plateau basalts. We infer a sub-volcanic arc petrogenesis for most Type III alpine-type peridotites. Type II alpine-type peridotites apparently reflect composite origins, such as the formation of an island-arc on ocean crust, resulting in large variations in the degree and provenance of melting over relatively short distances. The essential difference between Type I and Type III peridotites appears to be the presence or absence of diopside in the residue at the end of melting.Based on an examination of co-existing rock and spinel compositions in lavas, it appears that spinel is a sensitive indicator of melt composition and pressure of crystallization. The close similarity of spinel composition fields in genetically related basalts, dunites and peridotites at localities in the oceans and in ophiolite complexes indicates that its composition reflects the degree of melting in the mantle source region. Accordingly, we infer from the restricted range of spinel compositions in abyssal basalts that the degree of mantle melting beneath mid-ocean ridges is generally limited to that found in Type I alpine-type peridotites. It is apparent, therefore, that the phase boundary OL-EN-DI-SP +meltOL-EN-SP+melt has limited the degree of melting of the mantle beneath mid-ocean ridges. This was clearly not the case for many alpine-type peridotites, implying very different melting conditions in the mantle, probably involving the presence of water.     

Global travertine deposition modulated by oscillations in climate

Travertine deposits are important records of past fluid flow in the Earth's crust, and document fluid migration through both tectonic activity and changes in climate. While many studies hint at possible relationships between travertine formation and global climate, none have investigated these connections on a global scale. Here we compile 1649 published travertine ages from six continents to test the hypothesis that global and/or regional changes in climate regulate travertine deposition. Peaks in bedded travertine ages occur with main frequencies that correspond to 100‐kyr changes in global climate, where most peaks occur during glacial terminations or interglacial periods, including a large peak that coincides with the Early Holocene climatic optimum. Time–series analysis also suggests a possible connection with 41‐kyr obliquity cycles. At regional scales, many peaks also correspond with local times of high precipitation or wet conditions. This can be attributed to higher groundwater recharge rates, providing the necessary water to form travertine. Many bedded travertine‐depositing systems may therefore be water‐limiting and sufficient CO₂ may be present even during times of no travertine deposition. Exceptions to this conclusion are banded vein travertine deposits, which typically form during times of dry climate when water tables are low. Copyright © 2019 John Wiley & Sons, Ltd.     

Upper estimates for the element number of non-redundant antenna configurations on square and hexagonal grids

Estimates for the maximum number of elements of non-redundant configurations on integer square and hexagonal grids of given sizes are derived (a “non-redundant” configuration implies that the vector differences between its elements are all distinct). When projecting a large 2-D interferometer or a telescope, such an estimate can be used as a guide for evaluating the maximum possible number of antennas in a non-redundant configuration that can be arranged within a given area. The suggested estimates are empirical and based on the available data. They are obtained by reducing the problem to the linear case and by generalizing the method applied therein. Examples of applying the method are presented.     

Sustainable shellfish aquaculture in Saldanha Bay,South Africa

The carrying capacity for bivalve shellfish culture in Saldanha Bay, South Africa, was analysed through the application of the well-tested EcoWin ecological model, in order to simulate key ecosystem variables. The model was set up using: (i) oceanographic and water-quality data collected from Saldanha Bay, and (ii) culture-practice information provided by local shellfish farmers. EcoWin successfully reproduced key ecological processes, simulating an annual mean phytoplankton biomass of 7.5 µg Chl a l^–1 and an annual harvested shellfish biomass of about 3 000 tonnes (t) y^–1, in good agreement with reported yield. The maximum annual carrying capacity of Small Bay was estimated as 20 000 t live weight (LW) of oysters Crassostrea gigas, or alternatively 5 100 t LW of mussels Mytilus galloprovincialis, and for Big Bay as 100 000 t LW of oysters. Two production scenarios were investigated for Small Bay: a production of 4 000 t LW y^–1 of mussels, and the most profitable scenario for oysters of 19 700 t LW y^–1. The main conclusions of this work are: (i) in 2015–2016, both Small Bay and Big Bay were below their maximum production capacity; (ii) the current production of shellfish potentially removes 85% of the human nitrogen inputs; (iii) a maximum-production scenario in both Big Bay and Small Bay would result in phytoplankton depletion in the farmed area; (iv) increasing the production intensity in Big Bay would probably impact the existing cultures in Small Bay; and (v) the production in Small Bay could be increased, resulting in higher income for farmers.