Control‐volume mixed finite element methods

A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions.     

The foreshock sequence of haicheng earthquake and earthquake swarm—the use of foreshock sequences in earthquake prediction

We have compared the Haicheng foreshock sequence with several earthquake swarms which occurred in its neighborhood. The spatial distribution of the earthquakes is relatively concentrated. For the most part, the events occurred within a few kilometers of each other. The focal mechanisms are comparatively stable. However, there are several swarms in which the variations of focal mechanisms are quite obvious after the occurrence of the largest event of the sequence, which would allow it to be recognized as a swarm. However, there are also swarms whose focal mechanisms are no less stable throughout the sequence compared to the Haicheng foreshock sequence. This feature could thus not be used to identify a foreshock sequence. The temporal distributions of foreshocks and swarms are quite similar in some cases. This is again not a definite criterion for identifying foreshocks, but is worthy of further study. Thus, no definite criterion for identifying foreshock sequences has been found. However, some earthquake swarms may be recognized in their later stage.Finally, we introduced a magnitude sequence with gaps which can be used to see whether a large event is still forthcoming. This method (in conjunction with other methods) could be used in areas prone to large earthquakes, immediately before a large event, to improve the probability of predicting the occurrence of a large event. We also report that the temporal distribution of all the sequences showed a 12-hour recurrence pattern that corresponded with the earth tides, indicating that tidal forces might be influencing foreshocks and earthquake swarm occurrence.     

Hydroxyl stretching in phyllosilicates at high pressures and temperatures: an infrared spectroscopic study

The hydroxyl stretching frequencies of four phyllosilicates have been measured at high pressures and temperatures using an externally heated diamond-anvil cell and synchrotron infrared spectroscopy. Spectra were measured up to 26, 31, 21 and 8 GPa at room temperature for samples of talc, pyrophyllite, muscovite and 10-Å phase, respectively. Spectra were also measured in the range 273–500 K at ambient pressure for all samples and at 8–9 GPa for talc and pyrophyllite. The frequency of the Mg₃OH band in talc increases with pressure due to the absence of hydrogen bonding. The different orientation of the hydroxyl group in pyrophyllite and muscovite leads to hydrogen bonding and a decrease in the frequency of the Al₂OH band with pressure. 10-Å phase is approximately equivalent to talc with the addition of interlayer H₂O. In a spectrum of a sample synthesised for 143 h, two hydroxyl stretching bands are clearly resolved on compression. One is the same as the Mg₃OH band in talc, indicating the presence of intra-layer hydroxyl in a talc-like environment with no hydrogen bonding. The other, which separates from the talc-like band at 1 GPa, is associated with intra-layer hydroxyl that is hydrogen bonded to interlayer H₂O. There are equivalent bands in high-pressure spectra of a sample of deuterated 10-Å phase, synthesised for 400 h. This sample shows a greater extent of hydrogen bonding at ambient pressure than the 143 h sample. For all of the phases studied, increasing temperature leads to a decrease in frequency for every hydroxyl stretching vibration, both at low and high pressures. The shifts in frequency with temperature are an order of magnitude greater than the shifts with pressure when normalised to previously measured structural parameters.     

Laboratory Measurement and Correction of Thermal Properties for Application to the Rock Mass

The requirement for quantitative thermal property data in assessing the performance of an engineered repository is identified. Experimental methods for the measurement of thermal conductivity, thermal expansivity and heat capacity are outlined and typical test results are presented for two rock-types of interest in relation to the Sellafield site. The controls on the thermal properties are discussed and a series of corrections from the rock material to the rock mass properties are presented. Finally, a correction scheme that accounts for the effect of scale is followed.     

The geochemical coherence of Pu and Nd and the 244Pu238U ratio of the early solar system

Two very different sets of ${}^{244}\text{Pu}{}^{238}\text{U}$ ratios have been reported for early solar system materials. One group of samples yields high (0.015–0.016) ratios (Podosek, 1970a, 1972; Drozd et al., 1977) and calculations based on another group of analyses yield low ratios (~0.004) (Marti et al., 1977). Recently measured partition coefficients for Pu and Sm are used to evaluate the data of Marti et al. and $\text{sol}{}^{244}\text{Pu}{}^{238}\text{U}$ ratios from other sources are also considered. A low $\text{sol}{}^{244}\text{Pu}{}^{238}\text{U}$ ratio (~0.005) is favored, and some implications of this low ratio to galactic nucleosynthesis and meteorite age dating are briefly discussed.     

Lake ecosystem dynamics and links to climate change inferred from a stable isotope and organic palaeorecord from a mountain lake in southwestern China (ca. 22.6–10.5 cal ka BP)

A detailed understanding of long-term climatic and environmental change in southwestern China is hampered by a lack of long-term regional palaeorecords. Organic analysis (%TOC, %TN, C/N ratios and δ¹³C values) of a sediment sequence from Lake Shudu, Yunnan Province (ca. 22.6–10.5 cal ka BP) indicates generally low aquatic palaeoproductivity rates over millennial timescales in response to cold, dry climatic conditions. However, the record is punctuated by two marked phases of increased aquatic productivity from ca. 17.7 to 17.1 cal ka BP and from ca. 11.9 to 10.5 cal ka BP. We hypothesise that these shifts reflect a marked, stepwise lacustrine response to Asian summer monsoon strengthening during the last deglaciation.     

Tests for size and shape dependency in deep-sea mixing

Bioturbation is generally not a strongly size-dependent or shape-dependent process in deep-sea sediments. Different sizes of ash (11–250 μm) deposited during a single eruption are usually mixed about equally, although in a few cores we found minor tendencies toward greater mixing of finer ash. Platey volcanic ash and spherical microtektites initially deposited within a few centimeters of each other in the sediment column are mixed with roughly comparable intensity. Radiocarbon dating of three sizes of CaCO₃ near the top of one box core produced ages within the analytical error; this also suggests no size dependency in the mixing process. Another core showed age differences opposite in sense to those that would be predicted if the finer fractions were more intensely mixed than the coarser fractions.The sizes examined (11–250 μm) encompass most of the sedimentary components used in paleoclimatic analyses of deep-sea cores. From these findings, we infer that significant artificial lead/lag offsets in paleoclimatic signals will not be created by size-selective or shape-selective mixing. Mixing models calibrated on one sediment size fraction should be generally applicable to other sizes and kinds of material.