@article{PornsawadBoeckmannPanitsupakamon2022, author = {Pornsawad, Pornsarp and B{\"o}ckmann, Christine and Panitsupakamon, Wannapa}, title = {The Levenberg-Marquardt regularization for the backward heat equation with fractional derivative}, series = {Electronic transactions on numerical analysis - ETNA}, volume = {57}, journal = {Electronic transactions on numerical analysis - ETNA}, publisher = {Kent State University}, address = {Kent}, isbn = {978-3-7001-8258-0}, issn = {1068-9613}, doi = {10.1553/etna_vol57s67}, pages = {67 -- 79}, year = {2022}, abstract = {The backward heat problem with time-fractional derivative in Caputo's sense is studied. The inverse problem is severely ill-posed in the case when the fractional order is close to unity. A Levenberg-Marquardt method with a new a posteriori stopping rule is investigated. We show that optimal order can be obtained for the proposed method under a H{\"o}lder-type source condition. Numerical examples for one and two dimensions are provided.}, language = {en} } @article{StachanowNeumannBlankensteinetal.2022, author = {Stachanow, Viktoria and Neumann, Uta and Blankenstein, Oliver and Bindellini, Davide and Melin, Johanna and Ross, Richard and Whitaker, Martin J. J. and Huisinga, Wilhelm and Michelet, Robin and Kloft, Charlotte}, title = {Exploring dried blood spot cortisol concentrations as an alternative for monitoring pediatric adrenal insufficiency patients}, series = {Frontiers in pharmacology}, volume = {13}, journal = {Frontiers in pharmacology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {1663-9812}, doi = {10.3389/fphar.2022.819590}, pages = {8}, year = {2022}, abstract = {Congenital adrenal hyperplasia (CAH) is the most common form of adrenal insufficiency in childhood; it requires cortisol replacement therapy with hydrocortisone (HC, synthetic cortisol) from birth and therapy monitoring for successful treatment. In children, the less invasive dried blood spot (DBS) sampling with whole blood including red blood cells (RBCs) provides an advantageous alternative to plasma sampling. Potential differences in binding/association processes between plasma and DBS however need to be considered to correctly interpret DBS measurements for therapy monitoring. While capillary DBS samples would be used in clinical practice, venous cortisol DBS samples from children with adrenal insufficiency were analyzed due to data availability and to directly compare and thus understand potential differences between venous DBS and plasma. A previously published HC plasma pharmacokinetic (PK) model was extended by leveraging these DBS concentrations. In addition to previously characterized binding of cortisol to albumin (linear process) and corticosteroid-binding globulin (CBG; saturable process), DBS data enabled the characterization of a linear cortisol association with RBCs, and thereby providing a quantitative link between DBS and plasma cortisol concentrations. The ratio between the observed cortisol plasma and DBS concentrations varies highly from 2 to 8. Deterministic simulations of the different cortisol binding/association fractions demonstrated that with higher blood cortisol concentrations, saturation of cortisol binding to CBG was observed, leading to an increase in all other cortisol binding fractions. In conclusion, a mathematical PK model was developed which links DBS measurements to plasma exposure and thus allows for quantitative interpretation of measurements of DBS samples.}, language = {en} } @article{YenvonSpechtLinetal.2022, author = {Yen, Ming-Hsuan and von Specht, Sebastian and Lin, Yen-Yu and Cotton, Fabrice and Ma, Kuo-Fong}, title = {Within- and between-event variabilities of strong-velocity pulses of moderate earthquakes within dense seismic arrays}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {El Cerito, Calif.}, issn = {0037-1106}, doi = {10.1785/0120200376}, pages = {361 -- 380}, year = {2022}, abstract = {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.}, language = {en} } @article{PohleAdamBeumer2022, author = {Pohle, Jennifer and Adam, Timo and Beumer, Larissa}, title = {Flexible estimation of the state dwell-time distribution in hidden semi-Markov models}, series = {Computational statistics \& data analysis}, volume = {172}, journal = {Computational statistics \& data analysis}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0167-9473}, doi = {10.1016/j.csda.2022.107479}, pages = {15}, year = {2022}, abstract = {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.}, language = {en} } @article{PathirajaLeeuwen2022, author = {Pathiraja, Sahani Darschika and Leeuwen, Peter Jan van}, title = {Multiplicative Non-Gaussian model error estimation in data assimilation}, series = {Journal of advances in modeling earth systems : JAMES}, volume = {14}, journal = {Journal of advances in modeling earth systems : JAMES}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1942-2466}, doi = {10.1029/2021MS002564}, pages = {23}, year = {2022}, abstract = {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.}, language = {en} } @article{LilienkampvonSpechtWeatherilletal.2022, author = {Lilienkamp, Henning and von Specht, Sebastian and Weatherill, Graeme and Caire, Giuseppe and Cotton, Fabrice}, title = {Ground-Motion modeling as an image processing task}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120220008}, pages = {1565 -- 1582}, year = {2022}, abstract = {We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.}, language = {en} } @article{DubeBoeckmannRitter2022, author = {Dube, Jonas and B{\"o}ckmann, Christine and Ritter, Christoph}, title = {Lidar-Derived Aerosol Properties from Ny-{\AA}lesund, Svalbard during the MOSAiC Spring 2020}, series = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, volume = {14}, journal = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, number = {11}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs14112578}, pages = {17}, year = {2022}, abstract = {In this work, we present Raman lidar data (from a Nd:YAG operating at 355 nm, 532 nm and 1064 nm) from the international research village Ny-Alesund for the time period of January to April 2020 during the Arctic haze season of the MOSAiC winter. We present values of the aerosol backscatter, the lidar ratio and the backscatter Angstrom exponent, though the latter depends on wavelength. The aerosol polarization was generally below 2\%, indicating mostly spherical particles. We observed that events with high backscatter and high lidar ratio did not coincide. In fact, the highest lidar ratios (LR > 75 sr at 532 nm) were already found by January and may have been caused by hygroscopic growth, rather than by advection of more continental aerosol. Further, we performed an inversion of the lidar data to retrieve a refractive index and a size distribution of the aerosol. Our results suggest that in the free troposphere (above approximate to 2500 m) the aerosol size distribution is quite constant in time, with dominance of small particles with a modal radius well below 100 nm. On the contrary, below approximate to 2000 m in altitude, we frequently found gradients in aerosol backscatter and even size distribution, sometimes in accordance with gradients of wind speed, humidity or elevated temperature inversions, as if the aerosol was strongly modified by vertical displacement in what we call the "mechanical boundary layer". Finally, we present an indication that additional meteorological soundings during MOSAiC campaign did not necessarily improve the fidelity of air backtrajectories.}, language = {en} } @article{MalemShinitskiOjedaOpper2022, author = {Malem-Shinitski, Noa and Ojeda, Cesar and Opper, Manfred}, title = {Variational bayesian inference for nonlinear hawkes process with gaussian process self-effects}, series = {Entropy}, volume = {24}, journal = {Entropy}, number = {3}, publisher = {MDPI}, address = {Basel}, issn = {1099-4300}, doi = {10.3390/e24030356}, pages = {22}, year = {2022}, abstract = {Traditionally, Hawkes processes are used to model time-continuous point processes with history dependence. Here, we propose an extended model where the self-effects are of both excitatory and inhibitory types and follow a Gaussian Process. Whereas previous work either relies on a less flexible parameterization of the model, or requires a large amount of data, our formulation allows for both a flexible model and learning when data are scarce. We continue the line of work of Bayesian inference for Hawkes processes, and derive an inference algorithm by performing inference on an aggregated sum of Gaussian Processes. Approximate Bayesian inference is achieved via data augmentation, and we describe a mean-field variational inference approach to learn the model parameters. To demonstrate the flexibility of the model we apply our methodology on data from different domains and compare it to previously reported results.}, language = {en} } @article{BellingeriFrizPaychaetal.2022, author = {Bellingeri, Carlo and Friz, Peter and Paycha, Sylvie and Preiß, Rosa Lili Dora}, title = {Smooth rough paths, their geometry and algebraic renormalization}, series = {Vietnam journal of mathematics}, volume = {50}, journal = {Vietnam journal of mathematics}, number = {3}, publisher = {Springer}, address = {Singapore}, issn = {2305-221X}, doi = {10.1007/s10013-022-00570-7}, pages = {719 -- 761}, year = {2022}, abstract = {We introduce the class of "smooth rough paths" and study their main properties. Working in a smooth setting allows us to discard sewing arguments and focus on algebraic and geometric aspects. Specifically, a Maurer-Cartan perspective is the key to a purely algebraic form of Lyons' extension theorem, the renormalization of rough paths following up on [Bruned et al.: A rough path perspective on renormalization, J. Funct. Anal. 277(11), 2019], as well as a related notion of "sum of rough paths". We first develop our ideas in a geometric rough path setting, as this best resonates with recent works on signature varieties, as well as with the renormalization of geometric rough paths. We then explore extensions to the quasi-geometric and the more general Hopf algebraic setting.}, language = {en} } @article{DimitrovaKoppitz2022, author = {Dimitrova, Ilinka and Koppitz, J{\"o}rg}, title = {On relative ranks of the semigroup of orientation-preserving transformations on infinite chain with restricted range}, series = {Communications in algebra}, volume = {50}, journal = {Communications in algebra}, number = {5}, publisher = {Taylor \& Francis Group}, address = {Philadelphia}, issn = {0092-7872}, doi = {10.1080/00927872.2021.2000998}, pages = {2157 -- 2168}, year = {2022}, abstract = {Let X be an infinite linearly ordered set and let Y be a nonempty subset of X. We calculate the relative rank of the semigroup OP(X,Y) of all orientation-preserving transformations on X with restricted range Y modulo the semigroup O(X,Y) of all order-preserving transformations on X with restricted range Y. For Y = X, we characterize the relative generating sets of minimal size.}, language = {en} }