Spatial process and data models: Toward integration of agent-based models and GIS

The use of object-orientation for both spatial data and spatial process models facilitates their integration, which can allow exploration and explanation of spatial-temporal phenomena. In order to better understand how tight coupling might proceed and to evaluate the possible functional and efficiency gains from such a tight coupling, we identify four key relationships affecting how geographic data (fields and objects) and agent-based process models can interact: identity, causal, temporal and topological. We discuss approaches to implementing tight integration, focusing on a middleware approach that links existing GIS and ABM development platforms, and illustrate the need and approaches with example agent-based models.     

Geomasking sensitive health data and privacy protection: an evaluation using an E911 database

Geomasking is used to provide privacy protection for individual address information while maintaining spatial resolution for mapping purposes. Donut geomasking and other random perturbation geomasking algorithms rely on the assumption of a homogeneously distributed population to calculate displacement distances, leading to possible under-protection of individuals when this condition is not met. Using household data from 2007, we evaluated the performance of donut geomasking in Orange County, North Carolina. We calculated the estimated k-anonymity for every household based on the assumption of uniform household distribution. We then determined the actual k-anonymity by revealing household locations contained in the county E911 database. Census block groups in mixed-use areas with high population distribution heterogeneity were the most likely to have privacy protection below selected criteria. For heterogeneous populations, we suggest tripling the minimum displacement area in the donut to protect privacy with a less than 1% error rate.     

Natural Attenuation of Perchlorate in Denitrified Groundwater

Monitoring of a well‐defined septic system groundwater plume and groundwater discharging to two urban streams located in southern Ontario, Canada, provided evidence of natural attenuation of background low level (ng/L) perchlorate (ClO₄^?) under denitrifying conditions in the field. The septic system site at Long Point contains ClO₄^? from a mix of waste water, atmospheric deposition, and periodic use of fireworks, while the nitrate plume indicates active denitrification. Plume nitrate (NO₃^?‐N) concentrations of up to 103 mg/L declined with depth and downgradient of the tile bed due to denitrification and anammox activity, and the plume was almost completely denitrified beyond 35 m from the tile bed. The ClO₄^? natural attenuation occurs at the site only when NO₃^?‐N concentrations are <0.3 mg/L, after which ClO₄^? concentrations decline abruptly from 187 ± 202 to 11 ± 15 ng/L. A similar pattern between NO₃^?‐N and ClO₄^? was found in groundwater discharging to the two urban streams. These findings suggest that natural attenuation (i.e., biodegradation) of ClO₄^? may be commonplace in denitrified aquifers with appropriate electron donors present, and thus, should be considered as a remediation option for ClO₄^? contaminated groundwater.     

Effects of a Dual‐Pump Crude‐Oil Recovery System,Bemidji, Minnesota,USA

A crude‐oil spill occurred in 1979 when a pipeline burst near Bemidji, MN. In 1998, the pipeline company installed a dual‐pump recovery system designed to remove crude oil remaining in the subsurface at the site. The remediation from 1999 to 2003 resulted in removal of about 115,000 L of crude oil, representing between 36% and 41% of the volume of oil (280,000 to 316,000 L) estimated to be present in 1998. Effects of the 1999 to 2003 remediation on the dissolved plume were evaluated using measurements of oil thicknesses in wells plus measurements of dissolved oxygen in groundwater. Although the recovery system decreased oil thicknesses in the immediate vicinity of the remediation wells, average oil thicknesses measured in wells were largely unaffected. Dissolved‐oxygen measurements indicate that a secondary plume was caused by disposal of the pumped water in an upgradient infiltration gallery; this plume expanded rapidly immediately following the start of the remediation in 1999. The result was expansion of the anoxic zone of groundwater upgradient and beneath the existing natural attenuation plume. Oil‐phase recovery at this site was shown to be challenging, and considerable volumes of mobile and entrapped oil remain in the subsurface despite remediation efforts.     

Miocene silicic volcanism in southwestern Idaho: geochronology,geochemistry, and evolution of the central Snake River Plain

New ⁴⁰Ar-³⁹Ar geochronology, bulk rock geochemical data, and physical characteristics for representative stratigraphic sections of rhyolite ignimbrites and lavas from the west-central Snake River Plain (SRP) are combined to develop a coherent stratigraphic framework for Miocene silicic magmatism in this part of the Yellowstone ‘hotspot track’. The magmatic record differs from that in areas to the west and east with regard to its unusually large extrusive volume, broad lateral scale, and extended duration. We infer that the magmatic systems developed in response to large-scale and repeated injections of basaltic magma into the crust, resulting in significant reconstitution of large volumes of the crust, wide distribution of crustal melt zones, and complex feeder systems for individual eruptive events. Some eruptive episodes or ‘events’ appear to be contemporaneous with major normal faulting, and perhaps catastrophic crustal foundering, that may have triggered concurrent evacuations of separate silicic magma reservoirs. This behavior and cumulative time-composition relations are difficult to relate to simple caldera-style single-source feeder systems and imply complex temporal-spatial development of the silicic magma systems. Inferred volumes and timing of mafic magma inputs, as the driving energy source, require a significant component of lithospheric extension on NNW-trending Basin and Range style faults (i.e., roughly parallel to the SW–NE orientation of the eastern SRP). This is needed to accommodate basaltic inputs at crustal levels, and is likely to play a role in generation of those magmas. Anomalously high magma production in the SRP compared to that in adjacent areas (e.g., northern Basin and Range Province) may require additional sub-lithospheric processes. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. This paper constitutes part of a special issue dedicated to Bill Bonnichsen on the petrogenesis and volcanology of anorogenic rhyolites.     

Scalar advection schemes for ocean modelling on unstructured triangular grids

Several schemes for scalar advection on unstructured triangular grids are assessed for possible use in ocean modelling applications. Finite element, finite volume and finite volume–element approaches are evaluated. A series of tests, including a numerical order of convergence analysis, idealized rotating cone and cylinder experiments, and transport of a tracer through the Stommel Gyre representation of ocean basin-scale circulation, are carried out. Volume element Eulerian–Lagrangian and third-order Runge-Kutta discontinuous Galerkin schemes are recommended for use in tracer studies. Taylor–Galerkin and second-order Runge–Kutta discontinuous Galerkin are found to be robust and accurate second-order schemes. When positivity is required, a fluctuation redistribution scheme was found to be an easily implemented, accurate, and computationally efficient approach. Responsible editor: Phil Dyke     

Paleomagnetism of Late Cretaceous Poços de Caldas alkaline complex,Southern Brazil

The paleomagnetism of the Late Cretaceous Poços de Caldas alkaline complex (46.6°W, 21.9°S) was investigated through 42 oriented cores from seven sites. Six sites, reversed relative to the present magnetic field of the Earth, yield a pole at 127°W, 82°S (dp = 8°,dm = 13°). This pole is located close to other Late Cretaceous poles for South America obtained by Creer [1] from untreated paleomagnetic samples. The results are significantly different from those for the nearby Early Cretaceous Serra Geral basalt but close to the Triassic pole for South America. The polar wandering path for South America for the Mesozoic seems to be more complicated than anticipated. The available paleomagnetic information may not yet be precise enough to determine the time of opening of the Atlantic.