Modeling preferential flow in reactive barriers: implications for performance and design

Reactive barriers are passive and in situ ground water treatment systems. Heterogeneities in hydraulic conductivity (K) within the aquifer-reactive barrier system will result in higher flux rates, and reduced residence times, through portions of the barrier. These spatial variations in residence time will affect the treatment capacity of the barrier. A numerical flow model was used to evaluate the effects of spatial variations in K on preferential flow through barriers. The simulations indicate that the impact of heterogeneities in K will be a function of their location and distribution; the more localized the high K zone, the greater the preferential flow. The geometry of the reactive barrier will also strongly influence flow distribution. Aquifer heterogeneities will produce greater preferential flow in thinner barriers compared to thicker barriers. If the barrier K is heterogeneous, greater preferential flow will occur in thicker barriers. The K of the barrier will affect the flow distribution; decreasing the K of the barrier can result in more even distribution of flow. Results indicate that less variable flow will be attained utilizing thicker, homogeneous barriers. The addition of homogeneous zones to thinner barriers will be effective at redistributing flow only if installed immediately adjacent to both the up- and downgradient faces of the barrier.     

Sedimentation on the Ionian active margin (Hellenic arc) — Provenance of sediments and mechanisms of deposition

Lithological, sedimentological, mineralogical and faunal analyses of twenty cores from the western Peloponnesus margin and adjacent Matapan Trench permit the origin of the sediments and the processes of sedimentation to be defined in some detail.Except for the ash layers derived from the Ischian Province, the major part of the sediment is local in origin (Peloponnesus). Analyses of heavy, light and clay minerals enable three main provinces to be distinguished. Each of these provinces includes one or more slope basins and a portion of trench, without lateral connection between them and mixing of inputs.The diversity and complexity of sedimentary structures, the granulometric characters, the mixed fauna (molluscs and ostracods) and the variability of rates of sedimentation all testify to the prevalence of reworking phenomena and gravitative processes.These data confirm and elaborate on the model of sedimentation deduced from seismic studies in which sediments are successively trapped, released by the slope basins and transferred via channels and canyons towards the underlying trench basins. This “cascade feeding” appears characteristics of this subductive sedimentary environment.     

The ionospheric impact of the October 2003 storm event on Wide Area Augmentation System

The United States Federal Aviation Administrations (FAA) Wide-Area Augmentation System (WAAS) for civil aircraft navigation is focused primarily on the Conterminous United States (CONUS). Other Satellite-Based Augmentation Systems (SBAS) include the European Geostationary Navigation Overlay Service (EGNOS) and the Japanese Multi-transport Satellite-based Augmentation System (MSAS). Navigation using WAAS requires accurate calibration of ionospheric delays. To provide delay corrections for single frequency global positioning system (GPS) users, the wide-area differential GPS systems depend upon accurate determination of ionospheric total electron content (TEC) along radio links. Dual-frequency transmissions from GPS satellites have been used for many years to measure and map ionospheric TEC on regional and global scales. The October 2003 solar-terrestrial events are significant not only for their dramatic scale, but also for their unique phasing of solar irradiance and geomagnetic events. During 28 October, the solar X-ray and EUV irradiances were exceptionally high while the geomagnetic activity was relatively normal. Conversely, 29–31 October was geomagnetically active while solar irradiances were relatively low. These events had the most severe impact in recent history on the CONUS region and therefore had a significant effect on the WAAS performance. To help better understand the event and its impact on WAAS, we examine in detail the WAAS reference site (WRS) data consisting of triple redundant dual-frequency GPS receivers at 25 different locations within the US. To provide ground-truth, we take advantage of the three co-located GPS receivers at each WAAS reference site. To generate ground-truth and calibrate GPS receiver and transmitter inter-frequency biases, we process the GPS data using the Global Ionospheric Mapping (GIM) software developed at the Jet Propulsion Laboratory. This software allows us to compute calibrated high resolution observations of TEC. We found ionospheric range delays up to 35 m for the day-time CONUS during quiet conditions and up to 100 m during storm time conditions. For a quiet day, we obtained WAAS planar fit slant residuals less than 2 m (0.4 m root mean square (RMS)) and less than 25 m (3.4 m RMS) for the storm day. We also investigated ionospheric gradients, averaged over distances of a few hundred kilometers. The gradients were no larger than 0.5 m over 100 km for a quiet day. For the storm day, we found gradients at the 4 m level over 100 km. Similar level gradients are typically observed in the low-latitude region for quiet or storm conditions.     

THE NEW DECADE: DIGITS AND DIVERSITY

The President of the Photogrammetric Society, Mr. I Newton, chaired a panel discussion concerned with those developments in photogrammetry which may take place during the 1990s.     

Wet tropospheric effects on precise relative GPS height determination

Summary Considerable interest has been generated recently in the use of the Global Positioning System (GPS) for precise height determination. A major error source in these measurements is the propagation delay due to atmospheric water vapour. In order to achieve the high precisions required for such applications as absolute sea-level monitoring improvement of wet delay modelling is necessary. Results from a GPS campaign show a significant correlation (0.91) between the variability of the wet delay measured using a water vapour radiometer (WVR) at the Onsala site and the absolute value of the residual error in the height determination of a 134 km baseline from Onsala to Jönköping. This correlation indicates that the atmosphericvariability as inferred from the WVR data includes information on the quality of the GPS height estimate. During periods of high atmospheric activity, e.g., during the passage of a weather front, the use of a six-parameter gradient model reduces the spread for the vertical coordinate from 40 mm to 20 mm (with standard deviations of 17 mm and 9 mm respectively) over the 134 km baseline (less than 1 × 10^–7) using 8 hour data spans on 11 different days over a six month period.     

Spectrum estimation using maximum entropy and multiresolution considerations

Summary Given a sample autocovariance sequence of finite length for some observed random process, the spectrum estimation problem involves the extension of this sequence for the required Fourier transformation. The maximum entropy approach which is based on the optimal use of information contents, leads to a dual sequence of reflection coefficients with reciprocal spectrum of the process. The estimation of the maximum entropy spectrum implies results identical to those using autoregressive modeling in one dimension under appropriate white noise assumptions. In cases of a non-white noise component, the approach is generalized to an autoregressive-moving-average model. Recent developments in multiresolution analysis with spectral domain decompositions also offer possibilities of subband spectrum estimation for specific applications. Using a simulated data sequence with two close frequencies, the estimated spectrum from a two-level decomposition with autoregressive modeling shows better resolution than with conventional processing. Geodetic and geophysical applications are briefly indicated.     

The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.     

Response of GPS occultation signals to atmospheric gravity waves and retrieval of gravity wave parameters

We show that the amplitude of the Global Positioning System (GPS) signals in the radio occultation (RO) experiments is an indicator of the activity of the gravity waves (GW) in the atmosphere. The amplitude of the GPS RO signals is more sensitive to the atmospheric wave structures than is the phase. Early investigations used only the phase of the GPS occultation signals for statistical investigation of the GW activity in the height interval 10–40 km on a global scale. In this study, we use the polarization equations and Hilbert transform to find the 1-D GW radio image in the atmosphere by analyzing the amplitude of the RO signal. The radio image, also called the GW portrait, consists of the phase and amplitude of the GW as functions of height. We demonstrate the potential of this method using the amplitude data from GPS/Meteorology (GPS/MET) and satellite mission Challenge Mini-satellite Payload (CHAMP) RO events. The GW activity is nonuniformly distributed with the main contribution associated with the tropopause and the secondary maximums related to the GW breaking regions. Using our method we find the vertical profiles of the horizontal wind perturbations and its vertical gradient associated with the GW influence. The estimated values of the horizontal wind perturbations are in fairly good agreement with radiosonde data. The horizontal wind perturbations v(h) are ±1 to ±5 m s with vertical gradients dv/dh ±0.5 to ±15 m s km at height 10–40 km. The height dependence of the GW vertical wavelength was inferred through the differentiation of the GW phase. Analysis of this dependence using the dispersion relationship for the GW gives the estimation of the projection of the horizontal background wind velocity on the direction of the GW propagation. For the event considered, the magnitude of this projection changes between 1.5 and 10 m s at heights of 10–40 km. We conclude that the amplitude of the GPS occultation signals contain important information about the wave processes in the atmosphere on a global scale.     

Improved relativistic transformations in GPS

For GPS satellite clocks, a nominal (hardware) frequency offset and a conventional periodic relativistic correction derived as a dot product of the satellite position and velocity vectors, are used to compensate the relativistic effects. The conventional hardware clock rate offset of 38,575.008 ns/day corresponds to a nominal orbit semi-major axis of about 26,561,400 m. For some of the GPS satellites, the departures from the nominal semi-major axis can cause an apparent clock rate up to 10 ns/day. GPS orbit perturbations, together with the earth gravity field oblateness, which is largely responsible for the orbit perturbations, cause the standard GPS relativistic transformations to depart from the rigorous relativity transformation by up to 0.2 ns/day. In addition, the conventional periodic relativistic correction exhibits periodic errors with amplitudes of about 0.1 and 0.2 ns, with periods of about 6 h and 14 days, respectively. Using an analytical integration of the gravity oblateness term (J₂), a simple analytical approximation was derived for the apparent clock rate and the 6-h periodic errors of the standard GPS gravity correction. For daily linear representations of GPS satellite clocks, the improved relativistic formula was found to agree with the precise numerical integration of the GPS relativistic effects within about 0.015 ns. For most of the Block IIR satellites, the 6-h periodical errors of the GPS conventional relativistic correction are already detectable in the recent IGS final clock combinations.