Climate change impacts on U.S. Coastal and Marine Ecosystems

Increases in concentrations of greenhouse gases projected for the 21st century are expected to lead to increased mean global air and ocean temperatures. The National Assessment of Potential Consequences of Climate Variability and Change (NAST 2001) was based on a series of regional and sector assessments. This paper is a summary of the coastal and marine resources sector review of potential impacts on shorelines, estuaries, coastal wetlands, coral reefs, and ocean margin ecosystems. The assessment considered the impacts of several key drivers of climate change: sea level change; alterations in precipitation patterns and subsequent delivery of freshwater, nutrients, and sediment; increased ocean temperature; alterations in circulation patterns; changes in frequency and intensity of coastal storms; and increased levels of atmospheric CO₂. Increasing rates of sea-level rise and intensity and frequency of coastal storms and hurricanes over the next decades will increase threats to shorelines, wetlands, and coastal development. Estuarine productivity will change in response to alteration in the timing and amount of freshwater, nutrients, and sediment delivery. Higher water temperatures and changes in freshwater delivery will alter estuarine stratification, residence time, and eutrophication. Increased ocean temperatures are expected to increase coral bleaching and higher CO₂ levels may reduce coral calcification, making it more difficult for corals to recover from other disturbances, and inhibiting poleward shifts. Ocean warming is expected to cause poleward shifts in the ranges of many other organisms, including commercial species, and these shifts may have secondary effects on their predators and prey. Although these potential impacts of climate change and variability will vary from system to system, it is important to recognize that they will be superimposed upon, and in many cases intensify, other ecosystem stresses (pollution, harvesting, habitat destruction, invasive species, land and resource use, extreme natural events), which may lead to more significant consequences.     

A Comparison of Two-Phase Dynamic Upscaling Methods Based on Fluid Potentials

Although a number of methods for calculating dynamic pseudo-functions have been developed over the years, there is still a lack of understanding as to why a certain method will succeed in some cases but fail in others. In this paper, we describe the results of an assessment of several upscaling methods, namely the Kyte and Berry (KB) method, the Stone method, the Hewett and Archer (HA) method and the Transmissibility-Weighted (TW) method. We have analyzed the equations for deriving the methods and investigated the results of numerical simulations of gas displacing oil, in a variety of models to enable us to gain new insights into these, and related, upscaling methods. In particular, some novel observations on methods based on fluid potential are presented and the issue of using predicted fluid mobilities as a criterion of accuracy of an upscaling method is clarified.     

Inorganic speciation of dissolved elements in seawater: the influence of pH on concentration ratios

Assessments of inorganic elemental speciation in seawater span the past four decades. Experimentation, compilation and critical review of equilibrium data over the past forty years have, in particular, considerably improved our understanding of cation hydrolysis and the complexation of cations by carbonate ions in solution. Through experimental investigations and critical evaluation it is now known that more than forty elements have seawater speciation schemes that are strongly influenced by pH. In the present work, the speciation of the elements in seawater is summarized in a manner that highlights the significance of pH variations. For elements that have pH-dependent species concentration ratios, this work summarizes equilibrium data (S = 35, t = 25°C) that can be used to assess regions of dominance and relative species concentrations. Concentration ratios of complex species are expressed in the form log[A]/[B] = pH - C where brackets denote species concentrations in solution, A and B are species important at higher (A) and lower (B) solution pH, and C is a constant dependent on salinity, temperature and pressure. In the case of equilibria involving complex oxy-anions (MO _x (OH) _y ) or hydroxy complexes (M(OH) _n ), C is written as pK _n = -log K _n or pK _n * = -log K _n * respectively, where K _n and K _n * are equilibrium constants. For equilibria involving carbonate complexation, the constant C is written as pQ = -log(K ₂ ^l K _n [HCO₃ ^-]) where K ₂ ^l is the HCO₃ ^- dissociation constant, K _n is a cation complexation constant and [HCO₃ ^-] is approximated as 1.9 × 10^-3 molar. Equilibrium data expressed in this manner clearly show dominant species transitions, ranges of dominance, and relative concentrations at any pH.     

G-Probe-1 - An International Proficiency Test for Microprobe Laboratories - Report on Round 1 : February 2002 (TB-1 Basaltic Glass)

Results are presented for round one of a new international proficiency testing programme designed for microprobe laboratories involved in the routine analysis of silicate minerals. The sample used for this round was TB-1, a basaltic glass fused and prepared by the USGS. Thirty nine laboratories contributed data to this round, the majority of major element results being undertaken by EPMA and the majority of trace elements by LA-ICP-MS. Assigned values were derived from the median of results produced by nine selected laboratories that analysed powdered material by conventional ICP-MS, INAA and XRF techniques using bulk powders of the sample. Submitted microprobe results were evaluated using a target precision calculated using the Horwitz function, adopting the same criteria as those used for "applied" geochemistry laboratories in the companion GeoPT proficiency testing programme for laboratories involved in the routine bulk analysis of silicate rocks. An evaluation of results from participating microprobe laboratories indicated that overall, data were compatible with this precision function. A comparison between the performance of bulk and microprobe techniques used in the analysis of the basaltic glass showed remarkably good agreement, with significant bias only observed for the major oxide MgO.     

Freshwater inundation effects on emergent vegetation of a hypersaline salt marsh

The coastal marshlands of the Nueces estuary, Texas depend upon periodic freshwater inundation to support current community structure and promote further establishment and expansion of emergent halophytes. Decades of watershed modifications have dramatically decreased freshwater discharge into the upper estuary resulting in hypersaline and dry conditions. In an attempt to partially restore inflow, the U.S. Bureau of Reclamation excavated two overflow channels re-connecting the Nueces River to the marshlands. Freshwater-mediated (precipitation and inflow) changes in tidal creek and porewater salinity and emergent marsh vegetation were examined over a 5-yr period at three stations in the upper Nueces Marsh with the aid of a Geographical Information System (GIS). Two stations were potentially subjected to freshwater inflow through the channels, while one station experienced only precipitation. Decreased tidal creek and porewater salinity were significantly correlated with increased freshwater at all stations (R²=0.37 to 0.56), although porewater salinities remained hypersaline. GIS analyses indicated the most considerable vegetation change following freshwater inundation was increased cover of the annual succulentSalicornia bigelovii. Fall inundation allowed seed germination and rapid expansion of this species into previously bare areas during the subsequent winter and following spring. The station affected by both inflow and precipitation exhibited greaterS. bigelovii cover than the station affected solely by precipitation in both spring 1999 (58. 7% compared to 27.9%) and 2000 (48.6% compared to 1.9%). Percent cover of the perennialBatis maritima temporarily increased after periods of consistent rainfall. The response was short term, and cover quickly returned to pre-inundation conditions within 3 mo. Prolonged inundation led to longterm (>2yr) decreases in percent cover ofB. maritima. Our results suggest that the timing and quantity of freshwater inundation strongly dictate halophyte response to precipitation and inflow. Brief periods of freshwater inundation that occur at specific times of year alleviate stress and promote seed germination and growth, but extended soil saturation can act as a disturbance that has a negative impact on species adapted to hypersaline conditions.     

Integrated Two-Dimensional Modeling of Fluid Flow and Compaction in a Sedimentary Basin

In this paper we employ mixed finite elements and numerically study an integrated two-dimensional model of fluid flow and compaction in a sedimentary basin. This model describes a single phase incompressible flow in a two-dimensional section of a sedimentary basin with vertical compaction. At each time step, an iterative algorithm is used to solve this model. The determination of the grid movement is based on the mass conservation and movement of sediments in the basin, while the mixed method is utilized to solve the fluid flow over the moving grid. Numerical experiments are presented to verify this iterative algorithm and show representative solutions for the model under consideration.     

Unstructured, Control-Volume Distributed, Full-Tensor Finite-Volume Schemes with Flow Based Grids

Locally conservative flux-continuous, full-tensor, discretization schemes are presented for general unstructured grids. The schemes are control-volume distributed, where flow variables and rock properties are assigned to the polygonal control-volumes derived from the primal grid. A relationship between these finite volume schemes and the mixed finite element method is established. An extension for unstructured grids is described that leads to a general symmetric positive definite discretization matrix for both quadrilateral and triangular grids. A novel flow based gridding approach for unstructured mesh generation is also proposed for heterogeneous reservoir domains. Results computed with the flux continuous schemes on unstructured flow-based grids demonstrate the advantages of the methods.     

湖南印支期过铝质花岗岩的形成: 岩浆底侵与地壳加厚热效应的数值模拟

在合理构建华南印支期地质-物理模型的基础上, 利用FLAC软件, 模拟了该区印支期过铝质富钾花岗岩形成与基性岩浆底侵, 及陆壳变形叠置加厚两种动力学背景的可能联系. 模拟结果表明, 220 Ma±的基性岩浆底侵能导致地壳含水矿物相岩石的深熔, 但除非印支期存在大规模基性岩浆的底侵作用, 否则其热效应持续时间和热效应波及范围难以形成具大岩基规模的湖南印支期花岗岩. 在陆壳叠置加厚模型中, 地壳的叠置加厚可导致中下地壳界面温度升高到700℃以上, 引起片麻质岩石熔融, 当加厚因子达1.3, 白云母矿物脱水熔融产生的熔体达到熔体流动临界比例(≥20%), 从而形成花岗岩基; 结合印支期挤压逆冲推覆构造和同期基性火山岩极少出露的地质事实, 认为陆壳变形加厚可能是湖南印支期构造岩浆作用形成的主导因素.     

Improved NMO correction with a specific application to shallow-seismic data

The application of traditional NMO‐correction techniques in the processing of seismic data may result in severe distortion. This distortion is observed as a decrease in frequency content (NMO stretch) or even a total loss of data for steep velocity gradients. It is shown that the specific nature of the observed distortion is introduced by the particular implementation of the traditional NMO‐correction technique. It is also shown that other techniques that have been proposed in the literature suffer from similar artefacts. An alternative approach is suggested, based on the correction of tapered blocks of seismic data, followed by a coherence filter to compensate for the specific artefacts thus introduced. Correction of seismic amplitudes (e.g. geometrical spreading and attenuation) is implemented as an integral part of the NMO‐correction method thus allowing for both dynamic and kinematic reconstruction of interfering events. From the application of this method to synthetic and field data, it is concluded that the proposed technique may be particularly useful in the processing of shallow seismic (shear‐wave) data. Although not illustrated in this paper, it is emphasized that application of the proposed method is by no means restricted to shear waves; application to P‐wave data will prove equally useful, although the stretch effects involved may not be as severe as in the examples given in this paper.     

Analysis of the resolution function in seismic prestack depth imaging

We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions.