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.     

Ensemble simulation of land evapotranspiration in China based on a multi-forcing and multi-model approach

In order to reduce the uncertainty of offline land surface model(LSM) simulations of land evapotranspiration(ET), we used ensemble simulations based on three meteorological forcing datasets [Princeton, ITPCAS(Institute of Tibetan Plateau Research, Chinese Academy of Sciences), Qian] and four LSMs(BATS, VIC, CLM3.0 and CLM3.5), to explore the trends and spatiotemporal characteristics of ET, as well as the spatiotemporal pattern of ET in response to climate factors over mainland China during 1982–2007. The results showed that various simulations of each member and their arithmetic mean(Ens Mean) could capture the spatial distribution and seasonal pattern of ET sufficiently well, where they exhibited more significant spatial and seasonal variation in the ET compared with observation-based ET estimates(Obs MTE). For the mean annual ET, we found that the BATS forced by Princeton forcing overestimated the annual mean ET compared with Obs MTE for most of the basins in China, whereas the VIC forced by Princeton forcing showed underestimations. By contrast, the Ens Mean was closer to Obs MTE, although the results were underestimated over Southeast China. Furthermore, both the Obs MTE and Ens Mean exhibited a significant increasing trend during 1982–98; whereas after 1998, when the last big EI Ni ?no event occurred, the Ens Mean tended to decrease significantly between 1999 and 2007, although the change was not significant for Obs MTE. Changes in air temperature and shortwave radiation played key roles in the long-term variation in ET over the humid area of China, but precipitation mainly controlled the long-term variation in ET in arid and semi-arid areas of China.     

A modeling study of effective radiative forcing and climate response due to tropospheric ozone

This study simulates the effective radiative forcing(ERF) of tropospheric ozone from 1850 to 2013 and its effects on global climate using an aerosol–climate coupled model, BCC AGCM2.0.1 CUACE/Aero, in combination with OMI(Ozone Monitoring Instrument) satellite ozone data. According to the OMI observations, the global annual mean tropospheric column ozone(TCO) was 33.9 DU in 2013, and the largest TCO was distributed in the belts between 30°N and 45°N and at approximately 30°S; the annual mean TCO was higher in the Northern Hemisphere than that in the Southern Hemisphere;and in boreal summer and autumn, the global mean TCO was higher than in winter and spring. The simulated ERF due to the change in tropospheric ozone concentration from 1850 to 2013 was 0.46 W m~(-2), thereby causing an increase in the global annual mean surface temperature by 0.36℃, and precipitation by 0.02 mm d~(-1)(the increase of surface temperature had a significance level above 95%). The surface temperature was increased more obviously over the high latitudes in both hemispheres, with the maximum exceeding 1.4?C in Siberia. There were opposite changes in precipitation near the equator,with an increase of 0.5 mm d~(-1)near the Hawaiian Islands and a decrease of about-0.6 mm d~(-1)near the middle of the Indian Ocean.     

Incorporation of a dynamic root distribution into CLM4.5: Evaluation of carbon and water fluxes over the Amazon

Roots are responsible for the uptake of water and nutrients by plants and have the plasticity to dynamically respond to different environmental conditions. However, most land surface models currently prescribe rooting profiles as a function only of vegetation type, with no consideration of the surroundings. In this study, a dynamic rooting scheme, which describes root growth as a compromise between water and nitrogen availability, was incorporated into CLM4.5 with carbon–nitrogen(CN) interactions(CLM4.5-CN) to investigate the effects of a dynamic root distribution on eco-hydrological modeling. Two paired numerical simulations were conducted for the Tapajos National Forest km83(BRSa3) site and the Amazon, one using CLM4.5-CN without the dynamic rooting scheme and the other including the proposed scheme. Simulations for the BRSa3 site showed that inclusion of the dynamic rooting scheme increased the amplitudes and peak values of diurnal gross primary production(GPP) and latent heat flux(LE) for the dry season, and improved the carbon(C) and water cycle modeling by reducing the RMSE of GPP by 0.4 g C m~(-2)d~(-1), net ecosystem exchange by 1.96 g C m~(-2)d~(-1), LE by 5.0 W m~(-2), and soil moisture by 0.03 m~3m~(-3), at the seasonal scale, compared with eddy flux measurements, while having little impact during the wet season. For the Amazon, regional analysis also revealed that vegetation responses(including GPP and LE) to seasonal drought and the severe drought of 2005 were better captured with the dynamic rooting scheme incorporated.     

An algorithm to retrieve precipitation with synthetic aperture radar

This paper presents a new type of rainfall retrieval algorithm, called the model-oriented statistical and Volterra integration. It is a combination of the model-oriented statistical (MOS) and Volterra integral equation (VIE) approaches. The steps involved in this new algorithm can be briefly illustrated as follows. Firstly, information such as the start point and width of the rain is obtained through pre-analysis of the data received by synthetic aperture radar (SAR). Secondly, the VIE retrieval algorithm is employed over a short distance to obtain information on the shape of the rain. Finally, the rain rate can be calculated by using the MOS retrieval algorithm. Simulation results show that the proposed algorithm is effective and simple, and can lead to time savings of nearly 50% compared with MOS. An example of application of SAR data is also discussed, involving the retrieval of precipitation information over the South China Sea.     

贵州辐射雾的时空变化特征及其气象要素分析

选取2008—2012年贵州省3646站次的辐射雾天气过程,根据08时能见度大小将其划分为4级雾、3级雾、2级雾和1级雾四个等级,利用常规地面观测资料分析其时空变化特征,并根据辐射雾多发生区选取出修文、三穗和凤冈3个分布在贵州省高速干道上的站点,利用相应地面气象观测台站的逐日逐时实测气象要素资料,分析其不同等级的气象要素特征。结果表明:各等级辐射雾在10月—次年1月发生较多;空间分布呈"东多西少"格局,4级雾主要分布中心为修文和岑巩,其余等级辐射雾主要分布中心为三穗、凤冈、平塘、正安、松桃和锦屏;修文站的4级雾偏多与其平均相对湿度明显偏高、风速小于3 m/s和风向主要为偏东风相关性较高;三穗站和凤冈站的平均相对湿度与当日08时和前一日20时的气温差相对利于雾的生成,但其风速较小,垂直混合较弱,故其主要为较低等级雾。     

海河流域夏季降水与北半球环流异常的相关

使用1981—2013年6—8月的月平均海河流域站点降水资料以及NCEP/NCAR再分析资料,采用SVD方法分析了北半球高度场环流异常与海河流域夏季降水的关系。结果表明夏季海河流域的一致型降水异常主要受热带地区高度场异常的影响,而北半球中高纬的高度异常分布对海河流域东南部降水异常有重要作用。热带及中高纬地区的高度场异常与副热带西风急流和流域上空的水汽输送有一定的对应关系,急流强度和位置的变化以及水汽辐合异常造成海河流域降水异常,尤其中高纬高度异常的这种对应关系更加明显,并且两个地区的高度异常分布对应的海河流域异常水汽输送有不同来源。     

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.     

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.