Self-similar solutions for the blast wave problem in radiative gasdynamics,I

Similarity solutions, describing the flow of a perfect gas behind spherical shock waves, are investigated including the radiation heat flux. The shock is assumed to be propagating in a medium at rest. Shock radius varies exponentially with time and density is inversely proportional to fifth power of the shock radius immediately ahead of the shock front.     

High-frequency monitoring of hydrological and biogeochemical fluxes in forested catchments of southern Chile

The variability of rainfall-dependent streamflow at catchment scale modulates many ecosystem processes in wet temperate forests. Runoff in small mountain catchments is characterized by a quick response to rainfall pulses which affects biogeochemical fluxes to all downstream systems. In wet-temperate climates, water erosion is the most important natural factor driving downstream soil and nutrient losses from upland ecosystems. Most hydrochemical studies have focused on water flux measurements at hourly scales, along with weekly or monthly samples for water chemistry. Here, we assessed how water and element flows from broad-leaved, evergreen forested catchments in southwestern South America, are influenced by different successional stages, quantifying runoff, sediment transport and nutrient fluxes during hourly rainfall events of different intensities. Hydrograph comparisons among different successional stages indicated that forested catchments differed in their responses to high intensity rainfall, with greater runoff in areas covered by secondary forests (SF), compared to old-growth forest cover (OG) and dense scrub vegetation (CH). Further, throughfall water was greatly nutrient enriched for all forest types. Suspended sediment loads varied between successional stages. SF catchments exported 455 kg of sediments per ha, followed by OG with 91 kg/ha and CH with 14 kg/ha, corresponding to 11 rainfall events measured from December 2013 to April 2014. Total nitrogen (TN) and phosphorus (TP) concentrations in stream water also varied with rainfall intensity. In seven rainfall events sampled during the study period, CH catchments exported less nutrients (46 kg/ha TN and 7 kg/ha TP) than SF catchments (718 kg/ha TN and 107 kg/ha TP), while OG catchments exported intermediate sediment loads (201 kg/ha TN and 23 kg/ha TP). Further, we found significant effects of successional stage attributes (vegetation structure and soil physical properties) and catchment morphometry on runoff and sediment concentrations, and greater nutrients retention in OG and CH catchments. We conclude that in these southern hemisphere, broad-leaved evergreen temperate forests, hydrological processes are driven by multiple interacting phenomena, including climate, vegetation, soils, topography, and disturbance history.     

A numerical study of the nocturnal heat island over a medium-sized mid-latitude city (Columbus,Ohio)

A numerical investigation is conducted of the nocturnal heat island over Columbus, Ohio, a relatively flat mid-latitude city. Specifically, a cross-sectional steady-state numerical model to simulate the (thermal) structure of the nocturnal urban boundary layer is developed from a one-dimensional, time-dependent model due to Estoque. The model is applied to Columbus for special periods in September 1968 and March 1969 during which comparable experimental data were collected. The numerical simulations agree well with the observed data with respect to the detailed spatial form of the surface-based thermal boundary layer across the city. The use of the model for other metropolitan areas, for ascertaining the relative effect of city size and building geometry on the development of the thermal boundary layer, and for determining the effect of alternate land-use strategies on the thermal stratification are also discussed.Digest of portion of a dissertation submitted to the Pennsylvania State University in partial fulfillment of the requirements for the degree of Doctor of Philosophy.Presented in part at the Air Pollution, Turbulence and Diffusion Symposium, Las Cruces, New Mexico, 7–10 December 1971.On assignment from the National Oceanic and Atmospheric Administration (NOAA), U.S. Department of Commerce.     

A preliminary investigation of siliceous microfossil succession in late Quaternary sediments from Lake Baikal, Siberia

Siliceous microfossil assemblage succession was analyzed in a 100 m sediment core from Lake Baikal, Siberia. The core was recovered from the lake's central basin at a water depth of 365 m. Microfossil abundance varied greatly within the intervals sampled, ranging from samples devoid of siliceous microfossils to samples with up to 3.49 × 10¹¹ microfossils g^-1 sediment. Fluctuations in abundance appear to reflect trends in the marine ¹⁸O record, with peak microfossil levels generally representing climate optima. Microfossil taxa present in sampled intervals changed considerably with core depth. Within each sample a small number of endemic diatom species dominated the assemblage. Changes in dominant endemic taxa between sampled intervals ranged from extirpation of some taxa, to shifts in quantitative abundance. Differences in microfossil composition and the association of variations in abundance with climate fluctuations suggest rapid speciation in response to major climatic excursions.     

Helium isotopes in ferromanganese crusts from the central Pacific Ocean

Helium isotopes have been measured in samples of two ferromanganese crusts (VA13/2 and CD29-2) from the central Pacific Ocean. With the exception of the deepest part of crust CD29-2 the data can be explained by a mixture of implanted solar- and galactic cosmic ray-produced (GCR) He, in extraterrestrial grains, and radiogenic He in wind-borne continental dust grains. ⁴He concentrations are invariant and require retention of less than 12% of the in situ He produced since crust formation. Loss has occurred by recoil and diffusion. High ⁴He in CD29-2 samples older than 42 Ma are correlated with phosphatization and can be explained by retention of up to 12% of the in situ-produced ⁴He. ³He/⁴He of VA13/2 samples varies from 18.5 to 1852 R_a due almost entirely to variation in the extraterrestrial He contribution. The highest ³He/⁴He is comparable to the highest values measured in interplanetary dust particles (IDPs) and micrometeorites (MMs). Helium concentrations are orders of magnitude lower than in oceanic sediments reflecting the low trapping efficiency for in-falling terrestrial and extraterrestrial grains of Fe-Mn crusts. The extraterrestrial ³He concentration of the crusts rules out whole, undegassed 4-40 μm diameter IDPs as the host. Instead it requires that the extraterrestrial He inventory is carried by numerous particles with significantly lower He concentrations, and occasional high concentration GCR-He-bearing particles.     

Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Numerical modeling of three‐phase dissolution of underground cavities using a diffuse interface model

Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.     

A major lava tube system from Undara Volcano,North Queensland

The Undara Volcano erupted 0.19 m.y. ago and formed lava fields covering 1,500 km² with a volume of approximately 23 km³. One of the flows extended 160 km on a gradient that averaged only 0.3°. This great length was a result of very high effusion rates, favourable topography and lava tube efficiency. The Undara lavas are rather uniform hawaiites. Lava temperatures are estimated to have been somewhat less than 1175–1220°C and viscosities greater than 10 to 30 Pa s. Long, apparently single lava tubes are well preserved in many places. They are marked by depressions, caves and long level ridges. A system of lava tubes extends for perhaps more than 100 km. The features of the lava tubes are comparable with those described elsewhere. Aligned depressions associated with caves appear to have formed contemporaneously. Most are much wider than the caves and probably represent collapsed lava ponds. The lava tubes appear to have formed by roofing over of lava channels. Close to lava tubes, the rocks developed strongly oxidised characteristics, such as oxidised olivine phenocrysts, ferric clinopyroxene and extensively developed hematite. Differentiated lava forms drips in some caves and is also oxidised.     

Attenuation relations for strong seismic ground motion in the Himalayan region

Strong motion data from various regions of India have been used to study attenuation characteristics of horizontal peak acceleration and velocity. The strong ground motion data base considered in the present work consists of various earthquakes recorded in the northern part of India since 1986 with magnitudes 5.7 to 7.2. Using these data, relations for horizontal peak acceleration and velocity, which are $$\begin{gathered} log_{10} a = 1.14 + 0.31M + 0.65log_{10} R \hfill \\ log_{10} v = 0.571 + 0.41M + 0.768log_{10} R \hfill \\ \end{gathered} $$ have been proposed wherea is the peak horizontal acceleration in cm/sec²,v is the peak horizontal velocity in mm/sec,M is body wave magnitude, andR is the hypocentral distance in km. The proposed relations are in reasonable agreement with the small amount of strong ground motion data available for the northern part of India. The present results will be useful in estimating strong ground motion parameters and in the earthquake resistant design in the Himalayan region.     

New Method for Continuous Transmissivity Profiling in Fractured Rock

A new method is presented to search for hydraulically transmissive features in open boreholes in bedrock. A flexible borehole liner made of a watertight, nylon fabric is filled with water to create a constant driving head to evert (reverse of invert) the liner down the hole so that the liner pushes the borehole water out into transmissive fractures or other permeable features. The descent rate is governed by the bulk transmissivity of the remaining permeable features below the liner. Initially, the liner descent rate or velocity is a measure of transmissivity (T) of the entire hole. As the everting liner passes and seals each permeable feature, changes in the liner velocity indicate the position of each feature and an estimate of T using the Thiem equation for steady radial flow. This method has been performed in boreholes with diameters ranging from 96 to 330 mm. Profiling commonly takes a few hours in holes 200‐ to 300‐m long. After arrival of the liner at the bottom of the hole, the liner acts as a seal preventing borehole cross connection between transmissive features at different depths. Liner removal allows the hole to be used for other purposes. The T values determined using this method in a dolostone aquifer were found to be similar to the values from injection tests using conventional straddle packers. This method is not a replacement for straddle‐packer hydraulic testing of specific zones where greater accuracy is desired; however, it is effective and efficient for scanning entire holes for transmissive features.