Latest Pleistocene and Holocene glaciation of Baffin Island, Arctic Canada: key patterns and chronologies

Melting glaciers and ice caps on Baffin Island contribute roughly half of the sea-level rise from all ice in Arctic Canada, although they comprise only one-fourth of the total ice in the region. The uncertain future response of arctic glaciers and ice caps to climate change motivates the use of paleodata to evaluate the sensitivity of glaciers to past warm intervals and to constrain mechanisms that drive glacier change. We review the key patterns and chronologies of latest Pleistocene and Holocene glaciation on Baffin Island. The deglaciation by the Laurentide Ice Sheet occurred generally slowly and steadily throughout the Holocene to its present margin (Barnes Ice Cap) except for two periods of rapid retreat: An early interval 12 to 10 ka when outlet glaciers retreated rapidly through deep fiords and sounds, and a later interval 7 ka when ice over Foxe Basin collapsed. In coastal settings, alpine glaciers were smaller during the Younger Dryas period than during the Little Ice Age. At least some alpine glaciers apparently survived the early Holocene thermal maximum, which was several degrees warmer than today, although data on glacier extent during the early Holocene is extremely sparse. Following the early Holocene thermal maximum, glaciers advanced during Neoglaciation, beginning in some places as early as 6 ka, although most sites do not record near-Little Ice Age positions until 3.5 to 2.5 ka. Alpine glaciers reached their largest Holocene extents during the Little Ice Age, when temperatures were 1–1.5 °C cooler than during the late 20th century. Synchronous advances across Baffin Island throughout Neoglaciation indicate sub-Milankovitch controls on glaciation that could involve major volcanic eruptions and solar variability. Future work should further elucidate the state of glaciers and ice caps during the early Holocene thermal maximum and glacier response to climate forcing mechanisms.     

Radiation-induced defects in euclase: formation of O<Superscript>−</Superscript> hole and Ti<Superscript>3+</Superscript> electron centers

Irradiation techniques are often applied to gem minerals for color enhancement purposes. Natural green, blue and colorless specimens of rare gemological quality euclase, BeAlSiO₄(OH), from Brazil were irradiated with gamma rays in the dose range from 10 to 500 kGy. Although the colors of the different specimens were not strongly influenced, two different irradiation-induced paramagnetic defect centers were found by electron paramagnetic resonance (EPR). The first one is an O⁻ hole center interacting with one Al neighbor and the second is a Ti³⁺ electron center. The EPR angular rotation patterns of both irradiation-induced defects were measured and analyzed. The results suggest that O⁻ hole centers are formed by dissociation of the hydroxyl ions, similar as in topaz crystals. In euclase the OH⁻ ions interconnect distorted Al octahedra and Be tetrahedra in O₅ positions. During irradiation, the electrons are captured by titanium ions (Ti⁴⁺ + e⁻), leading to the formation of paramagnetic Ti³⁺ ions. From the EPR rotation patterns it is clear that these ions substitute for Al ions. The spin Hamiltonian parameters of the irradiation-induced defects are analyzed and compared to similar defect centers in other mineral specimens. Thermal annealing experiments show that the O⁻ hole centers and Ti³⁺ electron centers are directly connected through the radiation process.     

Provenance of the late Proterozoic to early Cambrian metaclastic sediments of the Sierra de San Luis (Eastern Sierras Pampeanas) and Cordillera Oriental, Argentina

Provenance studies have been performed utilising major and trace elements, Nd systematics, whole rock Pb–Pb isotopes and zircon U/Pb SHRIMP data on metasedimentary rocks of the Sierra de San Luis (Nogolí Metamorphic Complex, Pringles Metamorphic Complex, Conlara Metamorphic Complex and San Luis Formation) and the Puncoviscana Formation of the Cordillera Oriental. The goal was the characterisation of the different domains in the study area and to give insights to the location of the source rocks. An active continental margin setting with typical composition of the upper continental crust is depicted for all the complexes using major and trace elements. The Pringles Metamorphic Complex shows indications for crustal recycling, pointing to a bimodal provenance. Major volcanic input has to be rejected due to Th/Sc, Y/Ni and Cr/V ratios for all units. The εNd_{(540 Ma)} data is lower for the San Luis Formation and higher for the Conlara Metamorphic Complex, as compared to the other units, in which a good consistency is given. This is similar to the T_DM ages, where the metapsammitic samples of the San Luis Formation are slightly older. The spread of data is largest for the Pringles Metamorphic Complex, again implying two different sources. The whole rock ²⁰⁷Pb/²⁰⁶Pb isotopic data lies in between the South American and African sources, excluding Laurentian provenances. The whole rock Pb–Pb data is almost indistinguishable in the different investigated domains. Only the PMC shows slightly elevated ²⁰⁸Pb/²⁰⁴Pb values. Possible source rocks for the different domains could be the Quebrada Choja in the Central Arequipa–Antofalla domain, the Southern domain of the Arequipa–Antofalla basement, the Brazilian shield or southern Africa. Zircon SHRIMP data point to a connection between the Puncoviscana Formation and the Conlara Metamorphic Complex. Two maxima around 600 Ma and around 1000 Ma have been determined. The Nogolí Metamorphic Complex and the Pringles Metamorphic Complex show one peak of detrital zircons around 550 Ma, and only a few grains are older than 700 Ma. The detrital zircon ages for the San Luis Formation show age ranges between 590 and 550 Ma. A common basin can be assumed for the Conlara Metamorphic Complex and the Puncoviscana Formation, but the available data support different sources for the rest of the Complexes of the Sierra de San Luis. These share the diminished importance or the lack of the Grenvillian detrital peak, a common feature for the late Cambrian–early Ordovician basins of the Eastern Sierras Pampeanas, in contrast to the Sierras de Córdoba, the PVF and the Conlara Metamorphic Complex.     

Seismic Reflection Characteristics of Naturally-Induced Subsidence Affecting Transportation

High-resolution seismic reflections have been used effectively to investigate sinkholes formed from the dissolution of a bedded salt unit found throughout most of Central Kansas. Surface subsidence can have devastating effects on transportation structures. Roads, rails, bridges, and pipeliues can even be dramatically affected by minor ground instability. Areas susceptible to surface subsidence can put public safety at risk. Subsurface expressions significantly larger than surface depressions are consistently observed on seismic images recorded over sinkholes in Kansas. Until subsidence reaches the ground surface, failure appears to be controlled by compressional forces evidenced by faults with reverse orientation. Once a surface depression forms or dissolution of the salt slows or stops, subsidence structures are consistent with a tensional stress environment with prevalent normal faults. Detecting areas of rapid subsidence potential, prior to surface failure, is the ultimate goal of any geotechnical survey where the ground surface is susceptible to settling. Seismic reflection images have helped correlate active subsidence to dormant paleofeatures, project horizontal growth of active sinkholes based on subsurface structures, and appraise the risk of catastrophic failure.     

Quantitative Integration of High-Resolution Hydrogeophysical Data: A Novel Approach to Monte-Carlo-Type Conditional Stochastic Simulations and Implications for Hydrological Predictions

Geophysical techniques can help to bridge the inherent gap that exists with regard to spatial resolution and coverage for classical hydrological methods. This has led to the emergence of a new and rapidly growing research domain generally referred to as hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters, their inherent trade-off between resolution and range, as well as the notoriously site-specific nature of petrophysical parameter relations, the fundamental usefulness of multi-method surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative Integration of the resulting vast and diverse database into a unified model of the probed subsurface region that is consistent with all available measurements. To this end, we present a novel approach toward hydrogeophysical data integration based on a Monte-Carlo-type conditional stochastic simulation method that we consider to be particularly suitable for high-resolution local-scale studies. Monte Carlo techniques are flexible and versatile, allowing for accounting for a wide variety of data and constraints of differing resolution and hardness, and thus have the potential of providing, in a geostatistical sense, realistic models of the pertinent target parameter distributions. Compared to more conventional approaches, such as co-kriging or cluster analysis, our approach provides significant advancements in the way that larger-scale structural information contained in the hydrogeophysical data can be accounted for. After outlining the methodological background of our algorithm, we present the results of its application to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the detailed local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to a field dataset collected at the Boise Hydrogeophysical Research Site. Finally, we compare the performance of our data integration approach to that of more conventional methods with regard to the prediction of flow and transport phenomena in highly heterogeneous media and discuss the Implications arising.     

U-Pb SHRIMP ages of granitoides from the main batholith belt (North East Asia)

U-Pb SHRIMP-dating of zircons from twenty five intrusions representative of the Main granitoid batholith belt and associated dike swarms (Yano-Kolyma gold bearing province, North East of Asia) are mostly ∼150 ± 3 Ma (Kimmeridgian-Tithonian). Two less widespread impulses of magmatism dated at 160–155 Ma and 146–143 Ma representing the full range of ages present in the Main belt. Paleoproterozoic (∼1.8 Ga) inherited zircons were found in three intrusions from the south-western part of the belt where Precambrian crust of the North-Asia craton is inferred to underlie it. Published in Russian in Doklady Akademii Nauk, 2009, Vol. 426, No. 2, pp. 216–221. The article was translated by the authors.     

Estuarine Ecosystem Response Captured Using a Synoptic Climatology

Estuarine and coastal ecosystems respond strongly to proximate climate forcing. In this study, we present a regional, synoptic climatology as an approach to classify weather patterns that generate interannual variability in coastal and estuarine ecosystems. Synoptic climatology is a method that classifies sea level pressure data into distinct patterns representing common weather features for a specified region. A synoptic climatology was developed for the eastern United States and used to quantify surface conditions affecting Chesapeake Bay during wet and dry years. In a synthesis analysis, several mechanisms were identified that explained the link between weather patterns and ecosystem structure, principal among them is the delivery of freshwater to the Bay during spring. Wet and dry years were characterized by shifts in biogeography of the Chesapeake Bay. The shifts resulted from habitat changes and trophic interactions and included the timing and magnitude of the spring phytoplankton bloom, the distribution/abundance of mesozooplankton and gelatinous zooplankton, and juvenile indices of fish. Synoptic climatology resolved regional weather variability at a spatial scale not strongly controlled by larger-scale climate indices and explained ecosystem responses in Chesapeake Bay.     

A dynamical model with a new inversion technique applied to observations of Comet C/2000 WM1 (LINEAR)

Three continuum images of Comet C/2000 WM₁ (LINEAR) obtained on Nov 10, Nov 19, and Dec 03, 2001, are analyzed with the aid of a dynamical model, i.e. with a model that uses the size-dependent motion of dust grains under solar radiation pressure to determine the dust size distribution and its temporal change. The frames are photometrically calibrated in terms of the albedo filling factor product. On Nov 20.2 the Earth transited the orbital plane of the comet and an anti-tail was recognized in the image of Nov 19. For the determination of the particle fluxes describing the contribution of monodisperse particle shells to the cometary brightness the model uses a new regularization method employing Chebyshev polynomials of selected orders in emission time and particle size. It guarantees positiveness of the particle fluxes and imposes a varying degree of smoothness on their dependence on particle size and emission time. The particle emission velocities are still derived by trial and error. The dynamical model is described in detail. Results are presented for several low orders of the Chebyshev polynomials and are compared in order to understand the limitations imposed by the regularization process. The size distributions derived from the different observations do not always agree. This is particularly true for the earliest and most recent synchrones contributing to an image. In the observations of Nov 10 and Dec 03, i.e. excluding the anti-tail image, the integrated mass loss strongly decreases in the most recent time steps of the model although the comet is still approaching the Sun. This is interpreted as an artifact introduced by the overlap of the shells of large particle size emitted shortly before the observation. The model derives an increasing number of small particles released by the comet in the second half of November. This is at least in part considered as real and attributed to particle fragmentation occurring when the comet was at a heliocentric distance of about 1.4 AU.