Abrupt climate change and transient climates during the Paleogene: a marine perspective

Detailed investigations of high latitude sequences recently collected by the Ocean Drilling Program (ODP) indicate that periods of rapid climate change often culminated in brief transient climates, with more extreme conditions than subsequent long term climates. Two examples of such events have been identified in the Paleogene; the first in latest Paleocene time in the middle of a warming trend that began several million years earlier: the second in earliest Oligocene time near the end of a Middle Eocene to Late Oligocene global cooling trend. Superimposed on the earlier event was a sudden and extreme warming of both high latitude sea surface and deep ocean waters. Imbedded in the latter transition was an abrupt decline in high latitude temperatures and the brief appearance of a full size continental ice-sheet on Antarctica. In both cases the climate extremes were not stable, lasting for less than a few hundred thousand years, indicating a temporary or transient climate state. Geochemical and sedimentological evidence suggest that both Paleogene climate events were accompanied by reorganizations in ocean circulation, and major perturbations in marine productivity and the global carbon cycle. The Paleocene-Eocene thermal maximum was marked by reduced oceanic turnover and decreases in global delta 13C and in marine productivity, while the Early Oligocene glacial maximum was accompanied by intensification of deep ocean circulation and elevated delta 13C and productivity. It has been suggested that sudden changes in climate and/or ocean circulation might occur as a result of gradual forcing as certain physical thresholds are exceeded. We investigate the possibility that sudden reorganizations in ocean and/or atmosphere circulation during these abrupt transitions generated short-term positive feedbacks that briefly sustained these transient climatic states.     

A critical evaluation of multi‐proxy dendroclimatology in northern Finland

Twentieth‐century summer (July–August) temperatures in northern Finland are reconstructed using ring widths, maximum density and stable carbon isotope ratios (δ¹³C) of Scots pine tree rings, and using combinations of these proxies. Verification is based on the coefficient of determination (r²), reduction of error (RE) and coefficient of efficiency (CE) statistics. Of the individual proxies, δ¹³C performs best, followed by maximum density. Combining δ¹³C and maximum density strengthens the climate signal but adding ring widths leads to little improvement. Blue intensity, an inexpensive alternative to X‐ray densitometry, is shown to perform similarly. Multi‐proxy reconstruction of summer temperatures from a single site produces strong correlations with gridded climate data over most of northern Fennoscandia. Since relatively few trees are required (<15) the approach could be applied to long sub‐fossil chronologies where replication may be episodically low. Copyright © 2010 John Wiley & Sons, Ltd.     

Electron microprobe study of a ‘mysterite’-bearing inclusion from the Krymka LL-chondrite

The black inclusion from the Krymka LL3 chondrite previously found to contain ‘mysterite’ by Lewiset al. (1979) belongs to a hitherto unknown class of carbonaceous chondrites. Its olivine and pyroxene compositions. Fo 97–99 and En 96, respectively, are characteristic of carbonaceous chondrites and its plagioclase composition. An 100, is characteristic of C3's. It contains a peculiar group of Co-, Cr-rich metal grains whose compositions are similar, but not identical, to those in C2 chondrites and which also bear some similarities to those in Renazzo. Its weight ratios of total $\text{Fe}\text{SiO}{}_2$ and $\text{solSiO}{}_2\text{MgO}$ are 0.74 and 1.43, respectively, and its atomic ratio of $\text{Si}\text{Al}$ is 10.7, exactly the same as in carbonaceous chondrites. Its bulk chemical composition is very close to that of the Murchison C2 chondrite. The association of mysterite with a special type of carbonaceous chondrite material suggests that mysterite formed by low-temperature condensation in a different region of the nebula from other carbonaceous chondrites.     

The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Long-Term Trends of Nutrients and Phytoplankton in Chesapeake Bay

Climate effects on hydrology impart high variability to water-quality properties, including nutrient loadings, concentrations, and phytoplankton biomass as chlorophyll-a (chl-a), in estuarine and coastal ecosystems. Resolving long-term trends of these properties requires that we distinguish climate effects from secular changes reflecting anthropogenic eutrophication. Here, we test the hypothesis that strong climatic contrasts leading to irregular dry and wet periods contribute significantly to interannual variability of mean annual values of water-quality properties using in situ data for Chesapeake Bay. Climate effects are quantified using annual freshwater discharge from the Susquehanna River together with a synoptic climatology for the Chesapeake Bay region based on predominant sea-level pressure patterns. Time series of water-quality properties are analyzed using historical (1945–1983) and recent (1984–2012) data for the bay adjusted for climate effects on hydrology. Contemporary monitoring by the Chesapeake Bay Program (CBP) provides data for a period since mid-1984 that is significantly impacted by anthropogenic eutrophication, while historical data back to 1945 serve as historical context for a period prior to severe impairments. The generalized additive model (GAM) and the generalized additive mixed model (GAMM) are developed for nutrient loadings and concentrations (total nitrogen—TN, nitrate?+?nitrate—NO₂?+?NO₃) at the Susquehanna River and water-quality properties in the bay proper, including dissolved nutrients (NO₂?+?NO₃, orthophosphate—PO₄), chl-a, diffuse light attenuation coefficient (K _D (PAR)), and chl-a/TN. Each statistical model consists of a sum of nonlinear functions to generate flow-adjusted time series and compute long-term trends accounting for climate effects on hydrology. We present results identifying successive periods of (1) eutrophication ca. 1945–1980 characterized by approximately doubled TN and NO₂?+?NO₃ loadings, leading to increased chl-a and associated ecosystem impairments, and (2) modest decreases of TN and NO₂?+?NO₃ loadings from 1981 to 2012, signaling a partial reversal of nutrient over-enrichment. Comparison of our findings with long-term trends of water-quality properties for a variety of estuarine and coastal ecosystems around the world reveals that trends for Chesapeake Bay are weaker than for other systems subject to strenuous management efforts, suggesting that more aggressive actions than those undertaken to date will be required to counter anthropogenic eutrophication of this valuable resource.     

H2S fluxes from Mt. Etna, Stromboli, and Vulcano (Italy) and implications for the sulfur budget at volcanoes

We present here new measurements of sulfur dioxide and hydrogen sulfide emissions from Vulcano, Etna, and Stromboli (Italy), made by direct sampling at vents and by filter pack and ultraviolet spectroscopy in downwind plumes. Measurements at the F0 and FA fumaroles on Vulcano yielded SO₂/H₂S molar ratios of ≈0.38 and ≈1.4, respectively, from which we estimate an H₂S flux of 6 to 9 t · d⁻¹ for the summit crater. For Mt. Etna and Stromboli, we found SO₂/H₂S molar ratios of ≈20 and ≈15, respectively, which combined with SO₂ flux measurements, suggest H₂S emission rates of 50 to 113 t · d⁻¹ and 4 to 8 t · d⁻¹, respectively. We observe that “source” and plume SO₂/H₂S ratios at Vulcano are similar, suggesting that hydrogen sulfide is essentially inert on timescales of seconds to minutes. This finding has important implications for estimates of volcanic total sulfur budget at volcanoes since most existing measurements do not account for H₂S emission.     

The chemical compositions and surface structures of two roAp stars with lithium spots (HD 83368 and HD 60435)

The spectra of two roAp stars have been analyzed as part of a project to study lithium in magnetic Ap stars. Variability of the Li I 6708 Å resonance doublet and rare-earth lines was detected, which can be explained using an oblique rotator model with the lithium spots located at the magnetic poles. Synthetic spectra obtained at different rotational phases have yielded the first data on the atmospheric chemical compositions of these spotted stars. Using refined atomic data and the most complete line lists has enabled a detailed study of the spectra near the Li I 6708 Å line and computation of the Li I line profile taking into account the spotted distribution of the lithium over the stellar surfaces. The positions of two lithium spots and lithium abundances for each of the spots have been determined.     

Geological modelling of mineral deposits for prediction in mining

Accurate prediction of the shape, location, size and properties of the solid rock materials to be extracted during mining is essential for reliable technical and financial planning. This is achieved through geological modelling of the three-dimensional (3D) shape and properties of the materials present in mineral deposits, and the presentation of results in a form which is accessible to mine planning engineers. In recent years the application of interactive graphics software, offering 3D database handling, modelling and visualisation, has greatly enhanced the options available for predicting the subsurface limits and characteristics of mineral deposits. A review of conventional 3D geological interpretation methods, and the model struc- tures and modelling methods used in reserve estimation and mine planning software packages, illustrates the importance of such approaches in the modern mining industry. Despite the widespread introduction and acceptance of computer hardware and software in mining applications, in recent years, there has been little fundamental change in the way in which geology is used in orebody modelling for predictive purposes. Selected areas of current research, aimed at tackling issues such as the use of orientation data, quantification of morphological differences, incorporation of geological age relationships, multi-resolution models and the application of virtual reality hardware and software, are discussed.