Optimal, scalable forward models for computing gravity anomalies

We describe three approaches for computing a gravity signal from a density anomaly. The first approach consists of the classical "summation" technique, whilst the remaining two methods solve the Poisson problem for the gravitational potential using either a Finite Element (FE) discretization employing a multilevel preconditioner, or a Green's function evaluated with the Fast Multipole Method (FMM). The methods utilizing the PDE formulation described here differ from previously published approaches used in gravity modeling in that they are optimal, implying that both the memory and computational time required scale linearly with respect to the number of unknowns in the potential field. Additionally, all of the implementations presented here are developed such that the computations can be performed in a massively parallel, distributed memory computing environment. Through numerical experiments, we compare the methods on the basis of their discretization error, CPU time and parallel scalability. We demonstrate the parallel scalability of all these techniques by running forward models with up to $10^8$ voxels on 1000's of cores.Comment: 38 pages, 13 figures; accepted by Geophysical Journal Internationa

Manifestation of Strong Quadrupole Light-Molecule Interaction in the SER and SEHR Spectra of Pyrazine and Phenazine

It is demonstrated that explanation of appearance of strong forbidden lines in the SEHR spectra of pyrazine and phenazine, caused by totally symmetric vibrations and also other details of their SER and SEHR spectra can be made on the base of the dipole-quadrupole theory. The main point of this theory is conception of a strong quadrupole light molecule interaction, arising in surface fields strongly varying in space near a rough metal surface. Existence of the pointed lines is a good corroboration of existence of the strong quadrupole light-molecule interaction.Comment: 31 pages 6 figure

Dynamics of Arctic Permafrost Coasts in the 21st Century

Climate warming is particularly pronounced in the Arctic with temperatures rising twice as much as in the rest of the world. It seems natural that this warming has profound effects on the speed of erosion of Arctic coasts, since the majority consists of permafrost, composed of unlithified material and hold together by ice. Permafrost stores approximately 1307 Gt of carbon, which is almost 60 % more than currently being contained in the atmosphere. Understanding the main drivers and dynamics of permafrost coastal erosion is of global relevance, especially since floods and erosion are both projected to intensify. However, the assessment of the impacts of climate warming on Arctic coasts is impaired by little data availability. We reviewed relevant scientific literature on changing dynamics of Arctic coast, potential drivers of these changes and the impacts on the human and natural environment. We provide a comprehensive overview over the state of the art and share our thoughts on how we envision potential pathways of future Arctic coastal research. We found that the overwhelming majority of all studied Arctic coasts is erosive and that in most cases erosion rates per year are increasing, threatening coastal settlements, infrastructure, cultural sites and archaeological remains. The impacts on the natural environment are also manifold and reach from changing sediment fluxes which limit light availability in the water column to a higher input of carbon and nutrients into the nearshore zone with the potential to influence food chains

The Appearance and Modulation of Osteocyte Marker Expression during Calcification of Vascular Smooth Muscle Cells

Vascular calcification is an indicator of elevated cardiovascular risk. Vascular smooth muscle cells (VSMCs), the predominant cell type involved in medial vascular calcification, can undergo phenotypic transition to both osteoblastic and chondrocytic cells within a calcifying environment.In the present study, using in vitro VSMC calcification studies in conjunction with ex vivo analyses of a mouse model of medial calcification, we show that vascular calcification is also associated with the expression of osteocyte phenotype markers. As controls, the terminal differentiation of murine calvarial osteoblasts into osteocytes was induced in vitro in the presence of calcifying medium (containing ß-glycerophosphate and ascorbic acid), as determined by increased expression of the osteocyte markers DMP-1, E11 and sclerostin. Culture of murine aortic VSMCs under identical conditions confirmed that the calcification of these cells can also be induced in similar calcifying medium. Calcified VSMCs had increased alkaline phosphatase activity and PiT-1 expression, which are recognized markers of vascular calcification. Expression of DMP-1, E11 and sclerostin was up-regulated during VSMC calcification in vitro. Increased protein expression of E11, an early osteocyte marker, and sclerostin, expressed by more mature osteocytes was also observed in the calcified media of Enpp1(-/-) mouse aortic tissue.This study has demonstrated the up-regulation of key osteocytic molecules during the vascular calcification process. A fuller understanding of the functional role of osteocyte formation and specifically sclerostin and E11 expression in the vascular calcification process may identify novel potential therapeutic strategies for clinical intervention

Statistical mechanics of voting

Decision procedures aggregating the preferences of multiple agents can produce cycles and hence outcomes which have been described heuristically as `chaotic'. We make this description precise by constructing an explicit dynamical system from the agents' preferences and a voting rule. The dynamics form a one dimensional statistical mechanics model; this suggests the use of the topological entropy to quantify the complexity of the system. We formulate natural political/social questions about the expected complexity of a voting rule and degree of cohesion/diversity among agents in terms of random matrix models---ensembles of statistical mechanics models---and compute quantitative answers in some representative cases.Comment: 9 pages, plain TeX, 2 PostScript figures included with epsf.tex (ignore the under/overfull \vbox error messages

SARS Exposure and Emergency Department Workers

Of 193 emergency department workers exposed to severe acute respiratory syndrome (SARS), 9 (4.7%) were infected. Pneumonia developed in six workers, and assays showed anti-SARS immunoglobulin (Ig) M and IgG. The other three workers were IgM-positive and had lower IgG titers; in two, mild illness developed, and one remained asymptomatic

Mineralocorticoid receptor antagonist pretreatment to MINIMISE reperfusion injury after ST-elevation myocardial infarction (the MINIMISE STEMI Trial): rationale and study design.

Novel therapies capable of reducing myocardial infarct (MI) size when administered prior to reperfusion are required to prevent the onset of heart failure in ST-segment elevation myocardial infarction (STEMI) patients treated by primary percutaneous coronary intervention (PPCI). Experimental animal studies have demonstrated that mineralocorticoid receptor antagonist (MRA) therapy administered prior to reperfusion can reduce MI size, and MRA therapy prevents adverse left ventricular (LV) remodeling in post-MI patients with LV impairment. With these 2 benefits in mind, we hypothesize that initiating MRA therapy prior to PPCI, followed by 3 months of oral MRA therapy, will reduce MI size and prevent adverse LV remodeling in STEMI patients. The MINIMISE-STEMI trial is a prospective, randomized, double-blind, placebo-controlled trial that will recruit 150 STEMI patients from four centers in the United Kingdom. Patients will be randomized to receive either an intravenous bolus of MRA therapy (potassium canrenoate 200 mg) or matching placebo prior to PPCI, followed by oral spironolactone 50 mg once daily or matching placebo for 3 months. A cardiac magnetic resonance imaging scan will be performed within 1 week of PPCI and repeated at 3 months to assess MI size and LV remodeling. Enzymatic MI size will be estimated by the 48-hour area-under-the-curve serum cardiac enzymes. The primary endpoint of the study will be MI size on the 3-month cardiac magnetic resonance imaging scan. The MINIMISE STEMI trial will investigate whether early MRA therapy, initiated prior to reperfusion, can reduce MI size and prevent adverse post-MI LV remodeling.This research study was funded by the Rosetrees Trust, our local research networks, and the National Institute for Health Research University College London Hospitals Biomedical Research Centre (for a list of key study participants and committee members, see Supporting Information, Appendix, in the online version of this article). G.M.F. was supported by a research grant of the Swiss National Foundation and the SSMBS (Schweizerische Stiftung fur Medizinische und Biologische Stipendien). D.J.H. is supported by a BHF Senior Clinical Research ¨ Fellowship (FS/10/039/28270)

Rapid inversion of data from 2D resistivity surveys with electrode displacements

Resistivity monitoring surveys are used to detect temporal changes in the subsurface using repeated measurements over the same site. The positions of the electrodes are typically measured at the start of the survey program and possibly at occasional later times. In areas with unstable ground, such as landslide‐prone slopes, the positions of the electrodes can be displaced by ground movements. If this occurs at times when the positions of the electrodes are not directly measured, they have to be estimated. This can be done by interpolation or, as in recent developments, from the resistivity data using new inverse methods. The smoothness‐constrained least squares optimisation method can be modified to include the electrode positions as additional unknown parameters. The Jacobian matrices with the sensitivity of the apparent resistivity measurements to changes in the electrode positions are then required by the optimisation method. In this paper, a fast adjoint‐equation method is used to calculate the Jacobian matrices required by the least squares method to reduce the calculation time. In areas with large near‐surface resistivity contrasts, the inversion routine sometimes cannot accurately distinguish between electrode displacements and subsurface resistivity variations. To overcome this problem, the model for the initial time‐lapse dataset (with accurately known electrode positions) is used as the starting model for the inversion of the later‐time dataset. This greatly improves the accuracy of the estimated electrode positions compared to the use of a homogeneous half‐space starting model. In areas where the movement of the electrodes is expected to occur in a fixed direction, the method of transformations can be used to include this information as an additional constraint in the optimisation routine

Boundary line analysis of the effect of water-filled pore space on nitrous oxide emission from cores of arable soil

The boundary line has been proposed as a model of the effects of a variable on a biological response, when this variable might limit the response in only some of a set of observations. It is proposed that the upper boundary (in some circumstances the lower boundary) represents the response function of interest. Boundary-line analysis is a method for estimating this response function from data. The approach has been used to model the emission of N2O from soil in response to various soil properties. However, the methods that have been used to identify the boundary are based on somewhat ad hoc partitions of the data. A statistical model that we have presented previously has not been applied to this problem in soil science, and we do so here to represent how the water-filled pore space (WFPS) of the soil affects the rate of N2O emission. We derive a boundary-line response that can be shown to be a better model for the data than an unbounded alternative by statistical criteria. Furthermore, the fitted boundary-response model is consistent with past empirical observations and modelling studies with respect to both the WFPS at which the potential emission rate is largest and the measurement error for the emission rates themselves. We show how the fitted model might be used to interpret data on soil volumetric water content with respect to seasonal changes in potential emissions, and to compare potential emissions between soil series that have contrasting physical properties