Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

Nutrient pore water chemistry,Great Bay,New Hampshire: Benthic fluxes

Benthic fluxes of C, N, P and Si have been measured at two sites in Great Bay Estuary, New Hampshire. Higher fluxes of reactive phosphate, nitrate and reactive silicate were observed at the site where bioturbation is known to occur and the fluxes of NH₄ ⁺, PO₄ ^?3 and reactive silicate at this location were from 3 to 6 times higher than that calculated by simple pore water diffusion models. *** DIRECT SUPPORT *** A01BY019 00010     

Observational estimation of heat budgets on drifting ice and open water over the Arctic Ocean

It is of major scientific interests to determine the parameters of momentum, heat and vapor exchange in the planetary boundary layer in order to study the effects of ocean-ice-atmosphere interactions and their feedback mechanisms on global climate[1]. Lin…     

The geometry and petrogenesis of hydrothermal dolomites at Navan, Ireland

The dolomites at Navan, Ireland, formed in Courceyan peritidal and shallow-shelf limestones. The dolomite body has a plume-like geometry, cross-cutting both lithological boundaries and diagenetic barriers generated by sea-floor cementation and emergence. The dolomitizing fluids rose parallel to major faults to diffuse laterally through the succession, controlled by variations in permeability that reflect both facies variation on deposition and pre-dolomitization diagenesis. Cathodoluminescent zones reveal three principal stages of dolomite emplacement, separated by dissolution surfaces, with each stage reflecting several changes in the character of dolomitizing solutions. The predominance of dull zones indicates burial rather than surface conditions. The dolomites formed some time after burial in response to an areally limited hydrothermally-driven flow. Isotope values (σ¹⁸O of — 6σ6 to — 10.4%_δ and σ¹³C of — 0σ2 to +2σ5%_δ PDB), and fluid inclusion data, suggest that these fluids had compositions similar to those of Carboniferous seawater. However, these became hotter with time, with temperature increasing from 60 to 160δC. The Navan dolomites are closely associated with Europe's largest zinc-lead deposit. The distribution of the ores follows the same trend as that of the dolomites and paragenetic relationships indicate that dolomitization and mineralization were temporally and genetically related.     

School vulnerability to disaster: examination of school closure,demographic, and exposure factors in Hurricane Ike’s wind swath

Damage and destruction to schools from climate-related disasters can have significant and lasting impacts on curriculum and educational programs, educational attainment, and future income-earning potential of affected students. As such, assessing the potential impact of hazards is crucial to the ability of individuals, households, and communities to respond to natural disasters, extreme events, and economic crises. Yet, few studies have focused on assessing the vulnerability of schools in coastal regions of the USA. Using Hurricane Ike’s tropical storm wind swath in the State of Texas as our study area, we: (1) assessed the spatial distribution patterns of school closures and (2) tested the relationship between school closure and vulnerability factors (namely physical exposure and school demographics) using zero-inflated negative binomial regression models. The regression results show that higher probabilities of hurricane strikes, more urbanized school districts, and school districts located in coastal counties on the right side of Ike’s path have significant positive associations with an increase in the number of school closure days. Socioeconomic characteristics were not significantly associated with the number of days closed, with the exception of proportion of Hispanic youth in schools, a result which is not supported by the social vulnerability literature. At a practical level, understanding how hurricanes may adversely impact schools is important for developing appropriate preparedness, mitigation, recovery, and adaptation strategies. For example, school districts on the right side of the hurricane track can plan in advance for potential damage and destruction. The ability of a community to respond to future natural disasters, extreme events, and economic crises depends in part on mitigating these adverse effects.     

Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.