The anthracite smokers of eastern pennsylvania: PS2 (g)P_{S_2 (g)} -T stability diagram by TL analysis

Coal seams and culm banks associated with mine fires in the anthracite region of eastern Pennsylvania have been burning for decades. Many of the fires may have ignited by spontaneous combustion or by the burning of trash. Minerals associated with the combustion of anthracite form by the condensation of gas exhaled through surficial gas vents or anthracite smokers. A Pressure-Temperature (P-T) stability diagram is constructed for the condensation of orthorhombic sulfur from anthracite gas using Thermodynamic Loop Analysis (TL analysis). This method of analyzing chemical systems incorporates Kirchhoff's Law into a four step procedure structured around a closed thermodynamic cycle or thermodynamic loop. The four steps, referred to us The Four S S of Thermodynamic Loop Analysis, include: (1) Set Up—graphical characterization of the problem. (2) Sum—the application of thermodynamic principles. (3) Substitute—the use of materials data available from the literature, and (4) Solve—computation of one or more variables. The example presented demonstrates that thermodynamic loops can incorporate any number of polymorphic phase transformations. In addition, thermodynamic loop analysis is applicable to any geologic process involving the condensation of minerals from a gas. The stability diagram derived by TL analysis may have applicability in monitoring the release of sulfur gas into the atmosphere.     

Late Pleistocene sea-level oscillations (MIS 10–2) recorded in shallow marine and coastal plain sediments of the southern Wanganui Basin, New Zealand

The northern Wanganui Basin, New Zealand, is one of the key global sites for understanding marine cyclic sedimentation during the Quaternary. This paper presents the first evidence of marine cyclic sedimentation from its central-southern parts. Sedimentological, micropalaeontological and palynological analyses on a 280-m-deep borehole encountered units dating back to MIS 10. The sequence includes four marine cycles spanning MIS 9–5, which are overlain by terrestrial fluvial aggradation surfaces dating from MIS 4–2. Each marine unit represents a progressively shallowing depositional environment from the mid-shelf to coastal plain. This is overlain by a terrestrial sequence of lowstand fluvial terraces. Localized fault movements appear to have influenced the sedimentary character of the sequence during MIS 7a and 5e producing basement highs which provided protection to the shoreline. The cyclothems described in this paper now extend the already extensive, previously described record from MIS 17–10 to produce a combined eustatic record of Quaternary sea level change within the basin to MIS 5. They also provide an excellent example of the sedimentary response of a coastal basin to a progressive loss of sedimentation accommodation space.     

Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

Convolute laminations — a theoretical analysis: example of a Pennsylvanian sandstone

Data from an outcropping laminated interval were collected and analyzed to test the applicability of a theoretical model describing instability of layered systems. Rayleigh—Taylor wave perturbations result at the interface between fluids of contrasting density, viscosity, and thickness. In the special case where reverse density and viscosity interlaminations are developed, the deformation response produces a single wave with predictable amplitudes, wavelengths, and amplification rates.Physical measurements from both the outcropping section and modern sediments suggest the usefulness of the model for the interpretation of convolute laminations. Internal characteristics of the stratigraphic interval, and the developmental sequence of convoluted beds, are used to document the developmental history of these structures.     

Accretionary rims on inclusions in the Allende meteorite

Many inclusions in Allende, particularly those with irregular shapes, are surrounded by a sequence of thin layers which differ from one another in texture, mineralogy and mineral-chemistry. The layer underlying all others contains either: IA, pyroxene needles + olivine + clumps of hedenbergite and andradite; IB, olivine doughnuts; or IC, rectangular olivine crystals. The next layer outward, II, contains tiny (<5 μm) olivine plates and Layer III large (5–10 μm) olivine laths. The final layer, IV, occurs as clumps of andradite + hedenbergite surrounded by magnesium-rich pyroxene needles. It separates Layer III from the Allende matrix which is more poorly sorted and more sulfide-rich than Layer III. Nepheline and iron sulfide are common constituents of most layers, the latter being particularly fine-grained and abundant in Layer II. Although not every layer is present on every inclusion, the sequence of layers is constant. Evidence that the rims are accretionary aggregates includes the presence of highly disequilibrium mineral assemblages and the fact that they are highly porous masses consisting of many euhedral crystals with few intergrowths. In addition, the layers are thickest in topographic hollows on the surfaces of inclusions and the inner layers are absent or discontinuous beyond such irregularities, suggesting that the probability of accretion of crystals was low initially, except in pockets, and became greater later, after a soft cushion of accreted condensate crystals had already formed. Separation of assemblages of different mineralogy, mineral-chemistry and texture into different rim layers seems best explained by nebular models in which long, slow cooling histories allow differentiation during condensation by grain/gas separation processes.     

High-frequency radar observations of the June 2013 US East Coast meteotsunami

We report here on the observation and offline detection of the meteotsunami off the New Jersey coast on June 13, 2013, using coastal radar systems and tide gauges. This work extends the previous observations of tsunamis originating in Japan and Indonesia. The radars observed the meteotsunami 23 km offshore, 47 min before it arrived at the coast. Subsequent observations showed it moving onshore. The neighboring tide gauge height reading provides confirmation of the radar observations near the shore.     

The importance of record length in estimating the magnitude of climatic changes: an example using 175 years of lake ice-out dates in New England

Many studies have shown that lake ice-out (break-up) dates in the Northern Hemisphere are useful indicators of late winter/early spring climate change. Trends in lake ice-out dates in New England, USA, were analyzed for 25, 50, 75, 100, 125, 150, and 175 year periods ending in 2008. More than 100 years of ice-out data were available for 19 of the 28 lakes in this study. The magnitude of trends over time depends on the length of the period considered. For the recent 25-year period, there was a mix of earlier and later ice-out dates. Lake ice-outs during the last 50 years became earlier by 1.8 days/decade (median change for all lakes with adequate data). This is a much higher rate than for longer historical periods; ice-outs became earlier by 0.6 days/decade during the last 75 years, 0.4 days/decade during the last 100 years, and 0.6 days/decade during the last 125 years. The significance of trends was assessed under the assumption of serial independence of historical ice-out dates and under the assumption of short and long term persistence. Hypolimnion dissolved oxygen (DO) levels are an important factor in lake eutrophication and coldwater fish survival. Based on historical data available at three lakes, 32 to 46 % of the interannual variability of late summer hypolimnion DO levels was related to ice-out dates; earlier ice-outs were associated with lower DO levels.     

A theoretical investigation of bed-form shapes

The deformation of movable boundaries under the action of an applied turbulent shear stress is well known. The resulting bed forms often are highly organized and nearly two-dimensional, which makes them an intriguing focus of study considering that they are generated in both steady and oscillatory turbulent flows. Many past studies share a common approach in which an infinitesimal perturbation is prescribed and the resulting growth or decay patterns are examined. In this approach, the bed forms are usually sinusoidal and the perturbation analysis does not provide a theoretical prediction of equilibrium bed-form geometry. An alternative approach is suggested here in which the forcing terms (pressure and stress) are prescribed parametrically and the governing equations are solved for the flow velocity and the associated boundary deformation. Using a multilayered approach, in which the bottom boundary layer is divided into a discrete, yet, arbitrary number of finite layers, analytical solutions for the horizontal current and bed profile are derived. The derivations identify two nondimensional parameters, p₀/u₀² and ₀/kh₀u₀², which modulate the amplitude of the velocity fluctuations and boundary deformation. For the case of combined pressure and stress divergence anomalies, the magnitude of the front face and lee slopes exhibit an asymmetry that is consistent with observed bed forms in steady two-dimensional flows.Responsible Editor: Jens Kappenberg