A simulation study of the errors of omission and commission for GRACE RL01 gravity fields

This paper examines the influence that certain omission and commission errors can have on the gravity field models estimated from the initial release of data (RL01) from the Gravity Recovery And Recovery Experiment (GRACE) satellite mission. The effects of omission errors were analyzed by limiting the degree and order to which the GPS and K-band range-rate (KBR) measurement partials were extended in the solution process. The commission error studies focused on the impact of an imperfect mean reference gravity field model on the solution. Combinations of both of these error sources were also explored. The nature of these errors makes them difficult to distinguish from the true gravity signal, so the exploration of these error sources was performed using simulations; however, comparisons to real-data solutions are provided. The results show how each of the specific error sources investigated influences the gravity field solution. The simulations also show how all of the errors examined can be sufficiently mitigated through the appropriate choice of processing parameters.     

Jason-1 Precision Orbit Determination by Combining SLR and DORIS with GPS Tracking Data

With the implementation of the Jason-1 satellite altimeter mission, the goal of reaching the 1-cm level in orbit accuracy was set. To support the Precision Orbit Determination (POD) requirements, the Jason-1 spacecraft carries receivers for DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite) and GPS (Global Positioning System), as well as a retroreflector for SLR (Satellite Laser Ranging). The overall orbit accuracy for Jason will depend on the quality and the relative weighting of the available tracking data. In this study, the relative importance of the SLR, DORIS, and GPS tracking data is assessed along with the most effective parameterization for accounting for the unmodeled accelerations through the application of empirical accelerations. The optimal relative weighting for each type of tracking data was examined. It is demonstrated that GPS tracking alone is capable of supporting a radial orbit accuracy for Jason-1 at the 1-cm level, and that including SLR tracking provides additional benefits. It is also shown that the GRACE (Gravity Recovery and Climate Experiment) gravity model GGM01S provides a significant improvement in the orbit accuracy and reduction in the level of geographically correlated orbit errors.     

Air–Sea/Land Interaction in the Coastal Zone

Atmospheric turbulence measurements made at the U.S. Army Corps of Engineers Field Research Facility (FRF) located on the Atlantic coast near the town of Duck, North Carolina during the CASPER-East Program (October–November 2015) are used to study air–sea/land coupling in the FRF coastal zone. Turbulence and mean meteorological data were collected at multiple levels (up to four) on three towers deployed at different landward distances from the shoreline, with a fourth tower located at the end of a 560-m-long FRF pier. The data enable comparison of turbulent fluxes and other statistics, as well as investigations of surface-layer scaling for different footprints, including relatively smooth sea-surface conditions and aerodynamically rough dry inland areas. Both stable and unstable stratifications were observed. The drag coefficient and diurnal variation of the sensible heat flux are found to be indicators for disparate surface footprints. The drag coefficient over the land footprint is significantly greater, by as much as an order of magnitude, compared with that over the smooth sea-surface footprint. For onshore flow, the internal boundary layer in the coastal zone was either stable or (mostly) unstable, and varied dramatically at the land-surface discontinuity. The offshore flow of generally warm air over the cooler sea surface produced a stable internal boundary layer over the ocean surface downstream from the coast. While the coastal inhomogeneities violate the assumptions underlying Monin–Obukhov similarity theory (MOST), any deviations from MOST are less profound for the scaled standard deviations and the dissipation rate over both water and land, as well as for stable and unstable conditions. Observations, however, show a poor correspondence with MOST for the flux-profile relationships. Suitably-averaged, non-dimensional profiles of wind speed and temperature vary significantly among the different flux towers and observation levels, with high data scatter. Overall, the statistical dependence of the vertical gradients of scaled wind speed and temperature on the Monin–Obukhov stability parameter in the coastal area is weak, if not non-existent.     

The role of the Polochic Fault as part of the North American and Caribbean Plate boundary: Insights from the infill of the Lake Izabal Basin

The Lake Izabal Basin in Guatemala is a major pull-apart basin along the sinistral Polochic Fault, which is part of the North American and Caribbean plate boundary. The basin infill contains information about the tectonic and sedimentological processes that have imparted a significant control on its sedimentary section. The inception of the basin has been linked to the relative importance of the Polochic Fault in the tectonic history of the plate boundary; yet, its sedimentological record and its inception age have been poorly documented. This study integrates diverse datasets, including industry reports, well logs and reports, well cuttings, vintage seismic data, outcrop observations and geochronological data to constrain the initial infill and age of inception of the basin. The integrated data show that during the Oligocene–Miocene, a marine carbonate platform was established in the region which was later uplifted and eroded in the early Miocene. The fluvial–lacustrine deposits above this carbonate platform are part of the initial infill of the basin and are constrained with zircon weighted-mean ²⁰⁶Pb/²³⁸U ages of 12.060 ± 0.008 from a volcanic tuff ~30 m above the unconformity. Sandstone, mudstone and coal dominate the interval from 12 to 4 Ma, with an increase in conglomerate correlating to the uplift of the Mico Mountains and San Gil Hill at 4 Ma. Fault switch activity between the Polochic and Motagua faults has been hypothesized to explain total offset along the Polochic Fault and the geologic and geodetic slip rates along the two faults. The 12 Ma age determined for the initial infill of the basin confirms this hypothesis. Consequently, our study confirms that at ~12 Ma the Polochic Fault served as the main fault of the plate boundary with inferred slip rates ranging from 13 to 21 mm/yr with a strong possibility that the Polochic Fault was, at some point between 15 Ma and 7 Ma, the only active fault of the plate boundary. The results of this study show that tectonic records preserved in sediments of strike-slip basins improve the understanding of the relative significance of individual faults and the implications with respect to strain partitioning throughout its tectonic history.     

Exploring the Boundaries of Web Map Services: The Example of the Online Injury Atlas for Ontario

Injury from causes such as falls, traffic accidents, or violence is a major public health issue globally. Injury prevention research aims to identify vulnerable populations and places by analyzing the spatial patterns of demographic and socio‐economic risk factors associated with elevated injury rates. The stakeholders in injury prevention and control are often distributed across government and public health institutions, non‐profits, and even the private sector (e.g. insurance firms). While this situation calls for distributed, online research tools, their implementation may conflict with health data confidentiality and license limitations for socio‐economic data. In this article, we present the Online Injury Atlas for Ontario, which was designed with the explicit goal of making use of, and contributing to, the Canadian Geospatial Data Infrastructure. We propose a service‐based architecture that integrates publicly accessible map services with protected data layers. Thereby, we demonstrate the benefits of using spatial data infrastructures alongside private data at different levels of protection. In addition, we discuss the extensive data processing needs and specific cartographic design requirements of a Web atlas in the health and social sciences domain.     

Magmatic-vapor expansion and the formation of high-sulfidation gold deposits: Chemical controls on alteration and mineralization

Large bulk-tonnage high-sulfidation gold deposits, such as Yanacocha, Peru, are the surface expression of structurally-controlled lode gold deposits, such as El Indio, Chile. Both formed in active andesite–dacite volcanic terranes. Fluid inclusion, stable isotope and geologic data show that lode deposits formed within 1500 m of the paleo-surface as a consequence of the expansion of low-salinity, low-density magmatic vapor with very limited, if any, groundwater mixing. They are characterized by an initial ‘Sulfate’ Stage of advanced argillic wallrock alteration ± alunite commonly with intense silicification followed by a ‘Sulfide’ Stage — a succession of discrete sulfide–sulfosalt veins that may be ore grade in gold and silver. Fluid inclusions in quartz formed during wallrock alteration have homogenization temperatures between 100 and over 500 °C and preserve a record of a vapor-rich environment.Recent data for El Indio and similar deposits show that at the commencement of the Sulfide Stage, ‘condensation’ of Cu–As–S sulfosalt melts with trace concentrations of Sb, Te, Bi, Ag and Au occurred at > 600 °C following pyrite deposition. Euhedral quartz crystals were simultaneously deposited from the vapor phase during crystallization of the vapor-saturated melt occurs to Fe-tennantite with progressive non-equilibrium fractionation of heavy metals between melt-vapor and solid. Vugs containing a range of sulfides, sulfosalts and gold record the changing composition of the vapor.Published fluid inclusion and mineralogical data are reviewed in the context of geological relationships to establish boundary conditions through which to trace the expansion of magmatic vapor from source to surface and consequent alteration and mineralization. Initially heat loss from the vapor is high resulting in the formation of acid condensate permeating through the wallrock. This Sulfate Stage alteration effectively isolates the expansion of magmatic vapor in subsurface fracture arrays from any external contemporary hydrothermal activity. Subsequent fracturing is localized by the embrittled wallrock to provide high-permeability fracture arrays that constrain vapor expansion with minimization of heat loss. The Sulfide Stage vein sequence is then a consequence of destabilization of metal-vapor species in response to depressurization and decrease in vapor density.The geology, mineralogy, fluid inclusion and stable isotope data and geothermometry for high-sulfidation, bulk-tonnage and lode deposits are quite different from those for epithermal gold–silver deposits such as McLaughlin, California that formed near-surface in groundwater-dominated hydrothermal systems where magmatic fluid has been diluted to less than about 30%. High sulfidation gold deposits are better termed ‘Solfataric Gold Deposits’ to emphasize this distinction. The magmatic-vapor expansion hypothesis also applies to the phenomenology of acidic geothermal systems in active volcanic systems and equivalent magmatic-vapor discharges on the flanks of submarine volcanoes.     

A multi-criteria evaluation model of earthquake vulnerability in Victoria,British Columbia

Researchers have recently examined the geographic variability of the vulnerability of populations to earthquakes. These studies focus mainly on the complex modelling of geophysical processes or identification of socio-economically disadvantaged populations. However, no studies to date have integrated different components of vulnerability with metrics of travel distance to hospitals and trauma centres (systemic vulnerability). We argue that this previously unaccounted component is an important conceptual and practical aspect of earthquake vulnerability. Accordingly, this paper presents a multi-criteria model for combining physical, social, and systemic components, moving towards a more comprehensive assessment of vulnerability. An analytic hierarchy process is used to produce a place-specific index of social vulnerability, which we combine with soil liquefaction and amplification index and a road network model for access to hospitals and trauma services. Using a geographic information system, we implemented this model for the Greater Victoria region (483 km², 2011 population: 345,000) in British Columbia, Canada. Clustering of total vulnerability was found in outlying areas, highlighting the importance of access to trauma centres. We conclude by identifying challenges in measuring earthquake vulnerability and advocate integration of systemic vulnerability components in natural hazards research.     

River network travel time is correlated with dissolved organic matter composition in rivers of the contiguous United States

Most terrestrial allochthonous organic matter enters river networks through headwater streams during high flow events. In headwaters, allochthonous inputs are substantial and variable, but become less important in streams and rivers with larger watersheds. As allochthonous dissolved organic matter (DOM) moves downstream, the proportion of less aromatic organic matter with autochthonous characteristics increases. How environmental factors converge to control this transformation of DOM at a continental scale is less certain. We hypothesized that the amount of time water has spent travelling through surface waters of inland systems (streams, rivers, lakes, and reservoirs) is correlated to DOM composition. To test this hypothesis, we used established river network scaling relationships to predict relative river network flow-weighted travel time (FWTT) of water for 60 stream and river sites across the contiguous United States (3090 discrete samples over 10 water years). We estimated lentic contribution to travel times with upstream in-network lake and reservoir volume. DOM composition was quantified using ultraviolet and visible absorption and fluorescence spectroscopy. A combination of FWTT and lake and reservoir volume was the best overall predictor of DOM composition among models that also incorporated discharge, specific discharge, watershed area, and upstream channel length. DOM spectral slope ratio (R² = 0.77) and Freshness Index (R² = 0.78) increased and specific ultraviolet absorbance at 254 nm (R² = 0.68) and Humification Index (R² = 0.44) decreased across sites as a function of FWTT and upstream lake volume. This indicates autochthonous-like DOM becomes continually more dominant in waters with greater FWTT. We assert that river FWTT can be used as a metric of the continuum of DOM composition from headwaters to rivers. The nature of the changes to DOM composition detected suggest this continuum is driven by a combination of photo-oxidation, biological processes, hydrologically varying terrestrial subsidies, and aged groundwater inputs.     

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence Duffing's equation, the van der Pol equation and the oscillator with quadratic damping problem are solved using each method.