Petrogenesis of kyanite‐ and corundum‐bearing mafic granulite in a meta‐ophiolite,SE Turkey

Although ophiolitic rocks are abundant in Anatolia (Turkey), only in rare cases have they experienced high‐grade metamorphism. Even more uncommon, in Anatolia and elsewhere are high‐grade meta‐ophiolites that retain an oceanic lithosphere stratigraphy from upper crustal mafic volcanic rocks through lower crustal gabbro to mantle peridotite. The Berit meta‐ophiolite of SE Turkey exhibits both features: from structurally higher to lower levels, it consists of garnet amphibolite (metabasalt), granulite facies metagabbro (as lenses in amphibolite inferred to be retrogressed granulite) and metaperidotite (locally with metapyroxenite layers). Whole‐rock major and trace‐element data indicate a tholeiitic protolith that formed in a suprasubduction setting. This paper presents new results for the metamorphic P–T conditions and path of oceanic lower crustal rocks in the Berit meta‐ophiolite, and an evaluation of the tectonic processes that may drive granulite facies metamorphism of ophiolite gabbro. In the Do?an?ehir (Malatya, Turkey) region, granulite facies gabbroic rocks contain garnet (Grt)+clinopyroxene (Cpx)+plagioclase (Pl)+corundum (Crn)±orthopyroxene (Opx)±kyanite (Ky)±sapphirine (Spr)±rutile. Some exhibit symplectites consisting of Crn+Cpx, Ky+Cpx and/or coronas of garnet (outer shell) around a polygonal aggregate of clinopyroxene that in some cases surrounds a polygonal aggregate of orthopyroxene. Coronitic and non‐coronitic textures occur in proximity in mm‐ to cm‐scale layers; corona structures typically occur in plagioclase‐rich layers. Their formation is therefore related primarily to protolith type (troctolite v. gabbro) rather than P–T path. Phase diagrams calculated for a kyanite‐rich granulite, a plagioclase‐rich non‐coronitic granulite, and a plagioclase‐rich coronitic granulite (taking into account changes in effective bulk composition during texture development) predict peak conditions of ~800°C, 1.1–1.5 GPa; these conditions do not require invoking an unusually high geothermal gradient. In the coronitic metagabbro, reaction textures formed along the prograde path: Crn–Cpx symplectites grew at the expense of garnet, sapphirine and plagioclase. Peak conditions were followed by isobaric cooling of ~150°C. Hornblende–plagioclase thermometry results for host amphibolite (Hbl+Pl±Crn±Grt±relict Cpx) indicate retrograde conditions of 620–675°C and 0.5–0.8 GPa accompanied by infiltration of H₂O‐rich fluid. This anticlockwise P–T path differs from an isothermal decompression path previously proposed for these rocks based on the presence of symplectite. Metamorphism of the ophiolitic rocks was driven by closing of the southern Neotethys Ocean, as oceanic lithosphere was obducted (most SE Anatolian ophiolites) or underthrust (Berit meta‐ophiolite). This was followed by subduction of a continental margin, driving cooling of the Berit granulite after the thermal peak at depths of ~40 km.     

Reconstructing estuarine conditions: oyster shells as recorders of environmental change, Southwest Florida

Live-collected shells of the oyster, Crassostrea virginica, contain geochemical records of modern temperature and salinity, so records of prehistoric conditions may be obtained from subfossil shells. Restoration of channelized watersheds in Florida is receiving much attention, and plans for targeted watersheds require information about estuarine conditions before channelization. Lack of historical records necessitates alternative methods to understand pre-disturbance conditions. A ¹⁴C-calibrated, amino-acid geochronology based on racemization of glutamic acid yielded ages ranging from 190–1220 AD and from 1270–1860 AD for subfossil oysters from Blackwater River (near-natural watershed) and for Faka-Union Bay (channelized watershed), respectively. δ¹⁸O and δ¹³C values of subfossil shells from Blackwater River indicate salinity and summer temperatures similar to present. Winter temperatures recorded in shells from 190, 590, 720, and 1050 AD appear 1–5 °C colder than present winter temperatures, whereas the shell from 1220 AD records winter temperatures similar to modern winter temperatures. These temperature shifts may indicate change in climate or natural seasonal variation of winter temperature from year to year. Subfossils from Faka-Union Bay may reflect a complicated hydrology, which cannot be evaluated by isotopic compositions alone and demonstrates the need for development of independent elemental proxies for temperature and salinity. Decreases in δ¹³C from subfossil to modern shells may in part result from CO₂ added to the atmosphere from fossil fuel burning (the Suess effect). Subfossil δ¹³C that is >1‰ more positive than modern shells suggest a change in the dominant carbon sources from terrestrial C₄ or aquatic plants to C₃ plants (mangroves).     

Small scale study of groundwater flow in a fractured carbonate-rock aquifer at the St-Eustache quarry, Québec, Canada

A small-scale hydrogeological study was conducted in a fractured carbonate-rock aquifer located in a quarry to relate groundwater flow to the fracture network. The field study in the St-Eustache quarry, which integrates structural surveys, well logging and hydraulic tests, showed that the most important features that affect groundwater flow in the sedimentary aquifer are high hydraulic conductivity horizontal bedding planes. Vertical fractures are abundant in the quarry and throughout the region, but they have a minor effect on groundwater flow. To have a significant impact on the flow regime and lead to vertical groundwater flow, the permeability of all vertical joints need to be enhanced compared to what was generally observed at the site. Such an increase in permeability could potentially occur where dissolution and fracturing is more intense or at stress release locations such as near the surface in the quarry.     

Microplate motions in the hotspot reference frame

We have computed motions of the major plates (seven large plates and seven medium‐sized plates) and 38 microplates relative to the hotspot reference frame, and present velocities of these 52 plates. Moreover, using updated plate boundaries for the present, we have computed new geometrical factors for plates and microplates, useful for kinematic calculations and to obtain the net‐rotation of the lithosphere and plate velocities in the mean‐lithosphere reference frame. Instead of a continuum or gradational distribution of the plates by size, the plates clearly partition into three groups each having their own characteristics. For the seven large plates, rotation poles generally lie in high latitudes; the seven medium‐sized plates have rotation poles in a restricted equatorial area; the 38 small plates show the greatest scatter. Moreover subsets of the 52 plates reveal differing fractal behaviour: the large, middle and small groupings each have a characteristic fractal dimension, suggestive of microplate clustering. The highest angular velocities occur for some of the smallest plates, with the location of their rotation poles closeby. Terra Nova, 18, 276–281, 2006     

Forecasting vertical ground surface movement from shrinking/swelling soils with artificial neural networks

Artificial neural networks (ANNs) are used to estimate vertical ground surface movement when soils expand and contract due to changes in soil moisture content caused by changing climate conditions. Several counterpropagation ANN test cases were investigated to map climate data (i.e. temperature and rainfall) to vertical ground surface movement at field sites in Texas and Australia. Three of the four ANN test cases use a historical time series of climate data to forecast ground surface elevation relative to a specified datum. The fourth ANN test case predicts the rate of ground surface movement, and requires post‐processing of the predicted rates to calculate ground surface elevation relative to a specified datum. The counterpropagation network has demonstrated a successful mapping of temperature and rainfall data to vertical ground surface movement at a field site when it is trained with a subset of data from the same field site (test cases 1 and 2). The results of training an ANN on one field site and testing it on another field site (test cases 3 and 4) demonstrate the ability of the ANN to capture trends in vertical ground surface movement. When compared with the predictions from a physics‐based method (shrink test‐water content method) that requires measurements/estimates of changes in soil water content, the ANN‐based predictions (based on climatic changes) captured the trends in the field measurements of shrinking–swelling soil surface movements equally well. These findings are promising and merit further investigation with data from additional field sites. Copyright © 2007 John Wiley & Sons, Ltd.     

Interpretation of tectonic versus magmatic processes for resetting apatite fission track ages in the Niğde Massif, Turkey

A network of 10 OBS operating in continuous recording mode was deployed in June 2003 in the Galicia passive Margin, surrounding the sinking zone of the ‘Prestige’ oil tanker. The OBSs recorded very different signals, including air-gun shots, seismic events, whale vocalizations or noise associated to ships, as well as two sets of atypical seismic signals that are the object of this contribution and can be described as harmonic tremors and short-duration events. Harmonic tremors appear in most sites as intense monochromatic signals lasting for up to 4 h, with coherent polarization and a wide range of amplitudes. Even if they may be recorded at different sites within a time interval, their envelopes are not correlated, suggesting a source origin at the vicinity of each receiver. A time periodicity close to 6 h is observed, suggesting a correlation with the tidal variation. Most OBSs show also a large number of impulsive events of short-duration, between 1 and 2 s, a very regular decrease of amplitude in the signal coda and characteristic frequency spectrum with one or two narrow peaks. As for the tremors, those events are detected only in individual sites, suggesting very local source zones. We favor the interpretation of both seismic features as coming from the resonance of fluid-filled cracks induced by impulsive pressure transients, in agreement with the theoretical model of Chouet [Chouet, B. 1988. Resonance of a fluid-driven crack: radiation properties and implications for the source of long-period events and harmonic tremor. J. Geophys. Res., 93, 4375–4400.]. The observed harmonic tremors would be the result of sustained pressure fluctuations, probably related to the stress variations induced by the tidal change of the oceanic load, while the short duration events are interpreted as the impulse response to the tremor-generating system.     

A method for discovery and analysis of temporal patterns in complex event data

Pattern analysis techniques currently common within geography tend to focus either on characterizing patterns of spatial and/or temporal recurrence of a single event type (e.g., incidence of flu cases) or on comparing sequences of a limited number of event types where relationships between events are already represented in the data (e.g., movement patterns). The availability of large amounts of multivariate spatiotemporal data, however, requires new methods for pattern analysis. Here, we present a technique for finding associations among many different event types where the associations among these varying event types are not explicitly represented in the data or known in advance. This pattern discovery method, known as T-pattern analysis, was first developed within the field of psychology for the purpose of finding patterns in personal interactions. We have adapted and extended the T-pattern method to take the unique characteristics of geographic data into account and implemented it within a geovisualization toolkit for an integrated computational-geovisual environment we call STempo. To demonstrate how T-pattern analysis can be employed in geographic research for discovering patterns in complex spatiotemporal data, we describe a case study featuring events from news reports about Yemen during the Arab Spring of 2011–2012. Using supplementary data from the Global Database of Events, Language, and Tone, we briefly summarize and reference a separate validation study, then evaluate the scalability of the T-pattern approach. We conclude with ideas for further extensions of the T-pattern technique to increase its utility for spatiotemporal analysis.     

Creating a thinking space for the rising stars of Australian urban geography

The Institute of Australian Geographers Urban Study Group and Australian Geographer present the first essay in our Leading Insights/Thinking Space partnership. We discuss Dr Natalie Osborne's 2018 essay, flag Dr Marilu Melo's 2019 essay and talk about creating a thinking space for the rising stars of Australian urban geography.     

Modelling diatom responses to climate induced fluctuations in the moisture balance in continental Antarctic lakes

The water chemistry of lake systems on the edge of the Antarctic continent responds quickly to changes in the moisture balance. This is expressed as increasing salinity and decreasing lake water level during dry periods, and the opposite during wet periods. The diatom composition of the lakes also changes with these fluctuations in salinity and lake water depth. This is important, as their siliceous remains become incorporated into lake sediments and can provide long-term records of past salinity using transfer functions. In order to develop transfer functions, diatoms and water chemistry data were inter-calibrated from five different East Antarctic oases, namely the Larsemann Hills, the Bølingen Islands, the Vestfold Hills, the Rauer Islands and the Windmill Islands. Results indicate that salinity is the most important environmental variable explaining the variance in the diatom flora in East Antarctic lakes. In oligo- saline lakes the variance is mainly explained by lake water depth. This dataset was used to construct a weighted averaging transfer function for salinity in order to infer historical changes in the moisture balance. This model has a jack-knifed r² of 0.83 and a RMSEP of 0.31. The disadvantage of this transfer function is that salinity changes in oligo-saline lakes are reconstructed inaccurately due to the edge effect and due to the low species turnover along the salinity gradient at its lower end. In order to infer changes in the moisture balance in these lakes, a second transfer function using weighted averaging partial least squares (with two components) for depth was constructed. This model has a jack-knifed r² of 0.76 and a RMSEP of 0.22. Both transfer functions can be used to infer climate driven changes in the moisture balance in lake sediment cores from oligo-, hypo-, meso- and hyper-saline lakes in East Antarctic oases between 102–75°E. The transfer function for lake water depth is promising to track trends in the moisture balance of small freshwater lakes, where changes in shallow and deep-water sediments are readily reflected in changing diatom composition.     

An ensemble approach to space–time interpolation

The availability of spatial data on an unprecedented scale as well as advancements in analytical and visualization techniques gives researchers the opportunity to study complex problems over large urban and regional areas. Nevertheless, few individual data sets exist that provide both the requisite spatial and/or temporal observational frequency to truly facilitate detailed investigations. Some data are collected frequently over time but only at a few geographic locations (e.g., weather stations). Similarly, other data are collected with a high level of spatial resolution but not at regular or frequent time intervals (e.g., satellite data). The purpose of this article is to present an interpolation approach that leverages the relative temporal richness of one data set with the relative spatial richness of another to fill in the gaps. Because different interpolation techniques are more appropriate than others for specific types of data, we propose a space–time interpolation approach whereby two interpolation methods – one for the temporal and one for the spatial dimension – are used in tandem to increase the accuracy results.

We call our ensemble approach the space–time interpolation environment (STIE). The primary steps within this environment include a spatial interpolation processor, a temporal interpolation processor, and a calibration processor, which enforces phenomenon-related behavioral constraints. The specific interpolation techniques used within the STIE can be chosen on the basis of suitability for the data and application at hand. In this article, we first describe STIE conceptually including the data input requirements, output structure, details of the primary steps, and the mechanism for coordinating the data within those steps. We then describe a case study focusing on urban land cover in Phoenix, Arizona, using our working implementation. Our empirical results show that our approach increased the accuracy for estimating urban land cover better than a single interpolation technique.