Time series analysis of CO concentrations from an Eastern Kentucky coal fire

Coal fires in natural outcrops, in abandoned and active coal mines, and in coal and coal-refuse piles are responsible for the uncontrolled emissions of gases, including CO, CO₂, H₂S, hydrocarbons, and volatile aromatics. Typically, measurements of gases at a mine vent are made over a short time interval, perhaps no more than 10 min, including the time for replicate measurements. Such timing provides little information on longer-term variations in emissions, although comparisons of seasonal measurements suggest such changes do occur. To address this problem, we placed temperature and CO data loggers in coal fire vents to collect time series measurements for a period of up to three weeks. For one experiment, 11 days of data at one-minute intervals indicated that the CO emissions were generally in the 400–550 ppmv range. However near daily depressions in CO concentrations occurred and in some cases fell below 50 ppmv; followed by an increase to ~ 700–800 ppmv; in turn followed by a return to the ambient conditions. Data for a separate 21-day collection period in a different vent of the same fire exhibited similar trends, albeit at a higher CO concentration. The drop in CO concentration may be associated with a meteorologically-driven inhalation cycle of the fire, whereby air diluted the combustion generated CO. We propose this was followed by an increase in the intensity of the fire due to increased O₂ from the inhaled air, producing increased CO concentrations, before settling back to the ambient conditions.     

Comprehensive simulation of the middle atmospheric climate: some recent results

This study discusses the results of comprehensive time-dependent, three-dimensional numerical modelling of the circulation in the middle atmosphere obtained with the GFDL SKYHI troposphere-stratosphere-mesosphere general circulation model (GCM). The climate in a long control simulation with an intermediate resolution version (3° in horizontal) is briefly reviewed. While many aspects of the simulation are quite realistic, the focus in this study is on remaining first-order problems with the modelled middle atmospheric general circulation, notably the very cold high latitude temperatures in the Southern Hemisphere (SH) winter/spring, and the virtual absence of a quasi-biennial oscillation (QBO) in the tropical stratosphere. These problems are shared by other extant GCMs. It was noted that the SH cold pole problem is somewhat ameliorated with increasing horizontal resolution in the model. This suggests that improved resolution increases the vertical momentum fluxes from the explicitly resolved gravity waves in the model, a point confirmed by detailed analysis of the spectrum of vertical eddy momentum flux in the winter SH extratropics. This result inspired a series of experiments with the 3° SKYHI model modified by adding a prescribed zonally-symmetric zonal drag on the SH winter westerlies. The form of the imposed momentum source was based on the simple assumption that the mean flow drag produced by unresolved waves has a spatial distribution similar to that of the Eliassen-Palm flux divergence associated with explicitly resolved gravity waves. It was found that an appropriately-chosen drag confined to the top six model levels (above 0.35 mb) can lead to quite realistic simulations of the SH winter flow (including even the stationary wave fields) through August, but that problems still remain in the late-winter/springtime simulation. While the imposed momentum source was largely confined to the extratropics, it produced considerable improvement in the simulation of the equatorial semiannual oscillation, with both the easterly and westerly phases being somewhat more intense than in the control simulation. A separate experiment was conducted in which the SKYHI model was simplified so that it had no topography and so that the seasonal cycle was frozen in perpetual equinox conditions. These changes result in a model that has much reduced interhemispheric asymmetry. This model spontaneously produces a long period mean flow oscillation of considerable amplitude in the tropical upper stratopause. The implication of this result for the general issue of obtaining a QBO in comprehensive GCMs is discussed.     

Integration of the life of field seismic data with the reservoir model at the Valhall Field‡

The frequent time‐lapse observations from the life of field seismic system across the Valhall field provide a wealth of information. The responses from the production and injection wells can be observed through time‐shift and amplitude changes. These observations can be compared to modelled synthetic seismic responses from a reservoir simulation model of the Valhall Field. The observed differences between the observations and the modelling are used to update and improve the history match of the reservoir model. The uncertainty of the resulting model is reduced and a more confident prediction of future reservoir performance is provided. A workflow is presented to convert the reservoir model to a synthetic seismic response and compare the results to the observed time‐lapse responses for any time range and area of interest. Correlation based match quality factors are calculated to quantify the visual differences. This match quality factor allows us to quantitatively compare alternative reservoir models to help identify the parameters that best match the seismic observations. Three different case studies are shown where this workflow has helped to reduce the uncertainty range associated with specific reservoir parameters. By updating various reservoir model parameters we have been able to improve the match to the observations and thereby improve the overall reservoir model predictability. The examples show positive results in a range of different reservoir modelling issues, which indicates the flexibility of this workflow and the ability to have an impact in most reservoir modelling challenges.     

Nickel, Copper, Zinc and Cadmium Cycling with Manganese in Lake Vanda (Wright Valley, Antarctica)

Lake Vanda is a closed-basin, permanently ice-covered lake located in the Wright Valley of Antarctica. The lakes more important geochemical features include the fact that it is fed by a single glacial melt water stream for only 6–8weeks out of the year; that it has remained stratified for more than a millennium; and that, like other lakes in the region, it is remote from anthropogenic influence. These, together with the fact that it is among the least biologically productive lakes in the world, make it an ideal system for examining the transport, cycling and fate of trace metals in the aquatic environment. Like others before us, we view this lake as a natural geochemical laboratory, a flask in the desert. This paper presents the first set of closely spaced, vertical, profiles for dissolved and particulate Mn, Fe, Ni, Cu, Zn and Cd in the water column. Despite the absence of an outflow, metals in the fresh upper waters of the lake have extremely low concentrations, in the pico-molar to nano-molar range, and are partitioned largely into dissolved rather than particulate phases. Efficient metal scavenging by particles from these oxygen-rich waters is indicated. Significant increases in metal concentrations begin to appear at depth, between 57 and 60m, and these increases coincide with the onset of manganese oxide dissolution in oxic, but lower pH waters. Vertical profiles suggest that the entire suite of trace metals (Ni, Cu, Zn, and Cd) is being released from manganese oxide carrier phases. Thermodynamic analysis indicates that Mn3O4 (i.e., the mineral hausmannite) may be important in metal sequestration and recycling in the deeper waters of Lake Vanda. Manganese-reducing organisms reported by Bratina etal. (1998) are active in the zone of metal release and these could also contribute to the observed cycling.     

Application of airborne LiDAR data and airborne multispectral imagery to structural mapping of the upper section of the Troodos ophiolite, Cyprus

Structural maps are traditionally produced by mapping features such as faults, folds, fabrics, fractures and joints in the field. However, large map areas and the spatially limited ground perspective of the field geologist can potentially increase the likelihood that not all structural features will be identified within a given area. The ability to recognise and map both local and regional structural features using high-resolution remote sensing data provides an opportunity to complement field-based mapping to help generate more comprehensive structural maps. Nonetheless, vegetation cover can adversely affect the extraction of structural information from remotely sensed data as it can mask the appearance of subtle spectral and geomorphological features that correspond to geological structures. This study investigates the utility of airborne Light Detection And Ranging (LiDAR) data and airborne multispectral imagery for detailed structural mapping in vegetated ophiolitic rocks and sedimentary cover of a section of the northern Troodos ophiolite, Cyprus. Visual enhancement techniques were applied to a 4-m airborne LiDAR digital terrain model and 4-m airborne multispectral imagery to assist the generation of structural lineament maps. Despite widespread vegetation cover, dykes and faults were recognisable as lineaments in both data sets, and the predominant strike trends of lineaments in all resulting maps were found to be in agreement with field-based structural data. Interestingly, prior to fieldwork, most lineaments were assumed to be faults, but were ground-verified as dykes instead, emphasising the importance of ground-truthing. Dyke and fault trends documented in this study define a pervasive structural fabric in the upper Troodos ophiolite that reflects the original sea-floor spreading history in the Larnaca graben. This structural fabric has not previously been observed in such detail and is likely to be continuous in adjacent regions under sedimentary cover. This information may be useful to future exploration efforts in the region focused on identification of structurally controlled mineral and groundwater resources. Overall, our case study highlights the efficacy of airborne LiDAR data and airborne multispectral imagery for extracting detailed and accurate structural information in hard-rock terrain to help complement field-based mapping.     

Saturn’s tropospheric composition and clouds from Cassini/VIMS 4.6-5.1 μm nightside spectroscopy

The latitudinal variation of Saturn’s tropospheric composition (NH₃, PH₃ and AsH₃) and aerosol properties (cloud altitudes and opacities) are derived from Cassini/VIMS 4.6-5.1 μm thermal emission spectroscopy on the planet’s nightside (April 22, 2006). The gaseous and aerosol distributions are used to trace atmospheric circulation and chemistry within and below Saturn’s cloud decks (in the 1- to 4-bar region). Extensive testing of VIMS spectral models is used to assess and minimise the effects of degeneracies between retrieved variables and sensitivity to the choice of aerosol properties. Best fits indicate cloud opacity in two regimes: (a) a compact cloud deck centred in the 2.5-2.8 bar region, symmetric between the northern and southern hemispheres, with small-scale opacity variations responsible for numerous narrow light/dark axisymmetric lanes; and (b) a hemispherically asymmetric population of aerosols at pressures less than 1.4 bar (whose exact altitude and vertical structure is not constrained by nightside spectra) which is 1.5-2.0× more opaque in the summer hemisphere than in the north and shows an equatorial maximum between ±10° (planetocentric).Saturn’s NH₃ spatial variability shows significant enhancement by vertical advection within ±5° of the equator and in axisymmetric bands at 23-25°S and 42-47°N. The latter is consistent with extratropical upwelling in a dark band on the poleward side of the prograde jet at 41°N (planetocentric). PH₃ dominates the morphology of the VIMS spectrum, and high-altitude PH₃ at p < 1.3 bar has an equatorial maximum and a mid-latitude asymmetry (elevated in the summer hemisphere), whereas deep PH₃ is latitudinally-uniform with off-equatorial maxima near ±10°. The spatial distribution of AsH₃ shows similar off-equatorial maxima at ±7° with a global abundance of 2-3 ppb. VIMS appears to be sensitive to both (i) an upper tropospheric circulation (sensed by NH₃ and upper-tropospheric PH₃ and hazes) and (ii) a lower tropospheric circulation (sensed by deep PH₃, AsH₃ and the lower cloud deck).     

Probabilistic data integration for characterization of spatial distribution of residual LNAPL

The spatial distribution of residual light non-aqueous phase liquid (LNAPL) is an important factor in reactive solute transport modeling studies. There is great uncertainty associated with both the areal limits of LNAPL source zones and smaller scale variability within the areal limits. A statistical approach is proposed to construct a probabilistic model for the spatial distribution of residual NAPL and it is applied to a site characterized by ultra-violet-induced-cone-penetration testing (CPT–UVIF). The uncertainty in areal limits is explicitly addressed by a novel distance function (DF) approach. In modeling the small-scale variability within the areal limits, the CPT–UVIF data are used as primary source of information, while soil texture and distance to water table are treated as secondary data. Two widely used geostatistical techniques are applied for the data integration, namely sequential indicator simulation with locally varying means (SIS–LVM) and Bayesian updating (BU). A close match between the calibrated uncertainty band (UB) and the target probabilities shows the performance of the proposed DF technique in characterization of uncertainty in the areal limits. A cross-validation study also shows that the integration of the secondary data sources substantially improves the prediction of contaminated and uncontaminated locations and that the SIS–LVM algorithm gives a more accurate prediction of residual NAPL contamination. The proposed DF approach is useful in modeling the areal limits of the non-stationary continuous or categorical random variables, and in providing a prior probability map for source zone sizes to be used in Monte Carlo simulations of contaminant transport or Monte Carlo type inverse modeling studies.     

Use of suppressive subtractive hybridization and cDNA arrays to discover patterns of altered gene expression in the liver of dihydrotestosterone and 11-ketotestosterone exposed adult male largemouth bass (Micropterus salmoides)

In this study male largemouth bass (LMB) were exposed to the naturally occurring androgens, dihydrotestosterone (DHT) or 11-ketotestosterone (11-KT) in order to identify genes that are differentially regulated by these steroid hormones. Using subtractive hybridization on livers of fish treated with DHT against vehicle control, many novel LMB genes were cloned. These genes were added to our gene library and arrayed. Six genes were up-regulated and five were down-regulated by both androgens. But, each androgen also regulated specific genes. One gene that was identified as a potential androgen marker was spermidine-spermine-N(1)-acetyltransferase that was up-regulated by both androgens. Determining which genes are responsive to natural androgens will help to identify biochemical pathways that are impacted.     

Structure and subduction processes along the Oregon-Washington margin

Seismic reflection and refraction data off Washington and Oregon are used to determine the style of sediment deformation and to infer the physical properties of accreted sediments on the lower slope. Onshore-offshore seismic refraction data off Washington are used to determine the location of the trench, or where the plate bending starts.We find that off Washington the subduction zone is characterized by a trench whose physiographic expression is buried under several kilometers of sediments and is tens of kilometers landward of the lower slope, which is accreting seaward as the result of the offscraping of sediments.Seismic reflection data support previous observations that offscraping occurs along seaward and landward dipping thrust faults. Refraction data indicate that a sediment package thrust up along a seaward dipping fault (off Washington) was not measurably changed in velocity with respect to a Cascadia basin section. However a package uplifted by thrusting along a landward dipping fault (off Oregon) did have increased velocity. It is suggested that the increased velocities off Oregon could be the result of erosion and exposure of more deeply buried and compacted sediments, rather than the result of dewatering due to tectonic stress. Off Washington the sensitivity of velocity to porosity and resolution of the seismic method does not preclude dewatering due to tectonic stress, but it does limit the degree of dewatering.In the deeper parts of the lower slope section off Washington and Oregon velocities as high as 3 to 4 km/sec are found. Heat flow data indicate that the temperatures in this high velocity regime are greater than 100°C. It is hypothesized that lithification related to clay diagenesis may be partly responsible for the high velocities, rather than simply compaction. It also appears that the high velocity sediments are subducted while the unlithified low velocity sediments are offscraped.