Uncertainty as the overlap of alternate conditional distributions

An important task in modern geostatistics is the assessment and quantification of resource and reserve uncertainty. This uncertainty is valuable support information for many management decisions. Uncertainty at specific locations and uncertainty in the global resource is of interest. There are many different methods to build models of uncertainty, including Kriging, Cokriging, and Inverse Distance. Each method leads to different results. A method is proposed to combine local uncertainties predicted by different models to obtain a combined measure of uncertainty that combines good features of each alternative. The new estimator is the overlap of alternate conditional distributions.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Microwave emission from a sunspot

From the gyroresonance brightness temperature spectrum of a sunspot, one can determine the magnetic field strength by using the property that microwave brightness is limited above a frequency given by an integer-multiple of the gyrofrequency. In this paper, we use this idea to find the radial distribution of magnetic field at the coronal base of a sunspot in the active region, NOAA 4741. The gyroresonance brightness temperature spectra of this sunspot are obtained from multi-frequency interferometric observations made at the Owens Valley Radio Observatory at 24 frequencies in the range of 4.0–12.4 GHz with spatial resolution 2.2″–6.8″. The main results of present study are summarized as follows: first, by comparison of the coronal magnetic flux deduced from our microwave observation with the photospheric magnetic flux measured by KPNO magnetograms, we show that theo-mode emission must arise predominantly from the second harmonic of the gyrofrequency, while thex-mode arises from the third harmonic. Second, the radial distribution of magnetic fieldsB(r) at the coronal base of this spot (say, 2000–4000 km above the photosphere) can be adequately fitted by $$B(r) = 1420(1 \pm 0.080)\exp \left[ { - \left( {\frac{r}{{11.05''(1 \pm 0.014)}}} \right)^2 } \right]G,$$ wherer is the radial distance from the spot center at coronal base. Third, it is found that coronal magnetic fields originate mostly from the photospheric umbral region. Fourth, although the derived vertical variation of magnetic fields can be approximated roughly by a dipole model with dipole moment 1.6 × 10³⁰ erg G^?1 buried at 11000 km below the photosphere, the radial field distribution at coronal heights is found to be more confined than predicted by the dipole model.     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

A decade of change in the Skidaway River estuary I. Hydrography and nutrients

The Skidaway River estuary is a tidally-dominated subtropical estuary in the southeastern USA surrounded by extensiveSpartina salt marshes. Weekly smapling at high and low tide began in 1986 for hydrography, nutrients, chlorophylla, particulate matter, and microbial and plankton biomass and composition; hydrographic and nutrient data during 1986–1996 are reported here. Salinity varied inversely with river discharge and exhibited variability at all time scales but with no long-term trend. Water temperature typically ranged over 25°C and was without apparent long-term frend. Seasonal cycles in concentrations of NO₃, NH₄, PO₄, Si(OH)₄, and DON were observed, with annual maxima generally occurring in late summer. Superimposed on seasonal cycles, all five nutrients exhibited steady increases in minimum, mean, and maximum concentrations; mean concentrations increased c. 50–150% during the decade. Nutrient concentrations were highly correlated with water temperature over the ten-year period, but weakly related to salinity and discharge. Nutrients were strongly correlated with one another, and the relative ratios among inorganic nutrients showed little long-term trend. Correlations among temperature and nutrient concentrations exhibited considerable inter-annual variability. Major spikes in organic and inorganic nutrient concentrations coincided with significant rainfall events; concentrations increased hyperbolically with rainfall. Although pristine compared to more heavily impacted waterways primarily outside the region, residential development and population density have been increasing rapidly during the past 15–20 years. Land use is apparently altering nutrient loading over the long-term (months-years), and superimposed on this are stochastic meteorological events that accelerate these changes over the short term (days-weeks).     

Kosa (yellow sand) components in precipitation collected at central Japan

In east Asia, acidic gases derived from fossil fuel combustion have increased in the past decades. On the other hand, the Asian dust, also called Kosa (yellow sand) is transported following windstorms from arid lands in the Asian continent. Many researchers have been interested in the reaction between acidic aerosols and Kosa aerosols as well as the long-range transport of these emissions. To investigate the characteristics of chemical components in precipitation on a long-term basis over Japan, precipitation was sequentially collected from April 1984 to March 1997 at Kanazawa located near the coast of the Sea of Japan. Precipitation samples were collected at 1 mm intervals for the first 5 mm rainfall and all volume of rainwater after 6 mm for all precipitation events with an automatic wet only precipitation collector. According to the analyses of precipitation including Kosa aerosols during Kosa periods, the reaction in the air between Kosa and acidic components during the long-range transport was discussed.     

Getting ‘holistic’ about land rehabilitation in hot arid regions

The paper describes Western Australian examples and causes of land degradation. It outlines shortcomings in the methodologies used to rehabilitate these areas. From this a protocol is suggested for an ‘holistic’ approach to land rehabilitation.     

Fine structure of Pn velocity beneath Sichuan-Yunnan region

We use 23298 Pn arrival-time data from Chinese national and provincial earthquake bulletins to invert fine structure of Pn velocity and anisotropy at the top of the mantle beneath the Sichuan-Yunnan and its adjacent region. The results suggest that the Pn velocity in this region shows significant lateral variation; the Pn velocity varies from 7.7 to 8.3 km/s. The Pn-velocity variation correlates well with the tectonic activity and heat flow of the region. Low Pn velocity is observed in southwest Yunnan, Tengchong volcano area, and the Panxi tectonic area. These areas have very active seismicity and tectonic activity with high surface heat flow. On the other hand, high Pn velocity is observed in some stable regions, such as the central region of the Yangtze Platform; the most pronounced high velocity area is located in the Sichuan Basin, south of Chengdu. Pn anisotropy shows a complex pattern of regional deformation. The Pn fast direction shows a prominent clockwise rotation pattern from east of the Tibetan block to the Sichuan-Yunnan diamond block to southwest Yunnan, which may be related to southeastward escape of the Tibetan Plateau material due to the collision of the Indian Plate to the Eurasia Plate. Thus there appears to be strong correlation between the crustal deformation and the upper mantle structure in the region. The delay times of events and stations show that the crust thickness decreases from the Tibetan Plateau to eastern China, which is consistent with the results from deep seismic sounding.