N2 and M2 lunar tides: atmospheric resonance revisited

A numerical model has been used to calculate the atmospheric response to forcing at periods in the region of 12-13.5 h. The results show that the response is enhanced in the neighbourhood of 13 h. These results have been compared with lunar tidal analyses of mesospheric wind data and geomagnetic variations at a number of stations. It is found that the N₂ lunar tidal component (period 12.66 h) is significantly enhanced relative to the main lunar tidal component M₂ (period 12.42 h) in both types of data, compared with what would be expected from the gravitational tidal potential. This supports the predictions of the numerical model. An appreciable phase shift is also found in the experimental data between the N₂ and M₂ tides, agreeing in sense with what would be expected for a resonance at a period around 13 h.     

Numerical simulation of (T 2, T 1) 2D NMR and fluid responses

One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time (T ₁) and transverse relaxation time (T ₂) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times (T _W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings (T _E s) and T _W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.     

Melting relationships in CaO-CO2 and MgO-CO2 to 36 kilobars with comments on CO2 in the mantle

The melting curves of CaCO₃ and MgCO₃ have been extended to pressures of 36 kb by experiments in piston-cylinder apparatus. At 30 kb, the melting temperatures of calcite and magnesite are 1610°C and 1585°C, respectively. New data for the magnesite dissociation reaction permit the location of an invariant point for the assemblage magnesite + periclase + liquid + vapor near 26 kb-1550°C. New data are also presented for the calcite-aragonite transition at 800°C, 950°C and 1100°C. At pressures above 36–50 kb, calcite and magnesite melt at temperatures lower than the solidus of dry mantle peridotite. Natural and experimental evidence suggests that carbon dioxide in the Earth's mantle could be present in a variety of forms: (a) a free vapor phase, (b) vapor dissolved in silicate magma, (c) crystalline carbonate, (d) carbonatite liquid, (e) carbon-bearing silicate analogs, or (f) carbonato-silicates (such as scapolite, spurrite, tilleyite, and related compounds).     

Thermomagnetic analysis of meteorites, 2. C2 chondrites

Samples of all eighteen of the known C2 chondrites have been analyzed thermomagnetically. For eleven of these, initial Fe₃O₄ content is low (generally <1%) and theJ_s-T curves are irreversible. The heating curves show variable greater (up to 10 times) than it is initially. This behavior is attributed to the production of magnetite from a thermally unstable phase — apparently FeS. Four of the remaining seven C2 chondrites contain Fe₃O₄ as the only significant magnetic phase: initial magnetite contents range from 4 to 13%. The remaining three C2 chondrites contain iron or nickel-iron in addition to Fe₃O₄. These seven C2 chondrites show little evidence of the breakdown of a thermally unstable phase.     

2D Monte Carlo versus 2D Fuzzy Monte Carlo health risk assessment

Risk estimates can be calculated using crisp estimates of the exposure variables (i.e., contaminant concentration, contact rate, exposure frequency and duration, body weight, and averaging time). However, aggregate and cumulative exposure studies require a better understanding of exposure variables and uncertainty and variability associated with them. Probabilistic risk assessment (PRA) studies use probability distributions for one or more variables of the risk equation in order to quantitatively characterize variability and uncertainty. Two-dimensional Monte Carlo Analysis (2D MCA) is one of the advanced modeling approaches that may be used to conduct PRA studies. In this analysis the variables of the risk equation along with the parameters of these variables (for example mean and standard deviation for a normal distribution) are described in terms of probability density functions (PDFs). A variable described in this way is called a second order random variable. Significant data or considerable insight to uncertainty associated with these variables is necessary to develop the appropriate PDFs for these random parameters. Typically, available data and accuracy and reliability of such data are not sufficient for conducting a reliable 2D MCA. Thus, other theories and computational methods that propagate uncertainty and variability in exposure and health risk assessment are needed. One such theory is possibility analysis based on fuzzy set theory, which allows the utilization of incomplete information (incomplete information includes vague and imprecise information that is not sufficient to generate probability distributions for the parameters of the random variables of the risk equation) together with expert judgment. In this paper, as an alternative to 2D MCA, we are proposing a 2D Fuzzy Monte Carlo Analysis (2D FMCA) to overcome this difficulty. In this approach, instead of describing the parameters of PDFs used in defining the variables of the risk equation as random variables, we describe them as fuzzy numbers. This approach introduces new concepts and risk characterization methods. In this paper we provide a comparison of these two approaches relative to their computational requirements, data requirements and availability. For a hypothetical case, we also provide a comperative interpretation of the results generated.     

Kinetics of Methyldiethanolamine Mineralization by Using UV/H2O2 Process

The UV/H₂O₂ is one of the popular techniques in the advanced oxidation processes (AOPs) and has been applied in the wastewater treatment during recent two decades. UV exposure on the H₂O₂ generate highly reactive hydroxyl radicals (OH^?), which are used to degrade organic contaminants through oxidation processes in wastewater. This present study involves the estimation of hydroxyl radical rate constants of methyldiethanolamine (MDEA) mineralization at different temperatures by using UV/H₂O₂ in aqueous solution. Laboratory experiments have been conducted and the profile of MDEA mineralization has been established. The hydroxyl radical rate constants and the activation energy of mineralization process have been calculated. The estimated hydroxyl rate constants and the activation energy are in good agreement with those reported in the literature.     

Interference filter spectral imaging of twilight O+(2P-2D) emission

A spectral imager specifically designed to measure the O⁺(²P-²D) emission in the thermosphere during twilight has been constructed and tested in Toronto (43.8°N, 79.3°W), and found to show promise for long-term and campaign-mode operations. A modification of the mesopause oxygen rotational temperature imager (MORTI), it consists basically of a narrow-band interference filter (0.14 nm bandwidth) to separate wavelengths as a function of off-axis angle, a lens to focus the spectrum into a series of concentric rings, and a focal plane array (CCD) to record the spectral images in digital form. The instrument was built with two fields of view, one for the zenith and one for 20° above the horizon, movable to track the azimuth of the Sun, in order to provide appropriate data for inversion. Data gathered during June 1991 provided measurements of the column-integrated emission rate with a precision of about 3%. An atomic oxygen profile was deduced that showed good agreement with that predicted by the MSIS-90 model atmosphere. Geomagnetically induced variations of the O⁺ lines, calcium spectra resulting from meteor showers, and OH nightglow were also observed.     

Removal of Cu2+ and Pb2+ Ions Using CMC Based Thermoresponsive Nanocomposite Hydrogel

In this study, carboxymethylcellulose (CMC) based thermoresponsive nanocomposite hydrogel was synthesized for the removal of Cu²⁺ and Pb²⁺ ions from aqueous solutions. To prepare nanocomposite hydrogel, graft copolymerization of N‐isopropyl acrylamide (NIPAm) and acrylic acid (AA) onto CMC was carried out in Na‐montmorillonite (MMT)/water suspension media and ammonium persulfate (APS) used as initiator. The chemical structures of hydrogels were characterized by Fourier transform infrared (FT‐IR) and X‐ray diffraction spectroscopy (XRD). Lower critical solution temperature (LCST), pH responsivity, swelling, and deswelling properties of the hydrogels were also examined. In addition competitive and non‐competitive removal of Cu²⁺ and Pb²⁺ studies were carried out. According to heavy metal sorption studies results, removal capacities of nanocomposite hydrogel for both metal ions were found to be higher than those of pure hydrogel. The analyzed adsorption data showed that the adsorption process of Cu²⁺ and Pb²⁺ could be explained by pseudo‐second order kinetic model. Moreover, according to competitive sorption studies, it is found to be that both hydrogels are more selective to Cu²⁺ ion rather than Pb²⁺.     

The determination of deep temperatures by means of the CO-CO2-H2-H2O geothermometer: an example using fumaroles in the Campi Flegrei,Italy. A comment

Three points raised in the paper by Tedesco and Sabroux (1987) are dealt with. (1) The inconsistency between the water partial pressure calculated by Tedesco and Sabroux (1987) and saturation pressure is due to the improper use of the water-gas-shift reaction as a geothermometer. In fact Tedesco and Sabroux (1987) do not take into account the distribution of gas species between the coexisting vapour and liquid phases. (2) The depth of the “steam reservoir” is evaluated by Tedesco and Sabroux (1987) in too simplistic a way. This matter should be treated with greater care owing to the high social impact of any consideration on the Phlegraean Fields system. (3) The reliability of carbon monoxide determination at the concentration level encountered at Solfatara depends on the collection method rather than on the gas-chromatographic technique.     

(T2,D)二维核磁共振测井混合反演方法与参数影响分析

一维核磁共振(1D NMR)测井技术在流体识别中具有一定的局限性.二维核磁共振(2D NMR)测井能同时测量到多孔介质中横向弛豫时间(T₂)和扩散系数(D)等信息,利用这两个参数区分流体性质,较一维核磁共振测井技术具有明显的优越性.针对梯度场下的2D NMR测井弛豫机理和数学模型,提出了适用于求解大型稀疏矩阵方程的反演方法-基于非负最小二乘法(LSQR)和截断奇异值分解(TSVD)法的混合算法.为验证方法的有效性,先根据多回波观测模式合成回波串数据,然后再用混合反演算法进行反演,反演得到横向弛豫时间(T₂)和扩散系数(D),并构建T₂-D二维谱图.结果对比表明,该混合反演算法得到的T₂-D二维谱与流体模型一致性好,计算精度均比单一反演方法有较大改善,表明该混合反演方法可用于油气储集层2D NMR测井的反演和流体识别.此外,分别对油水同层和气水同层模型进行了正演模拟和反演实验, 系统考察了不同磁场梯度、不同回波间隔组合对反演效果的影响,为2D NMR参数设计提供依据.     

High-pressure decomposition of Ca2GeO4 and Ca2MnO4 with the K2NiF4-type structures

By using the diamond-anvil pressure cell coupled with laser heating, Ca₂GeO₄ in the K₂NiF₄-type structure has been found to decompose into the mixture Ca₃Ge₂O₇ plus CaO at pressures greater than 200 kbar and at about 1000°C, and the same type of structure for Ca₂MnO₄ has been found to decompose into the mixture CaMnO₃ (perovskite) plus CaO at pressures greater than 100 kbar and at about 1400°C. The decomposition product of Ca₃Ge₂O₇ is a new compound which is isostructural with Sr₃Ti₂O₇ and has the lattice parameters of a = 3.72 ± 0.01 and $\text{c = 19.32 ± 0.05}\text{A}\text{?}$ at room temperature and 1 bar pressure. The results of the study of Ca₂GeO₄ and Ca₂MnO₄ (both with the K₂NiF₄-type structure) strongly support the view that compounds possessing the K₂NiF₄-type structure are unstable relative to corresponding mixtures possessing the perovskite and rocksalt structures. It is concluded that, in the earth's mantle, the K₂NiF₄-type Ca₂SiO₄ would ultimately decompose into the mixture CaSiO₃ (perovskite) + CaO or would otherwise transform to other as-yet-unknown phase(s), and that the mixture of MgSiO₃ (perovskite) + MgO (the post-spinel phase of Mg₂SiO₄) would not adopt the K₂NiF₄-type structure.     

TTI介质的相速度研究

推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．     

2D Optimized Electrode Arrays for Borehole Resistivity Tomography and CO2 Sequestration Modelling

Modern optimization approaches for electrode configurations can significantly improve the resolution of 2.5D resistivity imaging surveys. This study presents a brief review of the 2.5D optimization approach, particularly for borehole–borehole surveys with applications for mapping virtual CO₂ plumes sequestrated in deep saline reservoir formations. The applied algorithm searches for arrays that maximize the spatial resolution of the survey among the comprehensive dataset of best possible spatial resolution (i.e. least temporal resolution). A main goal of this study is to increase the temporal resolution of ERT borehole–borehole surveys by selecting optimized electrode configurations in order to minimise the required data acquisition time while sustaining a high spatial resolution. The optimized dataset starts with a base set and is iteratively increased based on the model resolution matrix (R ) until the required number of data points is achieved. Among four different optimization methods, the compare R (CR) method of the best resolution is applied to directly calculate R for each new array added to the optimized dataset. Small optimized datasets generated by this technique are only <5% of their comprehensive sets but of an average resolution ratio (R _r) of >0.95 (i.e. almost the same resolution). With increasing the size of the optimized dataset (during its generation), the algorithm progressively enhances R _r values in the central interwell region (of low sensitivities and low resolution) far higher than in the near borehole region (of high sensitivities). Also the inverted tomogram reliability increases by increasing the optimized data size. Briefly, the optimized arrays improve the resolution in the interwell region which is commonly low in borehole–borehole ERT studies. The inverted output model is evaluated quantitatively using the model difference relative to the input model. The results reflect the common smearing effects and artefacts of varying degrees that overpredict volumes, underpredict magnitudes and blur boundaries of the target anomalies. This input model is a synthetic resistivity model that was used to generate synthetic (forward solution) data used during the inversion. Applications on synthetic CO₂ models show that the mapping resolution for optimized datasets is better than that for other highly resolving arrays of the same number of data points. Problems of smeared boundaries and thin layers are less visible in the optimized array than in the other highly resolving arrays.     

Variations of near-surface atmospheric CO_2 and H_2O concentrations during summer on Muztagata

1 Introduction As one of the most important greenhouse gases, atmospheric carbon dioxide (CO2) has increased in concentration rapidly since preindustrial times[1―4] and significantly contributes to the climate change[5] caused by anthropogenic emissions. Documentation of the global carbon cycle has been critical for under-standing the causative relationships between green-house gases and climate change. For example, the CO2 level in the atmosphere has been monitored at many sites worldwide …     

Chemical composition of boulder-2 rocks and soils,Apollo 17, station 2

Forty-two elements have been measured via INAA and RNAA in six samples of five rocks from a 2-m tan-gray boulder-2 breccia (South Massif), in four soils from the South Massif and a valley soil. The chemical composition of the four metaclastic rocks corresponds to “high alumina” (52% Pl) and medium-K KREEP-type rocks. Rock 72335,2 is a medium K anorthositic gabbro (74% Pl). Both the North and South Massifs appear to be medium-K KREEP in composition and thus may represent a single stratigraphic unit of the Serenitatis basin event. Four soils of the South Massif are identical in composition to medium-K KREEP; they are more feldspathic and lower in LIL trace elements relative to the boulder-2 rocks. The valley soil 75081 is like the 10084 soil; both soils are high in TiO₂ and both are deficient in KREEP. Th and U give a sharp distinction between the valley and highland soils. The South Massif rocks and soils contain siderophiles at the 2–4% Cl level and show an ancient meteoritic pattern. Five samples of the four rocks have Ir/Au ratios of ± 0.02 which we assign to the Serenitatis basis planetesimal. The valley soil at Camelot Crater has low siderophiles (1% Cl). Our systematic study of four shadowed and exposed soils does not support the labile hypothesis for Cd, In, Tl and Zn. We observe no volatile (atmophile) movement from the South Massif highland soils to the valley soil 75081. The volatiles Cs and Tl appear to have been fractionated in the boulder-2 rocks during cratering, brecciation and metamorphic processes. The uniform ratio of FeO/MnO = 80–85, observed for all previous mare and highland sites, also holds for the Taurus-Littrow site.     

反应HO2+O→OH+O2对OH夜气辉辐射的贡献

本文探讨了OH Meinel夜气辉辐射的光化学模式，给出OH(ν≤9)分子数密度分布的计算通式，重点研究了化学反应HO\-2+O→OH(ν≤6)+O2对ν≤6各振动能级上OH分子数密度分布以及(ν′—ν″)(ν′≤6)振动带气辉辐射的影响.结果表明，该反应对数密度的贡献随振动能级的减小而增大，对(ν′—ν″)振动带辐射的贡献随着较高振动能级ν′的减小而增大，以春分时为例，它可使第1振动能级上的OH分子最大数密度和（1—0）带的最大辐射率增加约33%，第6振动能级上的最大数密度和（6—ν″）带的最大辐射率增加约7%，（1—0）带的辐射强度增加约30%，（6—ν″）带的强度增加约11%.该反应使各振动能级上分子数密度的高度分布剖面以及各振动带体辐射率的高度分布剖面变宽，最大数密度、最大发射率所处的高度下降1 km左右.此外，该反应的影响程度随着原子氧密度的降低而增大，随着温度的升高而增大，并且在夏至时最大，在冬至时最小.     

Inversion of 2D spectral induced polarization imaging data

Laboratory measurements of various materials suggest that more information can be obtained by measuring the in‐phase and out‐of‐phase potentials at a number of frequencies. One common model used to describe the variation of the electrical properties with frequency is the Cole‐Cole model. Apart from the DC resistivity (ρ) and chargeability (m) parameters used in conventional induced‐polarization (IP) surveys, the Cole‐Cole model has two additional parameters, i.e. the time (τ) and relaxation (c) constants. Much research has been conducted on the use of the additional Cole‐Cole parameters to distinguish between different IP sources. Here, we propose a modified inversion method to recover the Cole‐Cole parameters from a 2D spectral IP (SIP) survey. In this method, an approximate inversion method is initially used to construct a non‐homogeneous starting model for the resistivity and chargeability values. The 2D model consists of a number of rectangular cells with constant resistivity (ρ), chargeability (m), time (τ) and relaxation (c) constant values in each cell. A regularized least‐squares optimization method is then used to recover the time and relaxation constant parameters as well as to refine the chargeability values in the 2D model. We present results from tests carried out with the proposed method for a synthetic data set as well as from a laboratory tank experiment.     

Two-component spectral analysis of geomagneticpi2 pulsations

Summary The method of two-component spectral analysis is described in two versions: as the so-called static analysis, in which the pulsation sample being studied is treated as a whole, and as the dynamic analysis, in which the given sample is processed piecewise in order to determine the time variations of various parameters. Both versions were applied to geomagnetic pi2 pulsations. The pi2 oscillations are resolved into several frequency components. The properties of these components are studied in dependence on the frequency and local time. The time variation of the amplitude (attenuation) and the polarization characteristics (precession) of pi2's are also studied.     

CO2-filled vesicles in mid-ocean basalt

Volatile-filled vesicles are present in minor amounts in all samples of mid-ocean basalt yet collected (and presumably erupted) down to depths of 4.8 km. When such vesicles are pierced in liquid under standard conditions, the volume expansion of the gas is 0.2 ± 0.05 times the eruption pressure in bars or 20 ± 5 times the eruption depth in km. Such expansion could be used as a measure of eruption depth.A variety of techniques: (1) vacuum crushing and gas chromatographic, freezing separation, and mass spectrographic analyses; (2) measurements of phase changes on a freezing microscope stage; (3) microscopic chemical and solubility observations; and (4) volume change measurements, all indicate that CO₂ comprises more than 95% by volume of the vesicle gas in several submarine basalt samples from the Atlantic and Pacific. The CO₂ held in vesicles is present in quantities about equal to or greater than that presumed to be dissolved in the glass (melt) and amounts to 400–900 ppm of the rock. The rigid temperature of the glass is 800–1000°C and increases for shallower samples. A sulfur gas was originally present in subordinate amounts in the vesicles, but has largely reacted with iron in the vesicle walls to produce sulfide spherules.     

True-amplitude migration of 2D synthetic data1

True-amplitude (TA) migration, which is a Kirchhoff-type modified weighted diffraction stack, recovers (possibly) complex angle-dependent reflection coefficients which are important for amplitude-versus-offset (AVO) inversion. The method can be implemented using existing prestack or post-stack Kirchhoff migration and fast Green's function computation programs. Here, it is applied to synthetic single-shot and constant-offset seismic data that include post-critical reflections (complex reflection coefficients) and caustics. Comparisons of the amplitudes of the TA migration image with theoretical reflection coefficients show that the (possibly complex) angle-dependent reflection coefficients are correctly estimated.