Historical simulation and twenty-first century prediction of oceanic CO2 sink and pH change

A global ocean carbon cycle model based on the ocean general circulation model POP and the improved biogeochemical model OCMIP-2 is employed to simulate carbon cycle processes under the historically observed atmospheric CO 2 concentration and different future scenarios (called Rep- resentative Concentration Pathways, or RCPs). The RCPs in this paper follow the design of Inter- governmental Panel on Climate Change (IPCC) for the Fifth Assessment Report (AR5). The model results show that the ocean absorbs CO 2 from atmosphere and the absorbability will continue in the 21st century under the four RCPs. The net air-sea CO 2 flux increased during the historical time and reached 1.87 Pg/a (calculated by carbon) in 2005; however, it would reach peak and then decrease in the 21st century. The ocean absorbs CO 2 mainly in the mid latitude, and releases CO 2 in the equator area. However, in the Antarctic Circumpolar Current (ACC) area the ocean would change from source to sink under the rising CO 2 concentration, including RCP4.5, RCP6.0, and RCP8.5. In 2100, the anthropogenic carbon would be transported to the 40 S in the Atlantic Ocean by the North Atlantic Deep Water (NADW), and also be transported to the north by the Antarctic Bottom Water (AABW) along the Antarctic continent in the Atlantic and Pacific oceans. The ocean pH value is also simulated by the model. The pH decreased by 0.1 after the industrial revolution, and would continue to decrease in the 21st century. For the highest concentration sce- nario of RCP8.5, the global averaged pH would decrease by 0.43 to reach 7.73 due to the absorption of CO 2 from atmosphere.     

Seismoelectric measurements in a porous quartz‐sand sample with anisotropic permeability

Seismoelectric coupling coefficients are difficult to predict theoretically because they depend on a large numbers of rock properties, including porosity, permeability, tortuosity, etc. The dependence of the coupling coefficient on rock properties such as permeability requires experimental data. In this study, we carry out a set of laboratory measurements to determine the dependence of seismoelectric coupling coefficient on permeability. We use both an artificial porous “sandstone” sample, with cracks, built using quartz‐sand and Berea sandstone samples. The artificial sample is a cube with 39% porosity. Its permeability levels are anisotropic: 14.7 D, 13.8 D, and 8.3 D in the x‐, y‐, and z‐directions, respectively. Seismoelectric measurements are performed in a water tank in the frequency range of 20 kHz–90 kHz. A piezoelectric P‐wave source is used to generate an acoustic wave that propagates through the sample from the three different (x, y, and z) directions. The amplitudes of the seismoelectric signal induced by the acoustic waves vary with the direction. The highest signal is in the direction of the highest permeability, and the lowest signal is in the direction of the lowest permeability. Since the porosity of the sample is constant, the results directly show the dependence of seismoelectric coefficients on permeability. Seismoelectric measurements with natural rocks are performed using Berea sandstone 500 and 100 samples. Because the Berea samples are nearly isotropic in permeability, the amplitudes of the seismoelectric signals induced in the different directions are the same within the measurement error. Because the permeability of Berea 500 is higher than that of Berea 100, the amplitude of the seismoelectric signals induced in Berea 500 is higher than those in Berea 100. To determine the relative contributions of porosity and permeability on seismoelectric conversion, we carried out an analysis, using Pride (1994) formulation and Kozeny–Carman relationship; the normalized amplitudes of seismoelectric coupling coefficients in three directions are calculated and compared with the experimental results. The results show that the seismoelectric conversion is related to permeability in the frequency range of measurements. This is an encouraging result since it opens the possibility of determining the permeability of a formation from seismoelectric measurements.     

会聚束电子衍射对葡萄石空间群的辅助测定

利用透射电子显微技术测定了湖北大冶铁山矿床中的微米级葡萄石空间群。选区电子衍射几何构图法分析表明,该矿物系斜方晶系,衍射消光符号为Ｐｎ－－晶胞参数为α＝０．４５８,ｂ＝０．５５５和ｃ＝１．８５３ｎｍ。（１００）和（３０１）方向的会聚束电子衍射花样研究表明,其空间群为Ｐｎｃｍ。     

Development of Climate and Earth System Models in China: Past Achievements and New CMIP6 Results

The Earth-Climate System Model (ECSM) is an important platform for multi-disciplinary and multi-sphere integration research, and its development is at the frontier of international geosciences, especially in the field of global change. The research and development (R&D) of ECSM in China began in the 1980s and have achieved great progress. In China, ECSMs are now mainly developed at the Chinese Academy of Sciences, ministries, and universities. Following a brief review of the development history of Chinese ECSMs, this paper summarized the technical characteristics of nine Chinese ECSMs participating in the Coupled Model Intercomparison Project Phase 6 and preliminarily assessed the basic performances of four Chinese models in simulating the global climate and the climate in East Asia. The projected changes of global precipitation and surface air temperature and the associated relationship with the equilibrium climate sensitivity under four shared socioeconomic path scenarios were also discussed. Finally, combined with the international situation, from the perspective of further improvement, eight directions were proposed for the future development of Chinese ECSMs.     

动电测井实验研究Ⅱ:伴随动电场和界面动电场

动电效应指孔隙介质中与固/液界面的双电层和孔隙流体渗流有关的弹性-电磁耦合现象,探索基于动电效应的勘探和测井新方法则是石油工业重点关注的研究方向之一.本文对流体饱和孔隙介质中的动电效应进行了实验测量研究,记录到了不同模型井中伴随声波的动电转换信号和界面动电转换信号,对比分析了这两种动电信号的产生条件及传播特性,验证了理论分析结果,并进一步研究了动电信号的分波成分及其与声波信号的关系,探讨了声源激发模式、电极接收方式及数据处理方法对动电信号的影响,为动电测井仪器设计奠定实验基础.     

聚晶金刚石复合片残余热应力的影响因素

残余热应力是影响聚晶金刚石复合片(PDC)性能好坏的最重要因素之一。考虑聚晶金刚石层(PCD)与硬质合金层厚度比以及PDC压制过程中烧结温度的波动对聚晶金刚石复合片残余热应力的影响,在ANSYS中建立PDC模型,运用热-结耦合法分析PDC的残余热应力。计算表明,随着PCD层与硬质合金层厚度比由0.067增加到0.333,PCD层表面中心的压应力由1.61 GPa降低到380 MPa,PCD层最大径向压应力由1.61 GPa降低1.03 GPa左右,而PCD层边缘靠近界面附近最大轴向拉应力逐渐增大;随着PDC压制过程中烧结温度由1000℃升高到1500℃,PCD层的最大径向压应力、最大轴向拉应力以及最大剪应力等均逐渐增大。认为,在研究PDC合成新工艺过程中,应在保证PDC使用寿命的前提下尽量降低PCD层与硬质合金层厚度比;必须尽量切断原材料以及人为操作对温度的影响。     

Effects of inorganic seed aerosols on the growth and chemical composition of secondary organic aerosol formed from OH-initiated oxidation of toluene

(NH₄)₂SO₄, CaCl₂, Na₂SiO₃ and NaNO₃ were selected as surrogates of inorganic seed aerosols of ambient atmosphere of Chinese urban areas, respectively, to study their effects on the formation of secondary organic aerosol (SOA) in the toluene/CH₃ONO/NO_x photooxidation system. The SMPS and aerosol laser time-of-flight mass spectrometer (ALTOFMS) was used to measure the aerodynamic size and chemical composition of individual SOA particles in real-time. Experimental results indicate that either the growth or products of SOA is affected by the presence of inorganic seed aerosol. Inorganic seed aerosols would promote growth rates of SOA formation at the start of the reaction and inhibits its formation rate with prolonging the reaction time. In the case of about 100 μg m^?3 seed aerosol load, the addition of Na₂SiO₃ induced a same growth rate of SOA formation as NaNO₃. The influence of four individual seed aerosols on the generation of SOA decreased in the order of CaCl₂ > (NH₄)₂SO₄ > NaNO₃, Na₂SiO₃. The presence of Na₂SiO₃ or NaNO₃ has no obvious effect on the growth rates of SOA formation, but it does increase the yield of organic acid and nitrogen-containing organic compounds, respectively. Besides the significantly effect on the growth rate of SOA formation, the presence of CaCl₂ or (NH₄)₂SO₄ can lead to the formation of high-molecular weight species which is found to be positively correlated with the hygroscopic behavior of seed aerosols. The CaCl₂ shows the strongest hygroscopic behavior among the four individual seed aerosols, and the most significant promotion effect on the formation of the high-molecular weight species. It is proposed that the SOA generation enhancement and high-molecular weight products are achieved by particle-phase heterogeneous reactions induced and catalyzed by the acidity of CaCl₂ and (NH₄)₂SO₄ seed aerosols.     

Prospective scenarios of the saltwater intrusion in an estuary under climate change context using Bayesian neural networks

In recent decades, saltwater intrusion over some low-lying coastal regions was deteriorated by rising sea-level and decreasing streamflow in the context of climate change. Though physically-based hydrodynamic models are the most detailed means to simulate salinity processes, they are commonly restricted by data insufficiency issues both in spatial resolution and temporal lasting. This motivates us to build a statistical model enable simulation and scenario analysis for coastal salinity change with limited observations. A Bayesian neural network (BNN) model is built hereby to simulate salinity. It offers more precise estimation compared with the conventional artificial neural network. Meanwhile, the model gives the uncertainty behaviors of the final salinity simulation which is not available for other methods. Future scenarios of salinity change are constructed and analyzed in different time periods on the basis of the validated BNN model. Results indicate that the water quality over lower Pearl River is degrading along with more significant uncertainties. Further analysis suggests that streamflow alteration has a more direct impact on salinity variations than the sea-level change does. The method allows a profound analysis of the potential influence on water quality degradation in coastal and low-lying regions in support of water management and adaptation toward global climate change.     

Arctic sea ice concentration and thickness data assimilation in the FIO-ESM climate forecast system

To improve the Arctic sea ice forecast skill of the First Institute of Oceanography-Earth System Model (FIO-ESM) climate forecast system, satellite-derived sea ice concentration and sea ice thickness from the Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS) are assimilated into this system, using the method of localized error subspace transform ensemble Kalman ?lter (LESTKF). Five-year (2014–2018) Arctic sea ice assimilation experiments and a 2-month near-real-time forecast in August 2018 were conducted to study the roles of ice data assimilation. Assimilation experiment results show that ice concentration assimilation can help to get better modeled ice concentration and ice extent. All the biases of ice concentration, ice cover, ice volume, and ice thickness can be reduced dramatically through ice concentration and thickness assimilation. The near-real-time forecast results indicate that ice data assimilation can improve the forecast skill significantly in the FIO-ESM climate forecast system. The forecasted Arctic integrated ice edge error is reduced by around 1/3 by sea ice data assimilation. Compared with the six near-real-time Arctic sea ice forecast results from the subseasonal-to-seasonal (S2S) Prediction Project, FIO-ESM climate forecast system with LESTKF ice data assimilation has relatively high Arctic sea ice forecast skill in 2018 summer sea ice forecast. Since sea ice thickness in the PIOMAS is updated in time, it is a good choice for data assimilation to improve sea ice prediction skills in the near-real-time Arctic sea ice seasonal prediction.     

Laser desorption/ionization mass spectrometric study of secondary organic aerosol formed from the photooxidation of aromatics

Five aromatic hydrocarbons – benzene, toluene, ethylbenzene, p-xylene and 1,2,4-trimethylbenzene – were selected to investigate the laser desorption/ionization mass spectra of secondary organic aerosols (SOA) resulting from OH-initiated photooxidation of aromatic compounds. The experiments were conducted by irradiating aromatic hydrocarbon/CH₃ONO/NO _X mixtures in a home-made smog chamber. The aerosol time-of-flight mass spectrometer (ATOFMS) was used to measure the aerodynamic size and chemical composition of individual secondary organic aerosol particles in real-time. Experimental results showed that aerosol created by aromatics photooxidation is predominantly in the form of fine particles, which have diameters less than 2.5 μm (i.e. PM2.5), and different aromatic hydrocarbons SOA mass spectra have eight same positive laser desorption/ionization mass spectra peaks: m/z = 18, 29, 43, 44, 46, 57, 67, 77. These mass spectra peaks may come from the fragment ions of the SOA products: oxo-carboxylic acids, aldehydes and ketones, nitrogenated organic compounds, furanoid and aromatic compounds. The possible reaction mechanisms leading to these products were also discussed.