Rock physics analysis for time-lapse seismic at Schiehallion Field, North Sea

Rock physics analysis plays a vital role in time‐lapse seismic interpretation because it provides the link between changes in rock and fluid properties and the resulting seismic data response. In this case study of the Schiehallion Field, we discuss a number of issues that commonly arise in rock physics analyses for time‐lapse studies. We show that:

• 1 Logarithmic fits of dry bulk (K_dry) and shear (G_dry) moduli vs. effective pressure (P_eff) are superior to polynomial fits.
• 2 2D surface fits of K_dry and G_dry over porosity (φ) and effective pressure using all the core data simultaneously are more useful and accurate than separate 1D fits over φ and P_eff for each individual core.
• 3 One average set (facies) of K_dry(φ, P_eff) and G_dry(φ, P_eff) can be chosen to represent adequately the entire Schiehallion reservoir.
• 4 Saturated velocities and densities modelled by fluid substitution of K_dry(φ, P_eff), G_dry(φ, P_eff) and the dry bulk density ρ_dry(φ) compare favourably with well‐log velocities and densities.
• 5 P‐ and S‐wave impedance values resulting from fluid substitution of K_dry(φ, P_eff), G_dry(φ, P_eff) and ρ_dry(φ) show that the largest impedance changes occur for high porosities and low effective pressures.
• 6 Uncertainties in K_dry(φ, P_eff) and G_dry(φ, P_eff) derived for individual cores can be used to generate error surfaces for these moduli that represent bounds for quantifying uncertainties in seismic modelling or pressure–saturation inversion.

     

Triggering of El Niño by westerly wind events in a coupled general circulation model

Two ten-members ensemble experiments using a coupled ocean-atmosphere general circulation model are performed to study the dynamical response to a strong westerly wind event (WWE) when the tropical Pacific has initial conditions favourable to the development of a warm event. In the reference ensemble (CREF), no wind perturbation is introduced, whereas a strong westerly wind event anomaly is introduced in boreal winter over the western Pacific in the perturbed ensemble (CWWE). Our results demonstrate that an intense WWE is capable of establishing the conditions under which a strong El Niño event can occur. First, it generates a strong downwelling Kelvin wave that generates a positive sea surface temperature (SST) anomaly in the central-eastern Pacific amplified through a coupled ocean-atmosphere interaction. This anomaly can be as large as 2.5°C 60 days after the WWE. Secondly, this WWE also initiates an eastward displacement of the warm-pool that promotes the occurrence of subsequent WWEs in the following months. These events reinforce the initial warming through the generation of additional Kelvin waves and generate intense surface jets at the eastern edge of the warm-pool that act to further shift warm waters eastward. The use of a ten-members ensemble however reveals substantial differences in the coupled response to a WWE. Whereas four members of CWWE ensemble develop into intense El Niño warming as described above, four others display a moderate warming and two remains in neutral conditions. This diversity between the members appears to be due to the internal atmospheric variability during and following the inserted WWE. In the four moderate warm cases, the warm-pool is initially shifted eastward following the inserted WWE, but the subsequent weak WWE activity (when compared to the strong warming cases) prevents to further shift the warm-pool eastwards. The seasonal strengthening of trade winds in June–July can therefore act to shift warm waters back into the western Pacific, reducing the central-eastern Pacific warming. This strong sensitivity of the coupled response to WWEs may therefore limit the predictability of El Niño events, as the high frequency wind variability over the warm pool region remains largely unpredictable even at short time lead.     

Neutral hydrogen absorption observations of the central region of NGC 5929

We present 0.15-arcsec (25-pc) resolution MERLIN observations of neutral hydrogen absorption detected towards the nuclear region of the type 2 Seyfert galaxy NGC 5929. Absorption is detected only towards the north-eastern radio component with a column density of (6.5 ± 0.6) × 10²¹ cm⁻². Based on comparison with an HST WFPC2 continuum image, we propose that the absorption is caused by a 1.5-arcsec structure of neutral gas and dust offset 0.3 arcsec south-east of the nucleus and running NE–SW. A separate cloud of dust is apparent 1.5 arcsec to the south-west of the nucleus in the HST image. A comparison of the centroid velocity (2358 ± 5 km s⁻¹) and full width at half-maximum (43 ± 6 km s⁻¹) of the absorbing gas with previous [O  III ] observations suggests that both the neutral and ionized gas are undergoing galactic rotation towards the observer in the north-east and away from the observer in the south-west. The main structure is consistent with an inclined ring of gas and dust encircling the active galactic nucleus (AGN); alternatively it may be a bar or inner spiral arm. We do not detect neutral hydrogen absorption or dust obscuration against the radio nucleus (column density < 3.1 × 10²¹ cm⁻²) expected by a torus of neutral gas and dust in unified models of AGNs for a type 2 Seyfert galaxy.     

Public health impacts of climate change in Washington State: projected mortality risks due to heat events and air pollution

Illness and mortality related to heat and worsening air quality are core public health concerns associated with climate change projections. We examined the historical relationship between age- and cause-specific mortality rates from 1980 through 2006 and heat events at the 99th percentile of humidex values in the historic period from January 1, 1970 to December 31, 2006 in the greater Seattle area (King, Pierce and Snohomish counties), Spokane County, the Tri-Cities (Benton and Franklin counties) and Yakima County; the relative risks of mortality during heat events were applied to population and climate projections for Washington State to calculate number of deaths above the baseline (1980–2006) expected during projected heat events in 2025, 2045 and 2085. Three different warming scenarios were used in the analysis. Relative risks for the greater Seattle area showed a significant dose-response relationship between heat event duration and daily mortality rates for non-traumatic deaths for persons ages 45 and above, typically peaking at four days of exposure to humidex values above the 99th percentile. The largest number of projected excess deaths in all years and scenarios for the Seattle region was found for age 65 and above. Under the middle warming scenario, this age group is expected to have 96, 148 and 266 excess deaths from all non-traumatic causes in 2025, 2045 and 2085, respectively. We also examined projected excess deaths due to ground-level ozone concentrations at mid century (2045–2054) in King and Spokane counties. Current (1997–2006) ozone measurements and mid-twenty-first century ozone projections were coupled with dose-response data from the scientific literature to produce estimates overall and cardiopulmonary mortality. Daily maximum 8-h ozone concentrations are forecasted to be 16–28% higher in the mid twenty-first century compared to the recent decade of 1997–2006. By mid-century in King County the non-traumatic mortality rate related to ozone was projected to increase from baseline (0.026 per 100,000; 95% confidence interval 0.013–0.038) to 0.033 (95% CI 0.017–0.049). For the same health outcome in Spokane County, the baseline period rate of 0.058 (95% CI 0.030–0.085) was estimated increase to 0.068 (95% CI 0.035–0.100) by mid-century. The cardiopulmonary death rate per 100,000 due to ozone was estimated to increase from 0.011 (95% CI 0.005–0.017) to 0.015 (0.007–0.022) in King County, and from 0.027 (95% CI 0.013–0.042) to 0.032 (95% CI 0.015–0.049) in Spokane County. Public health interventions aimed at protecting Washington’s population from excessive heat and increased ozone concentrations will become increasingly important for preventing deaths, especially among older adults. Furthermore, heat and air quality related illnesses that do not result in death, but are serious nevertheless, may be reduced by the same measures.     

Adsorption kinetics of CO2, CH4, and their equimolar mixture on coal from the Black Warrior Basin,West-Central Alabama

Laboratory experiments were conducted to investigate the adsorption kinetic behavior of pure and mixed gases (CO_2, CH₄, approximately equimolar CO₂ + CH₄ mixtures, and He) on a coal sample obtained from the Black Warrior Basin at the Littleton Mine (Twin Pine Coal Company), Jefferson County, west-central Alabama. The sample was from the Mary Lee coal zone of the Pottsville Formation (Lower Pennsylvanian). Experiments with three size fractions (45–150 µm, 1–2 mm, and 5–10 mm) of crushed coal were performed at 40 °C and 35 °C over a pressure range of 1.4–6.9 MPa to simulate coalbed methane reservoir conditions in the Black Warrior Basin and provide data relevant for enhanced coalbed methane recovery operations. The following key observations were made: (1) CO₂ adsorption on both dry and water-saturated coal is much more rapid than CH₄ adsorption; (2) water saturation decreases the rates of CO₂ and CH₄ adsorption on coal surfaces, but it appears to have minimal effects on the final magnitude of CO₂ or CH₄ adsorption if the coal is not previously exposed to CO₂; (3) retention of adsorbed CO₂ on coal surfaces is significant even with extreme pressure cycling; and (4) adsorption is significantly faster for the 45–150 μm size fraction compared to the two coarser fractions.     

The statistical properties of Λ cold dark matter halo formation

We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press–Schechter (EPS) formalism and the related galform Monte Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present-day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift, accretion rates, and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present-day mass distribution of haloes into which high-redshift progenitors of a specific mass become incorporated. We find that the EPS formalism and its Monte Carlo extension capture the qualitative behaviour of all these statistics, but as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.     

Potential impact of ash eruptions on dairy farms from a study of the effects on a farm in eastern Bay of Plenty,New Zealand; implications for hazard mitigation

This paper investigates the impact ash fall would have on dairy farming, based on a study of ‘Tulachard’, a dairy farming operation at Rerewhakaaitu, North Island, New Zealand. It includes analysis of the potential effects on the dairy shed and milking machine, electrical supply and distribution, water supply and distribution, tractors and other farm vehicles, farm buildings (haysheds, pump sheds, implement sheds, etc.), milk-tanker access to the farm and critical needs of dairy cows and farm to keep milking. One of the most vulnerable areas identified in the study was the cooling of milk at the milking shed, pending dairy tanker pick-up. The cooling system’s condenser is exposed to the atmosphere and falling ash would make it highly vulnerable. Laboratory testing with wet and dry ash was conducted to determine its resilience to ash ingestion. It was found to perform satisfactorily during dry testing, but during wet testing significant clogging/blocking of the condenser’s radiator occurred, dramatically reducing airflow through the condenser. Specific mitigation recommendations have been developed that include cleaning with compressed air and adapting farm management techniques to lessen usage of the condenser during an ash-fall event. Specific recommendations for management of dairy farm operation are given to mitigate the effects of an ash-fall event.

Spatial and temporal patterns of nutrient concentration and export in the tidal Hudson River

Spatial and temporal dynamics of N and P were examined in the tidal Hudson River between 1992 and 1996. For all seasons and at all locations in the river nutrient concentrations were generally quite high. TN averaged 60 μM and was above 50 μM in 75% of samples. TP averaged 1.7 μM and was above 1.2 μM in 75% of samples. NO₃ was the dominant form of N (60% of TN) while PO₄ comprised about 40% of TP. Seasonal and spatial variation in most N and P components was quite low but patterns were apparent. Seasonally, forms of N (TN, NO₃ and NH₄) and PO₄ showed opposite patterns. All N components showed summertime decreases, but PO₄ increased over the summer. Spatially, along the 200 km fresh to oligohaline stretch, N and P showed similar patterns—declining from upper to mid sections of the river but subsequently increasing in most down river, oligohaline stretches. The down river increase in nutrients is likely caused by a combination of sewage inputs and salinity-related geochemical release of P. A preliminary budget of the upper to the mid section of the river (a 100 km stretch) suggests that the decline in nutrient concentration in this section is due to the net retention of almost 2,000 mT N and 200 mT P per year or about 20% of the N and P input to this section of river. The retention in tidal rivers, like the Hudson, occurs immediately above the estuary and may, therefore, be relatively more significant than retention occurring higher in the watershed.     

Impacts of Methane on Carbon Dioxide Storage in Brine Formations

In the context of geological carbon sequestration (GCS), carbon dioxide (CO₂) is often injected into deep formations saturated with a brine that may contain dissolved light hydrocarbons, such as methane (CH₄). In this multicomponent multiphase displacement process, CO₂ competes with CH₄ in terms of dissolution, and CH₄ tends to exsolve from the aqueous into a gaseous phase. Because CH₄ has a lower viscosity than injected CO₂, CH₄ is swept up into a ‘bank’ of CH₄‐rich gas ahead of the CO₂ displacement front. On the one hand, this may provide a useful tracer signal of an approaching CO₂ front. On the other hand, the emergence of gaseous CH₄ is undesirable because it poses a leakage risk of a far more potent greenhouse gas than CO₂ if the cap rock is compromised. Open fractures or faults and wells could result in CH₄ contamination of overlying groundwater aquifers as well as surface emissions. We investigate this process through detailed numerical simulations for a large‐scale GCS pilot project (near Cranfield, Mississippi) for which a rich set of field data is available. An accurate cubic‐plus‐association equation‐of‐state is used to describe the non‐linear phase behavior of multiphase brine‐CH₄‐CO₂ mixtures, and breakthrough curves in two observation wells are used to constrain transport processes. Both field data and simulations indeed show the development of an extensive plume of CH₄‐rich (up to 90 mol%) gas as a consequence of CO₂ injection, with important implications for the risk assessment of future GCS projects.     

Investigating potential future changes in surface water flooding hazard and impact

Surface water flooding (SWF) is a recurrent hazard that affects lives and livelihoods. Climate change is projected to change the frequency of extreme rainfall events that can lead to SWF. Increasingly, data from Regional Climate Models (RCMs) are being used to investigate the potential water-related impacts of climate change; such assessments often focus on broad-scale fluvial flooding and the use of coarse resolution (>12 km) RCMs. However, high-resolution (<4 km) convection-permitting RCMs are now becoming available that allow impact assessments of more localised SWF to be made. At the same time, there has been an increasing demand for more robust and timely real-time forecast and alert information on SWF. In the UK, a real-time SWF Hazard Impact Model framework has been developed. The system uses 1-km gridded surface runoff estimates from a hydrological model to simulate the SWF hazard. These are linked to detailed inundation model outputs through an Impact Library to assess impacts on property, people, transport, and infrastructure for four severity levels. Here, a set of high-resolution (1.5 km and 12 km) RCM data has been used as input to a grid-based hydrological model over southern Britain to simulate Current (1996–2009) and Future (~2100s; RCP8.5) surface runoff. Counts of threshold-exceedance for surface runoff and precipitation (at 1-, 3- and 6-hr durations) are analysed. Results show that the percentage increases in surface runoff extremes, are less than those of precipitation extremes. The higher-resolution RCM simulates the largest percentage increases, which occur in winter, and the winter exceedance counts are greater than summer exceedance counts. For property impacts, the largest percentage increases are also in winter; however, it is the 12-km RCM output that leads to the largest percentage increase in impacts. The added-value of high-resolution climate model data for hydrological modelling is from capturing the more intense convective storms in surface runoff estimates.