About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

Deep-ocean vertical noise directionality

The structure of beam noise measured at the output of a vertical array in a range dependent ocean basin was investigated using the modified wide-angle parabolic equation (PE). Noise sources were distributed throughout the basin, and the field due to each noise source at an array located in the midbasin was calculated. The response of the array to the superposition of the noise sources was found by beamforming. An efficient and direct approach that superimposes the noise sources on the PE field as the field is marched toward the array was developed. Downslope calculations of the midbasin vertical directionality were made between 50 and 400 Hz with this technique. Use of a geoacoustic model shows that the bottom behaves as a low-pass filter     

Inverted echo sounder measurement of dynamic height through an ENSOcycle in the central equatorial Pacific

A four-year record from an inverted echo sounder deployed near Palmyra Island at 6°N in the central Pacific Ocean is compared with a simultaneous record of subsurface pressure from this island lagoon. A factor m, converting round-trip acoustic travel time to surface dynamic height relative to a deep pressure level, was estimated from the ratio of the spectra of the two records in the energetic synoptic oscillation band. Year-to-year variation in m was not statistically significant. For the overall record, m was found to be -70±8 dynamic m/s, where the error bounds represent a 90% confidence interval. This is consistent with first-baroclinic-mode excitation     

A stepwise method for determining the number of component distributions in a mixture

This paper considers the problem of estimatingm, the number of components in a finite mixture of distributions from a parametric family. A step-up procedure using the bootstrap method is proposed. Some properties of the procedure are illustrated with simulation studies. An example of the method, applied to orientation of beach clasts, is given.     

Three-dimensional linear solution for wave propagation with slopingbottom

This paper presents a three-dimensional analytic linear wave solution for surface gravity wave propagation over a sloping bottom that is valid for small, but realistic, slopes. The sloping-bottom linear model is compared to published laboratory data and to predictions of two-dimensional, constant-bottom nonlinear theories. The model is shown to describe the measured wave-height growth in the wave transformation region up to a limiting local Ursell number U_r of 0.35-1.0, depending on the wave type, although, as a linear model, it does not predict the harmonics observed in that range. For U_r<0.35, the harmonics can generally be neglected and the sloping-bottom linear theory agrees closely with both the published wave-height data and third-order Stokes nonlinear theory. As a three-dimensional linear model, superposition can be invoked to synthesize and relate wave structure in the transformation region to complex incident ocean spectra with both wind wave and swell components that arrive with a range of incidence angles. As such, the sloping-bottom linear model presented here should be a convenient useful tool for ocean modeling through a significant portion of the wave transformation region     

A time series analysis of sound propagation in a strongly multipathshallow water environment with an adiabatic normal mode approach

Measured time series were generated by small omnidirectional explosive sources in a shallow water area. A bottom-mounted hydrophone recorded sound signals that propagated over a sloping bottom. The time series in the 250-500 Hz band were analyzed with a broad-band adiabatic normal mode approach. The measured waveforms contain numerous bottom interacting multipaths that are complicated by the subbottom structure that contains high-velocity layers near the water-sediment interface. Several of the details of the geoacoustic structure and the depth of the water column at the receiver are inferred from comparisons of the measured data to simulated time series. The sensitivity of broad-band matched-field ambiguity surfaces in the range-depth plane for a single receiver to selected waveguide parameters is examined. A consistent analysis is made where the simulated time series are compared to the measured time series along with the single-receiver matched-field localization solutions for ranges out to 5 km. In this range interval, it was found that the peak cross-correlation between the measured and simulated time series varied between 0.84 and 0.69. The difference between the GPS range and the range obtained from the matched-field solution varied from 0 to 63 m. The geoacoustic structure obtained in the analysis consists of an 8-m low-velocity sediment layer over an 8-m high-velocity layer followed by a higher velocity, infinite half-space     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.     

Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997     

Depositional and post‐depositional history of warm stage deposits at Knocknacran,Co. Monaghan,Ireland: implications for preservation of Irish last interglacial deposits

Organic‐rich deposits, uncovered during overburden removal from mantled gypsum karst at Knocknacran opencast gypsum mine, Co. Monaghan, are the best candidate to date for a last interglacial record in Ireland. The two till and organic‐rich deposits (preserved at different quarry elevations) were emplaced on to a Tertiary dolerite surface during high‐energy flood events and subsequently folded and faulted by movement towards sinkholes in underlying gypsum. Uranium–thorium disequilibrium dating suggests that the organic‐rich deposits in the upper section were hydrologically isolated at ca. 41 ka and those in the lower section at ca. 86 ka. Interpretation of the pollen content, although tentative because of the depositional and post‐depositional history of the material, suggests that the organic material originated in a warm stage possibly warmer than the post‐Eemian interstadials. The unusual setting of preservation may indicate that in situ, last interglacial deposits have generally been removed by erosion in Ireland. Copyright © 2004 John Wiley & Sons, Ltd.