Electromagnetic induction studies in the Eyre Peninsula, South Australia

Magnetic field fluctuations have been recorded by an array of portable three-component magnetometers at 60 sites across the Eyre Peninsula in South Australia between December 1993 and March 1995. An additional 54 magnetometer data records, collected prior to 1989 and described by Milligan (1989) and Milligan, White & Chamalaun (1989), were included in the analysis. A major conductive feature in the crust, first noted by White & Milligan (1984) as the Eyre Peninsula Anomaly (EPA), is re-examined to assess its continuity to the north of the original arrays and to investigate its relationship with major tectonic features.
Magnetic-field time-series were converted to induction arrows in the frequency domain. These induction arrows were initially inverted using the minimum-structure 2-D Occam approach to estimate the electrical conductance of the crust. Following this, thin-sheet forward modelling was used to examine the relationship between the conductance and the dominant tectonic features. The principal results of the modelling are that a narrow conductive feature extends inland from the coast about 160 km before terminating, and the conductance is in the range 3000 to 10 000 S, which decreases inland.
A strong correlation exists between the electrical conductance of the Eyre Peninsula and Bouguer gravity anomalies, and in particular the EPA is coincident with a significant Bouguer gravity gradient. There is also good agreement between the locations of the foci of earthquakes of magnitude greater than 4.0 and the EPA. We believe that the anomaly is associated with a geological fracture in the Precambrian upper crust as a result of crustal extension prior to the rifting of Australia from Antarctica in the Jurassic (160 Ma).     

The Late Cretaceous palaeolatitude of the Neotethyan spreading axis in the eastern Mediterranean region

A compilation of available palaeomagnetic data from the Troodos (Cyprus) and Baër–Bassit (Syria) ophiolitic terranes of the eastern Mediterranean Tethyan orogenic belt is presented. The ophiolites represent fragments of oceanic lithosphere generated at a Neotethyan spreading axis in the Late Cretaceous, although debate continues over the tectonic setting of this spreading axis and its position within the eastern Mediterranean palaeogeography. Two types of model reconstructions have been proposed: Type 1—the ophiolites formed in a southerly Neotethyan basin by spreading above an oceanic subduction zone. The Baër–Bassit ophiolite was then emplaced a relatively short distance (tens of kilometers) southwards on to the Arabian continental margin, leaving the Troodos ophiolite isolated in an intra-oceanic setting to the west; and Type 2—the ophiolites formed in a northerly Neotethyan basin by spreading at a ‘normal’ oceanic ridge, with subsequent large-scale thrusting (hundreds of kilometers) to the south of emplaced ophiolites over microcontinental fragments to reach their present positions. Palaeomagnetic determination of the palaeolatitude of the Neotethyan spreading axis is, therefore, of considerable interest.Previous palaeomagnetic analyses have demonstrated the presence of significant, and in some cases extreme, relative tectonic rotations of a variety of origins in both ophiolites. To allow palaeomagnetic data from these rotated units to be combined, an inclination-only formulation of the palaeomagnetic tilt test is employed. This provides unequivocal evidence that both ophiolites retain pre-deformational remanent magnetizations, which are interpreted as original ocean-floor magnetizations acquired close to the time of crustal formation in the Late Cretaceous. The mean inclinations of 37.0±2.6° for the Troodos terrane and 41.1±3.4° for the Baër–Bassit terrane indicate respective palaeolatitudes for the spreading axes of 20.6°N±1.8° and 23.6°N±2.5°, consistent with a Late Cretaceous position between the Arabian and Eurasian margins. These data, together with a well-defined palaeolatitude of 25.5°N±4.5° for the eastern Pontides previously reported in the literature, provide constraints which must be incorporated in any successful tectonic reconstruction of the eastern Mediterranean Tethys. The implications of these constraints for Type 1 and 2 models are discussed using a series of plate tectonic cross-sections constructed along a line extending northwards from the Arabian continental margin. In the absence of palaeomagnetic data from Late Cretaceous rocks of the eastern Taurides, however, it is presently impossible to use these palaeolatitudinal constraints to resolve the root zone debate on a purely palaeomagnetic basis. Solutions which satisfy the constraints may be found for both types of model reconstruction. Additional, published field-based geological considerations, however, strongly support models in which the Troodos and Baër–Bassit (and other southerly) ophiolites were generated in a southern Neotethyan basin, rather than those involving generation in a northerly basin and subsequent large-scale thrust displacement to the south.     

Macro‐scale bed roughness of the siple coast ice streams in West Antarctica

Subglacial bed conditions are key to understanding ice stream behaviour and evolution, with bed roughness re?ecting substrate composition and ?ow resistance. Here we present an analysis of bed roughness in the Siple Coast region of West Antarctica from airborne radio‐echo sounding data. The ice streams are associated generally with low bed roughness values, which decrease downstream. The bed of the slow‐?owing Ice Stream C (～10 000 km²) is also characterized by being smooth at all scales (wavelengths ranging from 5 km to in excess of 40 km). Furthermore, the bed is smooth either side of Ice Stream C. This suggests the location of the ice stream is controlled by internal ice sheet dynamics rather than by bed morphology. If the ice stream were encouraged to migrate laterally, when active, there would be little resistance offered by the subglacial morphology. Other inter‐ice stream regions are rough, however, indicating a subglacial topographic in?uence on ice stream position. Bed roughness increases up‐?ow of ice streams, which, unless the bed is modi?ed, may limit the inland migration of these systems. Copyright © 2004 John Wiley & Sons, Ltd.     

Late Weichselian relative sea-level changes and ice sheet history in southeast Greenland

Relative sea-level (RSL) observations from the margins of the Greenland Ice Sheet (GIS) provide information regarding the timing and rate of deglaciation and constraints on geophysical models of ice sheet evolution. In this paper we present the first RSL record for the southeast sector of the GIS based on field observations completed close to Ammassalik. The local marine limit is c. 69 m above sea-level (asl) and is dated to c. 11 k cal. yrs BP (thousand calibrated years before present) and is a minimum date for ice free conditions at the study site. RSL fell to c. 24 m asl by 9.5 k cal. yrs BP and continued to fall at a decreasing rate to reach close to present by 6.5 k cal. yrs BP. Our chronology agrees with radiocarbon dates from offshore cores that indicate ice free conditions on the adjacent mid-shelf by 15 k cal. yrs BP. We compare the new RSL data with predictions generated using two recently published glaciological models of the GIS that differ in the amount and timing of ice loading and unloading over our study area. These two GIS models are coupled to the same Earth viscosity model and background (global) ice model to aid in the data-model comparison. Neither model provides a close fit to the RSL observations. Based on a preliminary sensitivity study using a suite of Earth viscosity models, we conclude that the poor data-model fit is most likely due to an underestimate of the local ice unloading. An improved fit could be achieved by delaying the retreat of a thicker ice sheet across the continental shelf. A thick ice sheet extending well onto the continental shelf is in agreement with other recent observations elsewhere in east and south Greenland.     

An empirical function to estimate the depths of linear hot sources: Laboratory modeling and field measurements of lava tubes

Estimating depths of buried lava tubes is important for determining the thermal budgets and effusion rates of basaltic volcanic systems. This research used a laboratory experiment scaled to a lava tube system to measure the 3D temperature field surrounding a hot viscous fluid flowing through a buried glass tube while varying conditions such as flow rate and temperature. The depth of the glass tube was changed for different experimental runs. Numerical techniques were applied to model the laboratory experiment. The surface thermal distributions from 166 thermal traverses, constrained to a depth to width ratio of 0.6 to 1.6, were analyzed to empirically derive a depth estimation function using regression techniques. This “Linear Anomaly Surface Transect (LAST)” depth function is a scaleable depth estimation technique which can be solved with thermal imaging data alone. The minimum temperature, maximum temperature and width of a Lorentzian distribution fit to a surface thermal transect, are the only inputs required for the LAST function to estimate the depths of the hot source. The input parameters were then applied to non-laboratory situations including the Kuhio lava tube system in Hawai’i. The LAST function produced depth estimates of ∼ 0.3 m for the Kuhio lava tube in Hawai’i, which did not agree with observations on the ground. This is the result of the complex composition and geometry of an actual lava tube where heat transfer is controlled by more than a simple fluid filling a tube, but also by convection of gasses and fluids in a partially filled passage. Though not effective for lava tubes at this time, the model provides promising results for simple cases applied to engineering and underground coal fires.     

Holocene floristic diversity and richness in northeast Norway revealed by sedimentary ancient DNA (sedaDNA) and pollen

We present a Holocene record of floristic diversity and environmental change for the central Varanger Peninsula, Finnmark, based on ancient DNA extracted from the sediments of a small lake (sedaDNA). The record covers the period c. 10 700 to 3300 cal. a BP and is complemented by pollen data. Measures of species richness, sample evenness and beta diversity were calculated based on sedaDNA sampling intervals and 1000‐year time windows. We identified 101 vascular plant and 17 bryophyte taxa, a high proportion (86%) of which are still growing within the region today. The high species richness (>60 taxa) observed in the Early Holocene, including representatives from all important plant functional groups, shows that modern shrub‐tundra communities, and much of their species complement, were in place as early as c. 10 700 cal. a BP. We infer that postglacial colonization of the area occurred prior to the full Holocene, during the Pleistocene‐Holocene transition, Younger Dryas stadial or earlier. Abundant DNA of the extra‐limital aquatic plant Callitriche hermaphroditica suggests it expanded its range northward between c. 10 200 and 9600 cal. a BP, when summers were warmer than present. High values of Pinus DNA occur throughout the record, but we cannot say with certainty if they represent prior local presence; however, pollen influx values >500 grains cm^?2 a^?1 between c. 8000 and 7300 cal. a BP strongly suggest the presence of pine woodland during this period. As the site lies beyond the modern tree limit of pine, it is likely that this expansion also reflects a response to warmer Early Holocene summers.     

Climate prediction for adaptation: Who needs what?

The role of short- and long-term climate predictions in determining the success of adaptation to climate change is investigated. A simple theoretical model that captures the relationship between adaptive performance, decision structure, and prediction accuracy at different temporal scales is developed, and its implications deduced. It is shown that users who face high adjustment costs (i.e. are inflexible) depend more heavily on accurate long-term predictions than those who are able to adjust their adaptation strategies at low cost. However, the constraints on the accuracy of long-term predictions in order for them to be ‘accurate enough’ to be recommended unreservedly are strongest for precisely those users who desire them most. There are thus intrinsic structural limitations to the utility of long-term predictions for adaptive decision making. In addition, the heterogeneous dependence of adaptive performance on predictions at different temporal scales gives rise to a resource allocation problem for investments aimed at improving predictions so as to maximize their social benefits. It is suggested that an explicit distributional framework that matches users’ needs with scientific feasibility is needed in order to guide such investments, and key issues any such framework will need to address are discussed.     

Galaxy shear estimation from stacked images

Statistics of the weak lensing of galaxies can be used to constrain cosmology if the galaxy shear can be estimated accurately. In general, this requires accurate modelling of unlensed galaxy shapes and the point spread function (PSF). I discuss suboptimal but potentially robust methods for estimating galaxy shear by stacking images such that the stacked-image distribution is closely Gaussian by the central limit theorem. The shear can then be determined by radial fitting, requiring only an accurate model of the PSF rather than also needing to model each galaxy accurately. When noise is significant, asymmetric errors in the centroid must be corrected, but the method may ultimately be able to give accurate unbiased results when there is a high galaxy density with constant shear. It provides a useful baseline for more optimal methods, and a test case for estimating biases, though the method is not directly applicable to realistic data. I test stacking methods on the simple toy simulations with constant PSF and shear provided by the GREAT08 project, on which most other existing methods perform significantly more poorly, and briefly discuss generalizations to more realistic cases. In the Appendix , I discuss a simple analytic galaxy population model where stacking gives optimal errors in a perfect ideal case.     

Ice sheet extent and early deglacial history of the southwestern sector of the Greenland Ice Sheet

The offshore and coastal geomorphology of southwest Greenland records evidence for the advance and decay of the Greenland Ice Sheet during the Last Glacial Maximum. Regional ice flow patterns in the vicinity of Sisimiut show an enlarged ice sheet that extended southwestwards on to the shelf, with an ice stream centred over Holsteinsborg dyb. High level periglacial terrain composed of blockfield and tors is dated to between 101 and 142 ka using ²⁶Al and ¹⁰Be cosmogenic exposure ages. These limit the maximum surface elevation of the Last Glacial Maximum ice sheet in this part of southwest Greenland to ca 750–810 m asl, and demonstrate that terrain above this level has been ice free since MIS 6. Last Glacial Maximum ice thickness on the coast of ca 700 m implies that the ice sheet reached the mid to outer continental shelf edge to form the Outer Hellefisk moraines. Exposure dates record ice surface thinning from 21.0 to 9.8 ka, with downwasting rates varying from 0.06 to 0.12 m yr⁻¹. This reflects strong surface ablation associated with increased air temperatures running up to the Bølling Interstadial (GIS1e) at ca 14 ka, and later marine calving under high sea levels. The relatively late retreat of the Itilleq ice stream inland of the present coastline is similar to the pattern observed at Jakobshavn Isbræ, located 250 km north in Disko Bugt, which also retreated from the continental shelf after ca 10 ka. We hypothesise that the ice streams of West Greenland persisted on the inner shelf until the early Holocene because of their considerable ice thickness and greater ice discharge compared with the adjacent ice sheet.     

Sulfur dioxide in the tropical marine boundary layer: dry deposition and heterogeneous oxidation observed during the Pacific Atmospheric Sulfur Experiment

Research flights with the National Center for Atmospheric Research (NCAR) C-130 airborne laboratory were conducted over the equatorial ocean during the Pacific Atmospheric Sulfur Experiment (PASE). The focused, repetitive flight plans provided a unique opportunity to explore the principal pathways of sulfur processing in remote marine environments in close detail. Fast airborne measurements of SO₂ using the Drexel University APIMS (Atmospheric Pressure Ionization Mass Spectrometer) instrument further provided unprecedented insight into the complete budget of this important sulfur gas. In general, turbulent mixing in the marine boundary layer (MBL) continuously depletes SO₂ due to the shallow convection of the tropical trade wind regime by venting the gas into the buffer layer (BuL) above. However, on nearly one-third of the flights a net import of SO₂ into the MBL from the BuL was observed. Concurrent measurements of the DMS budget allowed for a heterogeneous S(IV) oxidation rate to be inferred from the SO₂ budget residual. The average heterogeneous loss rate was found to be 0.05 h⁻¹, which taken in conjunction with the observed aerosol surface area distributions and O₃ levels indicates that the supermicron aerosols maintain a near neutral pH. The average dry deposition velocity of SO₂ was found to be 0.4 cm s⁻¹, about 30% lower than predicted by standard parameterizations. The yield of SO₂ from DMS oxidation was found to be near unity. The mission averages indicate that approximately 57% of the SO₂ in the MBL is lost to aerosols, 27% is subject to dry deposition, 7% is mixed into the BuL, and 10% is oxidized by OH.