The GPS Toolbox

The GPS Toolbox is dedicated to highlighting algorithms utilized by GPS engineers and scientists. If you have an interesting algorithm you would like to share with our readers or if you have a topic you would like to see covered in a future column, contact us at gps-toolbox@ngs.noaa.gov. To comment on the algorithms presented here, or to leave a request for an algorithm you may be looking for, visit our Web site (http://www.ngs.noaa.gov/gps-toolbox). ? 2000 John Wiley & Sons, Inc.     

A Hybrid Finite-Difference and Analytic Element Groundwater Model

Regional finite-difference models tend to have large cell sizes, often on the order of 1–2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW–MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models.     

Tectono-sedimentary evolution of Jurassic–Cretaceous diapiric structures: Miravete anticline,Maestrat Basin,Spain

Integration of extensive fieldwork, remote sensing mapping and 3D models from high-quality drone photographs relates tectonics and sedimentation to define the Jurassic–early Albian diapiric evolution of the N–S Miravete anticline, the NW-SE Castel de Cabra anticline and the NW-SE Cañada Vellida ridge in the Maestrat Basin (Iberian Ranges, Spain). The pre shortening diapiric structures are defined by well-exposed and unambiguous halokinetic geometries such as hooks and flaps, salt walls and collapse normal faults. These were developed on Triassic salt-bearing deposits, previously misinterpreted because they were hidden and overprinted by the Alpine shortening. The Miravete anticline grew during the Jurassic and Early Cretaceous and was rejuvenated during Cenozoic shortening. Its evolution is separated into four halokinetic stages, including the latest Alpine compression. Regionally, the well-exposed Castel de Cabra salt anticline and Cañada Vellida salt wall confirm the widespread Jurassic and Early Cretaceous diapiric evolution of the Maestrat Basin. The NE flank of the Cañada Vellida salt wall is characterized by hook patterns and by a 500-m-long thin Upper Jurassic carbonates defining an upturned flap, inferred as the roof of the salt wall before NE-directed salt extrusion. A regional E-W cross section through the Ababuj, Miravete and Cañada-Benatanduz anticlines shows typical geometries of salt-related rift basins, partly decoupled from basement faults. These structures could form a broader diapiric region still to be investigated. In this section, the Camarillas and Fortanete minibasins displayed well-developed bowl geometries at the onset of shortening. The most active period of diapiric growth in the Maestrat Basin occurred during the Early Cretaceous, which is also recorded in the Eastern Betics, Asturias and Basque-Cantabrian basins. This period coincides with the peak of eastward drift of the Iberian microplate, with speeds of 20 mm/year. The transtensional regime is interpreted to have played a role in diapiric development.     

Securing the Neoliberal City: Discourses of Creativity and Priority Neighborhoods in Toronto,Canada

This paper examines the convergence between the discourse of the creative city and the discourse on priority neighborhoods within urban policy imaginaries in Toronto, Canada. In particular, it examines the development of a number of arts programs targeted at low-income neighborhoods in the city. The twin objectives of these programs are to (a) foster creative and entrepreneurial subjectivities among “at-risk” youth, and (b) reduce the risk of violence that is presumed to be associated with youths living in poor neighborhoods. The paper analyzes how these two discourses are intertwined in a neoliberal politics devoted to enhancing regional quality-of-life competitiveness and to branding Toronto as an attractive creative-class destination. [Key words: creative city, poverty, low-income neighborhoods, neoliberalism, Toronto.]     

The thermodynamic influence of subgrid orography in a global climate model

Assessments of the impacts of climate change typically require information at scales of 10 km or less. Such a resolution in global climate simulations is unlikely for at least two decades. We have developed an alternative to explicit resolution that provides a framework for meeting the needs of climate change impact assessment much sooner. We have applied to a global climate model a physically based subgrid-scale treatment of the influence of orography on temperature, clouds, precipitation, and land surface hydrology. The treatment represents subgrid variations in surface elevation in terms of fractional area distributions of discrete elevation classes. For each class it calculates the height rise/descent of air parcels traveling through the grid cell, and applies the influence of the rise/descent to the temperature and humidity profiles of the elevation class. Cloud, radiative, and surface processes are calculated separately for each elevation class using the same physical parametrizations used by the model without the subgrid orography parametrization. The simulated climate fields for each elevation class can then be distributed in post-processing according to the spatial distribution of surface elevation within each grid cell. Parallel 10-year simulations with and without the subgrid treatment have been performed. The simulated temperature, precipitation and snow water are mapped to 2.5-minute (~5 km) resolution and compared with gridded analyses of station measurements. The simulation with the subgrid scheme produces a much more realistic distribution of snow water and significantly more realistic distributions of temperature and precipitation than the simulation without the subgrid scheme. Moreover, the 250-km grid cell means of most other fields are virtually unchanged by the subgrid scheme. This suggests that the tuning of the climate model without the subgrid scheme is also applicable to the model with the scheme.     

Evaluating techniques for sampling stream crayfish (Paranephrops planifrons)

We evaluated several capture and analysis techniques for estimating abundance and size structure of freshwater crayfish (Paranephrops planifrons) (koura) from a forested North Island, New Zealand stream to provide a methodological basis for future population studies. Direct observation at night and collecting with baited traps were not considered useful. A quadrat sampler was highly biased toward collecting small individuals. Handnetting at night and estimating abundances using the depletion method were not as efficient as handnetting on different dates and analysing by a mark‐recapture technique. Electrofishing was effective in collecting koura from different habitats and resulted in the highest abundance estimates, and mark‐recapture estimates appeared to be more precise than depletion estimates, especially if multiple recaptures were made. Handnetting captured more large crayfish relative to electrofishing or the quadrat sampler.     

Knowledge-Driven Prospectivity Mapping for Granite-Related Polymetallic Sn–F–(REE) mineralization,Bushveld Igneous Complex,South Africa

This paper presents mineral prospectivity mapping to identify potential new exploration ground for polymetallic Sn–F–REE mineralization associated with the Bushveld granites of the Bushveld Igneous Complex, South Africa. The Lebowa Granite Suite, commonly known as the Bushveld granites, is host to a continuum of polymetallic mineralization with a wide range of metal assemblages (Sn–Mo–W–Cu–Pb–Zn–As–Au–Ag–Fe–F–U–REE), ranging from a high-temperature to a low-temperature magmatic hydrothermal mineralizing environment. The prospectivity map was generated by fuzzy logic modeling and a selection of targeting criteria (or spatial proxies) based on a conceptual mineral system highlighting critical processes responsible for the formation of the polymetallic mineralization. The spatial proxies include proximity to differentiated granites (as heat and metal-rich fluid sources), Rb geochemical map (fluid-focusing mechanism such as fractionation process), principal component maps (PC 4 Y–Th and PC 14 Sn–W, fluid pathways for both high- and low-temperature mineralization) and proximity to roof rocks (traps for fluids). Logarithmic functions were used to rescale rasterized evidential maps into continuous fuzzy membership scores in a range of [0, 1]. The evidential maps were combined in two-staged integration matrix using fuzzy AND, OR and gamma operators to produce the granite-related polymetallic Sn–F–(REE) prospectivity map. The conceptual mineral system model and corresponding prospectivity model developed in this study yielded an encouraging result by delineating the known mineral deposits and occurrences of Sn–F–(REE) mineralization that were not used to assign weights to the evidential maps. The prospectivity model predicted, on average, 77% of the known mineral occurrences in the BIC (i.e., 56 of 73 Sn occurrences, 12 of 15 F occurrences and 6 of 8 REE occurrences). Based on this validation, 13 new targets were outlined in this study.     

Spatial and temporal variation in subtidal molluscan diversity amongst temperate estuarine habitats

Effective management of marine ecosystems is enhanced when detailed information on biodiversity is available. Key information to underpin management actions and conservation planning includes relationships between species assemblages and environmental gradients, and information on species distributions. We conducted a subtidal biodiversity assessment of surface‐dwelling subtidal molluscs in eight a priori defined habitat types using underwater visual censuses to quantitatively explore relationships between molluscan assemblages, and their correlation with benthic habitats and abiotic variables. In addition, variations in diversity were examined for two key habitat types (areas dominated by Dendronephthya australis and by filter feeders) over a period of 15 months to examine temporal change. We found that molluscs form distinct assemblages within subtidal habitats, but that assemblages within key habitats show inherent temporal variability. Regional (gamma) diversity of molluscs was found to result from a combination of: (i) within habitat alpha diversity, which increased with habitat complexity; (ii) between habitat beta diversity, with significant differences in molluscan assemblages amongst habitats with differing benthic growth, substrate type, and depth; and (iii) temporal beta diversity, with significant changes detected in molluscan assemblages over time. The results demonstrate how habitats and abiotic variables (principally depth and substrate type) combine to contribute to molluscan biodiversity in temperate estuaries, and illustrate the value of these factors as surrogates for surface‐dwelling subtidal molluscs in conservation planning.