
We investigate the lowdimensional structure of deterministic transformations
between random variables, i.e., transport maps between probability measures. In
the context of statistics and machine learning, these transformations can be
used to couple a tractable "reference" measure (e.g., a standard Gaussian) with
a target measure of interest. Direct simulation from the desired measure can
then be achieved by pushing forward reference samples through the map. Yet
characterizing such a mape.g., representing and evaluating itgrows
challenging in high dimensions. The central contribution of this paper is to
establish a link between the Markov properties of the target measure and the
existence of lowdimensional couplings, induced by transport maps that are
sparse and/or decomposable. Our analysis not only facilitates the construction
of transformations in highdimensional settings, but also suggests new
inference methodologies for continuous nonGaussian graphical models. For
instance, in the context of nonlinear statespace models, we describe new
variational algorithms for filtering, smoothing, and sequential parameter
inference. These algorithms can be understood as the natural
generalizationto the nonGaussian caseof the squareroot
RauchTungStriebel Gaussian smoother.

The accurate approximation of highdimensional functions is an essential task
in uncertainty quantification and many other fields. We propose a new function
approximation scheme based on a spectral extension of the tensortrain (TT)
decomposition. We first define a functional version of the TT decomposition and
analyze its properties. We obtain results on the convergence of the
decomposition, revealing links between the regularity of the function, the
dimension of the input space, and the TT ranks. We also show that the
regularity of the target function is preserved by the univariate functions
(i.e., the "cores") comprising the functional TT decomposition. This result
motivates an approximation scheme employing polynomial approximations of the
cores. For functions with appropriate regularity, the resulting
\textit{spectral tensortrain decomposition} combines the favorable
dimensionscaling of the TT decomposition with the spectral convergence rate of
polynomial approximations, yielding efficient and accurate surrogates for
highdimensional functions. To construct these decompositions, we use the
sampling algorithm \texttt{TTDMRGcross} to obtain the TT decomposition of
tensors resulting from suitable discretizations of the target function. We
assess the performance of the method on a range of numerical examples: a
modifed set of Genz functions with dimension up to $100$, and functions with
mixed Fourier modes or with local features. We observe significant improvements
in performance over an anisotropic adaptive Smolyak approach. The method is
also used to approximate the solution of an elliptic PDE with random input
data. The open source software and examples presented in this work are
available online.

We present an arbitraryorder spectral element method for generalpurpose
simulation of nonoverturning water waves, described by fully nonlinear
potential theory. The method can be viewed as a highorder extension of the
classical finite element method proposed by Cai et al (1998)
\cite{CaiEtAl1998}, although the numerical implementation differs greatly.
Features of the proposed spectral element method include: nodal Lagrange basis
functions, a general quadraturefree approach and gradient recovery using
global $L^2$ projections. The quartic nonlinear terms present in the Zakharov
form of the free surface conditions can cause severe aliasing problems and
consequently numerical instability for marginally resolved or very steep waves.
We show how the scheme can be stabilised through a combination of
overintegration of the Galerkin projections and a mild spectral filtering on a
per element basis. This effectively removes any aliasing driven instabilities
while retaining the highorder accuracy of the numerical scheme. The additional
computational cost of the overintegration is found insignificant compared to
the cost of solving the Laplace problem. The model is applied to several
benchmark cases in two dimensions. The results confirm the high order accuracy
of the model (exponential convergence), and demonstrate the potential for
accuracy and speedup. The results of numerical experiments are in excellent
agreement with both analytical and experimental results for strongly nonlinear
and irregular dispersive wave propagation. The benefit of using a highorder 
possibly adapted  spatial discretization for accurate water wave propagation
over long times and distances is particularly attractive for marine
hydrodynamics applications.

A major challenge in nextgeneration industrial applications is to improve
numerical analysis by quantifying uncertainties in predictions. In this work we
present a formulation of a fully nonlinear and dispersive potential flow water
wave model with random inputs for the probabilistic description of the
evolution of waves. The model is analyzed using random sampling techniques and
nonintrusive methods based on generalized Polynomial Chaos (PC). These methods
allow to accurately and efficiently estimate the probability distribution of
the solution and require only the computation of the solution in different
points in the parameter space, allowing for the reuse of existing simulation
software. The choice of the applied methods is driven by the number of
uncertain input parameters and by the fact that finding the solution of the
considered model is computationally intensive. We revisit experimental
benchmarks often used for validation of deterministic water wave models. Based
on numerical experiments and assumed uncertainties in boundary data, our
analysis reveals that some of the known discrepancies from deterministic
simulation in comparison with experimental measurements could be partially
explained by the variability in the model input. We finally present a synthetic
experiment studying the variance based sensitivity of the wave load on an
offshore structure to a number of input uncertainties. In the numerical
examples presented the PC methods have exhibited fast convergence, suggesting
that the problem is amenable to being analyzed with such methods.