Successive hydration and dehydration of high-P mafic granofels involving clinopyroxene–kyanite symplectites, Mt Daniel, Fiordland, New Zealand

Three texturally distinct symplectites occur in mafic granofels of the Arthur River Complex at MtDaniel, Fiordland, New Zealand. These include symplectic intergrowths of clinopyroxene and kyanite, described here for the first time. Pods of mafic granofels occur within the contact aureole of the Early Cretaceous Western Fiordland Orthogneiss batholith. The pods have cores formed entirely of garnet and clinopyroxene, and rims of pseudomorphous coarse‐grained symplectic intergrowths of hornblende and clinozoisite that reflect hydration at moderate to high‐P. These hornfelsic rocks are enveloped by a hornblende–clinozoisite gneissic foliation (S1). Narrow garnet reaction zones, in which hornblende and clinozoisite are replaced by garnet–clinopyroxene assemblages, developed adjacent to fractures and veins that cut S1. Fine‐grained symplectic intergrowths of (1) clinopyroxene and kyanite and (2) clinozoisite, quartz, kyanite and plagioclase form part of the garnet reaction zones and partially replace coarse‐grained S1 hornblende and clinozoisite. The development of the garnet reaction zones and symplectites was promoted by dehydration most probably following cooling of the contact aureole. Maps of oxide weight percent and cation proportions, calculated by performing matrix corrections on maps of X‐ray intensities, are used to study the microstructure of the symplectites.     

The Proterozoic P-T-t Evolution of the Kemp Land Coast, East Antarctica; Constraints from Si-saturated and Si-undersaturated Metapelites

Integrated metamorphic and geochronological data place new constraintson the metamorphic evolution of a Neoproterozoic orogen in eastAntarctica. Granulite-facies rocks from a 150 km stretch ofthe Kemp Land coast reflect peak conditions involving T 870–990°Cat P 7·4–10 kbar, with pressure increasing westwardtowards an Archaean craton. Electron microprobe-derived (Th+ U)–Pb monazite ages from metapelitic assemblages indicatethat the major mineral textures in these rocks developed duringthe c. 940 Ma Rayner Orogeny. Complex compositional zoning inmonazite suggests high-T recrystallization over c. 25 Myr. Diversityin metapelitic reaction textures reflects silica and ferromagnesiancontent: Si-saturated Fe-rich metapelites contain garnet thatis partially pseudomorphed by biotite and sillimanite, whereasSi-saturated Mg-rich metapelites and Si-undersaturated metapeliticpods have reaction microstructures involving cordierite enclosingorthopyroxene, garnet and/or sapphirine, cordierite + sapphirinesymplectites around sillimanite and coarse-grained orthopyroxene+ corundum separated by sapphirine coronae. Interpretationsbased on P–T pseudosections provide integrated bulk-rockconstraints and indicate a clockwise P–T–t pathcharacterized by a post-peak P–T trajectory with dP/dT 15–20 bar/ °C. This moderately sloped decompressive-coolingP–T path is in contrast to near-isothermal decompressionP–T paths commonly cited for this region of the RaynerComplex, with implications for the post-collisional tectonicresponse of the mid- to lower crust within this orogenic belt. KEY WORDS: electron microprobe monazite dating; granulite facies; Rayner Complex; sapphirine; THERMOCALCMinerals abbreviations: q, quartz; g, garnet; sill, sillimanite; ky, kyanite; opx, orthopyroxene; cd, cordierite; ksp, alkali feldspar; pl, plagioclase; bi, biotite; sp, spinel; ilm, ilmenite; mt, magnetite; ru, rutile; sa, sapphirine; cor, corundum; osm, osumilite; liq, silicate melt; mnz, monazite     

Tracing spiral density waves in M81

We use Spitzer IRAC 3.6 and 4.5 μm near-infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and Two-Micron All-Sky Survey K _s-band data to produce mass surface density maps of M81. The IRAC 3.6- and 4.5-μm data, whilst dominated by emission from old stellar populations, are corrected for small-scale contamination by young stars and polycyclic aromatic hydrocarbon emission. The I -band data are used to produce a mass surface density map by a   B − V   colour correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps, we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8-μm IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long-lived and allows the position of corotation to be determined.     

Changes in the benthic macroinvertebrate fauna of a large microtidal estuary following extreme modifications aimed at reducing eutrophication

An artificial channel was opened in 1994 between the microtidal Peel–Harvey Estuary and the Indian Ocean to increase tidal exchange and thus ameliorate the problems of eutrophication. Although this greatly reduced macroalgal and cyanobacterial growths and the amount of particulate organic matter, our data indicate that, contrary to managerial expectations, the benthic environment has deteriorated. Thus, although macroinvertebrate density has declined as predicted, taxonomic distinctness (Δ^*) has also declined and species composition has become more variable. Macroinvertebrate composition has also changed markedly at the species, family and even phylum levels. The Crustacea, the most sensitive of the major macrobenthic taxa to environmental stress, has become proportionally less abundant and speciose, whereas the Polychaeta, the least sensitive, was unique in showing the reverse trend. The benthos of the Peel–Harvey Estuary is thus apparently more stressed than previously, probably due to the multiple effects of a great increase in system use.     

Seismicity in the Rotorua and Kawerau geothermal systems,Taupo Volcanic Zone,New Zealand,based on improved velocity models and cross-correlation measurements

We analyze earthquakes occurring in and around the Rotorua and Kawerau geothermal systems, Taupo Volcanic Zone, New Zealand. The two data sets contain 504 and 1875 shallow (≤ 20 km deep) earthquakes, respectively, and span the 21 year period between 1984 and 2004. The arrival time data for these earthquakes are first used to calculate 1-D P- and S-wave seismic velocity models and accompanying station correction terms for both areas. In order to address the non-uniqueness of the joint hypocenter-velocity model estimation problem, we analyze suites of 1000 velocity models computed from random initial models. The final velocity models are well constrained, particularly at depths between 4 and 15 km, and consistent with the results obtained in previous seismic refraction studies of the central Taupo Volcanic Zone. Using a combination of cross-correlation-derived and catalog-based arrival times, we relocate subsets of the Rotorua and Kawerau data sets. In Rotorua, the relocated earthquakes cluster near the geothermally active parts of Rotorua City and beneath the Mount Ngongotaha rhyolite dome. Earthquake clusters and alignments reveal seismogenic structures in the mid-crust whose positions and geometries are consistent with previously published fault mechanisms and known near-surface faults. In Kawerau, the earthquakes within the geothermal field align along northeast-trending lineations, consistent with the predominant alignment of surface-mapped faults in the area.     

High-resolution precipitation and temperature downscaling for glacier models

The spatial resolution gap between global or regional climate models and the requirements for local impact studies motivates the need for climate downscaling. For impact studies that involve glacier modelling, the sparsity or complete absence of climate monitoring activities within the regions of interest presents a substantial additional challenge. Downscaling methods for this application must be independent of climate observations and cannot rely on tuning to station data. We present new, computationally-efficient methods for downscaling precipitation and temperature to the high spatial resolutions required to force mountain glacier models. Our precipitation downscaling is based on an existing linear theory for orographic precipitation, which we modify for large study regions by including moist air tracking. Temperature is downscaled using an interpolation scheme that reconstructs the vertical temperature structure to estimate surface temperatures from upper air data. Both methods are able to produce output on km to sub-km spatial resolution, yet do not require tuning to station measurements. By comparing our downscaled precipitation (1 km resolution) and temperature (200 m resolution) fields to station measurements in southern British Columbia, we evaluate their performance regionally and through the annual cycle. Precipitation is improved by as much as 30% (median relative error) over the input reanalysis data and temperature is reconstructed with a mean bias of 0.5°C at locations with high vertical relief. Both methods perform best in mountainous terrain, where glaciers tend to be concentrated.