Forest operations,tree species composition and decline in rainfall explain runoff changes in the Nacimiento experimental catchments,south central Chile

Few long-term studies have explored how intensively managed short rotation forest plantations interact with climate variability. We examine how prolonged severe drought and forest operations affect runoff in 11 experimental catchments on private corporate forest land near Nacimiento in south central Chile over the period 2008–2019. The catchments (7.7–414 ha) contain forest plantations of exotic fast-growing species (Pinus radiata, Eucalyptus spp.) at various stages of growth in a Mediterranean climate (mean long-term annual rainfall = 1381 mm). Since 2010, a drought, unprecedented in recent history, has reduced rainfall at Nacimiento by 20%, relative to the long-term mean. Pre-drought runoff ratios were <0.2 under 8-year-old Eucalyptus; >0.4 under 21-year-old Radiata pine and >0.8 where herbicide treatments had controlled vegetation for 2 years in 38% of the catchment area. Early in the study period, clearcutting of Radiata pine (85%–95% of catchment area) increased streamflow by 150 mm as compared with the year before harvest, while clearcutting and partial cuts of Eucalyptus did not increase streamflow. During 2008–2019, the combination of emerging drought and forestry treatments (replanting with Eucalyptus after clearcutting of Radiata pine and Eucalyptus) reduced streamflow by 400–500 mm, and regeneration of previously herbicide-treated vegetation combined with growth of Eucalyptus plantations reduced streamflow by 1125 mm (87% of mean annual precipitation 2010–2019). These results from one of the most comprehensive forest catchment studies in the world on private industrial forest land indicate that multiple decades of forest management have reduced deep soil moisture reservoirs. This effect has been exacerbated by drought and conversion from Radiata pine to Eucalyptus, apparently largely eliminating subsurface supply to streamflow. The findings reveal tradeoffs between wood production and water supply, provide lessons for adapting forest management to the projected future drier climate in Chile, and underscore the need for continued experimental work in managed forest plantations.     

湖南省似大地水准面2017模型及精度分析

利用近7万个湖南及邻近省份重力观测数据、502个GNSS/水准控制点及数字高程模型，以EIGEN-6C4全球重力场模型作参考重力场，采用顾及地球曲率影响的各类地形质量位及引力的第二类Helmert凝集法严密算法，利用高分辨率地形数据恢复甚短波扰动重力场，确定空间分辨率2′×2′的高精度湖南省似大地水准面模型（HNGG2017）。经外部检核，模型整体精度均优于±0.022 m。与历史模型相比，新模型在湖南北部常德汉寿、西南部永州江永等地区精度得到显著改善。     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

全球气候变化下南海诸岛保护优先区识别分析

全球变化下,珊瑚礁保护区是保护生物多样性、增强珊瑚礁对气候变暖抵抗力的有效方式,而维持珊瑚礁弹性是其核心内容。针对珊瑚礁最具有威胁性的热压力因子,基于南海1982—2009年卫星观测海表面温度（SST）数据和CMIP5加拿大地球系统模式CanESM2模型预估的2006—2100年南海SST数据构建热压力强度模型,从维持珊瑚礁弹性的角度识别IPCC RCP 4.5和RCP 8.5情景下南海诸岛保护优先区。结果表明：RCP 4.5和RCP 8.5情景下13%左右的南海诸岛珊瑚礁识别为保护优先区。根据热压力强度与珊瑚抵抗力及避难所关系,西沙群岛七连屿和晋卿岛近年观测与未来预估的热压力强度均比较低,在保障其服务功能的基础上建议实施完全保护;东沙群岛东沙环礁和中沙环礁排洪滩近年观测急性热压力强度较高但未来预估热压力强度较低,建议实施50%禁止利用保护;中沙群岛黄岩岛近年观测和未来预估的急性热压力强度均比较低,建议实施50%多用途保护。南沙群岛有14%左右的珊瑚礁识别为保护优先区,根据其热压力强度可实施30%~100%禁止利用保护或30%~50%多用途保护。RCP 4.5和RCP 8.5情景下的南海诸岛保护优先区及保护对策,可为维持珊瑚礁生态弹性及应对全球气候变化提供重要的参考价值。     

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.     

Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.     

Dynamical implications of Jupiter's tropospheric ammonia abundance

Groundbased radio observations indicate that Jupiter's ammonia is globally depleted from 0.6 bars to at least 4-6 bars relative to the deep abundance of ∼3 times solar, a fact that has so far defied explanation. The observations also indicate that (i) the depletion is greater in belts than zones, and (ii) the greatest depletion occurs within Jupiter's local 5-μm hot spots, which have recently been detected at radio wavelengths. Here, we first show that both the global depletion and its belt-zone variation can be explained by a simple model for the interaction of moist convection with Jupiter's cloud-layer circulation. If the global depletion is dynamical in origin, then important endmember models for the belt-zone circulation can be ruled out. Next, we show that the radio observations of Jupiter's 5-μm hot spots imply that the equatorial wave inferred to cause hot spots induces vertical parcel oscillation of a factor of ∼2 in pressure near the 2-bar level, which places important constraints on hot-spot dynamics. Finally, using spatially resolved radio maps, we demonstrate that low-latitude features exceeding ∼4000 km diameter, such as the equatorial plumes and large vortices, are also depleted in ammonia from 0.6 bars to at least 2 bars relative to the deep abundance of 3 times solar. If any low-latitude features exist that contain 3-times-solar ammonia up to the 0.6-bar ammonia condensation level, they must have diameters less than ∼4000 km.     

On the Line Formation in Stellar Magnetized Atmospheres

1 INTRODUCTIONIn the previous paper (on et al. 1999, hereafter Paper I), we investigated the wavelength-dependence of four colltribution functions (CFs) derived from dmerent formal solutions andreferring to different emergellt quantities in the unpolarized case. Because one cannot generallyassign a single formation region to the whole line band in a real stellar atmosphere, e.g., the solaratmosphere, instead, the line formation region can be defined as the layers deviating farthestfrom t…     

Propagation of Alfvén waves in multi-layer solar atmosphere model

In this paper, we idealize the actual solar atmosphere as a multi-isothermal-layer system so as to obtain the energy transmittance of the linear Alfvén wave that propagates through such a system in presence of a uniform oblique magnetic filed. The results indicate that the two-layer model is essentially different to the three-layer one. In the two-layer model, the temperature jump acts as a high pass filter. In the three-layer model, resonant transfer will take place and the transmittance undergoes oscillation as the trigonometric function terms dominate its behavior. For actual solar atmosphere, the result reveals that the lower parts of solar atmosphere are more suitable for those Alfvén waves with period of seconds to transfer their energy.     

Microwave observations of Saturn's rings: anisotropy in directly transmitted and scattered saturnian thermal emission

We present a new Very Large Array (VLA) image of Saturn, made from data taken in October 1998 at a wavelength of λ3.6 cm. The moderate ring opening angle (B≈15°) allows us to explore direct transmission of microwave photons through the A and C rings. We find a strong asymmetry of photons transmitted through the A ring, but not in the C ring, a new diagnostic of wake structure in the ring particles. We also find a weak asymmetry between east and west for the far side of the ansae. To facilitate quantitative comparison between dynamic models of the A ring and radio observations, we extend our Monte Carlo radiative transfer code (described in Dunn et al., 2002, Icarus 160, 132-160) to include idealized wakes. We show the idealized model can reproduce the properties of dynamic simulations in directly transmitted light. We examine the model behavior in directly transmitted and scattered light over a range of physical and geometric wake parameters. Finally, we present a wake model with a plausible set of physical parameters that quantitatively reproduces the observed intensity and asymmetry of the A ring both across the planet and in the ansae.