Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

双线偏振雷达资料在数值模式中的应用: 模拟器的构建

基于Rayleigh-Gans散射原理构造一个S波段双线偏振雷达模拟器。模拟器考虑了云冰、雪、雨和雹4种水凝物,以水凝物的混合比和数浓度以及水凝物粒子的轴比、相对介电常数、下落倾角为输入量,计算得到水平/垂直偏振反射率因子、差分反射率、比差分相位等偏振量。通过二维理想飑线模拟的试验结果表明,模拟器合理地再现了二维理想飑线系统成熟期的主要偏振特征:雹的反射率较高,差分反射率较低(又称Zdr洞);对流云区的Zdr柱;层云区的反射率和比差分相位的0℃层亮带特征;雨滴反射率与差分反射率因子的“雨线”统计特征。利用该模拟器建立了模式变量和偏振雷达观测的联系,有助于未来将偏振雷达观测资料应用于模式预报效果评估及对流尺度资料同化等方面的研究和应用。      

多普勒天气雷达双PRF径向速度资料质量控制

双PRF技术可用于扩展多普勒天气雷达不模糊速度的间隔,但是它也会导致径向速度资料受到离散奇异值的污染,在使用径向速度资料之前,需要对其进行质量控制。通过计算每个库的径向速度与其局地中值速度的差值,定量分析了CINRAD-SA、SC、CD和CC型业务多普勒天气雷达在双PRF模式下的径向速度数据质量,在此基础上提出了一种三步式的双PRF径向速度资料质量控制算法,即首先基于回波尺度大小对径向速度资料中的孤立噪声进行过滤,其次基于局地连续性原则对落在扩展不模糊速度间隔之外的模糊速度进行速度退模糊,最后基于已知的双PRF径向速度数据错误特征进行奇异值修正,并利用CINRAD-SA、SC、CD和CC型多普勒天气雷达的双PRF径向速度资料对该算法进行了测试。数据分析结果表明绝大部分点的速度都遵循连续性原则,与其局地中值速度的偏差较小或几乎没有偏差,而那些奇异值点的速度与其局地中值速度的偏差和高、低PRF的最大不模糊速度对应。对算法进行测试的结果表明该算法能有效地处理模糊速度并合理地修正速度奇异值。使用本算法对双PRF径向速度资料进行质量控制后,在保持原有数据空间分辨率不变的情况,径向速度数据质量得到显著提高。     

2014—2016年超强厄尔尼诺事件的气候影响

2014—2016年,赤道中东太平洋发生了一次超强厄尔尼诺事件,此次事件于2015年11月达到峰值,12月开始衰减。其峰值强度超过了1951年以来另外两次超强厄尔尼诺事件（1982/1983年和1997/1998年）的强度,成为了1951年以来最强的事件。截至2016年4月,事件已持续20个月,也成为了1951年以来持续时间最长的厄尔尼诺事件。在这次事件的发展过程中,热带太平洋至东亚副热带地区的大气环流表现出了显著的响应特征：赤道中东太平洋对流活动加强,异常上升运动发展,而赤道西太平洋对流活动受抑制,异常下沉运动控制;菲律宾附近异常反气旋生成并发展加强,西太平洋副热带高压强度偏强、西伸脊点异常偏西,尤其2015年冬季副热带高压强度为1980年以来最强。与此同时,2015年秋、冬季,我国长江以南大部降水偏多,尤其冬季华南地区（广东、广西、海南三省区）平均降水量达历史第一,较常年偏多1.6倍以上。近期,超强厄尔尼诺正处于衰减阶段,但是考虑到热带印度洋暖海温的“接力”作用,厄尔尼诺事件对2016年春、夏季我国气候异常的影响可能仍将持续。     

下扬子PTB界线深水相区粘土岩的火山成因研究及其对LPME的指示意义

全球,尤其是特提斯域二叠—三叠纪之交(PTB)剖面中普遍发育火山成因的粘土岩,对理解晚二叠世末生物大灭绝(LPME)的触发机制及相关基础地质问题具有重要意义。本文报道了在华南下扬子区新发现两条深水PTB剖面(皖南牛山和蔡村)的粘土岩工作,包括岩石学、矿物学和全岩地球化学等,填补了区域研究空白。研究结果发现,粘土岩主要由伊利石等粘土矿物,以及石英、岩浆锆石、长石等斑晶矿物组成。在地球化学上具有高K_2O、低Na_2O、相对富集大离子亲石元素(LILE)、相对亏损高场强元素(HFSE)等特征。据此,认为这些粘土岩为火山成因的钾质斑脱岩,其原岩可能为中酸性流纹英安岩,具有弧岩浆作用的源区背景。对比华南其他地区已发现的PTB界线粘土岩,发现它们成因类似,可能来源于古特提斯洋周缘陆陆碰撞,抑或是泛大洋俯冲潘吉亚大陆东缘(包括华南板块)所导致的长英质火山岩浆喷发,且具有多期多源性特点。PTB时期全球活跃的火山岩浆活动(包括镁铁质火山作用和华南地区火山灰所指征的长英质火山作用)可能是导致LPME的主要原因。     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

利用引力梯度数据计算径向梯度的优化方法

根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。     

体素模型的战场环境建模研究

针对基于面模型的战场环境建模难以展现战场环境要素的内部属性和规律等的不足,该文提出了利用体素模型实现战场环境建模的基本方法。该方法能够记录战场环境各要素的不同属性,反映其内部属性和规律,准确地模拟战场地理环境的变换过程,可实现战场环境各要素的统一建模。文章利用体素模型实现了对战场环境的认知与抽象,对战场环境体素数据模型进行了分类。最后以战场地形环境要素为例,利用战场环境基本体元实现了对战场地形环境的体素建模。     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.