Temporal variability in SO<Subscript>2</Subscript> exposures at Hawai’i Volcanoes National Park,USA

Sulfur dioxide measured at the Visitor’s Center of Hawai’i Volcanoes National Park, USA exhibits seasonal and diurnal patterns and is sensitive to rainfall, but the relationship between wind direction and SO₂ is too inconsistent to support deterministic predictions of hourly SO₂ based on hourly wind direction. Although SO₂ at the Visitor Center has usually been below regulatory levels, high SO₂ levels and adverse health effects remain a concern. This investigation identified patterns in hourly SO₂ based on wind direction, time of day, month, and rainfall occurrence using 4 years of hourly data. Empirical probabilities were investigated using a Bayesian approach. Winds from the volcanic vents were rare, but when they did occur SO₂ was elevated about half of the time. Conversely, half of the hours with elevated SO₂ occurred when volcanic vents were not directly upwind. Episodes of elevated SO₂ tended to occur during the months of November–March and between 8:00 a.m. and 5:00 p.m. Rainfall was associated with a marked reduction in SO₂ (29–81% depending on wind direction). Individuals that wish to avoid exposure to volcanic fumes can take these patterns into account.     

Automated placement of supplementary contour lines

Supplementary contour lines are placed between regular contour lines to visualize small but important forms that regular contour lines are unable to show. On topographic maps, typical forms are hillcrests, depressions, saddles, terraces, banks, and levees. No automated method for the selection of supplementary contour lines has been described so far. We document cartographic design principles for the selection of supplementary contour lines for topographic maps, and present an automated method for their placement. Results of the automated method are similar to manually placed supplementary contour lines. Our method helps map authors to create contour line maps that more effectively illustrate relevant small details in maps showing terrain elevation or other scalar fields.     

North Atlantic sea surface temperatures and their relation to the North Atlantic Oscillation during the last 230?years

August Sea Surface Temperatures (aSSTs) based on fossil diatom assemblages are generated with 2?year average resolution from a 230-year-long sediment core (Rapid 21-12B), from the Reykjanes Ridge in the subpolar North Atlantic. The results indicate a warming trend of ~0.5°C of the surface waters at the Reykjanes Ridge for the last 230?years. Superimposed on this warming trend there is a multidecadal to decadal aSST variability of up to 1°C. The interval from the 1770s to the 1830s represents the coldest period, whereas ~1860?C1880 represents the warmest period during the last 230?years. The last 25?years is characterized by a warming trend showing strong decadal aSST variability with several warm years, but also the coldest years since the 1820s. The time of these cold years in the mid-1970s, -1980s and -1990s correspond with the documented great salinity anomalies (GSA) in the North Atlantic suggesting increased fluxes of cold, low-salinity waters from the Arctic during the last decades. The aSST record and the August North Atlantic Oscillation (aNAO) index show similar multidecadal-scale variability indicating a close coupling between the oceanic and atmospheric patterns. The aSST record shows a negative correlation with the aNAO indicating cold aSST during the positive aNAO trend and vice versa. Results suggest that the wind driven variation in volume fluxes of the North Atlantic surface waters could be the major mechanism behind the observed relationship.     

Eigen analysis of tree-ring records: part 2, posing the eigen problem

The technique of expanding meteorological fields on eigenvectors of the field covariation matrix is popular. In this paper, we propose for the first time to use a mathematically similar technique to solve the main problem of dendrochronology: classifying variations in tree-ring records as either age- and microenvironment-dependent or climate-induced. Applying this technique to a sample of very long-lived Qilian junipers (Sabina przewalskii Kom.) from the Dulan region in western China, we demonstrate that the ring-width variations projected on the first eigenvector are age-dependent, but those projected on several of the first subsequent vectors are mainly climate-induced. In particular, the second and third projections capture multi-centennial climatic variations, and the variations projected on the fourth through seventh eigenvectors show periodic variations that are probably induced by the 178-year solar cycle. The projections on the smallest eigenvectors seem to be negligible.     

Models of the wave boundary layer

A general approach to model the structure of the wave boundary layer, based on the nonlinear Reynolds equations in a curvilinear system of coordinates, is described. Both spectral and numerical grid models are used. The energetic interactions between wind and wave in terms of Miles' parameter are studied. For waves outrunning or running against the wind, the range of the inverse flux of energy is found. For waves running slower than the wind, quadratic growth of is established. Vertical profiles of the wave momentum flux for different fetches are given. Following P. Janssen, a one-dimensional analytical model of the wave boundary layer is suggested. The effect of waves on the drag coefficient is analyzed.     

Validation and Verification of a Numerical Model for Tsunami Propagation and Runup

A robust numerical model to simulate propagation and runup of tsunami waves in the framework of non-linear shallow water theory is developed. The numerical code adopts a staggered leapfrog finite-difference scheme to solve the shallow water equations formulated for depth-averaged water fluxes in spherical coordinates. A temporal position of the shoreline is calculated using a free-surface moving boundary algorithm. For large scale problems, the developed algorithm is efficiently parallelized employing a domain decomposition technique. The developed numerical model is benchmarked in an exhaustive series of tests suggested by NOAA. We conducted analytical and laboratory benchmarking for the cases of solitary wave runup on simple beaches, runup of a solitary wave on a conically-shaped island, and the runup in the Monai Valley, Okushiri Island, Japan, during the 1993 Hokkaido-Nansei-Oki tsunami. In all conducted tests the calculated numerical solution is within an accuracy recommended by NOAA standards. We summarize results of numerical benchmarking of the model, its strengths and limits with regards to reproduction of fundamental features of coastal inundation, and also illustrate some possible improvements.     

Interaction of surface waves at very close wavenumbers

We show that interaction of two monochromatic waves at the water surface enters a different dynamic regime if their wavenumbers become very close. The study is conducted by means of a fully nonlinear wave model. In the course of evolution of the two waves, downshifting of the initial wave energy and growth of the first mode occur depending on wave steepness and dk/k. Behaviour of these features changes if dk/k?<?0.0025: both downshifting and growth rate become independent of dk/k, accompanied by rapid transfer of wave energy to large scales.     

Features of the long-term transformation of the Krasnodar reservoir,near the mouth of the Kuban River,Russia

The article considers the long-term(1941–2018) transformation of the Krasnodar valley reservoir, the largest in the North Caucasus. The main functions of the Krasnodar reservoir are irrigation of rice systems and flood protection of land in the Krasnodar reservoir region and the Republic of Adygea. According to topographic maps, Landsat satellite images(1974–2018) and field observations(2016–2018), four stages of transformation of the floodplain reservoir are identified. The selected stages are characterized by both natural causes(the transformation of the filling deltas into the extended deltas, etc.) and man-made causes(runoff diversions in the delta areas, etc.). The key factor of transformation is the formation of deltas of rivers flowing into the reservoir. Each of the selected stages, against the background of a gradual reduction in the area and volume of the reservoir, is characterized by the peculiarities of the formation of river deltas with the formation of genetically homogeneous sections of delta regions. During the period of operation of the reservoir, the delta of the main Kuban River moved up to 32.4 km and took away an area of 35.4 km~2 of the reservoir. During the formation of the deltas of the Kuban and Belaya rivers, a bridge was formed on the Krasnodar reservoir. The evolution of the delta regions led to the division of the reservoir into two autonomous reservoirs. The total area of the delta regions was 85.9 km~2 by 2018, i.e., 21% of the initial area of the reservoir. The transformation of the Krasnodar reservoir leads to a decrease in its regulated volume and gradual degradation.     

Post-glacial diatom and geochemical records of ecological status and water level changes of Lake Vorota,Western Beringia

Journal of Paleolimnology - In this study, radiocarbon-dated geochemical and diatom records from Lake Vorota located within the informal meteorological ‘pole of cold’ of the northern...