Sound Synthesis for Contact-Driven Musical Instruments via Discretisation of Hamilton's Equations

Physical modelling of musical instruments involves studying nonlinear interactions between parts of the instrument. These can pose several difficulties concerning the accuracy and stability of numerical algorithms. In particular, when the underlying forces are non-analytic functions of the phase-space variables, a stability proof can only be obtained in limited cases. An approach has been recently presented by the authors, leading to unconditionally stable simulations for lumped collision models. In that study, discretisation of Hamilton’s equations instead of the usual Newton’s equation of motion yields a numerical scheme that can be proven to be energy conserving. In this paper, the above approach is extended to collisions of distributed objects. Namely, the interaction of an ideal string with a flat barrier is considered. The problem is formulated within the Hamiltonian framework and subsequently discretised. The resulting nonlinearmatrix equation can be shown to possess a unique solution, that enables the update of the algorithm. Energy conservation and thus numerical stability follows in a way similar to the lumped collision model. The existence of an analytic description of this interaction allows the validation of the model’s accuracy. The proposed methodology can be used in sound synthesis applications involving musical instruments where collisions occur either in a confined (e.g. hammer-string interaction, mallet impact) or in a distributed region (e.g. string-bridge or reed-mouthpiece interaction).

Spectral impairment of enclosed higher order ambisonics

Ambisonics is spatial audio technique that attempts to recreate a physical sound field over as large an area as possible. Higher Order Ambisonic systems modelled with near field loudspeakers in free field as well as in a simulated room are investigated. The influence of reflections on the image quality is analysed objectively for both a studio-sized and large reproduction environment using the relative intensity of the reproduced sound field. The results of a simulated enclosed HOA system in the studio-sized room are compared to sound field measurements in the reproduced area.

A Spherical Array Approach for Simulation of Binaural Impulse Responses using the Finite Difference Time Domain Method

The finite difference time domain (FDTD) method has direct applications in musical instrument modeling, simulation of environmental acoustics, room acoustics and sound reproduction paradigms, all of which benefit from auralization. However, rendering binaural impulse responses from simulated
data is not straightforward to accomplish as the calculated pressure at FDTD grid nodes does not contain any directional information. This paper addresses this issue by introducing a spherical array to capture sound pressure on a finite difference grid, and decomposing it into a plane-wave density
function. Binaural impulse responses are then constructed in the spherical harmonics domain by combining the decomposed grid data with free field head-related transfer functions. The effects of designing a spherical array in a Cartesian grid are studied, and emphasis is given to the relationships
between array sampling and the spatial and spectral design parameters of several finite-difference
schemes.

Off-Centre Localisation Performance of Ambisonics and HOA For Large and Small Loudspeaker Array Radii

Ambisonics and Higher Order Ambisonics (HOA) are scalable spatial audio techniques that attempt to present a sound scene to listeners over as large an area as possible. A localisation experiment was carried out to investigate the performance of a first and third order system at three listening positions - one in the centre and two off-centre - using a 5 m radius loudspeaker array. The results are briefly presented and compared to those of an earlier experiment on a 2.2 m loudspeaker array. In both experiments the off-centre listeners were placed such that the ratio of distance from the centre to the array radius was constant in both experiments. The test used a reverse target-pointer adjustment method to determine the error, both signed and absolute, for each combination of listening position and system. The results for both arrays were found to be very similar, suggesting that the relative amplitude of the loudspeakers, which were the same in both cases, was more dominant for localisation than the arrival time differences, which differed between array sizes.

Rapid climate change did not cause population collapse at the end of the European Bronze Age

The impact of rapid climate change on contemporary human populations is of global concern. To contextualize our understanding of human responses to rapid climate change it is necessary to examine the archeological record during past climate transitions. One episode of abrupt climate change has been correlated with societal collapse at the end of the northwestern European Bronze Age. We apply new methods to interrogate archeological and paleoclimate data for this transition in Ireland at a higher level of precision than has previously been possible. We analyze archeological 14C dates to demonstrate dramatic population collapse and present high-precision proxy climate data, analyzed through Bayesian methods, to provide evidence for a rapid climatic transition at ca. 750 calibrated years B.C. Our results demonstrate that this climatic downturn did not initiate population collapse and highlight the nondeterministic nature of human responses to past climate change.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Air pollutant contamination and acidification of surface waters in the North York Moors, UK: Multi-proxy evidence from the sediments of a moorland pool

Despite the extensive geographical range of palaeolimnological studies designed to assess the extent of surface water acidification in the United Kingdom during the 1980s, little attention was paid to the status of surface waters in the North York Moors (NYM). In this paper, we present sediment core data from a moorland pool in the NYM that provide a record of air pollution contamination and surface water acidification. The 41-cm-long core was divided into three lithostratigraphic units. The lower two comprise peaty soils and peats, respectively, that date to between approximately 8080 and 6740 cal. BP. The uppermost unit comprises peaty lake muds dating from between approximately ad 1790 and the present day (ad 2006). The lower two units contain pollen dominated by forest taxa, whereas the uppermost unit contains pollen indicative of open landscape conditions similar to those of the present. Heavy metal, spheroidal carbonaceous particle, mineral magnetics and stable isotope analysis of the upper sediments show clear evidence of contamination by air pollutants derived from fossil-fuel combustion over the last c. 150years, and diatom analysis indicates that the naturally acidic pool became more acidic during the 20th century. We conclude that the exceptionally acidic surface waters of the pool at present (pH=c. 4.1) are the result of a long history of air pollution and not because of naturally acidic local conditions. We argue that the highly acidic surface waters elsewhere in the NYM are similarly acidified and that the lack of evidence of significant recovery from acidification, despite major reductions in the emissions of acidic gases that have taken place over the last c. 30years, indicates the continuing influence of pollutant sulphur stored in catchment peats, a legacy of over 150years of acid deposition.

Sea-level changes in Iceland and the influence of the North Atlantic Oscillation during the last half millennium

We present a new, diatom-based sea-level reconstruction for Iceland spanning the last -500 years, and investigate the possible mechanisms driving the sea-level changes. A sea-level reconstruction from near the Icelandic low pressure system is important as it can improve understanding of ocean-atmosphere forcing on North Atlantic sea-level variability over multi-decadal to centennial timescales. Our reconstruction is from Viarhólmi salt marsh in Snæfellsnes in western Iceland, a site from where we previously obtained a 2000-yr record based upon less precise sea-level indicators (salt-marsh foraminifera). The 20th century part of our record is corroborated by tide-gauge data from Reykjavik. Overall, the new reconstruction shows ca0.6m rise of relative sea level during the last four centuries, of which ca0.2m occurred during the 20th century. Low-amplitude and high-frequency sea-level variability is super-imposed on the pre-industrial long-term rising trend of 0.65m per 1000 years. Most of the relative sea-level rise occurred in three distinct periods: AD 1620-1650, AD 1780-1850 and AD 1950-2000, with maximum rates of ~3±2mm/yr during the latter two of these periods. Maximum rates were achieved at the end of large shifts (from negative to positive) of the winter North Atlantic Oscillation (NAO) Index as reconstructed from proxy data. Instrumental data demonstrate that a strong and sustained positive NAO (a deep Icelandic Low) generates setup on the west coast of Iceland resulting in rising sea levels. There is no strong evidence that the periods of rapid sea-level rise were caused by ocean mass changes, glacial isostatic adjustment or regional steric change. We suggest that wind forcing plays an important role in causing regional-scale coastal sea-level variability in the North Atlantic, not only on (multi-)annual timescales, but also on multi-decadal to centennial timescales.

Sediment accumulation rates in subarctic lakes: insights into age-depth modeling from 22 dated lake records from the Northwest Territories, Canada

Age-depth modeling using Bayesian statistics requires well-informed prior information about the behavior of sediment accumulation. Here we present average sediment accumulation rates (represented as deposition times, DT, in yr/cm) for lakes in an Arctic setting, and we examine the variability across space (intra- and inter-lake) and time (late Holocene). The dataset includes over 100 radiocarbon dates, primarily on bulk sediment, from 22 sediment cores obtained from 18 lakes spanning the boreal to tundra ecotone gradients in subarctic Canada. There are four to twenty-five radiocarbon dates per core, depending on the length and character of the sediment records. Deposition times were calculated at 100-year intervals from age-depth models constructed using the ‘classical’ age-depth modeling software Clam. Lakes in boreal settings have the most rapid accumulation (mean DT 20 ± 10 years), whereas lakes in tundra settings accumulate at moderate (mean DT 70 ± 10 years) to very slow rates, (>100 yr/cm). Many of the age-depth models demonstrate fluctuations in accumulation that coincide with lake evolution and post-glacial climate change. Ten of our sediment cores yielded sediments as old as c. 9,000 cal BP (BP = years before AD 1950). From between c. 9,000 cal BP and c. 6,000 cal BP, sediment accumulation was relatively rapid (DT of 20 to 60 yr/cm). Accumulation slowed between c. 5,500 and c. 4,000 cal BP as vegetation expanded northward in response to warming. A short period of rapid accumulation occurred near 1,200 cal BP at three lakes. Our research will help inform priors in Bayesian age modeling.

The Impact of Extension and Ecosystem Services on Smallholders’ Credit Constraint.

On formal credit markets, access to formal credit and reasonable credit terms of smallholder farmers
in rural sub-Saharan Africa is limited due to adverse selection. Financial institutions operating in
rural areas often cannot distinguish between borrowers (farmers) that are creditworthy and those that
are not, thus, allocate limited resource to agriculture to reduce credit risk. In the presence of limited business quality signaling by smallholder farmers, financial institutions shall demand for collateral and/or offer unfavorable contract terms. Moreover, agricultural productivity of rural sub-Saharan
Africa, dominated by subsistence or small-scale farmers, is also negatively impacted by the adverse
effect of climate change. A strategy that may make the farming practices of smallholder farmer’s
climate resilient and profitable may also improve smallholder farmer's access to formal credit. This
study investigates to what extent participating in ecosystem and extension services (EES) programs
signals business quality of smallholders, thus granting them credit accessibility. We collected data
on 210 smallholder farmers in 2013, comprising farmers that receive payments for ecosystem
services (PES) and farm management training from the International Small Group Tree Planting
Program (TIST) Kenya to test the aforementioned theory empirically. We use game theory,
particularly a screening and sorting model, to illustrate the prospects for farmers with EES to access
formal credit and to improve their credit terms given that they receive PES and banking services
training. Furthermore, the PES’ long term duration (10 – 30 years) generates stable cash-flow which
may be perceived as collateral substitute. Results suggest that smallholder farmers in the TIST
program were less likely to be credit constraint compared to non-TIST farmers. Distance to market,
education, livestock and farm income are factors that determine access to credit from microfinance
institutions in rural Kenya. Amongst farmers that have obtained loans, those keeping business records
enjoy more favorable formal credit conditions. These farmers were observed to pay ca. 5 percent less
interest rate in microfinance charges. For TIST farmers, this type of farm management practices may
be attributed to the banking services and other training they receive within the program. While the
availability of classical collateral (farmlands) and PES may reduce interest rate, the latter was found
to be statistically insignificant. This research underlines the importance of an effective extension
services in rural areas of developing countries and the need to improve gains from conservation
agriculture and ensuing PES. The benefits associated with EES and PES may encompass agricultural
financing.

Conserved histidine of metal transporter AtNRAMP1 is crucial for optimal plant growth under manganese deficiency at chilling temperatures

Manganese (Mn) is an essential nutrient required for plant growth, in particular in the process of photosynthesis. Plant performance is influenced by various environmental stresses including contrasting temperatures, light or nutrient deficiencies. The molecular responses of plants exposed to such stress factors in combination are largely unknown. 

Screening of 108 Arabidopsis thaliana (Arabidopsis) accessions for reduced photosynthetic performance at chilling temperatures was performed and one accession (Hog) was isolated. Using genetic and molecular approaches, the molecular basis of this particular response to temperature (GxE interaction) was identified. 

Hog showed an induction of a severe leaf chlorosis and impaired growth after transfer to lower temperatures. We demonstrated that this response was dependent on the nutrient content of the soil. Genetic mapping and complementation identified NRAMP1 as the causal gene. Chlorotic phenotype was associated with a histidine to tyrosine (H239Y) substitution in the allele of Hog NRAMP1. This led to lethality when Hog seedlings were directly grown at 4 degrees C. 

Chemical complementation and hydroponic culture experiments showed that Mn deficiency was the major cause of this GxE interaction. For the first time, the NRAMP-specific highly conserved histidine was shown to be crucial for plant performance.

Hydroclimatic shifts in northeast Thailand during the last two millennia - the record of Lake Pa Kho

The Southeast Asian mainland is located in the central path of the Asian summer monsoon, a region where paleoclimatic data are still sparse. Here we present a multi-proxy (TOC, C/N, δ¹³C, biogenic silica, and XRF elemental data) study of a 1.5m sediment/peat sequence from Lake Pa Kho, northeast Thailand, which is supported by 20 AMS ¹⁴C ages. Hydroclimatic reconstructions for Pa Kho suggest a strengthened summer monsoon between BC 170-AD 370, AD 800-960, and after AD 1450; and a weakening of the summer monsoon between AD 370-800, and AD 1300-1450. Increased run-off and a higher nutrient supply after AD 1700 can be linked to agricultural intensification and land-use changes in the region. This study fills an important gap in data coverage with respect to summer monsoon variability over Southeast Asia during the past 2000 years and enables the mean position of the Intertropical Convergence Zone (ITCZ) to be inferred based on comparisons with other regional studies. Intervals of strengthened/weaker summer monsoon rainfall suggest that the mean position of the ITCZ was located as far north as 35°N between BC 170-AD 370 and AD 800-960, whereas it likely did not reach above 17°N during the drought intervals of AD 370-800 and AD 1300-1450. The spatial pattern of rainfall variation seems to have changed after AD 1450, when the inferred moisture history for Pa Kho indicates a more southerly location of the mean position of the summer ITCZ.

Tephra from andesitic Shiveluch volcano, Kamchatka, NW Pacific: chronology of explosive eruptions and geochemical fingerprinting of volcanic glass

The ~16-ka-long record of explosive eruptions from Shiveluch volcano (Kamchatka, NW Pacific) is refined using geochemical fingerprinting of tephra and radiocarbon ages. Volcanic glass from 77 prominent Holocene tephras and four Late Glacial tephra packages was analyzed by electron microprobe. Eruption ages were estimated using 113 radiocarbon dates for proximal tephra sequence. These radiocarbon dates were combined with 76 dates for regional Kamchatka marker tephra layers into a single Bayesian framework taking into account the stratigraphic ordering within and between the sites. As a result, we report ~1,700 high-quality glass analyses from Late Glacial–Holocene Shiveluch eruptions of known ages. These define the magmatic evolution of the volcano and provide a reference for correlations with distal fall deposits. Shiveluch tephras represent two major types of magmas, which have been feeding the volcano during the Late Glacial–Holocene time: Baidarny basaltic andesites and Young Shiveluch andesites. Baidarny tephras erupted mostly during the Late Glacial time (~16–12.8 ka BP) but persisted into the Holocene as subordinate admixture to the prevailing Young Shiveluch andesitic tephras (~12.7 ka BP–present). Baidarny basaltic andesite tephras have trachyandesite and trachydacite (SiO₂ < 71.5 wt%) glasses. The Young Shiveluch andesite tephras have rhyolitic glasses (SiO₂ > 71.5 wt%). Strongly calc-alkaline medium-K characteristics of Shiveluch volcanic glasses along with moderate Cl, CaO and low P₂O₅ contents permit reliable discrimination of Shiveluch tephras from the majority of other large Holocene tephras of Kamchatka. The Young Shiveluch glasses exhibit wave-like variations in SiO₂ contents through time that may reflect alternating periods of high and low frequency/volume of magma supply to deep magma reservoirs beneath the volcano. The compositional variability of Shiveluch glass allows geochemical fingerprinting of individual Shiveluch tephra layers which along with age estimates facilitates their use as a dating tool in paleovolcanological, paleoseismological, paleoenvironmental and archeological studies. Electronic tables accompanying this work offer a tool for statistical correlation of unknown tephras with proximal Shiveluch units taking into account sectors of actual tephra dispersal, eruption size and expected age. Several examples illustrate the effectiveness of the new database. The data are used to assign a few previously enigmatic wide-spread tephras to particular Shiveluch eruptions. Our finding of Shiveluch tephras in sediment cores in the Bering Sea at a distance of ~600 km from the source permits re-assessment of the maximum dispersal distances for Shiveluch tephras and provides links between terrestrial and marine paleoenvironmental records.