On the long-term stability of terrestrial reference frame solutions based on Kalman filtering

The Global Geodetic Observing System requirement for the long-term stability of the International Terrestrial Reference Frame is 0.1 mm/year, motivated by rigorous sea level studies. Furthermore, high-quality station velocities are of great importance for the prediction of future station coordinates, which are fundamental for several geodetic applications. In this study, we investigate the performance of predictions from very long baseline interferometry (VLBI) terrestrial reference frames (TRFs) based on Kalman filtering. The predictions are computed by extrapolating the deterministic part of the coordinate model. As observational data, we used over 4000 VLBI sessions between 1980 and the middle of 2016. In order to study the predictions, we computed VLBI TRF solutions only from the data until the end of 2013. The period of 2014 until 2016.5 was used to validate the predictions of the TRF solutions against the measured VLBI station coordinates. To assess the quality, we computed average WRMS values from the coordinate differences as well as from estimated Helmert transformation parameters, in particular, the scale. We found that the results significantly depend on the level of process noise used in the filter. While larger values of process noise allow the TRF station coordinates to more closely follow the input data (decrease in WRMS of about 45%), the TRF predictions exhibit larger deviations from the VLBI station coordinates after 2014 (WRMS increase of about 15%). On the other hand, lower levels of process noise improve the predictions, making them more similar to those of solutions without process noise. Furthermore, our investigations show that additionally estimating annual signals in the coordinates does not significantly impact the results. Finally, we computed TRF solutions mimicking a potential real-time TRF and found significant improvements over the other investigated solutions, all of which rely on extrapolating the coordinate model for their predictions, with WRMS reductions of almost 50%.     

A multidisciplinary-based conceptual model of a fractured sedimentary bedrock aquitard: improved prediction of aquitard integrity

A hydrogeologic conceptual model that improves understanding of variability in aquitard integrity is presented for a fractured sedimentary bedrock unit in the Cambrian-Ordovician aquifer system of midcontinent North America. The model is derived from multiple studies on the siliciclastic St. Lawrence Formation and adjacent strata across a range of scales and geologic conditions. These studies employed multidisciplinary techniques including borehole flowmeter logging, high-resolution depth-discrete multilevel well monitoring, fracture stratigraphy, fluorescent dye tracing, and three-dimensional (3D) distribution of anthropogenic tracers regionally. The paper documents a bulk aquitard that is highly anisotropic because of poor connectivity of vertical fractures across matrix with low permeability, but with ubiquitous bed parallel partings. The partings provide high bulk horizontal hydraulic conductivity, analogous to aquifers in the system, while multiple preferential termination horizons of vertical fractures serve as discrete low vertical hydraulic conductivity intervals inhibiting vertical flow. The aquitard has substantial variability in its ability to protect underlying groundwater from contamination. Across widespread areas where the aquitard is deeply buried by younger bedrock, preferential termination horizons provide for high aquitard integrity (i.e. protection). Protection is diminished close to incised valleys where stress release and weathering has enhanced secondary pore development, including better connection of fractures across these horizons. These conditions, along with higher hydraulic head gradients in the same areas and more complex 3D flow where the aquitard is variably incised, allow for more substantial transport to deeper aquifers. The conceptual model likely applies to other fractured sedimentary bedrock aquitards within and outside of this region.     

Quantitative progradation dynamics and stratigraphic architecture of ancient shallow‐marine clinoform sets: a new method and its application to the Upper Jurassic Sognefjord Formation,Troll Field,offshore Norway

This article presents a new numerical inversion method to estimate progradation rates in ancient shallow‐marine clinoform sets, which is then used to refine the tectono‐stratigraphic and depositional model for the Upper Jurassic Sognefjord Formation reservoir in the super‐giant Troll Field, offshore Norway. The Sognefjord Formation is a 10–200‐m thick, coarse‐grained clastic wedge, that was deposited in ca. 6 Myr by a fully marine, westward‐prograding, subaqueous delta system sourced from the Norwegian mainland. The formation comprises four, 10–60‐m thick, westerly dipping, regressive clinoform sets, which are mapped for several tens of kilometres along strike. Near‐horizontal trajectories are observed in each clinoform set, and the sets are stacked vertically. Clinoform age and progradation rates are constrained by: (i) regionally correlatable bioevents, tied to seismically mapped clinoforms and clinoform set boundaries that intersect wells, (ii) exponential age–depth interpolations between bioevent‐dated surfaces and a distinctive foreset‐to‐bottomset facies transition within each well, and (iii) distances between wells along seismic transects that are oriented perpendicular to the clinoform strike and tied to well‐based stratigraphic correlations. Our results indicate a fall in progradation rate (from 170–500 to 10–65 km Myr^?1) and net sediment flux (from 6–14 to ≤1 km² Myr^?1) westwards towards the basin, which is synchronous with an overall rise in sediment accumulation rate (from 7–16 to 26–102 m Myr^?1). These variations are attributed to progradation of the subaqueous delta into progressively deeper waters, and a concomitant increase in the strength of alongshore currents that transported sediment out of the study area. Local spatial and temporal deviations from these overall trends are interpreted to reflect a subtle structural control on sedimentation. This method provides a tool to improve the predictive potential of sequence stratigraphic and clinoform trajectory analyses and offers a greater chronostratigraphic resolution than traditional approaches.     

Dark matter distribution in the Draco dwarf from velocity moments

If the observed relativistic plasma outflows in astrophysical jets are magnetically collimated and a single-component model is adopted, consisting of a wind-type outflow from a central object, then a problem arises with the inefficiency of magnetic self-collimation to collimate a sizeable portion of the mass and magnetic fluxes in the relativistic outflow from the central object. To solve this dilemma, we have applied the mechanism of magnetic collimation to a two-component model consisting of a relativistic wind-type outflow from a central source and a non-relativistic wind from a surrounding disc. By employing a numerical code for a direct numerical solution of the steady-state problem in the zone of super-fast magnetized flow, which allows us to perform a determination of the flow with shocks, it is shown that in this two-component model it is possible to collimate into cylindrical jets all the mass and magnetic fluxes that are available from the central source. In addition, it is shown that the collimation of the plasma in this system is usually accompanied by the formation of oblique shock fronts. The non-relativistic disc-wind not only plays the role of the jet collimator, but it also induces the formation of shocks as it collides with the initially radial inner relativistic wind and also as the outflow is reflected by the system axis. Another interesting feature of this process of magnetic collimation is a sequence of damped oscillations in the width of the jet.     

Numerical simulations of thermal convection in a rapidly rotating spherical shell cooled inhomogeneously from above

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell with and without inhomogeneous temperature anomalies on the top boundary have been carried out using a three-dimensional, time-dependent, spectral-transform code. The spherical shell of Boussinesq fluid has inner and outer radii the same as those of the Earth's liquid outer core. The Taylor number is 10⁷, the Prandtl number is 1, and the Rayleigh number R is 5R_c (R_c is the critical value of R for the onset of convection when the top boundary is isothermal and R is based on the spherically averaged temperature difference across the shell). The shell is heated from below and cooled from above; there is no internal heating. The lower boundary of the shell is isothermal and both boundaries are rigid and impermeable. Three cases are considered. In one, the upper boundary is isothermal while in the others, temperature anomalies with (l,m) = (3,2) and (6,4) are imposed on the top boundary. The spherically averaged temperature difference across the shell is the same in all three cases. The amplitudes of the imposed temperature anomalies are equal to one-half of the spherically averaged temperature difference across the shell. Convective structures are strongly controlled by both rotation and the imposed temperature anomalies suggesting that thermal inhomogeneities imposed by the mantle on the core have a significant influence on the motions inside the core. The imposed temperature anomaly locks the thermal perturbation structure in the outer part of the spherical shell onto the upper boundary and significantly modifies the velocity structure in the same region. However, the radial velocity structure in the outer part of the shell is different from the temperature perturbation structure. The influence of the imposed temperature anomaly decreases with depth in the shell. Thermal structure and velocity structure are similar and convective rolls are more columnar in the inner part of the shell where the effects of rotation are most dominant.     

Groundwater Monitoring on Guam: Management Responses to Recent Water Quality Violations

Guam has one of the finest limestone aquifers in the world. Located in the northern half of the island, this vital underground resource supplies the local community with about 80% of its drinking water needs. The majority of the island’s approximately 180,000 inhabitants live in the northern half of the island where significant economic growth and urban development have occurred over the last 30 years. The U.S. military has also occupied large tracts of land in this region for the past 60 years. The risks of groundwater contamination are, therefore, very real considering the population density in northern Guam and the rapid recharge rates to the underlying aquifer. Since April 1996, Guam Waterworks Authority (GWA) has regularly monitored the island’s drinking water resources for chemical and biological contaminants in accordance with the U.S. Safe Drinking Water Act requirements. This report summarizes and compares the groundwater data sets for two discrete time intervals: 1996 through 2001 and 2002 through 2007 for chemicals, and 1997(8) through 2002 and 2003 through 2007 for bacteria. Adaptive management strategies that evolved to deal with contaminant violations over this time frame are critically discussed. Overall, 95% of all GWA’s drinking water production wells were considered to be in relatively good shape from a chemical standpoint. Identified chemicals of concern were chlordane, ethylene dibromide, perchloroethylene, trichloroethylene, and nitrate. The microbiological integrity of Guam’s groundwater was, by comparison, less impressive prior to 2003 because of numerous wastewater spills, leaks, and overflows. Improved management practices introduced shortly afterward have significantly reduced the incidence of Escherichia coli contamination in recent years.