Array-based iterative measurements of SmKS travel times and their constraints on outermost core structure

Vigorous convection in Earth's outer core led to the suggestion that it is chemically homogeneous. However, there is increasing seismic evidence for structural complexities close to the outer core's upper and lower boundaries. Both body waves and normal mode data have been used to estimate a P wave velocity, V_p, at the top of the outer core (the E’ layer), which is lower than that in the Preliminary Reference Earth Model. However, these low V_p models do not agree on the form of this velocity anomaly. One reason for this is the difficulty in retrieving and measuring SmKS arrival times. To address this issue, we propose a novel approach using data from seismic arrays to iteratively measure SmKS-SKKS-differential travel times. This approach extracts individual SmKS signal from mixed waveforms of the SmKS series, allowing us to reliably measure differential travel times. We successfully use this method to measure SmKS time delays from earthquakes in the Fiji‐Tonga and Vanuatu subduction zones. SmKS time delays are measured by waveform cross correlation between SmKS and SKKS, and the cross‐correlation coefficient allows us to access measurement quality. We also apply this iterative scheme to synthetic SmKS seismograms to investigate the 3‐D mantle structure's effects. The mantle structure corrections are not negligible for our data, and neglecting them could bias the V_p estimation of uppermost outer core. After mantle structure corrections, we can still see substantial time delays of S3KS, S4KS, and S5KS, supporting a low V_p at the top of Earth's outer core

Seismic probes of solar interior magnetic structure

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfv\'{e}n waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.Comment: 4 pages, 4 figures, accepted Physical Review Letter

Hard x-ray photon-in-photon-out spectroscopy with lifetime resolution – of XAS, XES, RIXSS and HERFD

Spectroscopic techniques that aim to resolve the electronic configuration and local coordination of a central atom by detecting inner-shell radiative decays following photoexcitation using hard X-rays are presented. The experimental setup requires an X-ray spectrometer based on perfect crystal Bragg optics. The possibilities arising from non-resonant (X-Ray Emission Spectroscopy - XES) and resonant excitation (Resonant Inelastic X-Ray Scattering Spectroscopy – RIXSS, High-Energy-Resolution Fluorescence Detected (HERFD) XAS) are discussed when the instrumental energy broadenings of the primary (beamline) monochromator and the crystal spectrometer for x-ray emission detection are on the order of the core hole lifetimes of the intermediate and final electronic states. The small energy bandwidth in the emission detection yields line-sharpened absorption features. In transition metal compounds, electron-electron interactions as well as orbital splittings and fractional population can be revealed. Combination with EXAFS spectroscopy enables to extent the k-range beyond unwanted absorption edges in the sample that limit the EXAFS range in conventional absorption spectroscopy

An axisymmetric time-domain spectral-element method for full-wave simulations: Application to ocean acoustics

The numerical simulation of acoustic waves in complex 3D media is a key topic in many branches of science, from exploration geophysics to non-destructive testing and medical imaging. With the drastic increase in computing capabilities this field has dramatically grown in the last twenty years. However many 3D computations, especially at high frequency and/or long range, are still far beyond current reach and force researchers to resort to approximations, for example by working in 2D (plane strain) or by using a paraxial approximation. This article presents and validates a numerical technique based on an axisymmetric formulation of a spectral finite-element method in the time domain for heterogeneous fluid-solid media. Taking advantage of axisymmetry enables the study of relevant 3D configurations at a very moderate computational cost. The axisymmetric spectral-element formulation is first introduced, and validation tests are then performed. A typical application of interest in ocean acoustics showing upslope propagation above a dipping viscoelastic ocean bottom is then presented. The method correctly models backscattered waves and explains the transmission losses discrepancies pointed out in Jensen et al. (2007). Finally, a realistic application to a double seamount problem is considered.Comment: Added a reference, and fixed a typo (cylindrical versus spherical

Anelasticity across seismic to tidal timescales: a self-consistent approach

In a pioneering study, Wahr & Bergen developed the widely adopted, pseudo-normal mode framework for predicting the impact of anelastic effects on the Earth's body tides. Lau $\textit{et al.}$ have recently derived an extended normal mode treatment of the problem (as well as a minor variant of the theory known as the direct solution method) that makes full use of theoretical developments in free oscillation seismology spanning the last quarter century and that avoids a series of assumptions and approximations adopted in the traditional theory for predicting anelastic effects. There are two noteworthy differences between these two theories: (1) the traditional theory only considers perturbations to the eigenmodes of an elastic Earth, whereas the new theory augments this set of modes to include the relaxation modes that arise in anelastic behaviour; and (2) the traditional theory approximates the complex perturbation to the tidal Love number as a scaled version of the complex perturbation to the elastic moduli, whereas the new theory computes the full complex perturbation to each eigenmode. In this study, we highlight the above differences using a series of synthetic calculations, and demonstrate that the traditional theory can introduce significant error in predictions of the complex perturbation to the Love numbers due to anelasticity and the related predictions of tidal lag angles. For the simplified Earth models we adopt, the computed lag angles differ by ∼20 per cent. The assumptions in the traditional theory have important implications for previous studies that use model predictions to correct observables for body tide signals or that analyse observations of body tide deformation to infer mantle anelastic structure. Finally, we also highlight the fundamental difference between apparent attenuation (i.e. attenuation inferred from observations or predicted using the above theories) and intrinsic attenuation (i.e. the material property investigated through experiments), where both are often expressed in terms of lag angles or $\textit{Q}$$^{-1}$. In particular, we demonstrate the potentially significant (factor of two or more) bias introduced in estimates of $\textit{Q}$$^{-1}$ and its frequency dependence in studies that have treated $\textit{Q}$$^{-1}$ determined from tidal phase lags or measured experimentally as being equal. The observed or theoretically predicted lag angle (or apparent $\textit{Q}$$^{-1}$) differs from the intrinsic, material property due to inertia, self-gravity and effects associated with the energy budget. By accounting for these differences we derive, for a special case, an expression that accurately maps apparent attenuation predicted using the extended normal mode formalism of Lau $\textit{et al.}$ into intrinsic attenuation. The theory allows for more generalized mappings which may be used to robustly connect observations and predictions of tidal lag angles to results from laboratory experiments of mantle materials.This work was supported by NSF EAR-1464024, NSF EAR-1215061, and Harvard University

Rainfall threshold for hillslope outflow: an emergent property of flow pathway connectivity

Nonlinear relations between rain input and hillslope outflow are common observations in hillslope hydrology field studies. In this paper we use percolation theory to model the threshold relationship between rainfall amount and outflow and show that this nonlinear relationship may arise from simple linear processes at the smaller scale. When the rainfall amount exceeds a threshold value, the underlying elements become connected and water flows out of the base of the hillslope. The percolation approach shows how random variations in storage capacity and connectivity at the small spatial scale cause a threshold relationship between rainstorm amount and hillslope outflow. <br><br> As a test case, we applied percolation theory to the well characterized experimental hillslope at the Panola Mountain Research Watershed. Analysing the measured rainstorm events and the subsurface stormflow with percolation theory, we could determine the effect of bedrock permeability, spatial distribution of soil properties and initial water content within the hillslope. The measured variation in the relationship between rainstorm amount and subsurface flow could be reproduced by modelling the initial moisture deficit, the loss of free water to the bedrock, the limited size of the system and the connectivity that is a function of bedrock topography and existence of macropores. The values of the model parameters were in agreement with measured values of soil depth distribution and water saturation

A network analysis to compare biomarker profiles in patients with and without diabetes mellitus in acute heart failure

Aims: It is unclear whether distinct pathophysiological processes are present among patients with acute heart failure (AHF), with and without diabetes. Network analysis of biomarkers may identify correlative associations that reflect different pathophysiological pathways. Methods and results: We analysed a panel of 48 circulating biomarkers measured within 24 h of admission for AHF in a subset of patients enrolled in the PROTECT trial. In patients with and without diabetes, we performed a network analysis to identify correlations between measured biomarkers. Compared with patients without diabetes (n = 1111), those with diabetes (n = 922) had a higher prevalence of ischaemic heart disease and traditional coronary risk factors. After multivariable adjustment, patients with and without diabetes had significantly different levels of biomarkers across a spectrum of pathophysiological domains, including inflammation (TNFR-1a, periostin), cardiomyocyte stretch (BNP), angiogenesis (VEGFR, angiogenin), and renal function (NGAL, KIM-1) (adjusted P-value &lt;0.05). Among patients with diabetes, network analysis revealed that periostin strongly clustered with C-reactive protein and interleukin-6. Furthermore, renal markers (creatinine and NGAL) closely associated with potassium and glucose. These findings were not seen among patients without diabetes. Conclusion: Patients with AHF and diabetes, compared with those without diabetes, have distinct biomarker profiles. Network analysis suggests that cardiac remodelling, inflammation, and fibrosis are closely associated with each other in patients with diabetes. Furthermore, potassium levels may be sensitive to changes in renal function as reflected by the strong renal–potassium–glucose correlation. These findings were not seen among patients without diabetes and may suggest distinct pathophysiological processes among AHF patients with diabetes

Computing exact P-values for DNA motifs

Motivation: Many heuristic algorithms have been designed to approximate P-values of DNA motifs described by position weight matrices, for evaluating their statistical significance. They often significantly deviate from the true P-value by orders of magnitude. Exact P-value computation is needed for ranking the motifs. Furthermore, surprisingly, the complexity of the problem is unknown. Results: We show the problem to be NP-hard, and present MotifRank, software based on dynamic programming, to calculate exact P-values of motifs. We define the exact P-value on a general and more precise model. Asymptotically, MotifRank is faster than the best exact P-value computing algorithm, and is in fact practical. Our experiments clearly demonstrate that MotifRank significantly improves the accuracy of existing approximation algorithms

Annual survey air quality 2001

The air quality for 2001 is found to strongly resemble that for 2000. The decreasing trend in concentrations seen in the nineties, is still valid on the whole. Concerning the meteorological influences, 2001 can be considered as a normal year. Widespread exceedances occurred for ozone and PM10, and in the deposition of potential acid and total nitrogen. In the urban environment, exceedances of the annual mean limit value for NO2 occurred along road length of about 2000 km. Calculations for benzene and benzo[a]pyrene showed that exceedances were incidental. This annual air quality survey provides an overview of the air quality and the load placed on soils and surface water by atmospheric deposition in the Netherlands, on the basis of measurements and model calculations. The report describes global, photochemical, acidifying, particulate and local air pollution. Special attention has been paid to the Dutch air quality in relation to the new European legislation on air quality.Als we de luchtkwaliteit over 2001 vergelijken met die van 2000 blijken deze in grote lijnen met elkaar overeen te komen. De trendmatige daling van concentraties heeft zich in grote lijnen doorgezet. 2001 was wat meteorologische invloed betreft een normaal jaar. Normoverschrijdingen op landelijke schaal hebben zich voorgedaan van ozon, zwevende deeltjes (PM10), depositie van potentieel zuur en stikstof. In de stedelijke omgeving is de jaargemiddelde norm NO2 overschreden langs een weglengte van circa 2000 km. Voor benzeen en benzo[a]pyreen geven berekeningen aan dat incidenteel overschrijdingen van grenswaarden kunnen zijn voorgekomen. In het jaaroverzicht luchtkwaliteit 2001 wordt op basis van metingen en modelberekeningen een samenvattend beeld gegeven van de luchtkwaliteit en de belasting van bodem en oppervlaktewater door atmosferische depositie in Nederland over 2001. Het rapport beschrijft de onderwerpen: mondiale, fotochemische, verzurende, deeltjesvormige en lokale luchtverontreiniging. Als extra onderwerp is een hoofdstuk opgenomen over de Nederlandse luchtkwaliteit in relatie tot de nieuwe Europese wetgeving voor luchtkwaliteit