Implementation of an interior point source in the ultra weak variational formulation through source extraction

The Ultra Weak Variational Formulation (UWVF) is a powerful numerical method for the approximation of acoustic, elastic and electromagnetic waves in the time-harmonic regime. The use of Trefftz-type basis functions incorporates the known wave-like behaviour of the solution in the discrete space, allowing large reductions in the required number of degrees of freedom for a given accuracy, when compared to standard finite element methods. However, the UWVF is not well disposed to the accurate approximation of singular sources in the interior of the computational domain. We propose an adjustment to the UWVF for seismic imaging applications, which we call the Source Extraction UWVF. Differing fields are solved for in subdomains around the source, and matched on the inter-domain boundaries. Numerical results are presented for a domain of constant wavenumber and for a domain of varying sound speed in a model used for seismic imaging.

Enhanced MALDI MS sensitivity by weak base additives and glycerol sample coating

The concept of rationally designing MALDI matrices has been extended to the next “whole sample” level. These studies have revealed some unexpected and exploitable insights in improving MALDI sensitivity. It is shown that (i) additives which only provide additional laser energy absorption are best to be avoided; (ii) the addition of proton donors in the form of protonated weak bases can be highly beneficial; (iii) the addition of glycerol for coating crystalline samples is highly recommended. Overall, analytical sensitivity has been significantly increased compared to the current “gold” standards in MALDI MS, and new insights into the mechanisms and processes of MALDI have been gained.

Beyond the environmentalist's paradox and the debate on weak versus strong sustainability

Environmentalists generally argue that ecological damage will (eventually) lead to declines in human well-being. From this perspective, the recent introduction of the “environmentalist’s paradox” in BioScience by Raudsepp-Hearne and colleagues (2010) is particularly significant. In essence, Raudsepp-Hearne and colleagues (2010) claimed that although ecosystem services have been degraded, human well-being—paradoxically—has increased. In this article, we show that this debate is in fact rooted in a broader discussion on weak sustainability versus strong sustainability(the substitutability of human-made capital for natural capital). We warn against the reductive nature of focusing only on a stock–flow framework in which a natural-capital stock produces ecosystem services. Concretely, we recommend a holistic approach in which the complexity, irreversibility, uncertainty, and ethical predicaments intrinsic to the natural environment and its connections to humanity are also considered.

Long-term variations in the open solar flux and possible links to Earth’s climate

Recent paleoclimate studies provide strong evidence for an association between cosmogenic isotope production and Earth’s climate throughout the holecene. These isotopes are generated by the bombardment of Earth’s atmosphere by galactic cosmic rays, the fluxes of which vary in approximately inverse proportion to the total open magnetic flux of the Sun. This paper discusses how results from the Ulysses spacecraft allow us to quantify the open solar flux from observations of near-Earth interplanetary space and to study its long-term variations using the homogeneous record of geomagnetic activity. A study of the results and of their accuracy is presented. The two proposed mechanisms that could lead to the open solar flux being a good proxy for solar-induced climate change are discussed: the first is the modulation of the production of some types of cloud by the air ions produced by cosmic rays; the second is a variation in the total or spectral solar irradiance, in association with changes in the open flux. Some implications for our understanding of anthropogenic climate change are discussed.

EISCAT observations of unusual flows in the morning sector associated with weak substorm activity

A discussion is given of plasma flows in the dawn and nightside high-latitude ionospheric regions during substorms occurring on a contracted auroral oval, as observed using the EISCAT CP-4-A experiment. Supporting data from the PACE radar, Greenland magnetometer chain, SAMNET magnetometers and geostationary satellites are compared to the EISCAT observations. On 4 October 1989 a weak substorm with initial expansion phase onset signatures at 0030 UT, resulted in the convection reversal boundary observed by EISCAT (at \sim0415 MLT) contracting rapidly poleward, causing a band of elevated ionospheric ion temperatures and a localised plasma density depletion. This polar cap contraction event is shown to be associated with various substorm signatures; Pi2 pulsations at mid-latitudes, magnetic bays in the midnight sector and particle injections at geosynchronous orbit. A similar event was observed on the following day around 0230 UT (\sim0515 MLT) with the unusual and significant difference that two convection reversals were observed, both contracting poleward. We show that this feature is not an ionospheric signature of two active reconnection neutral lines as predicted by the near-Earth neutral model before the plasmoid is “pinched off”, and present two alternative explanations in terms of (1) viscous and lobe circulation cells and (2) polar cap contraction during northward IMF. The voltage associated with the anti-sunward flow between the reversals reaches a maximum of 13 kV during the substorm expansion phase. This suggests it to be associated with the polar cap contraction and caused by the reconnection of open flux in the geomagnetic tail which has mimicked “viscous-like” momentum transfer across the magnetopause.

Weak intermolecular interactions in an ionically bound molecular adsorbate: cyclopentadienyl=Cu(111)

The dissociative adsorption of cyclopentadiene (C5H6) on Cu(111) yields a cyclopentadienyl (Cp) species with strongly anionic characteristics. The Cp potential energy surface and frictional coupling to the substrate are determined from measurements of dynamics of the molecule together with density functional calculations. The molecule is shown to occupy degenerate threefold adsorption sites and molecular motion is characterized by a low diffusional energy barrier of 40 +/- 3 meV with strong frictional dissipation. Repulsive dipole-dipole interactions are not detected despite charge transfer from substrate to adsorbate.

Transition from suppressed to active convection modulated by a weak-temperature gradient derived large-scale circulation

Numerical simulations are performed to assess the influence of the large-scale circulation on the transition from suppressed to active convection. As a model tool, we used a coupled-column model. It consists of two cloud-resolving models which are fully coupled via a large-scale circulation which is derived from the requirement that the instantaneous domain-mean potential temperature profiles of the two columns remain close to each other. This is known as the weak-temperature gradient approach. The simulations of the transition are initialized from coupled-column simulations over non-uniform surface forcing and the transition is forced within the dry column by changing the local and/or remote surface forcings to uniform surface forcing across the columns. As the strength of the circulation is reduced to zero, moisture is recharged into the dry column and a transition to active convection occurs once the column is sufficiently moistened to sustain deep convection. Direct effects of changing surface forcing occur over the first few days only. Afterward, it is the evolution of the large-scale circulation which systematically modulates the transition. Its contributions are approximately equally divided between the heating and moistening effects. A transition time is defined to summarize the evolution from suppressed to active convection. It is the time when the rain rate within the dry column is halfway to the mean value obtained at equilibrium over uniform surface forcing. The transition time is around twice as long for a transition that is forced remotely compared to a transition that is forced locally. Simulations in which both local and remote surface forcings are changed produce intermediate transition times.

Isotopic examination of links between diet, social differentiation and DISH at the post-medieval Carmelite Friary of Aalst, Belgium

Stable isotope ratios (δ13C and δ15N) were measured in human burials from the post-medieval (16th–18th c. AD) Carmelite friary burial grounds at Aalst, a town in Flanders, Belgium. Dietary patterns of 39 adult individuals were analyzed, from a mixed monastic and lay population buried in three different locations, reflecting groups with differing social status. The data show significant variation in the consumption of perhaps meat, but certainly also marine protein between females and males. This result represents a remarkable continuity with medieval dietary patterns, suggesting that the social and economic changes of the early modern period had a limited effect on everyday life. When both sexes were examined together, individuals buried in the cloister garth consumed significantly less marine protein compared to people buried in the church, likely reflecting social stratification. No statistical differences were observed between isotopic values from the church and the cloister alley, suggesting a similarly diverse diet of the monastic part of the buried population and that of the richer lay population. Finally, the hypothesis that diffuse idiopathic skeletal hyperostosis (DISH) is linked to a diet rich in animal protein was tested. No systematic or statistically significant differences between pathological and non-pathological bones from the same individuals affected with DISH were observed, and no statistical differences were found between individuals with DISH and individuals without DISH

A variational approach to retrieve rain rate by combining information from rain gauges, radars, and microwave links

Accurate and reliable rain rate estimates are important for various hydrometeorological applications. Consequently, rain sensors of different types have been deployed in many regions. In this work, measurements from different instruments, namely, rain gauge, weather radar, and microwave link, are combined for the first time to estimate with greater accuracy the spatial distribution and intensity of rainfall. The objective is to retrieve the rain rate that is consistent with all these measurements while incorporating the uncertainty associated with the different sources of information. Assuming the problem is not strongly nonlinear, a variational approach is implemented and the Gauss–Newton method is used to minimize the cost function containing proper error estimates from all sensors. Furthermore, the method can be flexibly adapted to additional data sources. The proposed approach is tested using data from 14 rain gauges and 14 operational microwave links located in the Zürich area (Switzerland) to correct the prior rain rate provided by the operational radar rain product from the Swiss meteorological service (MeteoSwiss). A cross-validation approach demonstrates the improvement of rain rate estimates when assimilating rain gauge and microwave link information.

Elephants in pyjamas: testing the weak central coherence account of autism spectrum disorders using a syntactic disambiguation task

According to the weak central coherence (CC) account individuals with autism spectrum disorders (ASD) exhibit enhanced local processing and weak part-whole integration. CC was investigated in the verbal domain. Adolescents, recruited using a 2 (ASD status) by 2 (language impairment status) design, completed an aural forced choice comprehension task involving syntactically ambiguous sentences. Half the picture targets depicted the least plausible interpretation, resulting in longer RTs across groups. These were assumed to reflect local processing. There was no ASD by plausibility interaction and consequently little evidence for weak CC in the verbal domain when conceptualised as enhanced local processing. Furthermore, there was little evidence that the processing of syntactically ambiguous sentences differed as a function of ASD or language-impairment status.

Is children's reading “good enough”? Links between online processing and comprehension as children read syntactically ambiguous sentences

We monitored 8- and 10-year-old children’s eye movements as they read sentences containing a temporary syntactic ambiguity to obtain a detailed record of their online processing. Children showed the classic garden-path effect in online processing. Their reading was disrupted following disambiguation, relative to control sentences containing a comma to block the ambiguity, although the disruption occurred somewhat later than would be expected for mature readers. We also asked children questions to probe their comprehension of the syntactic ambiguity offline. They made more errors following ambiguous sentences than following control sentences, demonstrating that the initial incorrect parse of the garden-path sentence influenced offline comprehension. These findings are consistent with “good enough” processing effects seen in adults. While faster reading times and more regressions were generally associated with better comprehension, spending longer reading the question predicted comprehension success specifically in the ambiguous condition. This suggests that reading the question prompted children to reconstruct the sentence and engage in some form of processing, which in turn increased the likelihood of comprehension success. Older children were more sensitive to the syntactic function of commas, and, overall, they were faster and more accurate than younger children.

Conditioning of the weak-constraint variational data assimilation problem for numerical weather prediction

4-Dimensional Variational Data Assimilation (4DVAR) assimilates observations through the minimisation of a least-squares objective function, which is constrained by the model flow. We refer to 4DVAR as strong-constraint 4DVAR (sc4DVAR) in this thesis as it assumes the model is perfect. Relaxing this assumption gives rise to weak-constraint 4DVAR (wc4DVAR), leading to a different minimisation problem with more degrees of freedom. We consider two wc4DVAR formulations in this thesis, the model error formulation and state estimation formulation. The 4DVAR objective function is traditionally solved using gradient-based iterative methods. The principle method used in Numerical Weather Prediction today is the Gauss-Newton approach. This method introduces a linearised `inner-loop' objective function, which upon convergence, updates the solution of the non-linear `outer-loop' objective function. This requires many evaluations of the objective function and its gradient, which emphasises the importance of the Hessian. The eigenvalues and eigenvectors of the Hessian provide insight into the degree of convexity of the objective function, while also indicating the difficulty one may encounter while iterative solving 4DVAR. The condition number of the Hessian is an appropriate measure for the sensitivity of the problem to input data. The condition number can also indicate the rate of convergence and solution accuracy of the minimisation algorithm. This thesis investigates the sensitivity of the solution process minimising both wc4DVAR objective functions to the internal assimilation parameters composing the problem. We gain insight into these sensitivities by bounding the condition number of the Hessians of both objective functions. We also precondition the model error objective function and show improved convergence. We show that both formulations' sensitivities are related to error variance balance, assimilation window length and correlation length-scales using the bounds. We further demonstrate this through numerical experiments on the condition number and data assimilation experiments using linear and non-linear chaotic toy models.

Diagnosing links between atmospheric moisture and extreme daily precipitation over the UK

Atmospheric moisture characteristics associated with the heaviest 1% of daily rainfall events affecting regions of the British Isles are analysed over the period 1997–2008. A blended satellite/rain-gauge data set (GPCP-1DD) and regionally averaged daily rain-gauge observations (HadUKP) are combined with the ERA Interim reanalysis. These are compared with simulations from the HadGEM2-A climate model which applied observed sea surface temperature and realistic radiative forcings. Median extreme daily rainfall across the identified events and locations is larger for GPCP (32 mm day−1) than HadUKP and the simulations (∼25 mm day−1). The heaviest observed and simulated daily rainfall events are associated with increased specific humidity and horizontal transport of moisture (median 850 hPa specific humidity of ∼6 g kg−1 and vapour transport of ∼150 g kg−1 m s−1 for both observed and simulated events). Extreme daily rainfall events are less common during spring and summer across much of the British Isles, but in the south east region, they contribute up to 60% of the total number of distinct extreme daily rainfall events during these months. Compared to winter events, the summer events over south east Britain are associated with a greater magnitude and more southerly location of moisture maxima and less spatially extensive regions of enhanced moisture transport. This contrasting dependence of extreme daily rainfall on moisture characteristics implies a range of driving mechanisms that depend upon location and season. Higher spatial and temporal resolution data are required to explore these processes further, which is vital in assessing future projected changes in rainfall and associated flooding.