  The conditional extremes framework allows for event-based stochastic modeling
of dependent extremes, and has recently been extended to spatial and
spatio-temporal settings. After standardizing the marginal distributions and
applying an appropriate linear normalization, certain non-stationary Gaussian
processes can be used as asymptotically-motivated models for the process
conditioned on threshold exceedances at a fixed reference location and time. In
this work, we adapt existing conditional extremes models to allow for the
handling of large spatial datasets. This involves specifying the model for
spatial observations at $d$ locations in terms of a latent $m\ll d$ dimensional
Gaussian model, whose structure is specified by a Gaussian Markov random field.
We perform Bayesian inference for such models for datasets containing thousands
of observation locations using the integrated nested Laplace approximation, or
INLA. We explain how constraints on the spatial and spatio-temporal Gaussian
processes, arising from the conditioning mechanism, can be implemented through
the latent variable approach without losing the computationally convenient
Markov property. We discuss tools for the comparison of models via their
posterior distributions, and illustrate the flexibility of the approach with
gridded Red Sea surface temperature data at over $6,000$ observed locations.
Posterior sampling is exploited to study the probability distribution of
cluster functionals of spatial and spatio-temporal extreme episodes.
