Application of a Bayesian model to infer the contribution of coalbed natural gas produced water to the Powder River,Wyoming and Montana

The Powder River Basin (PRB) of Wyoming and Montana contains significant coal and coal bed natural gas (CBNG) resources. CBNG extraction requires the production of large volumes of water, much of which is discharged into existing drainages. Compared to surface waters, the CBNG produced water is high in sodium relative to calcium and magnesium, elevating the sodium adsorption ratio (SAR). To mitigate the possible impact this produced water may have on the quality of surface water used for irrigation, the State of Montana passed water anti‐degradation legislation, which could affect CBNG production in Wyoming. In this study, we sought to determine the proportion of CBNG produced water discharged to tributaries that reaches the Powder River by implementing a four end‐member mixing model within a Bayesian statistical framework. The model accounts for the ⁸⁷Sr/⁸⁶Sr, δ¹³C_DIC, [Sr] and [DIC] of CBNG produced water and surface water interacting with the three primary lithologies exposed in the PRB. The model estimates the relative contribution of the end members to the river water, while incorporating uncertainty associated with measurement and process error. Model results confirm that both of the tributaries associated with high CBNG activity are mostly composed of CBNG produced water (70–100%). The model indicates that up to 50% of the Powder River is composed of CBNG produced water downstream from the CBNG tributaries, decreasing with distance by dilution from non‐CBNG impacted tributaries from the point sources to ~10–20% at the Montana border. This amount of CBNG produced water does not significantly affect the SAR or electrical conductivity of the Powder River in Montana. Copyright © 2013 John Wiley & Sons, Ltd.     

Radon in soil anomaly observed at Los Azufres geothermal field,Michoacan: a possible precursor of the 1985 Mexico earthquake (Ms = 8.1)

A radon anomaly in a nuclear track detector placed on a fault was detected prior to the destructive (Ms = 8.1) Mexican earthquake of 19 September 1985. The fault is a structural feature of the geothermal field Los Azufres, 260 km NE of the epicentral area. Since no other phenomenon can be found as a possible cause for the radon anomaly, it is proposed that preseismic regional stress perturbations may produce changes in the fluid transport patterns at the fault, suggesting that radon measurements in similar conditions could provide a good precursor in high seismic risk areas.     

Fractal analysis of long-range paleoclimatic data: Oxygen isotope record of pacific core V28-239

R/S analysis of the oxygen isotope curve of Pacific core V28-239 yields a fractal dimension of 1.22. This value is considered to characterize global climatic change over the last 2 million years as expressed by changing O¹⁸ ratios and confirms that climatic variations are characterized by long-term persistence. The fractal dimension of 1.22 compares favorably with the approximate fractal dimension of 1.26 for annual precipitation records for nine major cities in the United States. Although the precipitation and oxygen isotope data are measured in different physical units and recorded at different time scales, fractal analysis allows for a mathematical comparison of the two phenomena. Additionally, since the fractal dimensions of the oxygen isotope and precipitation records are similar, it is implied that such fractal dimensions are characteristic of climate change over the spectral range of 10 to 10⁶ years. Given this temperature curves based on fractal parameters of long-term O¹⁸ data could be constructed which would allow examination of characteristics of temperature variation over tens and hundreds of years. Such studies may allow the establishment of limits on natural temperature variation and document the persistence of temperature trends through time. If these trends and limits can be resolved, long-range climatic prediction is feasible.     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Explorations into the formal structure of drainage basins

Every basin of higher than first order is drained by a channel network composed of two subnetworks. Their basins are separated by a drainage divide line, called the basin divider, which is the primary organizing feature of the main basin. Each basin of magnitude n contains n – 1 subnetworks of higher order, and is therefore organized by a set of n – 1 dividers. The dividers and the basin boundary are interconnected in a graph called the divider network of the basin; in graph-theoretic terms this network forms a tree and has the same magnitude and link numbers as the channel network draining the basin. While the subbasins and subnetworks of a drainage basin form a nesting hierarchy, the corresponding dividers do not; indeed, any two dividers share at most one node in common, and whether they do so is independent of whether the corresponding subbasins are nesting or disjoint. However, the dividers of nesting basins are linked by recursive relationships which permit the derivation of a set of algebraic equations; these equations relate the dividers of a basin to other basin components; for example, their combined length is equal to half the length of all first-order basin boundaries minus the length of the main basin boundary. The second part of the paper explores the dependence of the divider length on other basin parameters. The expected length, as predicted by the assumption of topological randomness, is clearly rejected by the data. An alternative approach (regression) is based on the observed magnitudes of the subbasins separated by each divider, and is reasonably successful in estimating divider length. The last section introduces the concept of the standardized basin defined by a boundary length of unity; the estimated lengths of the basin divider and the basin boundary permit an approximate reconstruction of the idealized basin shape and the location of the divider in it.     

Coastal Massachusetts pond development: edaphic,climatic, and sea level impacts since deglaciation

Kettle ponds in the Cape Cod National Seashore in southeastern Massachusetts differ in their evolution due to depth of the original ice block, the clay content of outwash in their drainage basins, and their siting in relation to geomorphic changes caused by sea-level rise, barrier beach formation, and saltmarsh development. Stratigraphic records of microfossil, carbon isotope, and sediment changes also document late-glacial and Holocene climatic changes.The ponds are separated into 3 groups, each of which follow different development scenarios. Group I ponds date from the late-glacial. They formed in clay-rich outwash, have perched aquifers and continuous lake sediment deposition. The earliest pollen and macrofossil assemblages in Group I pond sediments suggest tundra and spruce-willow parklands before 12 000 yr B.P., boreal forest between 12 000 and 10 500 yr B.P., bog/heath initiation and expansion during the Younger Dryas between 11 000 and 10 000 yr B.P., northern conifer forest between 10 500 and 9500 yr B.P., and establishment of the Cape oak and pitch pine barrens vegetation after 9500 yr B.P. Sedimentation rate changes suggest lowered freshwater levels between 9000 and 5000 yr B.P. caused by decreased precipitation on the Atlantic Coastal Plain. Lake sediment deposition began in the middle Holocene in Group II ponds which formed in clay-poor outwash. These ponds date from about 6000-5000 yr B.P. In these ponds sediment deposition began as sea level rose and the freshwater lens intersected the dry basins. The basal radiocarbon dates of these ponds and stable carbon isotope analyses of the pond sediments suggest a sea-level curve for Cape Cod Bay. Holocene topographic changes in upland and the landscape surrounding the ponds is reconstructed for this coastal area.Group III ponds in the late Holocene landscape of the Provincelands dunes originated as interdunal bogs about 1000 yr B.P. and became ponds more recently as water-levels increased. Peat formation in the Provincelands reflects climatic changes evident on both sides of the Atlantic region.This is the 8th in a series of papers published in this special AMQUA issue. These papers were presented at the 1994 meeting of the American Quaternary Association held 19–22 June, 1994, at the University of Minnesota, Minneapolis, Minnesota, USA. Dr Linda C. K. Shane served as guest editor for these papers.     

Influence of multipurpose trees on agricultural crops in arid regions of Haryana, India

The effects of multipurpose trees on the productivity of agricultural crops were studied in the arid regions of Haryana. In one of the experiments, wheat was grown on irrigated farms having scattered trees ofDalbergia sissoo, Azadirachta indica, Prosopis cinerariaandAcacia nilotica. Data of wheat yield for each tree species at different distances (1, 3, 5, and 7 m) and four directions (East, West, North and South) from the tree bases and control (no trees) were collected. Results indicate thatAzadirachta indicaandProsopis cinerariadid not produce any significant difference in the wheat yield while the other two species (Dalbergia sissooandAcacia nilotica) gave a reduction in yield.A. niloticahad a more significant and prominent effect and a reduction of 40 to 60 % wheat yield was observed.Dalbergia sissooreduced yield by 4 to 30 % but the reduction was only up to a distance of 3 m. In general, the impact of trees on productivity was observed up to 3 m distance and there was little, if any, impact up to 5 m and almost no impact at 7 m.In another experiment,D. sissoo(18-years-old) was raised as windbreak and productivity of cotton was assessed. The tree belt reduced the wind speed by 15 to 45 %, depending on season and wind speed. The influence of trees on cotton yield was observed to be negative up to 2H (36 m), where H is tree height, i.e. 18 m. Cotton productivity was observed to be maximum between 2H and 5H. Depending upon orientation of tree belt, the increase in cotton yield was found to be 4 to 10 %.     

海城震区深地震测深资料A/D转换方法及其效果的研究

介绍深地震测深模拟磁带资料在CROMEMCO型微机上实现A/D转换的方法及其效果。在实现A/D转换研究中,进行了时标信号采集方式、六响信号处理、模拟回放电路改进等方法的试验,从而提高了A/D转换精度与效果。经对海城震区实际资料分析可知:该系统A/D转换精度达到了规范要求。     

Chronology of Sanbagawa metamorphism

By collating age data based on the fossil age of the protoliths, radiometric dating of the metamorphic minerals, and sedimentary records of erosion at the earth's surface, the history of the Sanbagawa metamorphism can be summarized as follows. (1) The pre-metamorphic sedimentary rocks (Carboniferous-Jurassic + Early Cretaceous?) became mixed and formed a thickened packet in the vicinity of an ancient trench through a variety of subduction-related tectono-sedimentary processes, probably in Early Cretaceous time (c., 130-120 Ma). (2) The subducted protoliths underwent progressive metamorphism reaching a maximum depth of c. 30 km in late Early Cretaceous time (c. 116 ± 10 Ma). (3) The high-P/T metamorphic rocks began to rise toward the surface (during the interval 110-50 Ma) with minimum estimates for the average cooling rate around 9-12°C/Ma and an average uplift rate around 0.4-0.5 mm/year. (4) Finally, at some stage after reaching the erosional surface, the high-P/T metamorphic rocks were covered unconformably by the middle Eocene (c. 50-42 Ma) Kuma Group. On the basis of the present chronological summary of the Sanbagawa metamorphism, the areal extent of the Sanbagawa metamorphism is also discussed with respect to the weakly metamorphosed subduction-accretion complex of the next tectonic belt to the south, the Northern Chichibu belt.