The Bear Brook Watershed in Maine: Multi-decadal whole-watershed experimental acidification

The Bear Brook Watershed in Maine (BBWM) is a long-term research site established to study the response of forest ecosystem function to environmental disturbances of chronic acidic deposition and ecosystem nitrogen enrichment. Starting in 1989, the West Bear (treated) watershed received bimonthly applications of ammonium sulfate [(NH₄)₂SO₄] fertilizer from above the canopy, whereas East Bear (reference) received ambient deposition. The treatments were stopped in 2016, marking the beginning of the recovery phase. Research at the site has focused on soils, streams, and vegetation. Here, we describe data collected over three decades at the BBWM—input and stream output nutrient fluxes, quantitative soil pits and soil chemistry, and soil temperature and moisture.     

2007-2017年中国沙尘气溶胶的三维分布特征及输送过程

利用最新发布的CALIPSO产品，构建了2007-2017年中国沙尘气溶胶的三维分布，并结合HYSPLIT-4模式和再分析数据，探讨了沙尘的三维输送过程。结果表明：中国的沙尘排放源区主要是塔克拉玛干沙漠和巴丹吉林沙漠，沙尘气溶胶出现频率分别为60%和35%。塔克拉玛干沙漠排放的沙尘主要（50%~70%）停留在源地0~6 000 m高度，少部分向东输送至甘肃和内蒙古；巴丹吉林沙漠排放的沙尘则主要向东输送。中国沙尘排放量在春季最大，向东输送最强；夏季，东亚夏季风限制了沙尘向东输送；秋季，沙尘排放减弱，输送强度和夏季相当；沙尘排放量在冬季最小，输送最弱。夏季，沙尘在输送过程中可被抬升至高度5 000 m以上，春季次之，秋、冬季的沙尘主要在低层大气输送。沙尘在向东输送的过程中被抬升并和当地人为污染物混合变为污染性沙尘，华北地区污染性沙尘出现频率高达30%；输送到海洋的沙尘也会与洋面上（0~3 000 m高度）的海盐气溶胶混合，出现频率约为10%。     

乌鲁木齐大气混合层厚度和稳定度与大气污染的关系

利用乌鲁木齐市4座10层100 m梯度气象塔2013年6月~2014年4月气象观测资料和7个环境监测站[WTBX]AQI[WTBZ]资料，计算并分析了大气混合层厚度和稳定度特征，探讨了大气混合层厚度和稳定度与污染的关系。结果表明：乌鲁木齐市混合层厚度夏季郊区高、城区低，冬季从南郊—城区—北郊随地势降低依次降低；夏季和冬季分别在1 559～1 772 m和526～1 156 m之间。地面至2 km以上每500 m高度间隔统计混合层厚度，500～1 000 m出现频率最多；月变化为6～9月基本在500 m以上，且每个高度区间其概率均超过10%，10月～次年2月1 500 m以上区间概率明显减小；日变化为中午13:00～16:00达到最高值，下午和傍晚迅速下降。白天较大的感热输送提供充足的热力条件，这也体现出白天以不稳定层结为主，夜间则以稳定层结为主。大气稳定度分类结果，夏季郊区和城区不稳定（A～C类）所占比例差不多，冬季北郊稳定（E、F类）所占比较最大、城区最弱。[WTBX]AQI指数冬季最大，从南郊—城区—北郊依次增大，这与采暖期污染物多、南郊比北郊地势高有利于扩散输送有关。总体来看，乌鲁木齐大气混合层厚度空间分布与气象要素、大气稳定度、地形等密切相关，对AQI[WTBZ]指数分布有重要影响，这对近地层大气污染状况预报有着重要的指导意义。     

Stokes漂流对全球海表面温度模拟的影响

Stokes漂流对海洋上混合层中的流场和温度场结构具有不可忽视的作用。本文基于WAVEWATCHⅢ海浪模式模拟的海浪要素计算得到Stokes漂流,将其引入SBPOM模式的动量方程中,从体积输运的角度研究Stokes漂流对全球海表面温度的影响。分析发现Stokes漂流与Stokes输运在全球呈现高纬度强于中低纬度的带状分布特征,且这种流动与输运对全球海表面温度具有降温作用,该降温作用的分布与全球Stokes输运强度相对应,高纬降温作用大于中低纬度,特别是南极绕极流海域平均降温明显大于其余海域,最大降温可达1.5℃,且全球月平均降温超过0.1℃。     

A modified numerical model for moisture-salt transport in unsaturated sandy soil under evaporation

Soil salinization, caused by salt migration and accumulation underneath the soil surface, will corrode structures. To analyze the moisture-salt migration and salt precipitation in soil under evaporation conditions, a mathematical model consisting of a series of theoretical equations is briefly presented. The filling effect of precipitated salts on tortuosity factor and evaporation rate are taken into account in relevant equations. Besides, a transition equation to link the solute transport equation before and after salt precipitation is proposed. Meanwhile, a new relative humidity equation deduced from Pitzer ions model is used to modify the vapor transport flux equation. The results show that the calculated values are in good agreement with the published experimental data, especially for the simulation of volume water content and evaporation rate of Toyoura sand, which confirm the reliability and applicability of the proposed model.     

First evidence of IMF control of Jovian magnetospheric boundary locations: Cassini and Galileo magnetic field measurements compared

The Cassini spacecraft, en route to Saturn, passed close to Jupiter while the Galileo spacecraft was completing its 28th and 29th orbits of Jupiter, thus offering a unique opportunity for direct study of the solar wind-Jovian interaction. Here evidence is given of response of the Jovian magnetopause and bow shock positions to changes of the north-south component of the solar wind magnetic field, a phenomenon long known to occur in equivalent circumstances at Earth. The period analyzed starts with the passage over Cassini of an interplanetary shock far upstream of Jupiter. The shock's arrival at Galileo on the dusk-flank of the magnetosphere caused Galileo to exit into the solar wind. Using inter-spacecraft timing based on the time delay established from the shock arrival at each spacecraft, we point out that Galileo's position with respect to the Jovian bow shock appears to correlate with changes in the disturbed north-south reversing field seen behind the shock. We specifically rule out the alternative of changes in the shape of the bow shock with rotations of the interplanetary magnetic field as the cause.     

Gamma-Hadron Separation using Čerenkov Photon Density Fluctuations

In the atmospheric Čerenkov technique γ-rays are detected against the abundant background produced by hadronic showers. In order to improve the signal to noise ratio of theexperiment, it is necessary to reject a significant fraction of hadronic showers. Traditional background rejection methods based on image shape parameters have been extensively used for the data from imaging telescopes. However, non-imaging Čerenkov telescopes have to develop very different means of statistically identifying and removing cosmic ray events. Some of the parameters, which could be potentially important for non-imaging arrays, are the temporal and spectral differences, the lateral distributions and density fluctuations of Čerenkov photons generated by γ-ray and hadron primaries. Here we study the differences in fluctuations of Čerenkov photon density in the light pool at the observation level from showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the PACT array has been used to evaluate the efficiency of the new technique. Various types of density fluctuations like the short range and medium range fluctuations as well as flatness parameter are studied. The estimated quality factors reflect the efficiencies with which the hadrons can be rejected from the data. Since some of these parameters are independent, the cuts may be applied in tandem and we demonstrate that the proton rejection efficiency of ∼90% can be achieved. Use of density fluctuations is particularly suited for wavefront sampling observations and it seems to be a good technique to improve the signal to noise ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

The effects of topographically-controlled thermal tides in the martian upper atmosphere as seen by the MGS accelerometer

Mars Global Surveyor accelerometer observations of the martian upper atmosphere revealed large variations in density with longitude during northern hemisphere spring at altitudes of 130-160 km, all latitudes, and mid-afternoon local solar times (LSTs). This zonal structure is due to tides from the surface. The zonal structure is stable on timescales of weeks, decays with increasing altitude above 130 km, and is dominated by wave-3 (average amplitude 22% of mean density) and wave-2 (18%) harmonics. The phases of these harmonics are constant with both altitude and latitude, though their amplitudes change significantly with latitude. Near the South Pole, the phase of the wave-2 harmonic changes by 90° with a change of half a martian solar day while the wave-3 phase stays constant, suggesting diurnal and semidiurnal behaviour, respectively. We use a simple application of classical tidal theory to identify the dominant tidal modes and obtain results consistent with those of General Circulation Models. Our method is less rigorous, but simpler, than the General Circulation Models and hence complements them. Topography has a strong influence on the zonal structure.     

Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Single star scidar: atmospheric parameters profiling using the simulated annealing algorithm

In this paper, we present a new method to estimate, for each turbulent layer labelled i , the horizontal wind speed   v ( h _i )  , the standard deviation of the horizontal wind speed fluctuations  σ _v ( h_i )  and the integrated value of   C ² _n   over the thickness  Δ h_i   of the turbulent layer   C ² _n ( h_i )Δ h_i   , where   h_i   is the altitude of the turbulent layer. These parameters are extracted from single star scintillation spatiotemporal cross-correlation functions of atmospheric speckles obtained within the generalized mode. This method is based on the simulated annealing algorithm to find the optimal solution required to solve the problem. Astrophysics parameters for adaptive optics are also calculated using   C ² _n ( h_i )  and   v ( h_i )  values. The results of other techniques support this new method.