A quadratic element method for evaluating groundwater flow by the inversion of surface tilt with application to the Tono Area, Japan

By modifying a previous method with constant elements, we developed a quadratic element method for more accurately estimating groundwater flow by the inversion of tilt data. In this method: (1) a region of groundwater flow is divided into quadratic elements in which the change in groundwater volume per unit volume of rock (Δv) and the Skempton coefficient (B) vary in a quadratic manner with the coordinates, (2) the values of Δv are set to zero at the boundaries of the region of groundwater flow and (3) the sum of the squared second derivatives of Δv is adopted as a constraining condition that is weighted and added to the sum of the squared errors in tilt. First, analyses were performed for a flow model to determine the accuracy of this method for estimating groundwater flow and also to clarify the effect of the assumed size of a region of groundwater flow. These analyses showed that the quadratic element method proposed in this study gives a much better estimation of Δv than the constant element method and that a large region of groundwater flow should be assumed, rather than a small region, since the values of Δv at points outside of the actual region of groundwater flow are estimated to be nearly zero when a large region is assumed while these values are greatly overestimated when an excessively small region is assumed. Finally, the quadratic element method was applied to the site of the Mizunami Underground Research Laboratory in the Tono area, Japan. Inverse analyses were performed for tilt data measured by four tiltmeters with a resolution of 10⁻⁹ radians during the excavation of two shafts under the assumption that the rock mass is an isotropic and homogeneous half- space. The results showed that the method proposed in this study reproduced the tilt data very accurately. Thus, the distribution of Δv was estimated without sacrificing the reproducibility of the tilt data. The contour maps of B(1 + ν)Δv (ν: Poisson’s ratio) showed that the heterogeneous flow of groundwater occurred at the site and that groundwater volume decreased mainly in the area surrounded by two faults. The latter result is consistent with the finding obtained by previous investigations that these faults have low permeability in the direction perpendicular to the strike and may act as a flow barrier.     

Determining the dilation factor in 4D monitoring of compacting reservoirs by rock-physics models

Hydrocarbon depletion and fluid injection cause compaction and stretching of the reservoir and overburden layers. 4D prestack seismic data can be used to detect these changes because compaction/stretching causes changes in traveltimes and seismic velocities. We show that, by using two different petro‐elastic models at varying effective pressures, a good approximation is to assume that the fractional changes in layer thickness, ΔL/L, and seismic velocity, Δv/v, are related by a linear function of ΔL/L. The slope of this function (the dilation factor, α= (Δv/v)/(ΔL/L) ) is negative and its absolute value generally decreases (shale, low porosity) or increases (sandstone, high porosity) with increasing layer thickness and decreasing effective pressure. The analysis is mainly performed for isotropic deformations. The dilation factor for uniaxial deformations is smaller in absolute value. The dilation factor, which can be calculated from time‐lapse data, can be used to predict reservoir compaction/stretching as a function of depth and surface subsidence.     

Diffraction ofSH-waves originating from a moving point source by a rigid quarter space

Summary ApplyingWiener-Hopf techniques the diffraction ofSH-waves originating from a moving point source due to the presence of a step change in the elevation of the surface of a homogeneous half-space have been obtained. An interesting feature of resonance associated with the motion of a load at a velocityv ₀ ^–1 approachingv ₀ ^–1 has been observed.     

Downslope flow models of a Bingham liquid: Implications for lava flows

It is widely recognized that lavas behave as Bingham liquids, which are characterized by a yield stress σ and a plastic viscosity η. We consider two models describing downslope flows of a Bingham liquid with different aspect ratios A (= flow height/flow width): model 1 with A 1 and model 2 with A ≈ 1. Sufficiently uphill with respect to the front, such flows can be considered as laminar and locally isothermal. For both models, we obtain analytically the steady-state solution of the Navier-Stokes equations and the constitutive equation for a Bingham liquid. We study the flow height and velocity as functions of flow rate, rheological parameters and ground slope. It is found that such flows remain in the Newtonian regime at low yield stresses (σ 10³dyne/cm²), but the transition to the Bingham regime also depends on flow rate and occurs at higher values of σ for higher flow rates: for instance, a high aspect ratio flow (model 2) is still very close to the Newtonian regime at σ = 10⁴ dyne/cm², if the flow rate is greater than 10⁵ g/s. In the Bingham regime, flow heights are generally greater and flow velocities are smaller than in the Newtonian regime; moreover, flow heights are independent of flow rate, so that a change in flow rate results exclusively in a velocity change. After assuming a specific temperature dependence of σ and η between the solidus and the liquidus temperatures of an ideal Bingham liquid (1000°C and 1200 °C respectively), flow heights and velocities are examined as functions of temperature along the flow. Several effects observed in lava flows are predicted by these models and allow a more quantitative insight into the behaviour of lava flows.     

Scattering attenuation and intrinsic absorption using uniform and depth dependent model – Application to full seismogram envelope recorded in Northern Chile

Two seismic wave attenuation factors, scatteringattenuation Q _s ^-1 and intrinsicabsorption Q _i ^-1 are measured using theMultiple Lapse Time Window (MLTW) analysis method forthree different frequency bands, 1–2, 2–4, and 4–8 Hz.Data from 54 temporally deployed seismic stationslocated in northern Chile are used. This methodcompares time integrated seismic wave energies withsynthetic coda wave envelopes for a multiple isotropicscattering model. In the present analysis, the waveenergy is assumed to decay with distance in proportionto1/GSF·exp(- (Q _s ^-1+Q _i ^-1)· r/v), where r, and v are the propagationdistance, angular frequency and S wave velocity,respectively, and GSF is the geometricalspreading factor. When spatial uniformity of Q _s ^-1, Q _i ^-1 and v isassumed, i.e. GSF = 4r ², theestimates of the reciprocal of the extinction length,L _e ^-1 (= (Q _s ^-1+Q _i ^-1)·/v), are 0.017,0.012 and 0.010 km^-1, and those of the seismicalbedo, B ₀ (= Q _s ^-1/ (Q _s ^-1+Q _i ^-1)), are 0.48, 0.40and 0.34 for 1–2, 2–4 and 4–8 Hz, respectively, whichindicates that scattering attenuation is comparable toor smaller than intrinsic absorption. When we assumea depth dependent velocity structure, we also findthat scattering attenuation is comparable to orsmaller than intrinsic absorption. However, since thequantitative estimates of scattering attenuationdepend on the assumed velocity structure (strength ofvelocity discontinuity and/or Moho depth), it isimportant to consider differences in velocitystructure models when comparing attenuation estimates.     

Hydrothermal silica chimney fields in the Galapagos Spreading Center at 86°W

Silica chimneys were discovered in 1985 at 86°W in the rift valley of the Galapagos Spreading Center at 2600 m depth (“Cauliflower Garden”). The inactive chimneys lack any sulfides and consist almost entirely of amorphous silica (up to 96 wt.% SiO₂, opal-A); Fe and Mn oxides are minor constituents. Oxygen isotope data show that formation of the silica chimneys took place at temperatures between 32°C (+29.9‰ δ¹⁸O) and 42°C (+27.8‰ δ¹⁸O).Th/Udating reveals a maximum age of 1440 ± 300y. Amorphous silica solubility relations indicate that the silica chimneys were formed by conductive cooling of pure hydrothermal fluids or by conductive cooling of a fluid/seawater mixture. Assuming equilibrium with quartz at 500 bars, initial fluid temperatures of more than 175°C (i.e., a concentration of > 182 ppm SiO₂) were required to achieve sufficient supersaturation for the deposition of amorphous silica at 40°C and 260 bars. If the silica chimneys originate from the same or a similar fluid as higher-temperature ( < 300°C) sulfide-silica precipitates found nearby (i.e., 2.5 km away), then subsurface deposition of sulfides may have occurred.     

Model of vibrational level populations of Herzberg states of oxygen molecules at heights of the lower thermosphere and mesosphere

The paper presents a model of the kinetics of electronically excited O₂(c¹Σ_u⁻,v), O₂(A′³Δ_u,v), O₂(A³Σ_u⁺,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O₂(A³Σ_u⁺,v) and O₂(A′³Δ_u,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.     

An anomalous subauroral red arc on 4 August, 1972: comparison of ISIS-2 satellite data with numerical calculations

This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N⁺₂, O⁺₂, O⁺(⁴S), O⁺(²D), O⁻(²P), and O⁺(²P) ions, and O(¹D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm⁻³ s⁻¹ was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N₂(v > 0) and O₂(v > 0) in the calculations of the O⁺(⁴S) loss rate improves the agreement between the calculated N_e and the data on 4 August, 1972. The N₂(v > 0) and O₂(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O⁺ + N₂ rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O⁺ + N₂ loss rate due to vibrationally excited O₂ produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O₂(a¹δ_g) and O₂(b¹σ_g⁺), and rotational excitation of O₂, and found that the effect of these cooling rates on T_e can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(³P) with thermal electrons with excitation of O(¹D), and the electron cooling rates by vibrational excitation of O₂ and N₂ are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N₂(B³π_g⁺), N₂(C³π_u), and N₂(A³σ_u⁺) at several vibrational levels, O(¹S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm.     

Radioactive heat production in rocks and its relation to other petrophysical parameters

Summary Radioactive heat productionA is a scalar and isotropic petrophysical property independent of in situ temperature and pressure. Its value is usually expressed in HGU units (1 HGU=10^–13 cal/cm³ sec) and depends on the amounts of uranium, thorium and potassium.A varies with rock type over several orders of magnitude and reflects the geochemical conditions during rock formation (magmatic differentiation, sedimentation or metamorphism).In order to assign realistic thermal parameters to deeper-seated rocks correlations with seismic velocity (which can be determined from the surface) have been looked for. In the range characteristic for crystalline rocks of the crust (5–8 km/sec)A is strongly correlated with density and compressional wave velocityv _p:A decreases with increasingv _p orp. From this relationship it is now possible to estimate heat production values for any particular layer of a crustal section from measured seismic velocities. Contrary to earlier belief there is, as shown by experimental determinations, no correlation between heat productionA and thermal conductivityK in igneous and metamorphic rocks. In sediments however, especially in sand/shale sequences, a correlation betweenK andA is most likely: increasing clay mineral content, characterized by increasingA, causes the decrease ofK in these rocks.Contribution No. 111, Institute of Geophysics, Swiss Federal Institute of Technology, Zurich, Switzerland.     

A fracture mechanics study of subcritical tensile cracking of quartz in wet environments

Load relaxation and cross-head displacement rate-change experiments have been used to establish log₁₀ stress intensity factor (K) versus log₁₀ crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K _Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m^–3/2 and 1.002±0.048 MN·m^–3/2, respectively.Subcritical crack growth at velocities from 10^–3 m·s^–1 to 10^–9 m·s^–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole^–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.     

The 1975 sub-terminal lavas, mount etna: a case history of the formation of a compound lava field

The 1975 sub-terminal activity was characterised by low effusion rates (0.3–0.5 m³ s⁻¹) and the formation of a compound lava field composed of many thousands of flow units. Several boccas were active simultaneously and effusion rates from individual boccas varied from about 10⁻⁴ to 0.25 m³s⁻¹. The morphology of lava flows was determined by effusion rate (E): aa flows with well-developed channels and levees formed when E > 2 × 10⁻³ m³ s⁻¹, small pahoehoe flows formed when 2 × 10⁻³ m³ s⁻¹ >E > 5 > 10⁻⁴ m³ s⁻¹ and pahoehoe toes formed when E < 5 × 10⁻⁴ m³ s⁻¹. There was very little variation with time in the effusion temperature, composition or phenocryst content of the lava.New boccas were commonly formed at the fronts of mature lava flows which had either ceased to flow or were moving slowly. These secondary boccas developed when fluid lava in the interior of mature aa flows either found a weakness in the flow front or was exposed by avalanching of the moving flow front. The resulting release of fluid lava was accompanied by either partial drainage of the mature flow or by the formation of a lava tube in the parent flow. The temperature of the lava forming the new bocca decreased with increasing distance from the source bocca (0.035°C m⁻¹). It is demonstrated from the rate of temperature decrease and from theoretical considerations that many of the Etna lavas still contained a substantial proportion of uncooled material in their interior as they came to rest. The formation of secondary boccas is postulated to be one reason why direct measurements of effusion rates tend, in general, to overestimate the total effusion rates of sub-terminal Etna lava fields.     

On the periodic variations of geomagnetic activity indices Ap and ap

Yearly averages of geomagnetic activity indices Ap for the years 1967–1984 are compared to the respective averages of v² · B_s, where v is the solar wind velocity and B_s is the southward interplanetary magnetic field (IMF) component. The correlation of both quantities is known to be rather good. Comparing the averages of Ap with v² and B_s separately we find that, during the declining phase of the solar cycle, v² and during the ascending phase B_s have more influence on Ap. According to this observation (using Fourier spectral analysis) the semiannual and 27 days, Ap variations for the years 1932–1993 were analysed separately for years before and after sunspot minima. Only those time-intervals before sunspot minima with a significant 27-day recurrent period of the IMF sector structure and those intervals after sunspot minima with a significant 28–28.5-day recurrent period of the sector structure were used. The averaged spectra of the two Ap data sets clearly show a period of 27 days before and a period of 28–29 days after sunspot minimum. Moreover, the phase of the average semiannual wave of Ap is significantly different for the two groups of data: the Ap variation maximizes near the equinoxes during the declining phase of the sunspot cycle and near the beginning of April and October during the ascending phase of the sunspot cycle, as predicted by the Russell-McPherron (R-M) mechanism. Analysing the daily variation of ap in an analogue manner, the same equinoctial and R-M mechanisms are seen, suggesting that during phases of the solar cycle, when ap depends more on the IMF-B_s component, the R-M mechanism is predominant, whereas during phases when ap increases as v increases the equinoctial mechanism is more likely to be effective.     

Unified modeling of flow liquefaction and cyclic mobility

Flow liquefaction and cyclic mobility are two different phenomena that are encountered during strong earthquakes. Flow liquefaction is associated with the contractive behavior of loose granular materials and cyclic mobility is associated with the dilative response of both loose and dense granular materials at low confining stresses. These two types of response pertain to the same material and therefore should be modeled in a consistent manner. Whether a soil in a given state exhibits contractive or dilative behavior is dependent on its dilatancy, d=d_v^p/d_q^p. A form of the dilatancy d=d(η,ψ,C) is proposed, where η=q/p is the stress ratio, and ψ and C denote the internal state variables and the intrinsic properties, respectively. It has been shown that such a state-dependent dilatancy is effective in describing both the contractive and dilative behavior. This allows the soil behavior associated with flow liquefaction and cyclic mobility during earthquakes to be modeled in a unified way.     

Assessing the seismic torsional vulnerability of elevated tanks with RC frame-type staging

Recurrence of torsional failure of elevated water tanks in past earthquakes (including 1952 Kern County and recent 1993 Killari earthquakes) has highlighted the importance of this problem. It is established that these structures may have amplified torsion-induced rotation if their torsional-to-lateral natural period ratio τ is close to 1 and amplified displacement of structural elements due to the coupled lateral-torsional vibration if τ is within the critical range 0.7<τ<1.25. The present study aims to estimate the range of variation of τ for usually constructed reinforced concrete elevated water tanks with frame-type stagings for assessing their torsional vulnerability. Closed-form expressions for torsional and lateral stiffness of tank stagings are derived and verified by standard finite element software. These expressions are used for studying the variation of τ for feasible ranges of influencing parameters. It is seen that a very large number of such tanks may have τ within the said critical range. Closed-form expressions for moments and shear forces of columns and beams under torsion and that under lateral force are also derived. It is also seen with the help of these expressions that the frame stagings of these tanks normally designed for seismic lateral force, may yield through formation of plastic hinges simultaneously in all columns instead of in beams if they are subjected to large rotational response for having τ possibly very close to 1. Such a pattern of yielding generally converts the whole system suddenly into a mechanism causing immediate collapse. Therefore, torsional coupling seems to be a potential cause of failure for these structures.     

Solar wind velocity distribution on the heliospheric current sheet during Carrington rotations 1787-1795

The solar wind velocity distribution in the heliosphere is best represented using a v-map, where velocity contours are plotted in heliographic latitude-longitude coordinates. It has already been established that low-speed regions of the solar wind on the source surface correspond to the maximum bright regions of the K-corona and the neutral line of the coronal magnetic field. In this analysis, v-maps on the source surface for Carrington rotations (CRs) 1787-1795, during 1987, have been prepared using the interplanetary scintillation measurements at Research Institute of Atmospherics (RIA), Nagoya Univ., Japan. These v-maps were then used to study the time evolution of the low-speed (\leq450 km s⁻¹) belt of the solar wind and to deduce the distribution of solar wind velocity on the heliospheric current sheet. The low-speed belt of the solar wind on the source surface was found to change from one CR to the next, implying a time evolution. Instead of a slow and systematic evolution, the pattern of distribution of solar wind changed dramatically at one particular solar rotation (CR 1792) and the distributions for the succeeding rotations were similar to this pattern. The low-speed region, in most cases, was found to be close to the solar equator and almost parallel to it. However, during some solar rotations, they were found to be organised in certain longitudes, leaving regions with longitudinal width greater than 30^○ free of low-speed solar wind, i.e. these regions were occupied by solar wind with velocities greater than 450 km s⁻¹. It is also noted from this study that the low-speed belt, in general, followed the neutral line of the coronal magnetic field, except in certain cases. The solar wind velocity on the heliospheric current sheet (HCS) varied in the range 300–585 km s⁻¹ during the period of study, and the pattern of velocity distribution varied from rotation to rotation.     

Weiteres zur unmittelbaren Berechnung der durchschnittlichen Geschwindigkeit in der Reflexionsseismik

Zusammenfassung Die Anordnung der Schusspunkte I, II ...,A, B ... und der Geophone 1, 2, ...a, b ... nach Figur 1 ermöglicht — horizontalen Reflektor vorausgesetzt — bei demselben Reflexionspunkt eine gute Bestimmung der durchschnittlichen Geschwindigkeiit und des Reflektors durch Ausgleichung mit den Gleichungen (3)-(6). Die angegebenen Gleichungen können auch bei geneigtem Reflektor verwendet werden, da sogar bei 11° Neigung des Reflektors der dadurch verursachte Fehler inv undN unterhalb 0.5% bleibt. Der Reflexionspunkt wandert in diesem Falle allerdings (vgl. Figur 2) mit dem Betrag r im Sinne der Gleichung (17a) weiter.Kennt man die Neigung der Schnittgeraden in der Reflexionsebene, so kann man mit Hilfe der Gleichungen (12) und (14) die genaueren Werte vonv undN ermitteln.Zur Bestimmung des Neigungswinkels wird man vorteilhaft die Anordnung nach Figur 3 treffen, wo bei einem Schuss in I/A die Geophone in 1, 2 ...,a, b ... angeordnet sind. Aus den Messergebnissen können wir durch Ausgleichung nach den Gleichungen (23) und (24) bestimmen.Im Anhang werden Zahlenbeispiele mit praktischen Folgerungen angegeben.

---

Summary The system of shot points I, II, ...A, B ... and geophones 1, 2 ...a, b according to Figure 1 assures — assuming horizontal reflector — by identical reflection point an advantageous determination of the average velocity and of the reflector by adjustment with equations (3)-(6). These equations can also be used if the reflector dips, as the error caused even by a dip of 11° of the reflector inv andN does not exceed 0.5 percent. The reflection point moves, however, simultaneously (see Figure 2) with the quantily r according to equations (17a).If the dip of the intersection line in the reflection plane is known, the more precise values ofv andN can be computed with the aid of equations (12) and (14).To determine the dip , the system of Figure 3 is most convenient, where the geophons 1, 2 ...a, b ... are attached to a shot inI/A.In the appendix some numerical solutions are given and practical consequences drawn.

     

Mathematical modeling of PCB bioaccumulation in Perna viridis

In the present work, we built a mathematical model of polychlorinated biphenyl (PCB) bioaccumulation in Perna viridis, namely, a one-compartment model with a time dependent incorporation rate R (μg g⁻¹ lipid per ppb water per day), with positive substrate cooperativity as the underlying physical mechanism. The temporal change of the PCB concentration Q (μg g⁻¹ lipid) in the soft tissues of the mussel depends on the competition of the input rate RW and the output rate kQ, where W is the concentration of PCB in water (ppb water) and k is the elimination rate (per day). From our experimental data, k=0.181±0.017 d⁻¹. The critical concentration in water W_c for positive substrate cooperativity was found to be 2.4 ppb. Below W_c, R is a constant. For a water concentration of 0.5 ppb Aroclor 1254, R=24.0±2.4 μg g⁻¹ lipid ppb⁻¹ d⁻¹. Above W_c, positive substrate cooperativity comes into effect and R becomes a function of time and dependent on the concentration Q in a form R=γQ/(Q+δ). This is the case for a water concentration of 5 ppb Aroclor 1254, where γ=15.1 μg g⁻¹ lipid ppb⁻¹ d¹ and δ≈200 μg g⁻¹ lipid. From this model, the uptake is exponentially increasing when the PCB concentration in the mussel is small compared to 200 μg g⁻¹ lipid, and hyperbolically increasing when the concentration is large compared to 200 μg g⁻¹ lipid, which are consistent with the experimental data. The model is useful for understanding the true processes taking place during the bioaccumulation and for risk assessment with higher confidence. Future experimental data which challenge the present model are anticipated and in fact desirable for improvement and perfection of the model.     

The role of vibrationally excited oxygen and nitrogen in the ionosphere during the undisturbed and geomagnetic storm period of 6–12 April 1990

We present a comparison of the observed behavior of the F-region ionosphere over Millstone Hill during the geomagnetically quiet and storm periods of 6–12 April 1990 with numerical model calculations from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere. The major enhancement to the IZMIRAN model developed in this study is the use of a new loss rate of O⁺(⁴S) ions as a result of new high-temperature flowing afterglow measurements of the rate coefficients K₁ and K₂ for the reactions of O⁺(⁴S) with N₂ and O₂. The deviations from the Boltzmann distribution for the first five vibrational levels of O₂(v) were calculated, and the present study suggests that these deviations are not significant. It was found that the difference between the non-Boltzmann and Boltzmann distribution assumptions of O₂(v) and the difference between ion and neutral temperature can lead to an increase of up to about 3% or a decrease of up to about 4% of the calculated NmF2 as a result of a respective increase or a decrease in K₂. The IZMIRAN model reproduces major features of the data. We found that the inclusion of vibrationally excited N₂(v > 0) and O₂(v > 0) in the calculations improves the agreement between the calculated NmF2 and the data on 6, 9, and 10 April. However, both the daytime and nighttime densities are reproduced by the IZMIRAN model without the vibrationally excited nitrogen and oxygen on 8 and 11 April better than the IZMIRAN model with N₂(v > 0) and O₂(v > 0). This could be due to possible uncertainties in model neutral temperature and densities, EUV fluxes, rate coefficients, and the flow of ionization between the ionosphere and plasmasphere, and possible horizontal divergence of the flux of ionization above the station. Our calculations show that the increase in the O⁺ + N₂ rate factor due to N₂(v > 0) produces a 5–36% decrease in the calculated daytime peak density. The increase in the O⁺ + O₂ loss rate due to vibrationally excited O₂ produces 8–46% reductions in NmF2. The effects of vibrationally excited O₂ and N₂ on N_e and T_e are most pronounced during the daytime.     

The relationship between vertical eddy diffusion and buoyancy gradient in the deep sea

Vertical eddy diffusivities (K_v's) have been estimated at fourteen widely separated locations from fourteen²²²Rn profiles and two²²⁸Ra profiles measured near the ocean floor as part of the Atlantic and Pacific GEOSECS programs. They show an inverse proportionality to the local buoyancy gradient [(g/?)(??_pot/?z)] calculated from hydrographic measurements. The negative of the constant of proportionality is the buoyancy flux [?K_v(g/?)(??_pot/?z)] which has a mean of ?4 × 10^?6 cm²/sec³. Our results suggest that the buoyancy flux varies very little near the ocean floor. K_v's for the interior of the deep Pacific calculated from the relationship K_v = (4 × 10^?6cm²/sec³)/[(g/?)(??_pot/?z)] agree well with published estimates. K_v's calculated for the pycnocline are one to two orders of magnitude smaller than upper limits estimated from tritium and⁷Be distributions.Heat fluxes calculated with the model K_v's obtained from the²²²Rn profiles average 31 μcal cm^?2 sec^?1 in the Atlantic Ocean and 8 μcal cm^?2 sec^?1 in the Pacific Ocean.     

Permeability Evolution During Non-linear Viscous Creep of Calcite Rocks

Permeability, storage capacity and volumetric strain were measured in situ during deformation of hot-pressed calcite aggregates containing 10, 20, and 30 wt% quartz. Both isostatic and conventional triaxial loading conditions were used. The tests were performed at confining pressure of 300 MPa, pore pressures between 50 to 290 MPa, temperatures from 673 to 873 K and strain rates of 3 × 10⁻⁵ s⁻¹. Argon gas was used as the pore fluid. The initial porosities of the starting samples varied from 5% to 9%, with higher porosity correlated to higher quartz content. Microstructural observations after the experiment indicate two kinds of pores are present: 1) Angular, crack-like pores along boundaries between quartz grains or between quartz and calcite grains and 2) equant and tubular voids within the calcite matrix. Under isostatic loading conditions, the compaction rate covaries with porosity and increases with increasing effective pressure. Most of the permeability reduction induced during compaction is irreversible and probably owes to plastic processes. As has been found in previous studies on hot-pressed calcite aggregates, permeability, k, is nonlinearly related to porosity, ϕ. Over small changes in porosity, the two parameters are approximately related as k ∝ ϕⁿ. The exponent n strongly increases as porosity decreases to a finite value (from about 4 to 6% depending on quartz content), suggesting a porosity percolation threshold. When subjected to triaxial deformation, the calcite-quartz aggregates exhibit shear-enhanced compaction, but permeability does not decrease as rapidly as it does under isostatic conditions. During triaxial compaction the exponent n only varies between 2 and 3. Non-isostatic deformation seems to reduce the percolation threshold, and, in fact, enhances the permeability relative to that at the same porosity during isostatic compaction. Our data provide constraints on the governing parameters of the compaction theory which describes fluid flow through a viscous matrix, and may have important implications for expulsion of sedimentary fluids, for fluid flow during deformation and metamorphism, and melt extraction from partially molten rocks.