Trace element and Nd isotope analyses of apatite in granitoids and metamorphosed granitoids from the eastern Central Asian Orogenic Belt: Implications for petrogenesis and post-magmatic alteration

We present in situ trace element and Nd isotopic data of apatites from metamorphosed and metasomatized (i.e., altered) and unaltered granitoids in the Songnen and Jiamusi massifs in the eastern Central Asian Orogenic Belt, with the aim of fingerprinting granitoid petrogenesis, including both the magmatic and post-magmatic evolution processes. Apatites from altered granitoids (AG) and unaltered granitoids (UG) are characterized by distinct textures and geochemical compositions. Apatites from AG have irregular rim overgrowths and complex internal textures, along with low contents of rare earth elements (REEs), suggesting the re-precipitation of apatite during epidote crystallization and/or leaching of REEs from apatite by metasomatic fluids. ε_Nd(t) values of the these apatites are decoupled from zircon ε_Hf(t) values for most samples, which can be attributed to the higher mobility of Nd as compared to Sm in certain fluids. Apatites from UG are of igneous origin based on their homogeneous or concentric zoned textures and coupled Nd-Hf isotopic compositions. Trace element variations in igneous apatite are controlled primarily by the geochemical composition of the parental melt, fractional crystallization of other REE-bearing minerals, and changes in partition coefficients. Sr contents and Eu/Eu* values of apatites from UG correlate with whole-rock Sr and SiO₂ contents, highlighting the effects of plagioclase fractionation during magma evolution. Apatites from UG can be subdivided into four groups based on REE contents. Group 1 apatites have REE patterns similar to the host granitoids, but are slightly enriched in middle REEs, reflecting the influence of the parental melt composition and REE partitioning. Group 2 apatites exhibit strong light REE depletions, whereas Group 3 apatites are depleted in middle and heavy REEs, indicative of the crystallization of epidote-group minerals and hornblende before and/or during apatite crystallization, respectively. Group 4 apatites are depleted in heavy REEs, but enriched in Sr, which are features of adakites. Some unusual geochemical features of the apatites, including the REE patterns, Sr contents, Eu anomalies, and Nd isotopic compositions, indicate that inherited apatites are likely to retain the geochemical features of their parental magmas, and thus provide a record of small-scale crustal assimilation during magma evolution that is not evident from the whole-rock geochemistry.     

基于单三轴试验与直剪试验的层状软岩模拟

为探究模型试验中层状软岩模拟方法,依托九绵高速层状软岩隧道现场取样的室内试验结果,多次试验确定重晶石粉、石英砂、石膏粉、滑石粉和水的最优配比以模拟软岩基体,采用带孔薄膜模拟层理弱面的黏结作用,并通过直剪试验确定孔隙率,最后对不同层理角度与层理厚度的试样开展直剪与单、三轴试验以反映各向异性情况。结果表明：软岩基体模拟最优配比为0.55:0.15:0.07:0.06:0.17,重晶石粉对强度及破坏变形起决定作用,含量过低易压溃,含量过高产生上下贯通裂缝;采用30%孔洞率的薄膜模拟层理效果最佳;土样强度随层角呈U形变化,层厚减小（不小于2 cm）,强度减小。直剪试验中,45º层理剪切面朝层理方向倾斜,90º层理导致剪切面上下产生裂缝与碎裂,单三轴结果中,0º层理会产生小角度倾斜裂缝,45º层理产生垂直层理面的斜裂缝及二次破裂,90º层理产生顺层理的竖向劈裂;与现场结果对照后确定最优层厚为 3 cm,对层角与层厚的直剪与单三轴试验结果整体一致,揭示土样孔隙压密闭合−弹性−塑性胀裂破坏−徐变的过程。     

不同时期土地利用空间数据同化方法

针对不同时期土地利用空间数据之间的差异以及常用软件的不足,该文基于空间数据同化基本理论及策略,提出了FME与ArcGIS联合的不同时期土地利用空间数据同化方法。详细阐述了该方法的实施步骤及相关原理,并进行了初步的实验验证。实验结果表明,所提方法能够快速高效地消除不同时期土地利用数据之间在数据格式、数据语义、坐标系统等方面的差异,提高数据之间的协调一致性。该研究为不同时期土地利用空间数据的综合应用与分析提供了参考。     

土地变更调查之管理信息标注入库探讨

土地变更调查是土地利用和管理工作重要的基础性工作之一,用地管理信息标注入库是对土地利用现状变更信息的完善补充,是确保土地变更调查数据真实准确的重要措施。本文结合国家相关规程及技术要求,探讨总结了管理信息标注入库的一般工作流程及部分可能出现问题的解决方法。     

Characteristics of Fengyun-4A Satellite Atmospheric Motion Vectors and Their Impacts on Data Assimilation

The high observation efficiency, scanning speed and observation frequency of the Fengyun-4A (FY-4A) satellite indicates the progress of Chinese geostationary meteorological satellites. The characteristics of FY-4A atmospheric motion vectors (AMVs) derived from the high-level water vapor (WV-High) channel, mid-level water vapor (WV-Mid) channel, and infrared (IR) channel of FY-4A are analyzed, and their corresponding observation errors estimated. Then, the impacts of single-channel and multi-channel FY-4A AMVs on RMAPS-ST (the Rapid-refresh Multi-scale Analysis and Prediction System—Short Term) are evaluated based on one-month data assimilation cycling and forecasting experiments. Results show that the observation errors of FY-4A AMVs from the three channels have an explicit vertical structure. Results from the cycling experiments indicate that the assimilation of AMVs from WV-High produces more apparent improvement of the wind in the upper layer, while a more positive effect in the lower layer is achieved by the assimilation of AMVs from IR. Furthermore, the assimilation of AMVs from IR is more skillful for medium and moderate precipitation than from other channels owing to the good quality of data in the lower layer in the AMVs from IR. Assimilation of FY-4A AMVs from the three channels could combine the advantages of assimilation from each individual channel to improve the wind in the upper, middle and lower layers simultaneously.     

Impacts of Multigrid NLS-4DVar-based Doppler Radar Observation Assimilation on Numerical Simulations of Landfalling Typhoon Haikui (2012)

We applied the multigrid nonlinear least-squares four-dimensional variational assimilation(MG-NLS4DVar) method in data assimilation and prediction experiments for Typhoon Haikui(2012) using the Weather Research and Forecasting(WRF) model. Observation data included radial velocity(V_r) and reflectivity(Z) data from a single Doppler radar, quality controlled prior to assimilation. Typhoon prediction results were evaluated and compared between the NLS-4DVar and MG-NLS4DVar methods. Compared with a forecast that began with NCEP analysis data, our radar data assimilation results were clearly improved in terms of structure, intensity, track, and precipitation prediction for Typhoon Haikui(2012). The results showed that the assimilation accuracy of the NLS-4DVar method was similar to that of the MG-NLS4DVar method,but that the latter was more efficient. The assimilation of V_r alone and Z alone each improved predictions of typhoon intensity, track, and precipitation; however, the impacts of V_r data were significantly greater that those of Z data.Assimilation window-length sensitivity experiments showed that a 6-h assimilation window with 30-min assimilation intervals produced slightly better results than either a 3-h assimilation window with 15-min assimilation intervals or a 1-h assimilation window with 6-min assimilation intervals.     

GRAPES区域扰动预报模式动力框架设计及检验

设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法（GCR）求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。      

雷达资料同化对一次浙江初春罕见飑线过程数值模拟的影响分析

以2018年3月4日影响浙江的初春罕见飑线过程为例,利用WRF模式与GSI-3DVar同化系统开展了雷达资料同化对重大强对流天气的影响研究,分析了雷达资料同化对此次飑线过程的模拟改进作用和可能影响过程,对比了雷达反射率因子和雷达径向速度的同化效果,探讨了不同数量和位置的雷达资料同化对模拟效果的影响。结果表明:雷达资料(尤其是反射率因子)同化有效改善了飑线边界层特征的模拟,从而改进了模式对飑线过程中组合反射率因子、降水、大风等发展演变的模拟效果。雷达反射率因子同化通过直接调整水凝物质含量,修正风暴中降水模拟及由此引起的蒸发冷却,形成更接近实况的强冷池,进而产生了比雷达径向速度同化更大的正贡献。相对于同化非关键区域的雷达资料,同化飑线过程上游关键区域的雷达观测所包含的重要中小尺度信息,对飑线过程模拟效果提升更为重要。     

“Pilot chimney” on the waterloo campus

A brief survey is given of the sources, the methods of injection into the atmosphere, and the mechanisms of deposition into bodies of water for various atmospheric pollutants which contribute significantly to the pollution of large bodies of water. Specific examples of the relative importance of atmospheric sources are included where data are available. Programs of the Canadian Atmospheric Environment Service in this field of research are outlined.     

Assimilating Best Track Minimum Sea Level Pressure Data Together with Doppler Radar Data Using an Ensemble Kalman Filter for Hurricane Ike (2008) at a Cloud-Resolving Resolution

Extending an earlier study, the best track minimum sea level pressure (MSLP) data are assimilated for landfalling Hurricane Ike (2008) using an ensemble Kalman filter (EnKF), in addition to data from two coastal ground-based Doppler radars, at a 4-km grid spacing. Treated as a sea level pressure observation, the MSLP assimilation by the EnKF enhances the hurricane warm core structure and results in a stronger and deeper analyzed vortex than that in the GFS (Global Forecast System) analysis; it also improves the subsequent 18-h hurricane intensity and track forecasts. With a 2-h total assimilation window length, the assimilation of MSLP data interpolated to 10-min intervals results in more balanced analyses with smaller subsequent forecast error growth and better intensity and track forecasts than when the data are assimilated every 60 minutes. Radar data are always assimilated at 10-min intervals. For both intensity and track forecasts, assimilating MSLP only outperforms assimilating radar reflectivity (Z) only. For intensity forecast, assimilating MSLP at 10-min intervals outperforms radar radial wind (Vr) data (assimilated at 10-min intervals), but assimilating MSLP at 60-min intervals fails to beat Vr data. For track forecast, MSLP assimilation has a slightly (noticeably) larger positive impact than Vr(Z) data. When Vr or Z is combined with MSLP, both intensity and track forecasts are improved more than the assimilation of individual observation type. When the total assimilation window length is reduced to 1h or less, the assimilation of MSLP alone even at 10-min intervals produces poorer 18-h intensity forecasts than assimilating Vr only, indicating that many assimilation cycles are needed to establish balanced analyses when MSLP data alone are assimilated; this is due to the very limited pieces of information that MSLP data provide.