Application of a vanishing, quasi-sigma, vertical coordinate for simulation of high-speed, deep currents over the Sigsbee Escarpment in the Gulf of Mexico

Recent observations over the Sigsbee Escarpment in the Gulf of Mexico have revealed extremely energetic deep currents (near 1 m s⁻¹), which are trapped along the escarpment. Both scientific interest and engineering needs demand dynamical understanding of these extreme events, and can benefit from a numerical model designed to complement observational and theoretical investigations in this region of complicated topography. The primary objective of this study is to develop a modeling methodology capable of simulating these physical processes and apply the model to the Sigsbee Escarpment region. The very steep slope of the Sigsbee Escarpment (0.05–0.1) limits the application of ocean models with traditional terrain-following (sigma) vertical coordinates, which may represent the very complicated topography in the region adequately, can result in large truncation errors during calculation of the horizontal pressure gradient. A new vertical coordinate system, termed a vanishing quasi-sigma coordinate, is implemented in the Navy Coastal Ocean Model for application to the Sigsbee Escarpment region. Vertical coordinate surfaces for this grid have noticeably gentler slopes than a traditional sigma grid, while still following the terrain near the ocean bottom. The new vertical grid is tested with a suite of numerical experiments and compared to a classical sigma-layer model. The numerical error is substantially reduced in the model with the new vertical grid. A one-year, realistic, numerical simulation is performed to simulate strong, deep currents over the Escarpment using a very-high-resolution nested modeling approach. The model results are analyzed to demonstrate that the deep-ocean currents in the simulation replicate the prominent dynamical features of the observed intense currents in the region.     

Surveying Upstate NY Well Water for Pesticide Contamination: Cayuga and Orange Counties

Private wells in Cayuga and Orange counties in New York were sampled to determine the occurrence of pesticide contamination of groundwater in areas where significant pesticide use coincides with shallow or otherwise vulnerable groundwater. Well selection was based on local groundwater knowledge, risk modeling, aerial photo assessments, and pesticide application database mapping. Single timepoint samples from 40 wells in each county were subjected to 93‐compound chromatographic scans. All samples were nondetects (reporting limits ≤1 μg/L), thus no wells from either county exceeded any of 15 state groundwater standards or guidance values. More sensitive enzyme‐linked immunosorbent assays (ELISA) found two wells with quantifiable atrazine in each county (0.1–0.3 μg/L), one well with quantifiable diazinon (0.1 μg/L) in Orange County, and one well with quantifiable alachlor (0.2 μg/L) in Cayuga County. Trace detections (<0.1 μg/L) in Cayuga County included atrazine (five wells), metolachlor (six wells), and alachlor (one well), including three wells with multiple detections. All 12 Cayuga County wells with ELISA detections had either corn/grain or corn/forage rotations as primary surrounding land uses (although 20 other wells with the same land uses had no detections) and all quantified detections and most trace detections occurred in wells up to 9‐m deep. Orange County trace (<0.1 μg/L) ELISA detections (atrazine three wells, diazinon one well, and metolachlor five wells) and quantified detections were only generally associated with agricultural land uses. Finding acceptable drinking water quality in areas of vulnerable groundwater suggests that water quality in less vulnerable areas will also be good.     

The peterson limestone — early cretaceous lacustrine carbonate deposition in western wyoming and southeastern idaho

The lacustrine Peterson Limestone of western Wyoming and southeastern Idaho comprises six lithofacies throughout its 20,000 km² aerial extent. These are: (1) calcareous sandstone and shale, (2) red nodular limestone, (3) pink sandy micrite, (4) biomicrite, (5) graded silty micrite, and (6) limestone conglomerate. The first two represent floodplain deposition and paleosols, whereas the remaining are shallow nearshore and deeper lacustrine sediments.This sequence was developed in a large fresh, hardwater lake surrounded by fluvial systems and associated flood plains in a warm temperate climate. Well-oxidized sandy terrigenous rocks, together with calcareous paleosol nodules, indicate that flood-plain deposition both preceded and was concurrent with lacustrine carbonate deposition. Micrite and biomicrite formed in deeper parts of the basin while sandy and silty carbonate accumulated in shallower lake-margin areas. Less-calcareous shale units which are interbedded with deeper-water carbonate were deposited either during rapid basin subsidence and deepening of the lake center or during periods of slower carbonate precipitation. Turbidity currents and subaqueous debris flows generated along steeper lake margins, resulted in the deposition of rhythmic layers of graded silty micrite and diamictic limestone conglomerate in the deepest part of the basin. The carbonate-rich sediments comprising these two lithofacies were originally deposited on shallow lake-margin benches and subsequently were transported downslope toward the lake center.Comparison with other carbonate-precipitating lacustrine systems indicates that this lake was not like modern playas. Although no known modern lacustrine system is precisely like Lake Peterson, the flora, fauna, composition, and distribution of facies within modern temperate-region lakes most closely resemble those of the Peterson Limestone.     

An application of spatial filtering techniques for land cover mapping using TM‐images

Land cover classification using satellite imagery is commonly based on spectral information in the individual pixels. The information in neighbouring pixels is ignored. Spatial filtering techniques using information present in neighbouring pixels may however, contribute significantly to an improvement of the classification. In this study different methods of spatial filtering are applied to a part of a TM‐scene of Kenya to assess their relative reliability. The study area is characterized by extended, relatively homogeneous areas of eucalyptus forests and tea estates and by fragmentated areas of agricultural land use. Spectral information was combined with the results of different spatial filtering methods and then classified. The spatial filtering techniques applied were texture calculation by means of variance, “median minus original” filtering and fractal dimension computations using several sizes of templates. The obtained classification accuracy of several image combinations is compared using the percentage correctly classified and using an overall accuracy measure: the Kappa coefficient. It is concluded that in this case the spatial filtering techniques only slightly improve the classification. From the applied filtering methods texture calculation by means of variance yielded the best results.     

A skeletal assemblage classification system for non-tropical carbonate deposits based on New Zealand Cenozoic limestones

Cenozoic limestones in New Zealand are mainly skeletal grainstones and packstones formed under non-tropical climatic conditions in open marine shelf or ramp environments. Following petrographic analysis of the nature and abundance of the skeletal components in nearly 500 samples of these limestones, a complete linkage cluster analysis identified seven major skeletal assemblages that may be regarded as subdivisions of the single foramol skeletal association defined by Lees and Buller (1972) for temperate-region carbonate deposits. The seven assemblages are given contracted names, as follows: (a) BARNAMOL = barnacle/bivalve-dominated; (b) BIMOL = bivalve-dominated; (c) BRYOMOL = bryozoan/bivalve-dominated; (d) ECHINOFOR = echinoderm/benthic foraminiferal-dominated; (e) NANNOFOR = nannofossil/planktonic foraminiferal-dominated; (f) RHODALGAL = calcareous red algal-dominated; and (g) RHODECHFOR = calcareous red algal/echinoderm/benthic foraminiferal-dominated. A composite triangular classification diagram has been devised for naming the skeletal assemblage of an unknown sample on the basis of its three main skeletal components. The diagram successfully characterises more than 85% of the New Zealand Cenozoic limestone samples and also appears to be applicable for the skeletal assemblage designation of many overseas examples of non-tropical carbonate deposits. Limitations relate mainly to locally common skeletal types (e.g. serpulids, brachiopods) that are presently not incorporated into the New Zealand-based scheme. The general ecological preferences of the main skeletal contributors in each of the seven skeletal assemblages form a basis for relating the assemblages to broad shelf habitats. Consequently, as well as the benefits of providing a more consistent skeletal assemblage terminology for comparative studies between different workers, the scheme can assist with the paleoenvironmental interpretation of non-tropical skeletal carbonate facies.     

Formation of refractory inclusions by evaporation of condensate precursors

Berman’s (1983) activity-composition model for CaO-MgO-Al₂O₃-SiO₂ liquids is used to calculate the change in bulk chemical and isotopic composition during simultaneous cooling, evaporation, and crystallization of droplets having the compositions of reasonable condensate precursors of Types A and B refractory inclusions in CV3 chondrites. The degree of evaporation of MgO and SiO₂, calculated to be faithfully recorded in chemical and isotopic zoning of individual melilite crystals, is directly proportional to evaporation rate, which is a sensitive function of P_H2, and inversely proportional to the droplet radius and cooling rate. When the precursors are partially melted in pure hydrogen at peak temperatures in the vicinity of the initial crystallization temperature of melilite, their bulk chemical compositions evolve into the composition fields of refractory inclusions, mass-fractionated isotopic compositions of Mg, Si, and O are produced that are in the range of the isotopic compositions of natural inclusions, and melilite zoning profiles result that are similar to those observed in real inclusions. For droplets of radius 0.25 cm evaporating at P_H2 = 10⁻⁶ bar, precursors containing 8 to 13 wt.% MgO and 20 to 23% SiO₂ evolve into objects similar to compact Type A inclusions at cooling rates of 2 to 12 K/h, depending on the precise starting composition. Precursors containing 13 to 14 wt.% MgO and 23 to 26% SiO₂ evolve into objects with the characteristics of Type B1 inclusions at cooling rates of 1.5 to 3 K/h. The relatively SiO₂-poor members of the Type B2 group can be produced from precursors containing 14 to 16 wt.% MgO and 27 to 33% SiO₂ at cooling rates of <1 K/h. Type B2’s containing 27 to 35 wt.% SiO₂ and <12% MgO require precursors with higher SiO₂/MgO ratios at MgO > 15% than are found on any condensation curve. The characteristics of fluffy Type A inclusions, including their reversely zoned melilite, can only be understood in the context of this model if they contain relict melilite.     

Friction characteristics of Cd-rich carbonate films on calcite surfaces: implications for compositional differentiation at the nanometer scale

Lateral Force Microscopy (LFM) studies were carried out on cleaved calcite sections in contact with solutions supersaturated with respect to otavite (CdCO₃) or calcite-otavite solid solutions (SS) as a means to examine the potential for future application of LFM as a nanometer-scale mineral surface composition mapping technique. Layer-by-layer growth of surface films took place either by step advancement or by a surface nucleation and step advancement mechanisms. Friction vs. applied load data acquired on the films and the calcite substrate were successfully fitted to the Johnson Kendall Roberts (JKR) model for single asperity contacts. Following this model, friction differences between film and substrate at low loads were dictated by differences in adhesion, whereas at higher load they reflect differences in contact shear strength. In most experiments at fixed load, the film showed higher friction than the calcite surface, but the friction-load dependence for the different surfaces revealed that at low loads (0–40 nN), a calcian otavite film has lower friction than calcite; a result that is contrary to earlier LFM reports of the same system. Multilayer films of calcian-otavite displayed increasing friction with film thickness, consistent with the expectation that the film surface composition will become increasingly Cd-rich with increasing thickness. Both load- and thickness-dependence trends support the hypothesis that the contact shear strength correlates with the hydration enthalpy of the surface ions, thereby imparting friction sensitivity in the LFM to mineral-water interface composition.     

Behaviour of particle-laden flows into the ocean: experimental simulation and geological implications

The behaviour of subaerial particle-laden gravity currents (e.g. pyroclastic flows, lahars, debris flows, sediment-bearing floods and jökulhlaups) flowing into the sea has been simulated with analogue experiments. Flows of either saline solution, simple suspensions of silicon carbide (SiC) in water or complex suspensions of SiC and plastic particles in methanol were released down a slope into a tank of water. The excess momentum between subaerial and subaqueous flow is dissipated by a surface wave. At relatively low density contrasts between the tank water and the saline or simple suspensions, the flow mixture enters the water and forms a turbulent cloud involving extensive entrainment of water. The cloud then collapses gravitationally to form an underwater gravity current, which progresses along the tank floor. At higher density contrasts, the subaerial flow develops directly into a subaqueous flow. The flow slows and thickens in response to the reduced density contrast, which is driving motion, and then continues in the typical gravity current manner. Complex suspensions become dense flows along the tank floor or buoyant flows along the water surface, if the mixtures are sufficiently denser or lighter than water respectively. Flows of initially intermediate density are strongly influenced by the internal stratification of the subaerial flow. Material from the particulate-depleted upper sections of the subaerial flow becomes a buoyant gravity current along the water surface, whereas material from the particulate-enriched lower sections forms a dense flow along the tank floor. Sedimentation from the dense flow results in a reduction in bulk density until the mixture attains buoyancy, lifts off and becomes a secondary buoyant flow along the water surface. Jökulhlaups, lahars and debris flows are typically much denser than seawater and, thus, will usually form dense flows along the seabed. After sufficient sedimentation, the freshwater particulate mixture can lift off to form a buoyant flow at the sea surface, leading to a decoupling of the fine and coarse particles. Flood waters with low particulate concentrations (<2%) may form buoyant flows immediately upon entering the ocean. Subaerial pyroclastic flows develop a pronounced internal stratification during subaerial run-out and, thus, a flow-splitting behaviour is probable, which agrees with evidence for sea surface and underwater flows from historic eruptions of Krakatau and Mont Pelée. A pyroclastic flow with a bulk density closer to that of sea water may form a turbulent cloud, resulting in the deposition of much of the pyroclasts close to the shore. Dense subaqueous pyroclastic flows will eventually lift off and form secondary buoyant flows, either before or after the transformation to a water-supported nature.