Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

Hydrothermal alteration and tectonic evolution of an intermediate- to fast-spreading back-arc oceanic crust: Late Ordovician Solund-Stavfjord ophiolite, western Norway

Abstract The Solund‐Stavfjord ophiolite complex (SSOC) in western Norway represents a remnant of the Late Ordovician oceanic lithosphere, which developed in an intermediate‐ to fast‐spreading Caledonian back‐arc basin. The internal architecture and magmatic features of its crustal component suggest that the SSOC has a complex, multistage sea floor spreading history in a supra‐subduction zone environment. The youngest crustal section associated with the propagating rift tectonics consists of a relatively complete ophiolite pseudostratigraphy, including basaltic volcanic rocks, a transition zone between the sheeted dyke complex and the extrusive sequence, sheeted dykes, and high‐level isotropic gabbros. Large‐scale variations in major and trace element distributions indicate significant remobilization far beyond that which would result from magmatic processes, as a result of the hydrothermal alteration of crustal rocks. Whereas K₂O is strongly enriched in volcanic rocks of the extrusive sequence, Cu and Zn show the largest enrichment in the dyke complex near the dyke–volcanic transition zone or within this transition zone. The δ¹⁸O values of the whole‐rock samples show a general depletion structurally downwards in the ophiolite, with the largest and smallest variations observed in volcanic rocks and the transition zone, respectively. δ¹⁸O values of epidote–quartz mineral pairs indicate 260–290°C for volcanic rocks, 420°C for the transition zone, 280–345°C for the sheeted dyke complex and 290–475°C for the gabbros. The ⁸⁷Sr/⁸⁶Sr isotope ratios show the widest range and highest values in the extrusive rocks (0.70316–0.70495), and generally the lowest values and the narrowest range in the sheeted dyke complex (0.70338–0.70377). The minimum water/rock ratios calculated show the largest variations in volcanic rocks and gabbros (approximately 0–14), and generally the lowest values and range in the sheeted dyke complex (approximately 1–3). The δD values of epidote (?1 to ?12‰), together with the δ¹⁸O calculated for Ordovician seawater, are similar to those of present‐day seawater. Volcanic rocks experienced both cold and warm water circulation, resulting in the observed K₂O‐enrichment and the largest scatter in the δ¹⁸O values. As a result of metal leaching in the hot reaction zone above a magma chamber, Zn is strongly depleted in the gabbros but enriched in the sheeted dyke complex because of precipitation from upwelling of discharged hydrothermal fluids. The present study demonstrates that the near intact effect of ocean floor hydrothermal activity is preserved in the upper part of the SSOC crust, despite the influence of regional lower greenschist facies metamorphism.     

Seal capacity estimation from subsurface pore pressures

A cap rock's capacity to seal hydrocarbons depends on its wettability and the sizes of the pore throats within the interconnected pore system that the leaking hydrocarbons must penetrate. These critical pore throat sizes are often poorly constrained in hydrocarbon exploration, partly because measurements of pore throat sizes have not been performed, and partly because pore throat measurements on a few individual samples in the cap rock may not be representative for the seal capacity of the top seal as a whole. To the contrary, the presence of formation overpressure can normally be estimated in drilled exploration targets. The presence of overpressure in reservoirs testifies to small pore throats in the cap rocks, as large pore throats will result in sufficiently high cap rock permeability to bleed off the overpressure. We suggest a stepwise procedure that enables quantification of top seal capacities of overpressured traps, based on subsurface pressure information. This procedure starts with the estimation of cap rock permeabilities, which are consistent with observed overpressure gradients across the top seals. Knowledge of burial histories is essential for these estimations. Relationships between pore throat size and permeability from laboratory experiments are then applied to estimate critical pore throat diameters in cap rocks. These critical pore throat diameters, combined with information of the physical properties of the pore fluids, are then used to calculate membrane seal capacity of cap rocks. Estimates of top seal capacity based on this procedure are rather sensitive to the vertical fluid velocity, but they are also to some extent sensitive to inaccuracies of the pore throat/permeability relationship, overpressure gradient, interfacial tensions between pore fluids, hydrocarbon density and water viscosity values. Despite these uncertainties, applications of the above‐mentioned procedure demonstrated that the mere presence of reservoir overpressures testifies to sufficient membrane seal capacity of cap rocks for most geological histories. Exempt from this statement are basins with rapid and substantial sediment compaction in the recent past.     

Speckle interferometric observations of Pluto and Charon

We report speckle interferometric observations of Pluto and its moon (1978 P1) Charon obtained on 5 June 1980 with a single 1.8-m mirror of the Multiple Mirror Telescope. Our observations yield a separation of 0″.31 (±0″.05) between Pluto and Charon at position angle 285° (±7°) for JD 2444395.75. This result and other direct observations indicate an adjustment of +4.0 hr to the orbital epoch of R. S. Harrington and J. W. Christy [Astron.J.86, 442–443 (1981)]. Our observation, which represents the first resolution of the system near minimum separation, also suggests that the inclination of the orbit to the plane of the sky should be increased by 3°; this will delay the onset of the predicted eclipsee season by one apparition to 1984 or 1985. Our data are consistent with Pluto diameter 0″.14 (±0″.02) = 3000 (±400) km and Charon diameter 0″..05 (±0″.03) = 1100 (±600) km.     

Chemical and structural relations of epitaxial xenotime and zircon substratum in sedimentary and hydrothermal environments: a TEM study

Xenotime overgrowths on detrital zircon in siliciclastic sediments have been reported in numerous studies. However, in natural samples, solid solution of zircon and xenotime is limited to near-end-member compositions. In order to characterize the interface region between both minerals and to draw inferences on the growth mechanisms of authigenic xenotime, we studied xenotime overgrowths on detrital zircon grains from two Phanerozoic sandstone samples with contrasting post-depositional histories. In one sample, the small (≤10 μm), pyramidal xenotime overgrowths are of diagenetic origin and grew without major discontinuity on the detrital zircon grain. The second sample shows up to >50-μm-wide, porous and inclusion-rich, hydrothermal xenotime overgrowths on detrital zircon, whereas the transition zone between both minerals is accompanied by large pore volume. Chemical compositions of the xenotime precipitates from the two samples differ particularly in Y, REE, Th and Sc concentrations, whereas high MREE availability in the diagenetic sample and the presence of Sc in the hydrothermal sample, respectively, appear to have promoted xenotime growth. Transmission electron microscopy on electron-transparent foils cut from the interface region shows that both the diagenetic xenotime and the hydrothermal xenotime are crystalline and grew in optical and crystallographic continuity to their detrital zircon substrata. Only a narrow transition zone (≤90 nm—diagenetic sample, 200–300 nm—hydrothermal sample) between zircon and xenotime is in part made up of nanometre-scale crystalline domains that are slightly distorted and may have formed from dissolution–re-precipitation processes at the zircon rim along with precipitation from the respective fluid.     

Strain coupling between upper mantle and lower crust: natural example from the Běstvina granulite body, Bohemian Massif

Strain patterns within mantle rocks and surrounding coarse‐grained felsic granulites from the Kutná Hora Crystalline Complex in the Variscan Bohemian Massif have been studied in order to assess their strain coupling. The studied rock association occurs within low‐strain domains surrounded by fine‐grained granulite and migmatite. The Doubrava peridotite contains closely spaced and steeply dipping layers of garnet clinopyroxenite, which are parallel to the NE–SW‐striking, high‐temperature foliation in nearby granulites, while the Úhrov peridotite lacks such layering. The Spa?ice eclogite is not associated with peridotite and shows upright folds of alternating coarse‐ and fine‐grained varieties bearing NE–SW‐striking axial planes. Electron back‐scattered diffraction measurements revealed full strain coupling between clinopyroxenites and coarse‐grained granulites in the S₁ fabric that is superposed on the S₀ fabric preserved in peridotites. The B‐type olivine lattice preferred orientation (LPO) characterizes the S₀ fabric in peridotites and its reworking is strongly controlled by the presence of macroscopic clinopyroxenite layering. The S₁ in clinopyroxenites and coarse‐grained granulites is associated with the LS‐type clinopyroxene LPO and prism <c> slip in quartz respectively. While the S₁ fabric in these rock types is accompanied invariably by a sub‐vertical stretching lineation, the S₁ fabric developed in reworked Úhrov peridotite is associated with strongly planar axial (010) type of olivine LPO. The peridotites with the S₀ fabric are interpreted to be relicts of a fore‐arc mantle wedge hydrated to a various extent above the Saxothuringian subduction zone. The prograde metamorphism recorded in peridotites and eclogites occurred presumably during mantle wedge flow and was reaching UHP conditions. Strain coupling in the S₁ fabric between clinopyroxenites and granulites at Doubrava and upright folding of eclogites at Spa?ice document a link between tectonic and magmatic processes during orogenic thickening, coeval with intrusions of the arc‐related calcalkaline suites of the Central Bohemian Plutonic Complex (c. 360–345 Ma). Juxtaposition of peridotites and granulites could be explained by a rheological heterogeneity connected to the development of clinopyroxenite layering in the upper mantle and a previously published model of a lithospheric‐scale transpressional arc system. It invokes vertical shearing along NE–SW trending, sub‐vertical foliations in the upper mantle that could have led to an emplacement of mantle bodies into the granulitized, orogenic root in the sub‐arc region. Clearly, such a transpressional arc system could represent an important pathway for an emplacement of deep‐seated rocks in the orogenic lower crust.     

Estimation of regional meteorological and hydrological drought characteristics: a case study for Denmark

Information on regional drought characteristics provides critical information for adequate water resource management. This study introduces a method to calculate the probability of a specific area to be affected by a drought of a given severity and demonstrates its potential for calculating both meteorological and hydrological drought characteristics. The method is demonstrated using Denmark as a case study. The calculation procedure was applied to monthly precipitation and streamflow series separately, which were linearly transformed by the Empirical Orthogonal Functions (EOF) method. Denmark was divided into 260 grid-cells of 14×17 km, and the monthly mean and the EOF-weight coefficients were interpolated by kriging. The frequency distributions of the first two (streamflow) or three (precipitation) amplitude functions were then derived. By performing Monte Carlo simulations, amplitude functions corresponding to 1000 years of data were generated. Based on these simulated functions as well as interpolated mean and weight coefficients, long time series of precipitation and streamflow were simulated for each grid-cell. The probability distribution functions of the area covered by a drought and the drought deficit volumes were then derived and combined to produce drought severity-area-frequency curves. These curves allowed an estimation of the probability of an area of a certain extent to have a drought of a given severity, and thereby return periods could be assigned to historical drought events. A comparison of drought characteristics showed that streamflow droughts are less homogeneous over the region, less frequent and last for longer time periods than precipitation droughts.     

Searching for planets by differential astrometry with large telescopes

Traditional astrometric methods are limited in accuracy by the atmosphere in a way that does not show much improvement with increased telescope aperture. However, there is the potential for very high accuracy with large telescopes if advantage can be taken of these factors: First, the differential atmospheric distortion of images of closely adjacent stars is less with larger aperture; second, the diffraction limit is sharper, and third, photon statistics are improved. In this paper we analyze and give experimental tests of techniques that could be applied to the detection of planets with the mass of Jupiter or Uranus, if they are present in nearby binary star systems.The atmospheric perturbation of the relative position of the energy centroids measured in short exposure images of binary stars depends on the effective height of the turbulent distortion. For a 4-meter telescope, the error in centroid determination of a 4-arcsec binary can be as small as 20 milliarcsec (mas) in a single 20-millisecond (msec) exposure. The relative position measured by cross-correlation of short exposure speckle images, as suggested by McAlister (1977b), may give even higher accuracy. In this case, Roddier (Roddieret al., 1980) has shown that the atmospheric error depends on the thickness rather than the height of the layers that make the dominant contribution to the turbulence. Through Monte Carlo analysis we show that on occasions when the turbulence arises largely in a thin layer, a single 20-msec exposure of a 4-arcsec binary taken with a 4-m aperture can yield an astrometric accuracy of order 0.5 mas.We report on experiments made at the Steward Observatory 2.3-m telescope which achieved accuracies corresponding to 1.7 mas in a 2.24-arcsec binary and 16.1 mas in a 6.0-arcsec binary with only 15 and 18 specklegram pairs respectively. We plan to use the 6.5-m converted MMT to obtain much higher performance, between 4.0 mas and 0.40 masper independent specklegram pair, depending upon atmospheric conditions, for binaries of 4-arcsec separation. By cycling rapidly through perhaps 100 binaries, thus calibrating systematic errors through the average change in binary separation, Jupiter-mass planets may be detectable with small but regular access to the telescope.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

Climate Change Impact on Hydrological Extremes: Preliminary Results from the Polish-Norwegian Project

This paper presents the background, objectives, and preliminary outcomes from the first year of activities of the Polish–Norwegian project CHIHE (Climate Change Impact on Hydrological Extremes). The project aims to estimate the influence of climate changes on extreme river flows (low and high) and to evaluate the impact on the frequency of occurrence of hydrological extremes. Eight “twinned” catchments in Poland and Norway serve as case studies. We present the procedures of the catchment selection applied in Norway and Poland and a database consisting of near-natural ten Polish and eight Norwegian catchments constructed for the purpose of climate impact assessment. Climate projections for selected catchments are described and compared with observations of temperature and precipitation available for the reference period. Future changes based on those projections are analysed and assessed for two periods, the near future (2021–2050) and the far-future (2071–2100). The results indicate increases in precipitation and temperature in the periods and regions studied both in Poland and Norway.