Brown‒Resnick processes are max-stable processes that are associated to Gaussian processes. Their simulation is often based on the corresponding spectral representation which is not unique. We study to what extent simulation accuracy and efficiency can be improved by minimizing the maximal variance of the underlying Gaussian process. Such a minimization is a difficult mathematical problem that also depends on the geometry of the simulation domain. We extend Matheron's (1974) seminal contribution in two directions: (i) making his description of a minimal maximal variance explicit for convex variograms on symmetric domains, and (ii) proving that the same strategy also reduces the maximal variance for a huge class of nonconvex variograms representable through a Bernstein function. A simulation study confirms that our noncostly modification can lead to substantial improvements among Gaussian representations. We also compare it with three other established algorithms.

Extreme value theory for random vectors and stochastic processes with continuous trajectories is usually formulated for random objects where the univariate marginal distributions are identical. In the spirit of Sklar's theorem from copula theory, such marginal standardization is carried out by the pointwise probability integral transform. Certain situations, however, call for stochastic models whose trajectories are not continuous but merely upper semicontinuous (USC). Unfortunately, the pointwise application of the probability integral transform to a USC process does not, in general, preserve the upper semicontinuity of the trajectories. In this paper we give sufficient conditions to enable marginal standardization of USC processes and we state a partial extension of Sklar's theorem for USC processes. We specialize the results to max-stable processes whose marginal distributions and normalizing sequences are allowed to vary with the coordinate.

We propose a method to generate a max-stable process in C[0, 1] from a max-stable random vector in Rd by generalizing the max-linear model established by Wang and Stoev (2011). For this purpose, an interpolation technique that preserves max-stability is proposed. It turns out that if the random vector follows some finite-dimensional distribution of some initial max-stable process, the approximating processes converge uniformly to the original process and the pointwise mean-squared error can be represented in a closed form. The obtained results carry over to the case of generalized Pareto processes. The introduced method enables the reconstruction of the initial process only from a finite set of observation points and, thus, a reasonable prediction of max-stable processes in space becomes possible. A possible extension to arbitrary dimensions is outlined.

In this paper we provide the basis for new methods of inference for max-stable processes ξ on general spaces that admit a certain incremental representation, which, in important cases, has a much simpler structure than the max-stable process itself. A corresponding peaks-over-threshold approach will incorporate all single events that are extreme in some sense and will therefore rely on a substantially larger amount of data in comparison to estimation procedures based on block maxima. Conditioning a process η in the max-domain of attraction of ξ on being extremal, several convergence results for the increments of η are proved. In a similar way, the shape functions of mixed moving maxima (M3) processes can be extracted from suitably conditioned single events η. Connecting the two approaches, transformation formulae for processes that admit both an incremental and an M3 representation are identified.

Max-stable random fields play a central role in modeling extreme value phenomena. We obtain an explicit formula for the conditional probability in general max-linear models, which include a large class of max-stable random fields. As a consequence, we develop an algorithm for efficient and exact sampling from the conditional distributions. Our method provides a computational solution to the prediction problem for spectrally discrete max-stable random fields. This work offers new tools and a new perspective to many statistical inference problems for spatial extremes, arising, for example, in meteorology, geology, and environmental applications.

We develop classification results for max-stable processes, based on their spectral representations. The structure of max-linear isometries and minimal spectral representations play important roles. We propose a general classification strategy for measurable max-stable processes based on the notion of co-spectral functions. In particular, we discuss the spectrally continuous-discrete, the conservative-dissipative, and the positive-null decompositions. For stationary max-stable processes, the latter two decompositions arise from connections to nonsingular flows and are closely related to the classification of stationary sum-stable processes. The interplay between the introduced decompositions of max-stable processes is further explored. As an example, the Brown-Resnick stationary processes, driven by fractional Brownian motions, are shown to be dissipative.

In the characterization of multivariate extremal indices of multivariate stationary processes, multivariate maxima of moving maxima processes, or M4 processes for short, have been introduced by Smith and Weissman. Central to the introduction of M4 processes is that the extreme observations of multivariate stationary processes may be characterized in terms of a limiting max-stable process under quite general conditions, and that a max-stable process can be arbitrarily closely approximated by an M4 process. In this paper, we derive some additional basic probabilistic properties for a finite class of M4 processes, each of which contains finite-range clustered moving patterns, called signature patterns, when extreme events occur. We use these properties to construct statistical estimation schemes for model parameters.