Tree‐ring records unveil long‐term influence of the Pacific Decadal Oscillation on snowpack dynamics in the Stikine River basin,northern British Columbia

Long hydroclimate records are essential elements for the assessment and management of changing freshwater resources. These records are especially important in transboundary watersheds where international cooperation is required in the joint planning and management process of shared basins. Dendrochronological techniques were used to develop a multicentury record of April 1 snow water equivalent (SWE) for the Stikine River basin in northern British Columbia, Canada, from moisture‐sensitive white spruce (Picea glauca) tree rings. Explaining 43% of the instrumental SWE variability, to our knowledge, this research represents the first attempt to develop long‐term snowpack reconstructions in northern British Columbia. The results indicated that 15 extreme low April 1 SWE events occurred from 1789 to the beginning of the instrumental record in 1974. The reconstruction record also shows that the occurrence of hydrological extremes in the Stikine River basin is characterized by persistent below‐average periods in SWE consistent with phase shifts of the Pacific Decadal Oscillation (PDO). Spectral analyses indicate a very distinct in‐phase (positive) relationship between the multidecadal frequencies of variability (~40 years) extracted from the SWE tree‐ring reconstruction and other reconstructed winter and spring PDO indices. Comparison of the reconstructed SWE record with other tree‐ring‐derived PDO proxy records shows coherence at multidecadal frequencies of variability. The research has significant implications for regional watershed management by highlighting the hydrological response of the Stikine River basin to prior climate changes.     

On the damping of the bending waves in Saturn’s ring

We study numerically the motion of a single particle in the bending wave of finite thickness in Saturn’s ring. We include the forcing due to the planet, a moon, the coriolis force and the self gravity of the ring. In particular, we compute the variation of the velocity arising due to the variation of the amplitude and the phase of the epicyclic motion across the local vertical height of the ring. We suggest that the dissipation of energy due to the collision of ring particles in this shear layer damps out the bending wave of Saturn’s ring at the 5:3 vertical resonance of Mimas within a distance of 150 km from the site of its launching as is observed in Voyager data.     

Mapping trees in high resolution imagery across large areas using locally variable thresholds guided by medium resolution tree maps

Large area tree maps, important for environmental monitoring and natural resource management, are often based on medium resolution satellite imagery. These data have difficulty in detecting trees in fragmented woodlands, and have significant omission errors in modified agricultural areas. High resolution imagery can better detect these trees, however, as most high resolution imagery is not normalised it is difficult to automate a tree classification method over large areas. The method developed here used an existing medium resolution map derived from either Landsat or SPOT5 satellite imagery to guide the classification of the high resolution imagery. It selected a spatially-variable threshold on the green band, calculated based on the spatially-variable percentage of trees in the existing map of tree cover. The green band proved more consistent at classifying trees across different images than several common band combinations. The method was tested on 0.5 m resolution imagery from airborne digital sensor (ADS) imagery across New South Wales (NSW), Australia using both Landsat and SPOT5 derived tree maps to guide the threshold selection. Accuracy was assessed across 6 large image mosaics revealing a more accurate result when the more accurate tree map from SPOT5 imagery was used. The resulting maps achieved an overall accuracy with 95% confidence intervals of 93% (90–95%), while the overall accuracy of the previous SPOT5 tree map was 87% (86–89%). The method reduced omission errors by mapping more scattered trees, although it did increase commission errors caused by dark pixels from water, building shadows, topographic shadows, and some soils and crops. The method allows trees to be automatically mapped at 5 m resolution from high resolution imagery, provided a medium resolution tree map already exists.     

Aggregate impacts in Saturn's rings

Ring particle aggregates are formed in the outer parts of Saturn's main rings. We study how collisions between aggregates can lead to destruction or coalescence of these aggregates, with local N-body simulations taking into account the dissipative impacts and gravitational forces between particles. Impacts of aggregates with different mass ratios are studied, as well as aggregates that consist of particles with different physical properties. We find that the outcome of the collision is very sensitive to the shape of the aggregate, in the sense that more elongated aggregates are more prone to be destroyed. We were interested in testing the accretion criterion Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171] used in their F ring simulations, according to which accretion requires that the masses of the colliding bodies differ at least by a factor of 100. We confirm that such a critical mass ratio exists. In particular, simulations indicate that the exact critical mass ratio depends on the internal density and elasticity of particles, and the distance from the planet. The zone of transition, defined by the distance where individual particles or small aggregates first start to stick on the larger aggregate, and by the distance where two similar sized aggregates on the average eventually coalesce is only about 5000 km wide, if fixed particle properties are used. The rotational state of the aggregates that form via aggregate collision rapidly reaches synchronous rotation, similarly to the aggregates that form via gradual growth.     

Runoff reconstruction for the Bailong River from tree rings back to AD 1601, reveals changing hydrological signals of China north–south transition zone

We present a runoff reconstruction for the Bailong River based on the composite chronology developed from four sampling sites of Pinus tabulaeformis in the China north–south transition zone. The runoff reconstruction, spanning 1601–2013 CE, was developed by calibrating tree-ring data with the instrumental runoff record. Runoff reconstruction accounted for 44.3% of the actual runoff variance during the common period 1958–2010 and provided a long-term perspective on hydrological change in the China north–south transition zone. In the past 413 years, high- and low-runoff years accounted for 15.50% and 15.98%, respectively. Of all the 17 extreme hydrological events, 14 of them are extremely high-flow years, and 17th century was the wettest period during the past 413 years. The preliminary analysis results show that there is a relationship between our runoff reconstruction and the Atlantic multidecadal oscillation at multidecadal scale. Since the 1990s, runoff in the China north–south transition zone has also seen a significant decrease as a result of dry trends of the source region.     

巴丹吉林南缘近两百年来年降水重建及初步分析

沙漠及其边缘地区生态环境脆弱,对气候变化敏感。但沙漠地区有限的森林资源限制了区域百年到千年尺度上的历史气候变化研究。利用采自巴丹吉林南缘的青海云杉年轮宽度资料,重建了区域近191 a(1815—2005年)来的年降水(前一年7月至当年8月的总降水量)变化序列。重建的相关系数是0.636,方差解释量为40.4%,调整自由度后的解释方差R2adj为0.392。重建结果稳定可靠。分析区域过去年降水变化结果可见,19世纪该区域干湿变化频繁,20世纪前半段主要以干旱为主,干湿转变较少。20世纪20年代的干旱事件在巴丹吉林南缘的干旱持续时间更长。周期分析的结果表明,区域年降水量变化有2 a、4 a、64 a等周期。     

A procedure to analyze nonlinear density waves in Saturn's rings using several occultation profiles

Cassini radio science experiments have provided multiple occultation optical depth profiles of Saturn's rings that can be used in combination to analyze density waves. This paper establishes an accurate procedure of inversion of the wave profiles to reconstruct the wave kinematic parameters as a function of semi-major axis, in the nonlinear regime. This procedure is established using simulated data in the presence of realistic noise perturbations, to control the reconstruction error. It is then applied to the Mimas 5:3 density wave. There are two important concepts at the basis of this procedure. The first one is that it uses the nonlinear representation of density waves, and the second one is that it relies on a combination of optical depth profiles instead of just one profile. A related method to analyze density waves was devised by Longaretti and Borderies [Longaretti, P.-Y., Borderies, N., 1986. Icarus 67, 211-223] to study the nonlinear density wave associated with the Mimas 5:3 resonance, but the single photopolarimetric profile provided limited constraints. Other studies of density waves analyzing Cassini data [Colwell, J.E., Esposito, L.W., 2007. Bull. Am. Astron. Soc. 39, 461; Tiscareno, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34] are based on the linear theory and find inconsistent results from profile to profile. Multiple cuts of the rings are helpful in a fundamental way to ensure the accuracy of the procedure by forcing consistency among the various optical depth profiles. By way of illustration we have applied our procedure to the Mimas 5:3 density wave. We were able to recover precisely the kinematic parameters from the radio experiment occultation data in most of the propagation region; a preliminary analysis of the pressure-corrected dispersion allowed us to determine new but still uncertain values for the opacity (K?0.02 cm²/g) and velocity dispersion of (c₀?0.6 cm/s) in the wave region. Our procedure constitutes the first step in our planned analysis of the density waves of Saturn's rings. It is very accurate and efficient in the far-wave region. However, improvements are required within the first wavelength. The ways in which this method can be used to establish diagnostics of ring physics are outlined.     

Structure of self-gravity wakes in Saturn's A ring as measured by Cassini CIRS

The CIRS infrared spectrometer onboard the Cassini spacecraft has scanned Saturn's A ring azimuthally from several viewing angles since its orbit insertion in 2004. A quadrupolar asymmetry has been detected in this ring at spacecraft elevations ranging between 16° to 37°. Its fractional amplitude decreases from 22% to 8% from 20° to 37° elevations. The patterns observed in two almost complete azimuthal scans at elevations 20° and 36° strongly favor the self-gravity wakes as the origin of the asymmetry. The elliptical, infinite cylinder model of Hedman et al. [Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D., Buratti, B.J., Baines, K.H., Brown, R.H., Clark, R.N., 2007. Astron. J. 133, 2624-2629] can reproduce the CIRS observations well. Such wakes are found to have an average height-to-spacing ratio H/λ=0.1607±0.0002, a width-over-spacing W/λ=0.3833±0.0008. Gaps between wakes, which are filled with particles, have an optical depth τG=0.1231±0.0005. The wakes mean pitch angle ΦW is 70.70°±0.07°, relative to the radial direction. The comparison of ground-based visible data with CIRS observations constrains the A ring to be a monolayer. For a surface mass density of 40 g cm⁻² [Tiscarino, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34], the expected spacing of wakes is λ≈60 m. Their height and width would then be H≈10 m and W≈24 m, values that match the maximum size of particles in this ring as determined from ground-based stellar occultations [French, R.G., Nicholson, P.D., 2000. Icarus 145, 502-523].