To move or not to move: an applied study of veterinary clinic location in Wellington,New Zealand

This paper demonstrates the utility, in an applied context, of a standard Teitz and Bart heuristic for solving location-allocation problems. The example discussed is the location of veterinary clinics in Wellington, New Zealand. A distinct ‘service-poor’ area is identified in southern Wellington suggesting that the existing system of clinics, especially the largest multi-veterinary practice, is vulnerable to intrusion by competitors. A variety of strategies, involving new clinics and clinic relocation are analysed to ascertain alternative courses of action for both existing clinic operators and competitors. The appropriateness of this type of analysis as an aid to decision-making is reflected in the locational response of one practice to protect its share of the market for veterinary services.     

Pattern and chronology of glacial Lake Peace shorelines and implications for isostacy and ice‐sheet configuration in northeastern British Columbia,Canada

Recognition of positions of glacial lakes along the margin of continental ice sheets is critical in reconstructing ice configuration during deglaciation. Advances in remote sensing technology (e.g. LiDAR) have enabled the generation of accurate digital‐elevation models (DEMs) that reveal unprecedented geomorphic detail. Combined with geographical information systems, these tools have considerably advanced the mapping and correlation of geomorphic features such as relict shorelines. Shorelines of glacial Lake Peace (GLP) developed between the Laurentide and Cordilleran ice sheets in northeastern British Columbia and northwestern Alberta. Shoreline mapping from high resolution DEMs produced more than 55 500 elevation data points from 3231 shorelines, enabling the identification of four major phases of GLP: Phase I (altitude 960–990 m a.s.l.); Phase II (890–915 m a.s.l.); Phase III (810–865 m a.s.l.); and Phase IV (724–733 m a.s.l.). The timing of Phase II of GLP is estimated by two optical ages of <16.0±2.5 and 14.2±0.5 ka BP. Extensive mapping of the shorelines allows for measuring of glacial isostatic adjustment as ice retreated. Shorelines currently dip to the northeast at around 0.4–0.5 m km^?1. This slope reflects the asynchronous retreat of the Cordilleran (CIS) and Laurentide (LIS) ice sheets. The relative uplift in the southwest of the study area within the Rocky Mountains and foothills suggests that the Late Wisconsinan (MIS 2) CIS persisted in the foothill after the LIS lost mass and retreated, or that the Late Wisconsinan CIS was very thick and caused deep crustal loading, which resulted in more uplift in the southwest before reaching equilibrium during, or shortly after deglaciation.     

Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids

We consider the exterior Neumann problem of the Laplacian with boundary condition on spheroids. We propose to use spherical radial basis functions in the solution of the boundary integral equation arising from the Dirichlet-to-Neumann map. Our meshless approach with radial basis functions is particularly suitable for handling scattered satellite data. We also propose a preconditioning technique based on an overlapping domain decomposition method to deal with ill-conditioned matrices arising from the approximation problem.     

Drained cavity expansions in soils of finite radial extent subjected to two boundary conditions

The problem of drained cavity expansion in soils of finite radial extent is investigated. Cylindrical and spherical cavities expanded from zero radius subjected to either constant stress or zero displacement at the finite boundary are considered. The generalised analytical solution procedure presented enables more advanced constitutive models to be implemented than have been possible in previous studies. Results generated for Sydney quartz sand highlight substantial differences between cavity limit pressures for the two boundary conditions and for boundaries of finite and infinite radial extent. This is of significance in accounting for boundary effects when interpreting cone penetration tests conducted in calibration chambers. Copyright © 2011 John Wiley & Sons, Ltd.     

Understanding the fate of iron in a modern temperate estuary: Leirárvogur,Iceland

Fluvial dissolved Fe concentrations decrease upon mixing with seawater, resulting in the formation of Fe-floccules. However, a clear understanding of the fate of these floccules has yet to be established. Assessing how tidal processes affect the formation of Fe-colloids in the Leirárvogur estuary, SW Iceland, is an important step in understanding the formation and potential deposition of estuarine Fe-rich minerals within this estuarine system. The Leirárvogur estuary drains predominately Fe-rich basalt, increasing the likelihood of detecting changes in Fe-phases. Fluvial waters and local lake waters that drain into the estuary were compared and the effects of seasonal changes were considered, in an attempt to understand how varying end-members and external factors play a role in Fe-rich mineral formation. Aqueous and colloidal Fe concentrations were found to be greater towards the head of the Leirárvogur estuary, suggesting that potential Fe-rich minerals and complexes are forming at sites of fluvial input. Increasing suspended colloidal Fe towards the estuary mouth suggests that Fe-colloids are readily transported seaward.     

Landscape controls on long‐term runoff in subhumid heterogeneous Boreal Plains catchments

We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km²) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP^?1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.     

Cyclic strength and deformation of crushable carbonate sand

Cyclic triaxial tests have been carried out on a skeletal carbonate sand from the west coast of Eire. Results are presented for undrained cyclic shear tests on samples with 80% of relative density consolidated under both isotropic and anisotropic conditions. Failure was defined as a 5% double amplitude cyclic strain and a 5% peak axial strain for stress-reversal and non-reversal stress conditions, respectively. Using this definition the cyclic strength for isotropically consolidated samples was affected by the confining pressure although the angular platey nature of the sand resulted in higher cyclic strengths than for a comparable silica sand. For anisotropically consolidated samples the cyclic strength increased with increasing initial static shear stress while on the other hand the cyclic strength normalised in the usual way with respect to the initial confining pressure decreased as this pressure increased.     

Interaction between the controls on fluvial system development: tectonics,climate, base level and river capture – Rio Alias,Southeast Spain

Fluvial systems in uplifting terrain respond to tectonic, climatic, eustatic and local base‐level controls modified by specific local factors, such as river capture. The Rio Alias in southeast Spain is an ephemeral, transverse‐to‐structure fluvial system. The river drains two interconnected Neogene sedimentary basins, the Sorbas and Almeria basins, and crosses two major geological structures, the Sierras de Alhamilla/Cabrera and the Carboneras Fault Zone. Regional epeirogenic uplift resulted in sustained fluvial incision during the Quaternary, punctuated by major climatically driven periods of aggradation and dissection, which created a suite of five river terraces. The river terrace sequence was radically modified in the late Pleistocene by a major river capture (itself a response to regional tectonics), localized tectonic activity and eustatic base‐level change. The Rio Alias is defined by four reaches; within each the climatically‐generated, region‐wide, fluvial response was modified by tectonics, base‐level change or river capture to varying degrees. In the upper part of the basin (Lucainena reach), climate was the dominant control on river development, with limited modification of the sequence by uplift of the Sierra Alhamilla and local drainage reorganization by a local river capture. Downstream of the Sierra Alhamilla in the Polopus reach, the climatic signal is dominant, but its expression is radically modified by the response to a major river capture whereby the Alias system lost up to 70% of its pre‐capture drainage area. In the reach adjacent to the Carboneras Fault Zone (Argamason reach), modification of the terrace sequence by local tectonic activity and a resultant local base‐level fall led to a major local incisional event (propagating c. 3–4 km upstream from the area of tectonic disturbance). At the seaward end of the system (El Saltador reach) Quaternary sea‐level changes modified the patterns of erosion and incision and have resulted in steep incisional terrace profiles. The signals generated by regional tectonics and the Quaternary climate change can be identified throughout the basin but those generated by ongoing local tectonics, river capture and sea‐level change are spatially restricted and define the four reaches. The connectivity of the system from the headwaters to the coast decreased through time as incision progressed, resulting in changes in local coupling characteristics. Copyright © 2012 John Wiley & Sons, Ltd.     

Metasomatic overprinting by juvenile igneous fluids,Igdlerfigsalik, South Greenland

The late Igdlerfigsalik centre, part of the Igaliko nepheline syenite complex in the Gardar Province, S Greenland, is bounded to the N and E by compositionally similar, earlier syenites forming the early Igdlerfigsalik and South Qôroq centres. In a circa 1 km wide zone parallel to its contact with late Igdlerfigsalik, the South Qôroq centre shows recrystallisation and alteration of mafic phases. South Qôroq therefore comprises two zones, termed the “unaltered” and the “recrystallised” zones. A study of the biotites from the rocks of the area of the present study has demonstrated variations in biotite chemistry, notably variations in elements inferred to reside on the octahedral and hydroxyl sites. Samples were chosen to provide a representative selection from each centre, avoiding pegmatites and late-stage veins. Electron probe microanalyses demonstrate little Cl in these biotites. Fluorine contents of some biotites can be demonstrated to vary with Fe/(Fe+Mg) and ∑(octahedral Al and Ti). Statistical analysis of the data sets as a whole, however, does not demonstrate simple relationships between fluorine content and these parameters and more complex crystallographic controls over fluorine uptake are inferred. Despite these variations, the relative fluorine contents characteristic of each centre can be determined, which, it is believed, reflect the characteristic F activity of fluids associated with each period of magmatism. In South Qôroq, the unaltered zone has a distinctive F-content in biotite, reflecting the original F-content of fluids associated with South Qôroq magmatism. However in the recrystallised zone, the F-contents are distinct, more closely resembling F-contents characteristic of late Igdlerfigsalik rocks. Reaction between a fluorine-rich, younger fluid and older fluorine-poor rocks is inferred. Despite evidence in the micas for fluid-rock interaction, whole-rock chemical changes are irregular in nature and appear not to conform to simple trends. Mass-balance equations are therefore unable to predict the levels of fluid-rock interaction. Similar exchange between late-stage fluids from younger centres and pre-existing rocks may be commonplace in igneous plutonic environments, and fluorine in micas may provide a relatively sensitive guide to fluid-rock interactions.     

Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.