Researcher profile

Philip Jonathan

Philip Jonathan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

A penalised piecewise-linear model for non-stationary extreme value analysis of peaks over threshold

Metocean extremes often vary systematically with covariates such as direction and season. In this work, we present non-stationary models for the size and rate of occurrence of peaks over threshold of metocean variables with respect to one- or two-dimensional covariates. The variation of model parameters with covariate is described using a piecewise-linear function in one or two dimensions defined with respect to pre-specified node locations on the covariate domain. Parameter roughness is regulated to provide optimal predictive performance, assessed using cross-validation, within a penalised likelihood framework for inference. Parameter uncertainty is quantified using bootstrap resampling. The models are used to estimate extremes of storm peak significant wave height with respect to direction and season for a site in the northern North Sea. A covariate representation based on a triangulation of the direction-season domain with six nodes gives good predictive performance. The penalised piecewise-linear framework provides a flexible representation of covariate effects at reasonable computational cost.

preprint2022arXiv

Enhanced monitoring of atmospheric methane from space over the Permian basin with hierarchical Bayesian inference

Methane is a strong greenhouse gas, with a higher radiative forcing per unit mass and shorter atmospheric lifetime than carbon dioxide. The remote sensing of methane in regions of industrial activity is a key step toward the accurate monitoring of emissions that drive climate change. Whilst the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinal-5P satellite is capable of providing daily global measurement of methane columns, data are often compromised by cloud cover. Here, we develop a statistical model which uses nitrogen dioxide concentration data from TROPOMI to efficiently predict values of methane columns, expanding the average daily spatial coverage of observations of the Permian basin from 16% to 88% in the year 2019. The addition of predicted methane abundances at locations where direct observations are not available will support inversion methods for estimating methane emission rates at shorter timescales than is currently possible.

preprint2022arXiv

Extremal Characteristics of Conditional Models

Conditionally specified models are often used to describe complex multivariate data. Such models assume implicit structures on the extremes. So far, no methodology exists for calculating extremal characteristics of conditional models since the copula and marginals are not expressed in closed forms. We consider bivariate conditional models that specify the distribution of $X$ and the distribution of $Y$ conditional on $X$. We provide tools to quantify implicit assumptions on the extremes of this class of models. In particular, these tools allow us to approximate the distribution of the tail of $Y$ and the coefficient of asymptotic independence $η$ in closed forms. We apply these methods to a widely used conditional model for wave height and wave period. Moreover, we introduce a new condition on the parameter space for the conditional extremes model of Heffernan and Tawn (2004), and prove that the conditional extremes model does not capture $η$, when $η<1$.

preprint2022arXiv

Multivariate spatial conditional extremes for extreme ocean environments

The joint extremal spatial dependence of wind speed and significant wave height in the North East Atlantic is quantified using Metop satellite scatterometer and hindcast observations for the period 2007-2018, and a multivariate spatial conditional extremes (MSCE) model, ultimately motivated by the work of Heffernan and Tawn (2004). The analysis involves (a) registering individual satellite swaths and corresponding hindcast data onto a template transect (running approximately north-east to south-west, between the British Isles and Iceland), (b) non-stationary directional-seasonal marginal extreme value analysis at a set of registration locations on the transect, (c) transformation from physical to standard Laplace scale using the fitted marginal model, (d) estimation of the MSCE model on the set of registration locations, and assessment of quality of model fit. A joint model is estimated for three spatial quantities: Metop wind speed, hindcast wind speed and hindcast significant wave height. Results suggest that, when conditioning on extreme Metop wind speed, extremal spatial dependence for all three quantities decays over approximately 600-800 km.

preprint2022arXiv

Statistical estimation of spatial wave extremes for tropical cyclones from small data samples: validation of the STM-E approach using long-term synthetic cyclone data for the Caribbean Sea

Occurrences of tropical cyclones at a location are rare, and for many locations, only short periods of observations or hindcasts are available. Hence, estimation of return values (corresponding to a period considerably longer than that for which data is available) for cyclone-induced significant wave height (SWH) from small samples is challenging. The STM-E (space-time maximum and exposure) model was developed to provide reduced bias in estimates of return values compared to competitor approaches in such situations, and realistic estimates of return value uncertainty. STM-E exploits data from a spatial neighbourhood satisfying certain conditions, rather than data from a single location, for return value estimation. This article provides critical assessment of the STM-E model for tropical cyclones in the Caribbean Sea near Guadeloupe for which a large database of synthetic cyclones is available, corresponding to more than 3,000 years of observation. Results indicate that STM-E yields values for the 500-year return value of SWH and its variability, estimated from 200 years of cyclone data, consistent with direct empirical estimates obtained by sampling 500 years of data from the full synthetic cyclone database; similar results were found for estimation of the 100-year return value from samples corresponding to approximately 50 years of data. In general, STM-E also provides reduced bias and more realistic uncertainty estimates for return values relative to single location analysis.

preprint2020arXiv

Modelling non-stationary extremes of storm severity: a tale of two approaches

Models for extreme values accommodating non-stationarity have been amply studied and evaluated from a parametric perspective. Whilst these models are flexible, in the sense that many parametrizations can be explored, they assume an asymptotic distribution as the proper fit to observations from the tail. This paper provides a holistic approach to the modelling of non-stationary extreme events by iterating between parametric and semi-parametric approaches, thus providing an automatic procedure to estimate a moving threshold with respect to a periodic covariate in circular data. By exploiting advantages and mitigating pitfalls of each approach, a unified framework is provided as the backbone for estimating extreme quantiles, including that of the $T$-year level and finite right endpoint, which seeks to optimize bias-variance trade-off. To this end, two tuning parameters related to the spread of peaks over threshold are introduced. We provide guidance for applying the methodology to the directional modelling of hindcast storm peak significant wave heights recorded in the North Sea. Although the theoretical underpinning for adaptation of well-known estimators in statistics of extremes to circular data is given in some detail, the derivation of their asymptotic properties lays beyond the scope of this paper. A bootstrap technique is implemented for obtaining direction-driven confidence bounds in such a way as to account for the relevant boundary restrictions with minimal sensitivity to initial point. This provides a template for other applications where the analysis of directional extremes is of importance.

preprint2018arXiv

Statistics of extreme ocean environments: Non-stationary inference for directionality and other covariate effects

Numerous approaches are proposed in the literature for non-stationarity marginal extreme value inference, including different model parameterisations with respect to covariate, and different inference schemes. The objective of this article is to compare some of these procedures critically. We generate sample realisations from generalised Pareto distributions, the parameters of which are smooth functions of a single smooth periodic covariate, specified to reflect the characteristics of actual samples from the tail of the distribution of significant wave height with direction, considered in the literature in the recent past. We estimate extreme values models (a) using Constant, Fourier, B-spline and Gaussian Process parameterisations for the functional forms of generalised Pareto shape and (adjusted) scale with respect to covariate and (b) maximum likelihood and Bayesian inference procedures. We evaluate the relative quality of inferences by estimating return value distributions for the response corresponding to a time period of $10 \times$ the (assumed) period of the original sample, and compare estimated return values distributions with the truth using Kullback-Leibler, Cramer-von Mises and Kolmogorov-Smirnov statistics. We find that Spline and Gaussian Process parameterisations estimated by Markov chain Monte Carlo inference using the mMALA algorithm, perform equally well in terms of quality of inference and computational efficiency, and generally perform better than alternatives in those respects.

preprint2013arXiv

Bayes linear variance structure learning for inspection of large scale physical systems

Modelling of inspection data for large scale physical systems is critical to assessment of their integrity. We present a general method for inference about system state and associated model variance structure from spatially distributed time series which are typically short, irregular, incomplete and not directly observable. Bayes linear analysis simplifies parameter estimation and avoids often-unrealistic distributional assumptions. Second-order exchangeability judgements facilitate variance learning for sparse inspection time-series. The model is applied to inspection data for minimum wall thickness from corroding pipe-work networks on a full-scale offshore platform, and shown to give materially different forecasts of remnant life compared to an equivalent model neglecting variance learning.

preprint2012arXiv

Evaluating environmental joint extremes for the offshore industry

Understanding extreme ocean environments and their interaction with fixed and floating structures is critical for the design of offshore and coastal facilities. The joint effect of various ocean variables on extreme responses of offshore structures is fundamental in determining the design loads. For example, it is known that mean values of wave periods tend to increase with increasing storm intensity, and a floating system responds in a complex way to both variables. However, specification of joint extremes in design criteria has often been somewhat \textit{ad hoc}, being based on fairly arbitrary combinations of extremes of variables estimated independently. Such approaches are even outlined in design guidelines. Mathematically more consistent estimates of the joint occurrence of extreme environmental variables fall into two camps in the offshore industry -- response-based and response-independent. Both are outlined here, with emphasis on response-independent methods, particularly those based on the conditional extremes model recently introduced by Heffernan and Tawn (2004) which has a solid theoretical motivation. Several applications using the new methods are presented.