Late glacial and Holocene environmental change in the Lake Baikal region documented by oxygen isotopes from diatom silica

We investigate late glacial and Holocene climate change recorded in Lake Baikal using the oxygen isotope composition of diatom silica (δ¹⁸O_DIAT). Evaporation from the lake is minor, and the temperature fractionations of δ¹⁸O are unable to explain variations in the δ¹⁸O_DIAT record alone. Isotopically, low meltwater input from glaciers may have some influence on δ¹⁸O_DIAT, but the assumed periods of climatic warming and wastage do not coincide with large shifts in δ¹⁸O_DIAT. There is a gradual oxygen isotope lowering from 27.0‰ to 20.6‰ over the late glacial, while, during the Holocene, δ¹⁸O_DIAT values return to relatively high values. Previous studies of the modern oxygen and hydrogen isotope composition of Lake Baikal's inputs reveal that fluvial input to the lake's North Basin are isotopically lower than fluvial input from South Basin rivers. This north–south gradient of river δ¹⁸O and δD is mainly due to the greater input from isotopically low winter precipitation in the north and isotopically higher summer precipitation in the south. As a result, the δ¹⁸O_DIAT record from Lake Baikal can at least in part be explained by varying input from these sources related to seasonal changes in precipitation. Changes in atmospheric conditions may have a role in altering seasonality and the distribution of precipitation over Lake Baikal's catchment. A feedback mechanism is well known linking higher Eurasian spring snow cover extent (ESSC) to the development of anticyclonic conditions and low precipitation the following summer in the areas south of Lake Baikal. A simultaneous increase in the importance of depleted water (snowmelt) input from the north and decreased enriched summer precipitation in the south is needed to explain depletions in δ¹⁸O of lake water and subsequently δ¹⁸O_DIAT during colder periods. The opposite of this situation is required to enrich lake water during warmer periods. The analysis of δ¹⁸O from diatom silica is a useful proxy for environmental change, especially in lakes, like Lake Baikal, where carbonates are absent or diluted. However, analysis must be based on near pure diatom samples as even trace amounts of silt can have a dominating effect on δ¹⁸O_DIAT values.     

Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r∼118,800 km, 4.5 m at r∼120,700 km) and the inner A Ring (10-15 m for r<127,000 km).     

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21 cm line from the earliest phases of star and galaxy formation in the Universe. This 21 cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.     

The glacilacustrine sedimentary environment of Bowser Lake in the northern Coast Mountains of British Columbia, Canada

Bowser Lake, a fiord lake in the northern Coast Mountains of British Columbia, contains a thick Holocene fill consisting mainly of silt and clay varves. These sediments were carried into the lake by proglacial Bowser River which drains a high-energy, heavily glacierized basin. Sedimentation in the lake is controlled by seasonal snow and ice melt, by autumn rainstorms, and by rare, but very large jökulhlaups from glacier-dammed lakes in the upper Bowser River basin which complicate environmental inferences from the sedimentary record. Sediment is dispersed through the deep western part of the lake by energetic turbidity currents. The turbidity currents apparently do not overtop a sill that separates the western basin from much shallower areas to the east. Large amounts of silt and clay are deposited from suspension in the eastern part of the lake, but sediment accumulation rates there are much lower than to the west. Several strong acoustic reflectors punctuate the varved fill in the western basin; these may be thick or relatively coarse beds deposited during jökulhlaups or exceptionally large storms. The contemporary sediment yield to Bowser Lake, estimated from sediments in the lake, is about 360 t km^-2a^-1. This is a relatively high value, but it is less than yields insome other, similar montane basins with extensive snow and ice cover.The most likely explanation for the difference is that large amounts of sediment have been, and continue to be, stored on the Bowser delta andin small proglacial lakes.     

Storm rainfall conditions for floods and debris flows from recently burned areas in southwestern Colorado and southern California

Debris flows generated during rain storms on recently burned areas have destroyed lives and property throughout the Western U.S. Field evidence indicate that unlike landslide-triggered debris flows, these events have no identifiable initiation source and can occur with little or no antecedent moisture. Using rain gage and response data from five fires in Colorado and southern California, we document the rainfall conditions that have triggered post-fire debris flows and develop empirical rainfall intensity–duration thresholds for the occurrence of debris flows and floods following wildfires in these settings. This information can provide guidance for warning systems and planning for emergency response in similar settings.Debris flows were produced from 25 recently burned basins in Colorado in response to 13 short-duration, high-intensity convective storms. Debris flows were triggered after as little as six to 10 min of storm rainfall. About 80% of the storms that generated debris flows lasted less than 3 h, with most of the rain falling in less than 1 h. The storms triggering debris flows ranged in average intensity between 1.0 and 32.0 mm/h, and had recurrence intervals of two years or less. Threshold rainfall conditions for floods and debris flows sufficiently large to pose threats to life and property from recently burned areas in south-central, and southwestern, Colorado are defined by: I = 6.5D^− 0.7 and I = 9.5D^− 0.7, respectively, where I = rainfall intensity (in mm/h) and D = duration (in hours).Debris flows were generated from 68 recently burned areas in southern California in response to long-duration frontal storms. The flows occurred after as little as two hours, and up to 16 h, of low-intensity (2–10 mm/h) rainfall. The storms lasted between 5.5 and 33 h, with average intensities between 1.3 and 20.4 mm/h, and had recurrence intervals of two years or less. Threshold rainfall conditions for life- and property-threatening floods and debris flows during the first winter season following fires in Ventura County, and in the San Bernardino, San Gabriel and San Jacinto Mountains of southern California are defined by I = 12.5D^−0.4, and I = 7.2D^−0.4, respectively. A threshold defined for flood and debris-flow conditions following a year of vegetative recovery and sediment removal for the San Bernardino, San Gabriel and San Jacinto Mountains of I = 14.0D^−0.5 is approximately 25 mm/h higher than that developed for the first year following fires.The thresholds defined here are significantly lower than most identified for unburned settings, perhaps because of the difference between extremely rapid, runoff-dominated processes acting in burned areas and longer-term, infiltration-dominated processes on unburned hillslopes.     

High spatial resolution hyperspectral mapping of in-stream habitats, depths, and woody debris in mountain streams

This article evaluates the potential of 1-m resolution, 128-band hyperspectral imagery for mapping in-stream habitats, depths, and woody debris in third- to fifth-order streams in the northern Yellowstone region. Maximum likelihood supervised classification using principal component images provided overall classification accuracies for in-stream habitats (glides, riffles, pools, and eddy drop zones) ranging from 69% for third-order streams to 86% for fifth-order streams. This scale dependency of classification accuracy was probably driven by the greater proportion of transitional boundary areas in the smaller streams. Multiple regressions of measured depths (y) versus principal component scores (x₁, x₂,…, x_n) generated R² values ranging from 67% for high-gradient riffles to 99% for glides in a fifth-order reach. R² values were lower in third-order reaches, ranging from 28% for runs and glides to 94% for pools. The less accurate depth estimates obtained for smaller streams probably resulted from the relative increase in the number of mixed pixels, where a wide range of depths and surface turbulence occurred within a single pixel. Matched filter (MF) mapping of woody debris generated overall accuracies of 83% in the fifth-order Lamar River. Accuracy figures for the in-stream habitat and wood mapping may have been misleadingly low because the fine-resolution imagery captured fine-scale variations not mapped by field teams, which in turn generated false “misclassifications” when the image and field maps were compared.The use of high spatial resolution hyperspectral (HSRH) imagery for stream mapping is limited by the need for clear water to measure depth, by any tree cover obscuring the stream, and by the limited availability of airborne hyperspectral sensors. Nonetheless, the high accuracies achieved in northern Yellowstone streams indicate that HSRH imagery can be a powerful tool for watershed-wide mapping, monitoring, and modeling of streams.     

Dependence of normal albedo of lunar regions on wavelength

Normal albedos of thirteen lunar regions are determined at five wavelengths between 4000 and 8000 Å. The deduced values agree with those of the previous investigators except those given by Gehrelset al. (1964) and Van Diggelen (1965). They increase monotonically with the wavelength. At the shortest wavelength, a slight enhancement is detected for most of the regions observed.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Numerical simulations of the tropical air-sea planetary boundary layer

A one-dimensional numerical model of the planetary boundary layer was used to investigate thermal and kinetic energy budgets. The simulation experiments were based on two sets of data. The first set was based on a ‘typical’ June with climatological data extracted for the oceanic region slightly northeast of Barbados. The second set used data from the third phase of project BOMEX, for approximately the same area and time of year as the first set. Comparison with observations of three simulated elements (viz., sea surface temperature and wind and humidity at 6 m) which are important in determining the near-interface energy transports shows that:

1. the model is capable of realistic simulations of both ‘typical’ conditions, and conditions for a specific four-day period;
2. the model is capable of realistically simulating the differences between prevailing values of these parameters in the two cases (‘typical’ and specific four-day period).

The simulated interface fluxes are those of incoming and outgoing short- and long-wave radiation; transmitted radiation at -0.5 m in the ocean, sensible heat transfer into the ocean and air, and latent heat flux of evaporation. Comparison with observational analyses shows that the diurnal variations in net radiation and heat storage in the mixed layer are realistically simulated. The simulated values of evaporation are consistent with other estimates for both ‘typical’ conditions and specific conditions during this four-day period. The rate of heat storage varies between +51 and -37 percent of the diurnal maximum incoming radiation, and the evaporation varies between +16% and -13% of this term. The non-dimensional transfer coefficients (C _D, C_T, C_q) computed from the model show general agreement with the coefficients calculated from observations in the simulated region (Pondet al., 1971). The simulated vertical profiles of temperature are in general agreement with observed profiles, except in the uppermost portions of the atmospheric boundary layer where deviations of approximately 1.5C occur. Simulated vertical profiles of wind speed are generally consistent with observed profiles, with the largest deviations appearing to be of the order of 0.5 m s^-1. Simulated vertical profiles of the eddy fluxes of sensible heat, water vapor, and momentum are generally consistent with Bunker's (1970) aircraft-based measurements of these quantities. The time averages of these simulated profiles show regular decreases with height, while simulated profiles for specific hours of the day show intermediate maxima and minima, which are also seen in the measured profiles. The vertically integrated kinetic energy budgets of the modelled atmospheric layer are presented through the four terms of the kinetic energy budget, viz., the upper and the lower boundary drags, dissipation, and potential-to-kinetic conversion. The dominant terms in the atmospheric energy budgets are the production and dissipation terms, with kinetic energy being exported both to the overlying atmospheric layer and to the underlying oceanic layer at rates of about 2 to 6% of the production, respectively. Comparisons between the climatological and BOMEX simulations are presented. The vertically integrated humidity budgets are presented for the two simulation experiments. Under ‘typical’ conditions, the humidity budget reveals an upper boundary flux of about +29% of the lower boundary flux with the vertically integrated advective flux being -59% of the lower flux. For the specific four-day simulation, the upper boundary flux and advection are about +28 and -70%, respectively, of the lower boundary flux.