TY - JOUR A1 - Keller, Matthias A1 - Liu, Shiping A1 - Peyerimhoff, Norbert T1 - A note on eigenvalue bounds for non-compact manifolds JF - Mathematische Nachrichten N2 - In this article we prove upper bounds for the Laplace eigenvalues lambda(k) below the essential spectrum for strictly negatively curved Cartan-Hadamard manifolds. Our bound is given in terms of k(2) and specific geometric data of the manifold. This applies also to the particular case of non-compact manifolds whose sectional curvature tends to -infinity, where no essential spectrum is present due to a theorem of Donnelly/Li. The result stands in clear contrast to Laplacians on graphs where such a bound fails to be true in general. KW - Cheeger inequality KW - eigenvalues KW - Laplacian KW - negative curvature KW - Riemannian manifold Y1 - 2021 U6 - https://doi.org/10.1002/mana.201900209 SN - 0025-584X SN - 1522-2616 VL - 294 IS - 6 SP - 1134 EP - 1139 PB - Wiley-VCH CY - Weinheim ER - TY - JOUR A1 - Saynisch-Wagner, Jan A1 - Bärenzung, Julien A1 - Hornschild, Aaron A1 - Irrgang, Christopher A1 - Thomas, Maik T1 - Tide-induced magnetic signals and their errors derived from CHAMP and Swarm satellite magnetometer observations JF - Earth, planets and space : EPS N2 - Satellite-measured tidal magnetic signals are of growing importance. These fields are mainly used to infer Earth's mantle conductivity, but also to derive changes in the oceanic heat content. We present a new Kalman filter-based method to derive tidal magnetic fields from satellite magnetometers: KALMAG. The method's advantage is that it allows to study a precisely estimated posterior error covariance matrix. We present the results of a simultaneous estimation of the magnetic signals of 8 major tides from 17 years of Swarm and CHAMP data. For the first time, robustly derived posterior error distributions are reported along with the reported tidal magnetic fields. The results are compared to other estimates that are either based on numerical forward models or on satellite inversions of the same data. For all comparisons, maximal differences and the corresponding globally averaged RMSE are reported. We found that the inter-product differences are comparable with the KALMAG-based errors only in a global mean sense. Here, all approaches give values of the same order, e.g., 0.09 nT-0.14 nT for M2. Locally, the KALMAG posterior errors are up to one order smaller than the inter-product differences, e.g., 0.12 nT vs. 0.96 nT for M2. KW - Tides KW - Electromagnetic induction KW - Error covariance KW - Satellite magnetometer observations Y1 - 2021 U6 - https://doi.org/10.1186/s40623-021-01557-3 SN - 1880-5981 VL - 73 IS - 1 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Wormell, Caroline L. A1 - Reich, Sebastian T1 - Spectral convergence of diffusion maps BT - Improved error bounds and an alternative normalization JF - SIAM journal on numerical analysis / Society for Industrial and Applied Mathematics N2 - Diffusion maps is a manifold learning algorithm widely used for dimensionality reduction. Using a sample from a distribution, it approximates the eigenvalues and eigenfunctions of associated Laplace-Beltrami operators. Theoretical bounds on the approximation error are, however, generally much weaker than the rates that are seen in practice. This paper uses new approaches to improve the error bounds in the model case where the distribution is supported on a hypertorus. For the data sampling (variance) component of the error we make spatially localized compact embedding estimates on certain Hardy spaces; we study the deterministic (bias) component as a perturbation of the Laplace-Beltrami operator's associated PDE and apply relevant spectral stability results. Using these approaches, we match long-standing pointwise error bounds for both the spectral data and the norm convergence of the operator discretization. We also introduce an alternative normalization for diffusion maps based on Sinkhorn weights. This normalization approximates a Langevin diffusion on the sample and yields a symmetric operator approximation. We prove that it has better convergence compared with the standard normalization on flat domains, and we present a highly efficient rigorous algorithm to compute the Sinkhorn weights. KW - diffusion maps KW - graph Laplacian KW - Sinkhorn problem KW - kernel methods Y1 - 2021 U6 - https://doi.org/10.1137/20M1344093 SN - 0036-1429 SN - 1095-7170 VL - 59 IS - 3 SP - 1687 EP - 1734 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - JOUR A1 - Chang, Der-Chen A1 - Khalil, Sara A1 - Schulze, Bert-Wolfgang T1 - Analysis on regular corner spaces JF - The journal of geometric analysis N2 - We establish a new approach of treating elliptic boundary value problems (BVPs) on manifolds with boundary and regular corners, up to singularity order 2. Ellipticity and parametrices are obtained in terms of symbols taking values in algebras of BVPs on manifolds of corresponding lower singularity orders. Those refer to Boutet de Monvel's calculus of operators with the transmission property, see Boutet de Monvel (Acta Math 126:11-51, 1971) for the case of smooth boundary. On corner configuration operators act in spaces with multiple weights. We mainly study the case of upper left entries in the respective 2 x 2 operator block-matrices of such a calculus. Green operators in the sense of Boutet de Monvel (Acta Math 126:11-51, 1971) analogously appear in singular cases, and they are complemented by contributions of Mellin type. We formulate a result on ellipticity and the Fredholm property in weighted corner spaces, with parametrices of analogous kind. KW - Boutet de Monvel's calculus KW - Pseudo-differential operators KW - Singular cones KW - Mellin symbols with values in the edge calculus KW - Parametrices of elliptic operators KW - Kegel space Y1 - 2021 U6 - https://doi.org/10.1007/s12220-021-00614-3 SN - 1050-6926 SN - 1559-002X VL - 31 IS - 9 SP - 9199 EP - 9240 PB - Springer CY - New York ER - TY - JOUR A1 - Yen, Ming-Hsuan A1 - von Specht, Sebastian A1 - Lin, Yen-Yu A1 - Cotton, Fabrice A1 - Ma, Kuo-Fong T1 - Within- and between-event variabilities of strong-velocity pulses of moderate earthquakes within dense seismic arrays JF - Bulletin of the Seismological Society of America N2 - Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering. However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood. In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects. We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes. We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations. Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect. We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances. Y1 - 2021 U6 - https://doi.org/10.1785/0120200376 SN - 0037-1106 SN - 1943-3573 VL - 112 IS - 1 SP - 361 EP - 380 PB - Seismological Society of America CY - El Cerito, Calif. ER - TY - JOUR A1 - Rodríguez Zuluaga, Juan A1 - Stolle, Claudia A1 - Yamazaki, Yosuke A1 - Xiong, Chao A1 - England, Scott L. T1 - A synoptic-scale wavelike structure in the nighttime equatorial ionization anomaly JF - Earth and Space Science : ESS N2 - Both ground- and satellite-based airglow imaging have significantly contributed to understanding the low-latitude ionosphere, especially the morphology and dynamics of the equatorial ionization anomaly (EIA). The NASA Global-scale Observations of the Limb and Disk (GOLD) mission focuses on far-ultraviolet airglow images from a geostationary orbit at 47.5 degrees W. This region is of particular interest at low magnetic latitudes because of the high magnetic declination (i.e., about -20 degrees) and proximity of the South Atlantic magnetic anomaly. In this study, we characterize an exciting feature of the nighttime EIA using GOLD observations from October 5, 2018 to June 30, 2020. It consists of a wavelike structure of a few thousand kilometers seen as poleward and equatorward displacements of the EIA-crests. Initial analyses show that the synoptic-scale structure is symmetric about the dip equator and appears nearly stationary with time over the night. In quasi-dipole coordinates, maxima poleward displacements of the EIA-crests are seen at about +/- 12 degrees latitude and around 20 and 60 degrees longitude (i.e., in geographic longitude at the dip equator, about 53 degrees W and 14 degrees W). The wavelike structure presents typical zonal wavelengths of about 6.7 x 10(3) km and 3.3 x 10(3) km. The structure's occurrence and wavelength are highly variable on a day-to-day basis with no apparent dependence on geomagnetic activity. In addition, a cluster or quasi-periodic wave train of equatorial plasma depletions (EPDs) is often detected within the synoptic-scale structure. We further outline the difference in observing these EPDs from FUV images and in situ measurements during a GOLD and Swarm mission conjunction. KW - equatorial ionization anomaly KW - equatorial ionosphere KW - equatorial plasma bubbles KW - wave structure KW - forcing from below Y1 - 2021 U6 - https://doi.org/10.1029/2020EA001529 SN - 2333-5084 VL - 8 IS - 2 PB - American Geophysical Union CY - Malden, Mass. ER - TY - JOUR A1 - Eshghi, Nasim A1 - Mach, Thomas A1 - Reichel, Lothar T1 - New matrix function approximations and quadrature rules based on the Arnoldi process JF - Journal of computational and applied mathematics N2 - The Arnoldi process can be applied to inexpensively approximate matrix functions of the form f (A)v and matrix functionals of the form v*(f (A))*g(A)v, where A is a large square non-Hermitian matrix, v is a vector, and the superscript * denotes transposition and complex conjugation. Here f and g are analytic functions that are defined in suitable regions in the complex plane. This paper reviews available approximation methods and describes new ones that provide higher accuracy for essentially the same computational effort by exploiting available, but generally not used, moment information. Numerical experiments show that in some cases the modifications of the Arnoldi decompositions proposed can improve the accuracy of v*(f (A))*g(A)v about as much as performing an additional step of the Arnoldi process. KW - Arnoldi process KW - Matrix function approximation KW - Quadrature rule Y1 - 2021 U6 - https://doi.org/10.1016/j.cam.2021.113442 SN - 0377-0427 SN - 1879-1778 VL - 391 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Pohle, Jennifer A1 - Adam, Timo A1 - Beumer, Larissa T1 - Flexible estimation of the state dwell-time distribution in hidden semi-Markov models JF - Computational statistics & data analysis N2 - Hidden semi-Markov models generalise hidden Markov models by explicitly modelling the time spent in a given state, the so-called dwell time, using some distribution defined on the natural numbers. While the (shifted) Poisson and negative binomial distribution provide natural choices for such distributions, in practice, parametric distributions can lack the flexibility to adequately model the dwell times. To overcome this problem, a penalised maximum likelihood approach is proposed that allows for a flexible and data-driven estimation of the dwell-time distributions without the need to make any distributional assumption. This approach is suitable for direct modelling purposes or as an exploratory tool to investigate the latent state dynamics. The feasibility and potential of the suggested approach is illustrated in a simulation study and by modelling muskox movements in northeast Greenland using GPS tracking data. The proposed method is implemented in the R-package PHSMM which is available on CRAN. KW - Penalized likelihood KW - Smoothing KW - Time series KW - Animal movement modeling Y1 - 2022 U6 - https://doi.org/10.1016/j.csda.2022.107479 SN - 0167-9473 SN - 1872-7352 VL - 172 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Pathiraja, Sahani Darschika A1 - Leeuwen, Peter Jan van T1 - Multiplicative Non-Gaussian model error estimation in data assimilation JF - Journal of advances in modeling earth systems : JAMES N2 - Model uncertainty quantification is an essential component of effective data assimilation. Model errors associated with sub-grid scale processes are often represented through stochastic parameterizations of the unresolved process. Many existing Stochastic Parameterization schemes are only applicable when knowledge of the true sub-grid scale process or full observations of the coarse scale process are available, which is typically not the case in real applications. We present a methodology for estimating the statistics of sub-grid scale processes for the more realistic case that only partial observations of the coarse scale process are available. Model error realizations are estimated over a training period by minimizing their conditional sum of squared deviations given some informative covariates (e.g., state of the system), constrained by available observations and assuming that the observation errors are smaller than the model errors. From these realizations a conditional probability distribution of additive model errors given these covariates is obtained, allowing for complex non-Gaussian error structures. Random draws from this density are then used in actual ensemble data assimilation experiments. We demonstrate the efficacy of the approach through numerical experiments with the multi-scale Lorenz 96 system using both small and large time scale separations between slow (coarse scale) and fast (fine scale) variables. The resulting error estimates and forecasts obtained with this new method are superior to those from two existing methods. KW - model uncertainty KW - non-Gaussian KW - data-driven KW - uncertainty KW - quantification KW - Lorenz 96 KW - sub-grid scale Y1 - 2022 U6 - https://doi.org/10.1029/2021MS002564 SN - 1942-2466 VL - 14 IS - 4 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Hethey, Christoph Philipp A1 - Hartung, Niklas A1 - Wangorsch, Gaby A1 - Weisser, Karin A1 - Huisinga, Wilhelm T1 - Physiology-based toxicokinetic modelling of aluminium in rat and man JF - Archives of toxicology : official journal of EUROTOX N2 - A sufficient quantitative understanding of aluminium (Al) toxicokinetics (TK) in man is still lacking, although highly desirable for risk assessment of Al exposure. Baseline exposure and the risk of contamination severely limit the feasibility of TK studies administering the naturally occurring isotope Al-27, both in animals and man. These limitations are absent in studies with Al-26 as a tracer, but tissue data are limited to animal studies. A TK model capable of inter-species translation to make valid predictions of Al levels in humans-especially in toxicological relevant tissues like bone and brain-is urgently needed. Here, we present: (i) a curated dataset which comprises all eligible studies with single doses of Al-26 tracer administered as citrate or chloride salts orally and/or intravenously to rats and humans, including ultra-long-term kinetic profiles for plasma, blood, liver, spleen, muscle, bone, brain, kidney, and urine up to 150 weeks; and (ii) the development of a physiology-based (PB) model for Al TK after intravenous and oral administration of aqueous Al citrate and Al chloride solutions in rats and humans. Based on the comprehensive curated Al-26 dataset, we estimated substance-dependent parameters within a non-linear mixed-effect modelling context. The model fitted the heterogeneous Al-26 data very well and was successfully validated against datasets in rats and humans. The presented PBTK model for Al, based on the most extensive and diverse dataset of Al exposure to date, constitutes a major advancement in the field, thereby paving the way towards a more quantitative risk assessment in humans. KW - PBTK KW - Toxicokinetics KW - Al-26 KW - Aluminium Y1 - 2021 U6 - https://doi.org/10.1007/s00204-021-03107-y SN - 0340-5761 SN - 1432-0738 VL - 95 IS - 9 SP - 2977 EP - 3000 PB - Springer CY - Berlin ; Heidelberg ER -