The role of reaction affinity and secondary minerals in regulating chemical weathering rates at the Santa Cruz Soil Chronosequence, California

In order to explore the reasons for the apparent discrepancy between laboratory and field weathering rates and to determine the extent to which weathering rates are controlled by the approach to thermodynamic equilibrium, secondary mineral precipitation, and flow rates, a multicomponent reactive transport model (CrunchFlow) was used to interpret soil profile development and mineral precipitation and dissolution rates at the 226 ka Marine Terrace Chronosequence near Santa Cruz, CA. Aqueous compositions, fluid chemistry, transport, and mineral abundances are well characterized [White A. F., Schulz M. S., Vivit D. V., Blum A., Stonestrom D. A. and Anderson S. P. (2008) Chemical weathering of a Marine Terrace Chronosequence, Santa Cruz, California. I: interpreting the long-term controls on chemical weathering based on spatial and temporal element and mineral distributions. Geochim. Cosmochim. Acta72 (1), 36-68] and were used to constrain the reaction rates for the weathering and precipitating minerals in the reactive transport modeling. When primary mineral weathering rates are calculated with either of two experimentally determined rate constants, the nonlinear, parallel rate law formulation of Hellmann and Tisserand [Hellmann R. and Tisserand D. (2006) Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar. Geochim. Cosmochim. Acta70 (2), 364-383] or the aluminum inhibition model proposed by Oelkers et al. [Oelkers E. H., Schott J. and Devidal J. L. (1994) The effect of aluminum, pH, and chemical affinity on the rates of aluminosilicate dissolution reactions. Geochim. Cosmochim. Acta58 (9), 2011-2024], modeling results are consistent with field-scale observations when independently constrained clay precipitation rates are accounted for. Experimental and field rates, therefore, can be reconciled at the Santa Cruz site.Additionally, observed maximum clay abundances in the argillic horizons occur at the depth and time where the reaction fronts of the primary minerals overlap. The modeling indicates that the argillic horizon at Santa Cruz can be explained almost entirely by weathering of primary minerals and in situ clay precipitation accompanied by undersaturation of kaolinite at the top of the profile. The rate constant for kaolinite precipitation was also determined based on model simulations of mineral abundances and dissolved Al, SiO₂(aq) and pH in pore waters. Changes in the rate of kaolinite precipitation or the flow rate do not affect the gradient of the primary mineral weathering profiles, but instead control the rate of propagation of the primary mineral weathering fronts and thus total mass removed from the weathering profile. Our analysis suggests that secondary clay precipitation is as important as aqueous transport in governing the amount of dissolution that occurs within a profile because clay minerals exert a strong control over the reaction affinity of the dissolving primary minerals. The modeling also indicates that the weathering advance rate and the total mass of mineral dissolved is controlled by the thermodynamic saturation of the primary dissolving phases plagioclase and K-feldspar, as is evident from the difference in propagation rates of the reaction fronts for the two minerals despite their very similar kinetic rate laws.     

Evidence for high temperature and 18O‐enriched fluids in the Arab‐D of the Ghawar Field,Saudi Arabia

Using the clumped isotope method, the temperature of dolomite and calcite formation and the oxygen isotopic composition (δ¹⁸O_w) of the diagenetic fluids have been determined in a core taken from the Arab‐D of the Ghawar field, the largest oil reservoir in the world. These analyses show that while the dolomites and limestones throughout the major zones of the reservoir recrystallized at temperatures between ca 80°C and 100°C, the carbonates near the top of the reservoir formed at significantly lower temperatures (20 to 30°C). Although the δ¹⁸O values of the diagenetic fluids show large variations ranging from ca <0‰ to ca +8‰, the variations exhibit consistent downhole changes, with the highest values being associated with the portion of the reservoir with the highest permeability and porosity. Within the limestones, dolomites and dolomites associated with the zone of high permeability, there are statistically significant different trends between the δ¹⁸O_w values and recrystallization temperature. These relationships have different intercepts suggesting that fluids with varying δ¹⁸O_w values were involved in the formation of dolomite and limestone compared to the formation of dolomite associated with the zone of high permeability. These new data obtained using the clumped isotope technique show how dolomitization and recrystallization by deep‐seated brines with elevated δ¹⁸O_w values influence the δ¹⁸O values of carbonates, possibly leading to erroneous interpretations unless temperatures can be adequately constrained.     

基于深度学习卷积神经网络的地震波形自动分类与识别

发展高效、高精度、普适性强的自动波形拾取算法在地震大数据时代背景下显得越来越重要.波形自动拾取算法的主要挑战来自如何适应不同区域的不同类型地震事件的分类与筛选.本文针对地震事件-噪音分类这一问题, 使用13839个汶川地震余震事件建立数据集, 应用深度学习卷积神经网络(CNN)方法进行训练, 并用8900个新的汶川余震事件作为检测数据集, 其训练和检测准确率均达到95%以上.在对连续波形的检测中, CNN方法在精度和召回率上优于STA/LTA和Fbpicker传统方法, 并能找出大量人工挑选极易遗漏的微震事件.最后, 我们应用训练好的最优模型对选自全国台网的441个台站8天的连续波形数据进行了识别、到时挑取及与参考地震目录关联, CNN检出7016段波形, 用自动挑选算法拾取到1380对P, S到时, 并与540个地震目录事件成功关联, 对1级以上事件总体识别准确率为54%, 二级以上为80%, 证明了CNN模型具有泛化能力, 初步展示了CNN在发展兼具效率、精度、普适性算法, 实时地震监测等应用上具有巨大潜力.

     

Dynamic response characterization and simplified analysis methods for viscoelastic dampers considering heat transfer

Viscoelastic (VE) dampers are sensitive to temperature, excitation frequency, and strain level. As they dissipate the kinetic energy from earthquake or wind-induced structural vibrations, their temperature increases from the heat generated, consequently softening their VE materials and lowering their dynamic mechanical properties. Temperature increase can be significant for long-duration loading, but can be limited by heat conduction and convection which depend on damper configuration. The writers analytically explored such effect on the six different dampers by using their previously proposed three-dimensional finite-element analysis method. Results provided better understanding of how heat is generated within the VE material, conducted and stored in different damper parts, and dispersed to the surrounding air. These results also led to characterization of both local (e.g., temperatures, properties, and strain energy density) and global (e.g., hysteresis loops, and stiffness) behavior of VE dampers, and provided a framework for a new simplified one-dimensional (1D) modeling approach for time-history analysis. This new proposed 1D method greatly improves the computation time of the previously proposed long-duration method coupling fractional time-derivatives VE constitutive rule with 1D heat transfer analysis. Unlike the previous method, it idealizes uniform shear strain and VE material property distributions for computational efficiency, but still simulating non-uniform temperature distribution along the thickness direction of the VE material. Despite the approximations, it accurately predicts VE damper global responses.