Mechanical and transport constitutive models for fractures subject to dissolution and precipitation

Transient changes in the permeability of fractures in systems driven far‐from‐equilibrium are described in terms of proxy roles of stress, temperature and chemistry. The combined effects of stress and temperature are accommodated in the response of asperity bridges where mineral mass is mobilized from the bridge to the surrounding fluid. Mass balance within the fluid accommodates mineral mass either removed from the flow system by precipitation or advection, or augmented by either dissolution or advection. Where the system is hydraulically closed and initially at equilibrium, reduction in aperture driven by the effects of applied stresses and temperatures will be augmented by precipitation on the fracture walls. Where the system is open, the initial drop in aperture may continue, and accelerate, where the influent fluid is oversaturated with respect to the equilibrium mineral concentration within the fluid, or may reverse, if undersaturated. This simple zero‐dimensional model is capable of representing the intricate behavior observed in experiments where the feasibility of fracture sealing concurrent with net dissolution is observed. This zero‐order model is developed as a constitutive model capable of representing key aspects of changes in the transport parameters of the continuum response of fractured media to changes in stress, temperature and chemistry. Copyright © 2009 John Wiley & Sons, Ltd.     

Terrestrial – Ocean environmental change in the northwestern Pacific from the glacial times to Holocene

The Mixed Water Region (MWR) between the Oyashio Current and the Kuroshio Extension and surrounding terrestrial region have been affected by the global climatic change. The sedimentary core KR02-15 PC6 collected at 40°23.8871′N, 143°29.8663′E in the MWR provides detailed records of pollen and spores and alkenone sea surface temperature (SST) during the last 15 kyr, showing that both ocean and terrestrial environments have been appreciably linked. Abundance of Picea, a proxy for terrestrial cool climate, fluctuated, coinciding with the alkenone-SST change during the interstadial period (15.0–7.7 cal kyr B.P.). Alkenone-SST increased to the modern level at 8.0 cal kyr B.P. and remained almost constant until 2.0 kyr B.P. while the terrestrial warming still continued until 4.0 cal kyr B.P. The rainfall under terrestrial warm condition increased gradually from 7.9 cal kyr B.P. and rapidly from 2.0 cal kyr B.P. Both terrestrial and marine environments cooled down during the last 2 kyr possibly due to the intensified cold Oyashio Current.     

A numerical model simulating reactive transport and evolution of fracture permeability

A numerical model is presented to describe the evolution of fracture aperture (and related permeability) mediated by the competing chemical processes of pressure solution and free‐face dissolution/precipitation; pressure (dis)solution and precipitation effect net‐reduction in aperture and free‐face dissolution effects net‐increase. These processes are incorporated to examine coupled thermo‐hydro‐mechano‐chemo responses during a flow‐through experiment, and applied to reckon the effect of forced fluid injection within rock fractures at geothermal and petroleum sites. The model accommodates advection‐dominant transport systems by employing the Lagrangian–Eulerian method. This enables changes in aperture and solute concentration within a fracture to be followed with time for arbitrary driving effective stresses, fluid and rock temperatures, and fluid flow rates. This allows a systematic evaluation of evolving linked mechanical and chemical processes. Changes in fracture aperture and solute concentration tracked within a well‐constrained flow‐through test completed on a natural fracture in novaculite (Earth Planet. Sci. Lett. 2006, in press) are compared with the distributed parameter model. These results show relatively good agreement, excepting an enigmatic abrupt reduction in fracture aperture in the early experimental period, suggesting that other mechanisms such as mechanical creep and clogging induced by unanticipated local precipitation need to be quantified and incorporated. The model is applied to examine the evolution in fracture permeability for different inlet conditions, including localized (rather than distributed) injection. Predictions show the evolution of preferential flow paths driven by dissolution, and also define the sense of permeability evolution at field scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Observation of cosmic soft X-rays

Cosmic soft X-rays in the energy range between 0.14 and 7 keV were observed with thin polypropylene window proportional counters on board a sounding rocket. The field of view crossed the galactic plane in the Cygnus-Cassiopeia region at a large angle and reached the galactic latitudes of –55° and +30°. Referring also to the result with Be window counters, we obtained the energy spectrum of Cyg XR-2, the flux from the Cas A region and the distribution of the intensity of diffuse X-rays over the scanned region. The turn-over of the Cyg XR-2 spectrum at about 1 keV indicates that the distance of the Cyg XR-2 source lies between 600 and 800 pc, if the turn-over is due entirely to interstellar absorption. The flux from the Cas A region is obtained as 0.23±0.05 photons cm^–2 sec^–1 in the energy range between 1.1 and 4.1 keV. The intensity of diffuse soft X-rays depends on the galactic latitude more weakly than expected from the interstellar absorption of extragalactic X-rays and shows asymmetry with respect to the galactic equator, thus suggesting a contribution of galactic X-rays. The spectrum of extragalactic X-rays is approximately represented by a power lawE ^–1.8.     

Laboratory experiments and thermal calculations for the development of a next-generation glacier-ice exploration system: Development of an electro-thermal drilling device

A next-generation drilling system, equipped with a thermal drilling device, is proposed for glacier ice. The system is designed to penetrate glacier ice via melting of the ice and continuously analyze melt-water in a contamination-free sonde. This new type of drilling system is expected to provide analysis data in less time and at less cost than existing systems. Because of the limited number of parameters that can be measured, the proposed system will not take the place of conventional drilling systems that are used to obtain ice cores; however, it will provide a useful method for quickly and simply investigating glacier ice.An electro-thermal drilling device is one of the most important elements needed to develop the proposed system. To estimate the thermal supply required to reach a target depth in a reasonable time, laboratory experiments were conducted using ice blocks and a small sonde equipped solely with heaters. Thermal calculations were then performed under a limited range of conditions. The experiments were undertaken to investigate the effects of the shape and material of the drill head and heater temperature on the rate of penetration into the ice. Additional thermal calculations were then performed based on the experimental results.According to the simple thermal calculations, if the thermal loss that occurs while heat is transferred from the heater to ice (in melting the ice) is assumed to be 50%, the total thermal supply required for heaters in the sonde and cable is as follows: (i) 4.8 kW (sonde) plus 0 W (cable) to penetrate to 300 m depth over 10 days into temperate glacier ice for which the temperature is 0 °C at all depths and to maintain a water layer along 300 m of cable; (ii) 10 kW (sonde) plus 19–32 kW (cable) to penetrate to 1000 m depth over 1 month into cold glacier ice for which the temperature is −25 °C at the surface and 0 °C at 1000 m depth and to maintain a water layer along 1000 m of cable; and (iii) 19 kW (sonde) plus 140–235 kW (cable) to penetrate to 3000 m depth over 2 months into an ice sheet for which the temperature is −55 °C at the surface and 0 °C at 3000 m depth and to maintain a water layer along 3000 m of cable. The thermal supply required for the cable is strongly affected by the thickness of the water layer, cable diameter, and the horizontal distance from the ice wall at which the ice temperature was maintained at its initial temperature. A large thermal supply is required to heat 3000 m of cable in an ice sheet (scenario (iii) above), but penetration into glacier ice (scenarios (i) and (ii) above) could be realistic with the use of a currently employed generator.     

Interpretation of the GRACE-derived mass trend in Enderby Land, Antarctica

Monthly gravity solutions of the Gravity Recovery and Climate Experiment (GRACE) reveal three areas in Antarctica with striking interannual mass trends. The positive mass trend in Enderby Land, East Antarctica, is poorly understood because of uncertainties in the surface ice-sheet mass balance, post-glacial rebound (PGR), and processing of GRACE data. In this study, we compare the GRACE mass trends with values estimated from in situ snow-stake measurements, and Ice Cloud and land Elevation Satellite (ICESat) data. The mass trends estimated from ICESat data show a strong correlation with GRACE mass trends. In contrast, the snow-stake data show discrepancies with temporal variations in GRACE mass, especially in 2006. The discrepancies are probably associated with basal ice-sheet outflow, which is difficult to observe using snow stakes. We conclude that the bulk of the GRACE mass trend can be explained by snow accumulation and basal ice-sheet outflow.     

Quaternary paleolake formation and cataclysmic flooding along the upper Yenisei River

A suite of geomorphological and sedimentological features in the catchment of the upper Yenisei River in the Sayan mountains of southern Siberia testifies to the occurrence of cataclysmic floods that flowed down the river. Evidence of large-scale high-energy flood events includes: 1) gravel dunes, up to a few meters high and spaced 50 to 80 m apart, in the Kyzyl Basin 2) landforms such as hanging valleys and paleochannels and 3) flood sediments in a tributary valley. The origins of the Yenisei floods were likely diverse due to complex hydrological processes operating in the Sayan mountains. The possibilities include failures of multiple, variably impounded (ice, sedimentary, tectonic scarp, and lava flow dams) paleolakes in the two large intermontane basins of Darkhadyn Khotgor and Todza, and other minor basins, in the upper Yenisei River catchment. Dating techniques applied to the paleolakes in the Darkhadyn Khotgor and Todza basins revealed their formation during various periods in the middle–late Pleistocene and Holocene. Flooding from the Darkhadyn Khotgor appears to explain many of the inferred flood features, although contributions by flooding from other paleolake basins cannot be ruled out. Computer simulation of the flooding caused by a Darkhadyn Khotgor paleolake ice-dam failure indicates a probable peak discharge of  3.5 × 10⁶ m³ s^− 1, approximately one-fifth that of the floods that formed the Channeled Scabland in the U.S.A. Many of the outburst events probably occurred in the late Quaternary, but earlier floods could also have occurred.     

Assessing flood disaster impacts in agriculture under climate change in the river basins of Southeast Asia

Natural Hazards - This study focused on flood damage assessment for future floods under the impact of climate change. Four river basins of Southeast Asia were selected for the study. They included...     

Effects of Bathymetric Grid-Cell Sizes on Habitat Suitability Analysis of Cold-water Gorgonian Corals on Seamounts

We investigated the effects of the spatial resolution of bathymetry data on the predicted local distribution of large gorgonian corals on the southern Emperor Seamounts in the high seas of North Pacific. Terrain attributes were calculated from bathymetric survey data at six different grid-cell sizes (25–800 m) to investigate how the values of bottom terrain attributes and the prediction of MaxEnt habitat suitability models would vary with grid-cell size. The values of each terrain attribute showed different patterns among different grid-cell sizes, and MaxEnt demonstrated better performance at grid-cell size of 100 m or less than at larger sizes.