Quality Control and Evaluation of the Observed Daily Data in the North American Soil Moisture Database

The North American Soil Moisture Database (NASMD) was initiated in 2011 to assemble and homogenize in situ soil moisture measurements from 32 observational networks in the United States and Canada encompassing more than 1800 stations. Although statistical quality control (QC) procedures have been applied in the NASMD, the soil moisture content tends to be systematically underestimated by in situ sensors in frozen soils, and using a single maximum threshold (i.e., 0.6 m3 m-3) may not be sufficient for robust QC because of the diverse soil textures in North America. In this study, based on the in situ soil porosity and North American Land Data Assimilation System phase 2 (NLDAS-2) Noah soil temperature, the simple automated QC method is revised to supplement the existing QC approach. This revised QC method is first validated based on the assessment at 78 of the Soil Climate Analysis Network (SCAN) stations where the manually checked data are available, and is then applied to all stations in the NASMD to produce a more strict quality-controlled dataset. The results show that the revised automated QC procedure can flag the spurious and erroneous soil moisture measurements for the SCAN stations, especially for those located in high altitudes and latitudes. Relative to station measurements in the original NASMD, the quality-controlled data show a slightly better agreement with the manually checked soil moisture content. It should be noted that this quality-controlled dataset may be over-flagged for some valid soil moisture measurements due to potential errors of the soil temperature and soil porosity data, and validation in this study is limited by the availability of benchmark soil moisture data. The updated QC and additional validation will be desirable to boost confidence in the product when high-quality data become available in the future.     

Current Status of Global Ocean Reanalysis Datasets

With the development of the ocean satellite remote sensing technology, the reanalysis of past oceanic observations using modern data assimilation technique and the restructuring of the long-term and consistent gridded data products have made great progress. Such datasets provide us with the most primary research tools to identify the state and evolution of ocean, and understand the role of ocean in climate change and variability at different spatial-temporal scales. In this paper, the current research status in the global reanalysis datasets including some of international global ocean reanalysis projects and the corresponding reanalyzed products were systematically reviewed. In addition, the present status of the domestic research of ocean reanalysis datasets was briefly introduced. The validation of the reanalysis datasets and some quality problems represented by the reanalyzed products, and the Ocean Reanalyses Intercomparison Project were systematically reviewed. Moreover, the prospects of the studies of oceanic reanalysis in the future were also discussed in this paper.     

A new methodology for the extension of the impact of data assimilation on ocean wave prediction

It is a common fact that the majority of today's wave assimilation platforms have a limited, in time, ability of affecting the final wave prediction, especially that of long-period forecasting systems. This is mainly due to the fact that after “closing” the assimilation window, i.e., the time that the available observations are assimilated into the wave model, the latter continues to run without any external information. Therefore, if a systematic divergence from the observations occurs, only a limited portion of the forecasting period will be improved. A way of dealing with this drawback is proposed in this study: A combination of two different statistical tools—Kolmogorov–Zurbenko and Kalman filters—is employed so as to eliminate any systematic error of (a first run of) the wave model results. Then, the obtained forecasts are used as artificial observations that can be assimilated to a follow-up model simulation inside the forecasting period. The method was successfully applied to an open sea area (Pacific Ocean) for significant wave height forecasts using the wave model WAM and six different buoys as observational stations. The results were encouraging and led to the extension of the assimilation impact to the entire forecasting period as well as to a significant reduction of the forecast bias.     

WRF模式同化系统在“碧利斯”台风暴雨数值模拟中的应用

利用美国高分辨率中尺度模式WRF(weacherResearch and Forccast)式和、WRF三维同化系统(WRF 3DVAR),以2006年"碧利斯"台风低压引发的暴雨天气过程为例,通过控制试验和同化试验的对比分析,探讨了高空和地面实况资料同化对台风低压"碧利斯"暴雨过程分析和预报的影响.初步的结果显示,同化高空和地面实况资料后对模式的初始场有明显的改进、对暴雨过程的降水落区和强度有不同程度的正反馈,更接近实况的降水.     

Komatiites and nickel sulfide ores of the Black Swan area,Yilgarn Craton,Western Australia. 3: Komatiite geochemistry,and implications for ore forming processes

The Black Swan komatiite sequence is a package of dominantly olivine-rich cumulates with lesser volumes of spinifex textured rocks, interpreted as a section through an extensive komatiite lava flow field. The sequence hosts a number of nickel sulfide orebodies, including the Silver Swan massive shoot and the Cygnet and Black Swan disseminated orebodies. A large body of whole rock analyses on komatiitic rocks from the Black Swan area has been filtered for metasomatic effects. With the exception of mobile elements such as Ca and alkalis, most samples retain residual igneous geochemistry, and can be modelled predominantly by fractionation and accumulation of olivine. Whole rock MgO–FeO relationships imply a relatively restricted range of olivine compositions, more primitive than the olivine which would have been in equilibrium with the transporting komatiite lavas, and together with textural data indicate that much of the cumulus olivine in the sequence was transported. Flow top compositions show evidence for chromite saturation, but the cumulates are deficient in accumulated chromite. Chromite compositions are typical of those found in compound flow-facies komatiites, and are distinct from those in komatiitic dunite bodies. Incompatible trace element abundances show three superimposed influences: control by the relative proportion of olivine to liquid; a signature of crustal contamination and an overprint of metasomatic introduction of LREE, Zr and Th. This overprint is most evident in cumulates, and relatively insignificant in the spinifex rocks. Platinum and palladium behaved as incompatible elements and are negatively correlated with MgO. They show no evidence for wholesale depletion due to sulfide extraction, which was evidently restricted to specific lava tubes or pathways. The lack of correspondence between PGE depletion and contamination by siliceous material implies that contamination alone is insufficient to generate S-saturation and ore formation in the absence of sulfide in the assimilant. Contamination signatures in spinifex-textured rocks may be a guide to Ni-sulfide mineralisation, but are not entirely reliable in the absence of other evidence. The widespread vesicularity of the sequence may be attributable to assimilated water rather than to primary mantle-derived volatiles, and cannot be taken as evidence for primary volatile-rich magmas. The characteristic signature of the Black Swan Succession is the presence of highly localised disseminated sulfide within a sequence showing more widespread evidence for crustal contamination and interaction with its immediate substrate. This has important implications for the applicability of trace element geochemistry in exploration for komatiite-hosted nickel deposits.Electronic Supplementary Material Supplementary material is available in the online version of this article at Editorial handling: Peter Lightfoot     

Assimilation of SLA and SST data into an OGCM for the Indian Ocean

 Remotely sensed observations of sea-level anomaly and sea-surface temperature have been assimilated into an implementation of the Miami Isopycnic Coordinate Ocean Model (MICOM) for the Indian Ocean using the Ensemble Kalman Filter (EnKF). The system has been applied in a hindcast validation experiment to examine the properties of the assimilation scheme when used with a full ocean general circulation model and real observations. This work is considered as a first step towards an operational ocean monitoring and forecasting system for the Indian Ocean. The assimilation of real data has demonstrated that the sequential EnKF can efficiently control the model evolution in time. The use of data assimilation requires a significant amount of additional processing and computational resources. However, we have tried to justify the cost of using a sophisticated assimilation scheme by demonstrating strong regional and temporal dependencies of the covariance statistics, which include highly anisotropic and flow-dependent correlation functions. In particular, we observed a marked difference between error statistics in the equatorial region and at off-equatorial latitudes. We have also demonstrated how the assimilation of SLA and SST improves the model fields with respect to real observations. Independent in situ temperature profiles have been used to examine the impact of assimilating the remotely sensed observations. These intercomparisons have shown that the model temperature and salinity fields better resemble in situ observations in the assimilation experiment than in a model free-run case. On the other hand, it is also expected that assimilation of in situ profiles is needed to properly control the deep ocean circulation. Received: 8 January 2002 / Accepted: 8 April 2002     

Time scales and length scales in magma ?ow pathways and the origin ofmagmatic Ni-Cu-PGE ore deposits

Ore forming processes involve the redistribution of heat, mass and momentum by a wide range of processes operating at different time and length scales. The fastest process at any given length scale tends to be the dominant control. Applying this principle to the array of physical processes that operate within magma flow pathways leads to some key insights into the origins of magmatic Ni-Cu-PGE sulfide ore deposits. A high proportion of mineralised systems, including those in the super-giant Noril'sk-Talnakh camp, are formed in small conduit intrusions where assimilation of country rock has played a major role. Evidence of this process is reflected in the common association of sulfides with varitextured contaminated host rocks containing xenoliths in varying stages of assimilation. Direct incorporation of S-bearing country rock xenoliths is likely to be the dominant mechanism for generating sulfide liquids in this setting. However, the processes of melting or dissolving these xenoliths is relatively slow compared with magma flow rates and, depending on xenolith lithology and the composition of the carrier magma, slow compared with settling and accumulation rates. Chemical equilibration between sulfide droplets and silicate magma is slower still, as is the process of dissolving sulfide liquid into initially undersaturated silicate magmas. Much of the transport and deposition of sulfide in the carrier magmas may occur while sulfide is still incorporated in the xenoliths, accounting for the common association of magmatic sulfide-matrix ore breccias and contaminated "taxitic" host rocks. Effective upgrading of so-formed sulfide liquids would require repetitive recycling by processes such as reentrainment, back flow or gravity flow operating over the lifetime of the magma transport system as a whole. In contrast to mafic-hosted systems, komatiite-hosted ores only rarely show an association with externally-derived xenoliths, an observation which is partially due to the predominant formation of ores in lava flows rather than deep-seated intrusions, but also to the much shorter timescales of key component systems in hotter, less viscous magmas. Nonetheless, multiple cycles of deposition and entrainment are necessary to account for the metal contents of komatiite-hosted sulfides. More generally, the time and length scale approach introduced here may be of value in understanding other igneous processes as well as non-magmatic mineral systems.     

富闪深成岩的成因及其地球动力学意义

富闪深成岩是一套以角闪石为标志矿物的侵入杂岩,可由角闪石岩、角闪辉长岩、角闪闪长岩、(英云)闪长岩、黑云母花岗岩等不同酸度系列的岩石组成。在元素地球化学方面,富闪深成岩可呈低钾拉斑质、钙碱性或钾玄质,富集大离子亲石元素和轻稀土元素。起源于经俯冲作用改造的岩石圈地幔或软流圈地幔部分熔融的基性母岩浆通过分离结晶、同化-混染,或与壳源岩浆混合,形成富闪深成岩系列中不同的岩性端元。高水含量作为富闪岩浆体系的基本特征不仅塑造了杂岩体的标志性角闪石矿物学特征,也主导了岩浆的钙碱性演化趋势。富闪杂岩形成环境相对特殊,通常产于板块会聚终末期的洋脊-海沟交互场景和后碰撞阶段,并受到同时期区域构造作用控制。基于野外产出上的共生性和地球化学的相似性,富闪深成岩与高Ba-Sr花岗岩、晚太古代赞岐岩类存在密切的成因关联,对其进行类比研究,不仅有助于富闪深成岩成因解析,也可为探讨大陆地壳生长-分异机制及地球早期构造体制提供崭新视角。最后,作为显生宙克拉通早期破坏过程的重要岩浆记录,华北克拉通北缘发育的多期晚古生代富闪深成杂岩有望为揭开富闪深成岩成因和地球动力学奥秘提供重要的野外实验室。     

有云条件下HIRS/3资料的同化及其对“珍珠”台风的影响

卫星红外辐射资料汽接同化容易受到云的影响,大量卫星资料在有云情况下无法直接同化,其原因之一就是初始云结构难以获取.而云结构是同化系统观测算子的重要输入参数,缺乏云的结构信息,观测算子尤法模拟有云环境下的亮温.热带地区以及台风、暴雨等极端天气系统与云密不可分,要改进卫星辐射资料特别是红外资料直接同化的效果,需要获得较高精度的云三维结构.在GRAPES-3DVAR(Global and Regional Assimilation and Prediction Enhanced System,全球/区域同化预报系统)基础上,结合WRF(Weather Research and Forecasting Model,天气研究和预报模式)对"珍珠"台风云结构的模拟,增加模拟的云参数作为GRAPES-3DVAR同化系统中RTTOV模式(fast radiative transfer model for TOVS,快速辐射传输模式)的输入参数,改进了有云环境下卫星资料直接同化的观测算子,研究了红外辐射资料直接同化的方法,并分析其对"珍珠"台风预报的影响.针对NOAAl6(National Oceanic and Atmospheric Administration,美国国家海洋和大气管理局)的HIRS/3(High Resolution Infrared Radiation Sounder,高分辨率红外辐射探测仪3型)第4-8、10、13 19通道辐射亮温的同化试验表明,结合初始云结构信息改进同化系统观测算子,可以将红外辎射资料的观测信息同化到模式之中,改进台风温度、湿度结构,对短时降水预报有正面影响.     

内蒙古地区下垫面变化对土壤湿度数值模拟的影响

利用第二次全国土壤调查土壤质地数据（SNSS）和中国区域陆地覆盖资料（CLCV）将陆面过程模式CLM3.5（Community Land Model version 3.5）中基于联合国粮食农业组织发展的土壤质地数据（FAO）和MODIS卫星反演的陆地覆盖数据（MODIS）进行了替换,使用中国气象局陆面数据同化系统（CMA Land Data Assimilation System,CLDAS）大气强迫场资料,分别驱动基于同时改进土壤质地和陆地覆盖数据的CLM3.5（CLM-new）、基于只改进陆地覆盖数据的CLM3.5（CLM-clcv）、基于只改进土壤质地数据的CLM3.5（CLM-snss）和基于原始下垫面数据的CLM3.5（CLM-ctl）,对内蒙古地区2011～2013年土壤湿度的时空变化进行模拟试验,研究下垫面改进对CLM3.5模拟土壤湿度的影响。将四组模拟结果与46个土壤水分站点观测数据进行对比分析,结果表明:相对于控制试验,CLM-clcv、CLM-snss和CLM-new都能不同程度地改进土壤湿度模拟,其中CLM-clcv主要在呼伦贝尔改进明显,CLM-snss则在除呼伦贝尔以外的大部地区改进显著,CLM-ctl模拟的土壤湿度在各层上均系统性偏大,而CLM-new模拟土壤湿度最好地反映出内蒙古地区观测的土壤湿度的时空变化特征,显著改善了土壤湿度的模拟,体现在与观测值有着更高的相关系数和更小的平均偏差与均方根误差。