On the number of excursion sets of planar Gaussian fields

The Nazarov–Sodin constant describes the average number of nodal set components of smooth Gaussian fields on large scales. We generalise this to a functional describing the corresponding number of level set components for arbitrary levels. Using results from Morse theory, we express this functional as an integral over the level densities of different types of critical points, and as a result deduce the absolute continuity of the functional as the level varies. We further give upper and lower bounds showing that the functional is at least bimodal for certain isotropic fields, including the important special case of the random plane wave.


The Nazarov-Sodin constant
Let f : R 2 → R be a continuous stationary planar Gaussian field normalised to have zero mean and unit variance. The nodal set of f is the random set We assume throughout that κ is C 4+ , which ensures that almost surely f is C 2+ . Since κ is positive definite, continuous and κ(0) = 1, by Bochner's theorem there exists a probability measure ρ such that this is known as the spectral measure of f , and must be Hermitian (that is, ρ(−A) = ρ(A) for all Borel sets A). Since the distribution of f is uniquely determined by its covariance function (Kolmogorov's theorem), (1) shows that the distribution of f is uniquely determined by ρ. The geometric properties of N are of interest, in part, because in the case that f is a random eigenfunction of the Laplacian they relate to a significant conjecture in the physics literature: the Berry conjecture [4]. A summary of this conjecture and other research on this topic may be found in [17]. One of the main analytical results concerning this set, due to Nazarov and Sodin [13,15], states that the number of components of N in a large domain scales like the area of the domain. Specifically, if N R denotes the number of components of N inside the centred ball of radius R > 0, then provided f is ergodic, there exists a constant c L S = c L S (ρ) ≥ 0 such that N R /(π R 2 ) → c L S as R → ∞, where convergence occurs almost surely and in L 1 . Nazarov-Sodin also obtained analogous results in higher dimensions and for Gaussian ensembles on manifolds [13]. In the case that f is not ergodic, it has been shown [10] (under the additional assumption that ρ has compact support) that the expected number of nodal components, scaled by the area, still converges, i.e.
as R → ∞. Further, in [10] it was also shown that among fields with compactly supported spectral measures, the constant c L S varies continuously with ρ (in the weak- * topology).

The main results
The first contribution of this paper is to extend the results of Nazarov-Sodin and [10] to arbitrary levels. For u ∈ R 2 and R > 0 let B(u, R) be the ball of radius R centred at u and B(R) := B(0, R). Let N = {x ∈ R 2 : f (x) = } denote a level set of f and let N L S,R ( ) be the number of components of N contained in B(R) (i.e. those which intersect B(R) but not ∂ B(R)). We consider fields satisfying the following assumptions: Conditions 1. 1 The Gaussian field f satisfies: 1. For some ν > 0, f ∈ C 2+ν loc (R 2 ) almost surely; 2. ∇ 2 f (0) is a non-degenerate Gaussian vector (here, and later on, we treat ∇ 2 f as a vector of distinct partial derivatives ( f x x , f xy , f yy )); (0)) is a non-degenerate Gaussian vector.
We note that these assumptions are quite minimal. The first condition holds if the covariance function κ is C 4+ for some > 2ν, or equivalently when The second condition is equivalent to the support of ρ not being contained in the union of two lines through the origin. The third condition holds provided that the support of ρ is not too degenerate; in particular it holds if the support of ρ contains an open set. as R → ∞. The constant implied by the O(·) notation may depend on ρ but is independent of . If f is also ergodic, then N E S,R ( )/(π R 2 ) → c E S (ρ, ) almost surely and in L 1 .

Remark 1.4 Our use of the domain B(R) in the definitions of N L S,R ( ) and N E S,R ( )
is mainly for simplicity, and after minor modifications the proofs of Theorems 1.2 and 1.3 go through equally well for rescaled copies of any bounded reference domain Ω ⊂ R 2 , provided that Ω is convex and ∂Ω is piecewise smooth (and probably more generally as well). Moreover the limiting constants c L S and c E S do not depend on Ω (after replacing the scaling factor π R 2 by Area(Ω)R 2 ). The symmetry of f along with the observations in the previous remark immediately give the following corollary. Theorems 1.2 and 1.3 are, in isolation, only a modest improvement on previous results; they could be proven by slightly adapting the analysis in [10,13]. The main novelty of our work is to relate the functionals c L S and c E S to the density of critical points of f of different types and at different levels. To state the relationship, we shall require, in particular, a classification of the saddle points of the field into two types: Definition 1.7 Let x 0 be a saddle point of a function g : R 2 → R such that there are no other critical points at the same level as x 0 . We say that x 0 is an upper connected saddle if it is in the closure of only one component of {x ∈ R 2 : g(x) > g(x 0 )}. Similarly, x 0 is said to be lower connected if it is in the closure of only one component of {x ∈ R 2 : g(x) < g(x 0 )}.
We say that a Gaussian field f satisfying Conditions 1.1 is periodic if there exists x = 0 with κ(x) = 1 and is aperiodic otherwise. A Gaussian field which is aperiodic almost surely has no two critical points at the same level (see Lemma 2.4), and we will show that for such fields, all saddle points of f satisfy exactly one of the two conditions in Definition 1.7. In the case that f is periodic, we will require a more general definition for classifying saddle points as upper or lower connected, which is given in Sect. 3.
Previous work has shown that, for sufficiently regular isotropic 1 Gaussian fields, the expected number of local maxima, local minima or saddle points with value in a certain interval can be expressed as the integral of an explicit density function [5,6].
In Sect. 3 we prove the following version of this result, which applies to more general Gaussian fields and also isolates the upper and lower connected saddles, but does not explicitly identify the densities. This proposition uses Definition 3.1 for upper and lower connected saddle points (which coincides with Definition 1.7 for aperiodic fields). The main theorem of the paper gives an explicit expression for the functionals c L S and c E S in terms of the densities introduced in Proposition 1.8. As a result, we deduce the absolute continuity of these functionals as the level varies. In the case of spectral measures with compact support, we also show the joint continuity of these functionals with respect to both the level and the spectral measure.
For each f satisfying Conditions 1.1 we can associate a spectral measure ρ via (1). Let P c denote the collection of such measures that are supported in the closure of B (1). By rescaling the axes, the results we state for P c can be shown to hold for any spectral measures with compact support. Theorem 1.9 Let f satisfy Conditions 1.1 and let p m + , p m − , p s + , p s − denote the densities specified in Proposition 1.8. Then and and hence c L S and c E S are absolutely continuous in . In addition c L S and c E S are jointly continuous in (ρ, ) ∈ P c × R where P c is given the weak- * topology. Remark 1.10 Theorem 1.9 provides a new tool with which to analyse the Nazarov-Sodin constant. Since the densities p m + , p m − and p s are in principle known by the Kac-Rice formula, our result demonstrates that the study of the Nazarov-Sodin constant can be reduced to an analysis of the density p s − (or, equivalently, p s + = p s − p s − ), which may be an easier quantity to handle.
The number of excursion set components increases by one on passing through the lower connected saddle x 1 and is constant on passing through the upper connected saddle x 2 Remark 1.11 Theorems 1.2, 1.3 and 1.9 can be generalised to many examples of non-Gaussian stationary random fields, since our proof requires only that the field satisfies certain topological properties almost surely and is sufficiently regular to apply the Kac-Rice formula (see the more general version of our main results stated in Propositions 2.6 and 2.7 below). We also believe Theorems 1.2, 1.3 and 1.9 could be generalised to higher dimensions, with the analogues of (2) and (3) still valid once the saddle points defining p s + and p s − are replaced with critical points of index 1 and d − 1 respectively (using a more general definition for upper and lower connected saddles). However, since some of the topological arguments in the proof increase in complexity in higher dimensions, in the interest of simplicity we do not pursue this generalisation here.
Let us mention the key intuition behind the proof of Theorem 1.9. This theorem is based on a deterministic relationship between the excursion sets and critical points of sufficiently regular planar functions which is closely related to Morse theory. The excursion set of such a function above a level deforms continuously as the level increases, provided it does not pass through a critical point. In particular, there is no change in the number of components of the excursion set. When passing through a critical point, the topology of the excursion set changes in a way that is predicted by Morse theory and depends on the index of the critical point. For local maxima and local minima, the number of components of the excursion set changes in a consistent way. For saddle points, the change in the number of components is determined by whether it is upper connected or lower connected (see Fig. 1). Ultimately, Theorem 1.9 is a probabilisitic expression, in the setting of Gaussian fields, of this deterministic relationship between excursion set components and critical points of various types.
We also briefly discuss the assumptions required for our results (i.e. Conditions 1.1).
The assumption that f is almost surely C 2+ν is necessary to apply the topological arguments that we borrow from Morse theory. The non-degeneracy assumptions on ) are used with the Kac-Rice theorem (see Sect. 2) to show certain non-degeneracy properties of f , which are again necessary for our topological arguments.
It is natural, therefore, to ask whether our results still apply when these assumptions are weakened. While we suspect that this is true, we are not able to show it with our methods. On the other hand, for a certain special case of non-trivial degenerate field -namely, the case where the spectral measure ρ is supported on at most five points -we are able to give a complete description of c L S and c E S , which in particular shows that the main results hold also in this case (see Sect. 1.4 below).

Bounds on c LS and c ES
It is possible to bound the expected number of excursion sets or level sets of stationary Gaussian fields using local estimates (i.e. estimates which depend only on the derivatives of κ at the origin). Here we outline how these estimates apply to the functionals c E S and c L S and how they can be better characterised by making use of our results.

Corollary 1.12 Let f be a Gaussian field satisfying Conditions
where φ denotes the standard normal probability density function.
and not on the covariance of higher order derivatives of f .

Remark 1.
14 This result was proven for isotropic fields in [16] using a winding number calculation. We prove this result using Theorem 1.9 which simplifies the calculation in the non-isotropic case and also highlights an interesting identity. Specifically, the proof is as follows: substituting (3) into the left hand side of (4) and using the symmetries and p s − + p s + = p s , we see that where N m + ( ), N s ( ) and N m − ( ) are the number of critical points above level as defined in Proposition 1.8 for Ω = B(1/ √ π). Lemma 11.7.1 of [1] states that this alternating sum is precisely the right hand side of (4). The alternating sum of critical points of a function of different indices above a certain level can be used to calculate the Euler characteristic of the excursion set of the function above the level (see Chapter 9 of [1]). When working on finite subsets of the plane, boundary effects must be considered, but these become negligible as the area of the subset increases. Formally, if we let ϕ(A) denote the Euler characteristic of a set A, then Theorem 11.7.2 of [1] gives an expression for the expected Euler characteristic of an excursion set on a cube (including boundary effects) which implies that We now consider an upper bound on c L S . This is most easily formulated when f is an isotropic Gaussian field. In this case its covariance function may be expressed as κ(x) = K (|x|). It is shown in [5] that parameterise the critical point densities p m + , p m − and p s of isotropic fields in any dimension. It can also be shown that, for planar isotropic fields, if λ = √ 2 then κ(x) = J 0 ( 8/η 2 |x|), where J 0 is the 0-th Bessel function. When η 2 = 8, this particular field is known as the random plane wave (hereafter abbreviated to RPW), and is an object of great interest as it is the subject of the Berry conjecture [17].

Proposition 1.15 ([16]) Let f be an isotropic Gaussian field satisfying Conditions
where φ and Φ denote the standard normal probability density and cumulative density functions respectively. In particular, if λ 2 > 6e 2e+π ≈ 1.9 then by Corollary 1.12, c L S (0) < c L S (1) so that c L S is at least bimodal (that is, it has at least two local maxima).

Remark 1.16
This result is proven in [16] using the method of 'flip points'. Specifically, for any fixed direction u, the number of level set components at level in a finite region is bounded above by half the number of points t in the region such that f (t) = and ∂ u f (t) = 0 where ∂ u denotes the partial derivative in the direction u. The expected number of such points can be computed using the Kac-Rice formula. For isotropic fields, the choice of direction u is irrelevant. This method could also be applied to non-isotropic fields, and the bound could be optimised over the direction u, but we omit this for simplicity.
It is also possible to construct an upper bound on c L S using the inequality 0 ≤ p s − ≤ p s and the densities p m + , p s and p m − , which are explicitly known for isotropic fields. However this bound is larger than the bound in (6) at all levels . Figure 2 shows lower and upper bounds for c E S and c L S for the RPW based on (4), (6) and the equality c L S ( ) = c E S ( ) + c E S (− ) (recall that λ = √ 2 and η 2 = 8 in this case). Although these bounds are not particularly tight for close to 0, they quickly become accurate as | | increases. For ≥ 1 both the upper and lower bounds on c L S ( ) are within 5.1% of the true value (since the upper bound is within 5.1% of the lower bound) while for ≥ 1.5 both bounds are within 0.6% of the true value. For fields with lower values of λ, the percentage difference between the upper and lower bounds on c L S ( ) is a bit larger. Figure 3 shows the corresponding bounds on c E S and c L S for the Bargmann-Fock field: the centred, planar Gaussian field with covariance kernel κ(x) = exp(−|x| 2 /2), for which λ = 1 and η 2 = 2. In this case for ≥ 1 the upper and lower bounds on c L S ( ) are within 40% of the true value, and the corresponding accuracies for ≥ 1.5, ≥ 2 and ≥ 2.5 are 14%, 5% and 2% respectively. Remark 1.18 Corollary 1.12 states that c E S (ρ, ) > 0 for all > 0, and so c L S (ρ, ) > 0 for all = 0, however the positivity of c L S (ρ, 0) is not addressed by this method. In [13] Nazarov and Sodin provide sufficient conditions for the positivity of c L S (ρ, 0) in terms of the spectral measure which cover almost all non-degenerate cases.
Although the bounds in the previous two corollaries can be proven by local methods (i.e. without using our analysis of saddle points) there is no scope for improving these estimates using the same methods. However it may be possible to get tighter bounds on c L S and c E S through Theorem 1.9 and a better characterisation of p s − .
Using the explicitly known densities p m + , p m − and p s along with Theorem 1.9 allows us to derive some monotonicity properties of c E S and c L S .

A special class of degenerate fields
As mentioned above, our methods do not cover the case in which ρ is supported on the union of two lines through the origin. In this section we consider a certain special class of such fields for which we can give a more or less complete description of c L S and c E S .
If the spectral measure ρ is supported on a single line through the origin, it is well known that f is almost surely constant in one direction and so in particular has no compact level domains. In this case, all of our results hold trivially with c L S , c E S and the critical point densities identically equal to zero. The simplest non-trivial case of degenerate ρ are spectral measures which are supported on four or five points. In this case, we can compute c E S and c L S explicitly.
Recall that a random variable Y is said to be Rayleigh distributed with parameter

Proposition 1.20 Let f be the Gaussian field with spectral measure
where β, γ > 0, α = 1 − β − γ ≥ 0 and K , L ∈ R 2 are linearly independent. Then Furthermore there exist functions p m + , p s − , p m − , p s + and p s satisfying the conclusions of Proposition 1.8, and these are defined by where p Z denotes the probability density of a random variable Z . Therefore the equalities in the conclusion of Theorem 1.9 hold for f . Remark 1.21 A stationary centred continuous Gaussian field is ergodic if and only if its spectral measure has no atoms [12, Section 6.1], and so the fields considered in this proposition are not ergodic. The fact that N L S,R ( )/(π R 2 ) does not generally converge to c L S ( ) for these fields shows that the ergodicity requirement in the second part of Theorems 1.2 and 1.3 cannot be entirely relaxed.

Remark 1.22
It has previously been shown that if spt(ρ) has precisely four points then c L S (0) = 0 (see [10]). Proposition 1.20 shows that c L S (0) > 0 when spt(ρ) has five points, and that c L S ( ) > 0 for = 0 when spt(ρ) has four or five points.  Figure 4 suggests that c L S is bimodal whenever α = 0 (we know that it is at least bimodal, since c L S (0) = 0 and c L S ( ) > 0 for = 0 in this case), however for α > 0, c L S can be bimodal or unimodal. Combining this observation with Corollary 1.12 raises the interesting question of determining for which Gaussian fields c L S ( ) is bimodal. and these densities are explicitly known. In particular, it can be shown that p Y 1 +Y 2 is everywhere continuous, whereas p |Y 1 −Y 2 | has a jump discontinuity at zero (but is continuous elsewhere). So c E S is continuous except at zero. An analogous calculation permits the same conclusion for c L S . On the other hand, when α = 0, c E S and c L S are continuously differentiable everywhere. A major unresolved question arising from our work is to determine whether c L S and c E S are continuously differentiable more generally (outside degenerate cases).

Summary of the rest of the paper
The remainder of the paper is organised as follows. In Sect. 2 we state the main (completely deterministic) topological lemma underlying our results (Lemma 2.5) and combine this with Proposition 1.8 to prove Theorems 1.2, 1.3 and 1.9 and their

Proof of the main results
The first step in proving Theorems 1.2, 1.3 and 1.9 is to translate our assumptions on the Gaussian field f into topological properties; these properties differ slightly depending on whether or not f is periodic.
such that for all x ∈ R 2 and a ∈ R, g(x) = g(x +av). We say that g is doubly periodic if there exist two linearly independent vectors y, z ∈ R 2 such that for all x ∈ R 2 , g(x + y) = g(x) = g(x + z). We say that g is singly periodic if it is not doubly periodic and there exists y ∈ R 2 \{0} such that for all x ∈ R 2 , g(x) = g(x + y). If neither of these conditions holds we say that g is aperiodic.
We note that this definition applies to deterministic functions. We say that a random field f is doubly periodic if there exist linearly independent vectors y, z ∈ R 2 such that with probability one, We make an analogous definition for singly periodic fields, and we say that a random field is aperiodic if it is neither doubly periodic nor singly periodic. The Gaussian field with spectral measure ρ is doubly periodic if and only if there exists A ∈ G L(2, R) such

and similarly is singly periodic if and only if it is not doubly periodic and there exists
For our purposes, it is more natural to consider periodic functions as being defined on the torus or cylinder, and so we will specify the topological properties of Gaussian fields in terms of these domains.
For any doubly periodic function g : R 2 → R we can choose two linearly independent vectors y, z ∈ R 2 satisfying the conditions in Definition 2.1 with minimum distance to the origin, and g is then specified entirely by its values on the parallelogram We call this the associated parallelogram for g and call y, z periodic vectors for g (they are not unique since we could also choose −y, −z). By identifying P with the 2-dimensional torus T, we let g T be g defined on the torus. By the choice of y, z above, when g is a stationary random field there are no two distinct points in P on which g is almost surely equal.
If g is singly periodic then we can choose y ∈ R 2 \{0} satisfying the conditions in Definition 2.1 with minimum distance to the origin; we call this a periodic vector for g. By rotating the axes we may assume that g is periodic in the direction of the x-axis (i.e. y = (y 1 , 0)) so that g is specified entirely by its values on [0, y 1 ] × R. This strip can be identified with the infinite cylinder C : We will often work with the compact subset C(n) := [0, y 1 ] × [−n, n] under the same quotient relationship. We let g C and g C(n) denote g restricted to these surfaces under the previous identification.
We next introduce some elements of Morse theory, in particular the concept of a Morse function on a manifold; here we follow [8]. We emphasise that, unless explicitly stated, a manifold M may or may not have a boundary ∂ M. We also emphasise that we will only ever work with M being the plane R 2 , the closure of the disc B(R), the torus T, or the cylinders C and C(n), and so it is sufficient to bear these examples in mind. 3. If ∇ f is tangent to ∂ M at p and u is a unit vector normal to ∂ M at p, then and ∂ u denote the directional derivatives in the directions ∇ f ( p) and u respectively; 4. The critical points of f and f | ∂ M all occur at distinct levels.
If M is a manifold without boundary this simplifies to the requirement that all critical points of f are non-degenerate and occur at distinct levels.
We note that in the above definition f | ∂ M is viewed as a function on the (n − 1)dimensional manifold ∂ M so that a critical point of f | ∂ M is not necessarily a critical point of f .
Using the definition of a Morse function, we next introduce a set of assumptions that a stationary random field must satisfy in order for our main results to hold. In the lemma that immediately follows, we will claim in particular that these assumptions are satisfied for Gaussian fields satisfying Conditions 1.1, but we have chosen to isolate these assumptions to illustrate the limited role of Gaussianity in the proof of the main results.
We say that a point t ∈ ∂ M is a tangent point of f if t is a critical point of the restricted function f | ∂ M (the name comes from the fact that the level set { f = f (t)} will be tangent to ∂ M at t). We will also use the term 'tangent point' to describe a local extremum of the restriction of f to a finite union of line segments. Let N tang (A) denote the number of tangent points of f in A (where A is the boundary of a manifold or a finite union of line segments) and N crit (M) the number of critical points of f in M. We will only use these terms when A and M are compact, so that by the first point in the definition of a Morse function, the number of critical and tangent points will be finite.

Conditions 2.3 A stationary random field f satisfies the following:
1. If f is aperiodic, then for each R > 0, f | B(R) is almost surely Morse (taking a union over R ∈ N, this implies also that f is almost surely Morse on R 2 );

If f is doubly periodic, then f T is almost surely Morse;
3. If f is singly periodic, then for each n ∈ N, f C(n) is almost surely Morse (taking a union over N, this implies that f C is almost surely Morse); 4. If f is periodic and L ⊂ R 2 is a line segment of length R, then E(N tang (L)) = c θ R for a constant c θ > 0 that depends only on the direction of L; Proof All six conditions are established using variations of the Kac-Rice formula. This is an important tool in the study of Gaussian fields, and random functions more generally, which relates moments of the number of zeroes of a field to conditional moments of its gradient/derivative. See [14] for an introduction to the Kac-Rice formula or [1, Chapter 11] for a rigorous derivation.
(1)-(3). We first note that since f ∈ C 2 loc (R 2 ) almost surely by assumption, f | M ∈ C 2 loc (M) almost surely whenever M is B(R), T or C(n). Corollary 11.3.5 of [1] states that f | M has the first two properties of a Morse function if there exists a countable atlas 0)) are non-degenerate (where ∂ x denotes taking the derivative with respect to the first variable, ∂ x x denotes taking the second derivative with respect to the first variable and we make analogous definitions for ∂ y , ∂ yy and ∂ xy ). For each choice of M that we consider (i.e. B(R), T or C(n)) it is clear that we can choose a finite atlas of charts x i and that f | M • x −1 i will be a translation of f (restricted to some open set) for each such chart. Therefore the covariance kernel of f (i) is a restriction of κ, which is C 4+ by assumption, and all that remains is to show Gaussian vectors for any t ∈ R 2 . By assumption ∇ f (t) and ∇ 2 f (t) are non-degenerate Gaussian vectors for any t ∈ R 2 . It is a standard fact, for Gaussian fields with constant variance, that ∇ f (t) and ∇ 2 f (t) are independent for fixed t (see e.g. [1,Chapter 5]) and this proves the necessary non-degeneracy. (This part of the result does not depend on whether f is periodic. ) We note that f | B(R) has the third property of a Morse function provided that there is no θ ∈ [0, 2π ] such that where p = (R cos(θ ), R sin(θ )) (see Fig. 6).
If we assumed that f is C 3 , then Bulinskaya's lemma would imply that no such θ exists almost surely. Since we assume only that f is C 2+ν , we use a different argument which is adapted from the proof of [1,Lemma 11.2.11]. Since f ∈ C 2+ν loc (R 2 ) almost surely, g 2 is almost surely ν-Hölder continuous. Therefore for each δ > 0 we can find C δ , D δ > 0 such that, with probability at least 1 − δ, on ∂ B(R) the ν-Hölder norm of g 2 is at most C δ and each element of ∇ 2 f is bounded in absolute value by D δ ; we denote this event by A δ . Now let I m j be a collection of open intervals in R such that [0, 2π ] ⊂ ∪ j I m j for each m, j Length(I m j ) → 2π as m → ∞, and Length(I m j ) → 0 uniformly in j as m → ∞. Furthermore let θ m j be the midpoint of each I m j and define On the event A δ ∩ E m j we note that where p g 1 (θ m j ) is the probability density function of g 1 (θ m j ) and we have used the fact that g 1 (θ m j ) is independent of g 2 (θ m j ). Each g i (θ ) is a linear combination of elements of a non-degenerate Gaussian vector and so has non-zero variance, therefore on the compact region [0, 2π ] these variances are bounded away from zero. Since these are Gaussian fields, their univariate densities are bounded above by a constant times the inverse of their standard deviations, so we see that the density of each g i (θ m j ) is bounded above uniformly in m and j. We can therefore sum the expression on either side of (7) over j. Since Length(I m j ) converges to zero uniformly in j, so does the integrand in each term of the sum, and so we conclude that as m → ∞. Since this is true for any δ > 0 and P(A δ ) ≥ 1 − δ we conclude that there is almost surely no θ ∈ [0, 2π ] such that g 1 (θ ) = g 2 (θ ) = 0, and hence f almost surely has the third property of a Morse function.
The proof that f C(n) satisfies the third property of a Morse function is near identical. The argument above can be repeated for the boundary S 1 ×{−n, n}. Finally, f T trivially satisfies the third property of a Morse function.
Next we consider the fourth property of a Morse function. We first suppose f is aperiodic and define h 1 : Bulinskaya's lemma (Lemma 11.2.10 of [1]) states that P(h −1 1 (0) = ∅) = 1 provided that h 1 ∈ C 1 (T k ) almost surely and the univariate densities of h 1 are bounded uniformly in T k . By assumption h 1 ∈ C 1 loc (R 2 ) and so, since h 1 is Gaussian, we need only show that the determinant of the covariance matrix of Gaussian variable for all s = t, and so by Conditions 1.1 (in particular using stationarity of f ), the determinant of this covariance matrix is non-zero for all s = t. Since this determinant is continuous in (s, t) it is bounded away from zero on the compact set T k . Taking the countable union of T k for k ∈ N shows that almost surely there are no points To show that f has no two tangent points at the same level and no tangent points at the level of any critical point, we apply identical arguments to the following two functions: This completes the proof that f | B(R) is a Morse function when f is aperiodic. When f is doubly periodic, we can repeat these arguments with is the periodic vector of f . In both of these cases, the choice of domain ensures that f (t) − f (s) is non-degenerate, so that Conditions 1.1 can be used. This completes the proof of the first three points of the lemma.
(4). Let L be a line segment of length R and u a unit length vector parallel to L. Applying a Kac-Rice formula (for the particular form, see Corollary 11.2.2 of [1]) to f | L and using the stationarity of f along with the independence of ∇ f (t) and ∇ 2 f (t), we see that can be expressed as a linear combination of components of the non-degenerate Gaussian vector ∇ 2 f (0) and so has finite moments of all orders.
(5). Applying the same Kac-Rice formula to f , and using the same arguments as above, shows that and ∇ 2 f ( p) are independent, a simple calculation shows that (g (θ ), g (θ )) has non-degenerate distribution, so we can apply a Kac-Rice formula (specifically, Corollary 11.2.2 of [1]) to show that Since f is stationary The expectation on the right side of this equation is non-zero for each θ , since ∇ f (0) is non-degenerate, and by compactness is bounded away from zero uniformly in θ . Therefore, by considering the density of a Gaussian random variable, there exists Using this observation, along with the fact that ∇ f (0) is independent of ∇ 2 f (0), and substituting (8) into (9), we see that there exists a constant c 1 > 0 independent of θ and R such that We now introduce the deterministic relationship between level sets and critical points that is the foundation of our results. We denote the number of critical points of denotes local maxima, local minima, upper connected saddles and lower connected saddles respectively (for an aperiodic function that is Morse on R 2 , Definition 1.7 defines the upper/lower connected saddle points; the definitions in the general case will be given in Sect. 3). Let us also generalise our earlier notation: for a deterministic, planar function f , we denote the number of components of { f = } and { f ≥ } in B(R) by N L S,R ( ) and N E S,R ( ) respectively.
In the case of a singly periodic function f : R 2 → R we also need to define internal and external rectangular tilings of B(x, R) which have a finite number of horizontal and vertical line segments as boundary. Let (y 1 , 0) be the periodic vector of f , and let S(n 1 , n 2 ) = [n 1 y 1 , (n 1 +1)y 1 ]×[n 2 , n 2 +1] where n 1 , n 2 ∈ Z. We define B int (x, R) to be the union over n 1 , n 2 ∈ Z of all S(n 1 , n 2 ) contained in B(x, R). Similarly we define B ext (x, R) to be the union over n 1 , n 2 ∈ Z of all S(n 1 , n 2 ) which intersect B(x, R). Lemma 2.5 Let f : R 2 → R be a deterministic function. Suppose first that f is aperiodic and assume that f and f | B(R) are Morse. Then there exists c > 0, independent of f , such that for each ∈ R, and where max{|η ,R |, |γ ,R |} ≤ cN tang (∂ B(R)).
Suppose instead that f : R 2 → R is doubly periodic and assume that f T is Morse and f | ∂ P has a finite number of local extrema. Then there exists a constant c f > 0 depending only on P, such that, for each R > 0 and ∈ R, (10) and (11) hold with (12) replaced by Finally, suppose that f : R 2 → R is singly periodic with periodic vector (y, 0) and assume that f C(n) is Morse for each n ∈ Z and that N tang ({0} × [−2R, 2R]) < ∞. Then (10) and (11) still hold with (12) replaced by where c f and r f are constants that depend only on the periodic vector y = (y 1 , 0) of f (so in particular, they are independent of and R). We are now ready to state the core technical results of the paper, which are versions of Theorems 1.2, 1.3 and 1.9 that hold for arbitrary random fields (i.e. not necessarily Gaussian) using only the properties contained in Conditions 2.3.

Proposition 2.6 Let f be a stationary random field satisfying Conditions 2.3. For each
The constant implied by the O(·) notation may depend on the distribution of f but is independent of . If in addition f is ergodic, then

Proposition 2.7 Let f be a stationary random field satisfying Conditions
The constant implied by the O(·) notation may depend on the distribution of f but is independent of . If in addition f is ergodic, then

Remark 2.8
Since f is stationary, we could equivalently define c L S ( ) in Proposition 2.6 as for any R > 0 or make an analogous definition for c E S ( ) in Proposition 2.7.
Together with Lemma 2.4 and Proposition 1.8, these propositions imply Theorems 1.2 and 1.3, and also the first two parts of Theorems 1.9.
Proof (Proposition 2.6) As f is stationary, the expected number of critical points of f of a particular type in a domain is proportional to the area of the domain. Since the field satisfies Conditions 2.3 we can apply Lemma 2.5, and taking expectations yields Sending R → ∞ proves the first two statements of the proposition.
The remainder of the proof follows the general roadmap of the original derivation of the existence of the Nazarov-Sodin constant in [13]. Suppose that f is ergodic and let N (u) h,R denote the number of critical points of type h in B(u, R) with level in [ , ∞) or ( , ∞) for h = m + , s + and h = m − , s − respectively. The 'sandwich estimate' of [13] (Lemma 1) can be slightly altered to show that for any r ∈ (0, R) Applying Wiener's ergodic theorem (see [13, Section 6.1]) both of these integrals converge almost surely and in L 1 to the limit E N (0) h,1 /π . So in particular, N (0) h,R /(π R 2 ) has the same limit, and Applying Lemma 2.5 with the bound on E(|η ,R |) implied by Conditions 2.3 shows that as R → ∞. Combining these results completes the proof of L 1 convergence. We now extend our notation by defining N (u) L S,r to be the number of components of B(u, r ). The original 'sandwich estimate' of [13] states that we can rearrange this estimate as for some universal constant c 0 > 0. Applying Lemma 2.5 inside the integral term on the left hand side of (13) shows that (The upper bound here is a result of adding the upper bounds on the error terms in Lemma 2.5 in the cases that f is aperiodic, doubly periodic and singly periodic respectively, so that we can deal with all three cases at once.) From Wiener's ergodic theorem, the integral term within the absolute value signs will converge almost surely to c L S ( ) and the first integral term on the right hand side will converge almost surely to zero. Applying the same argument to the remaining integral terms shows that they will each converge to a constant, Φ 1 (r ), Φ 2 (r ), Φ 3 (r ) and Φ 4 (r ) respectively, where (0, r )) + E N tang (B ext (0, r )) , We now fix r > 0 and take the limit superior of (14) as R → ∞ to show that lim sup almost surely. Since the number of boundary points of f at level is deterministically bounded above by a constant times the number of tangent points, we see that Therefore we may take a countable sequence r n → ∞, such that (15) holds almost surely for each r n and the right hand side of (15) becomes arbitrarily small, to show that N  B(u, r ) and noting that the difference in these two terms is also bounded above by the number of tangent points of ∂ B(u, r ).
To complete the proof of Theorem 1.9, it remains only to show the joint continuity of c L S and c E S with respect to the level and spectral measure.
Proof (Theorem 1.9 -Joint continuity) By Prokhorov's theorem the weak- * topology on P c is metrisable and so we can use the sequential definition of continuity for c L S and c E S . Let (ρ n , n ) be a sequence in P c × R converging to (ρ, ). By the triangle inequality Theorem 1.3 of [10] states that the second term on the right hand side converges to 0 with n in the special case = 0. However the proof of this theorem can be repeated verbatim replacing the field f ρ with f ρ − to show that |c L S (ρ n , )− c L S (ρ, )| → 0. Now assume that n ≤ . By the first part of Theorem 1.9 we see that where N (n) crit [ n , ] denotes the number of critical points of f n (the Gaussian field with spectral measure ρ n ) in the circle of radius 1 with level in [ n , ]. By the Kac-Rice theorem (Corollary 11.2.2 of [1]) where we have used the independence of a field and its gradient at a point along with the Cauchy-Schwarz inequality. It is well known that the covariance structure of a stationary Gaussian field along with its derivatives at a point can be expressed in terms of the spectral measure (see [1,Chapter 5]). For example, with similar expressions for other derivatives and covariances. Since all spectral measures we are considering are supported on B(1), the definition of weak- * convergence implies that the covariance structure associated with each ρ n at the origin converges to that of ρ. Specifically, if f denotes the field with spectral measure ρ, then Applying this to (18), we see that det Var(∇ f n (0)) is bounded away from 0 in n. Similarly E |det ∇ 2 f n (0)| 2 is uniformly bounded above in n. Finally we note that Var( f n (0)) is uniformly bounded away from 0, so that choosing n sufficiently close to ensures that P ( f n (0) ∈ [ n , ]) is arbitrarily small. An identical argument works for n ≥ , and so combining these observations with (16) Proof (Corollary 1.19) Since f is isotropic we can use the explicitly derived critical point densities for local maxima, local minima and saddle points from [5]. We recall that these are parametrised in terms of λ ∈ (0, √ 2] and η 2 ∈ [0, ∞) as defined by (5).
where φ and Φ denote the standard normal probability density function and cumulative density function respectively. If λ = √ 2 then By Theorem 1.9, We denote the integrand above by I (x). In the case λ ∈ (0, √ 2), evaluating this expression using the standard Gaussian inequality 1 which is less than our upper bound for c L S ( + )− c L S ( ), and the first result follows. When λ = √ 2 and ≤ 0 we see that p m + ( ) = 0 and so by Theorem 1.9 Therefore c E S is weakly increasing in this case. When f is an aperiodic Gaussian field, by Lemma 2.4 it will be Morse on R 2 so we may take this as our Riemannian manifold and use the definition of lower/upper connected saddle points directly (we show in Sect. 4 that this definition coincides with that given in Sect. 1). If f is doubly periodic, recall that it is completely specified by its values on the associated parallelogram P which we identify with the torus. We then say that a saddle point of f is lower connected if it corresponds to a lower connected saddle of f T by the definition above (see Fig. 8). Similarly if f is singly periodic we say that a saddle point of f is lower connected if it corresponds to a lower connected saddle point of f C . Upper connected saddle points are defined analogously.
The following lemma shows that our above definitions partition the set of saddle points in all cases of interest. This result is proven in Sect. 4. We are now ready to prove Proposition 1.8.
Proof (Proposition 1.8) Let f satisfy Conditions 1.1. We fix a compact Ω ⊂ R 2 such that ∂Ω has finite Hausdorff-1 measure and consider N h ( ), the number of critical points of f | Ω of type h with value greater than , where h = m − , s, m + correspond to local minima, saddle points and local maxima respectively. If ( f (0), ∇ 2 f (0)) has a non-degenerate distribution, then by a Kac-Rice formula (specifically, Corollary 11.2.2 of [1], which requires that ∂Ω has finite Hausdorff-1 measure) where i is the index corresponding to h and φ denotes the standard Gaussian probability density. The third equality follows from the definition of conditioning on a Gaussian variable. If ( f (0), ∇ 2 f (0)) has a degenerate distribution then f (0) can be expressed as a linear combination of the elements of ∇ 2 f (0) almost surely. Substituting in this expression for f (0) allows us to apply Kac-Rice and the arguments above to derive (19) in this case too. Then, by Gaussian regression, is a continuous function of u (see, e.g., Proposition 1.2 of [2]). This proves the existence and continuity of the densities p m + , p m − and p s , and the fact that p m + (x) = p m − (−x) follows from the symmetry of f .
Since E (N s [0, , ]) = 0 p s (x)dx, we know that E(N s [0, ]) is absolutely continuous in . It is also clear that for any 1 is absolutely continuous in . Therefore there exists a function p s + : R → [0, ∞) such that Since E (N s [0, ∞)) is finite, the monotone convergence theorem shows that p s + ∈ L 1 (R). By symmetry of the Gaussian distribution, and the definition of lower connected saddles, this also shows the existence of p s − (x) = p s + (−x). The fact that p s − + p s + = p s follows from Lemma 3.2.

Topological lemmas
In this section we prove the deterministic Lemma 2.5 using topological arguments. We also establish Lemma 3.2 and Proposition 1.20 using similar methods, thereby completing the proof of all results in the paper.
To prove Lemmas 2.5 and 3.2 we require several results from Morse theory which we now introduce. We note that several aspects of this theory require us to work with compact manifolds, and so in the aperiodic and singly periodic cases we work with B(R) and C(n) rather than directly with R 2 and C; in the doubly periodic case, by contrast, we work directly with the torus T. We also emphasise that we only ever work with the five examples of manifolds just mentioned, so it is sufficient to have them in mind.
Recall that, for a Morse function f defined on a manifold M, the tangent points of f are defined as the critical points of f | ∂ M . We define a k-cell to be a copy of the closed unit disc in R k and temporarily denote this by B k . If Y is a topological space then we define the following operation to be 'attaching a k-cell to Y '. First we find a continuous function g : ∂ B k → Y , then we take the disjoint union Y B k and identify each point in ∂ B k with its image under g. By attaching a 0-cell, we simply mean taking the disjoint union of Y and a single point. These theorems are proven using methods very similar to those of the standard proofs for manifolds without boundary which can be found in almost any text on Morse theory.
We now work towards a proof of Lemma 2.5. First we will need another definition which, in the case of aperiodic or singly periodic fields, identifies a subset of the upper/lower connected saddle points which have unfavourable topological properties (for our purposes). Recall that if M is a manifold then A ⊂ M is said to be simple if A is compact, connected and every loop in A is M-contractible. In the case that M = R 2 or B(R) and A is an excursion set component, this is just the condition that A is bounded. We say that x 0 is an infinite-four-arm saddle point if A + has two unbounded components for all > 0 sufficiently small.

Remark 4.4
In the proof of Lemma 3.2, we will show that when M ∞ = R 2 (the most important case for our analysis), the level set at the level of an infinite-four-arm saddle takes the form shown in Fig. 9b. This corresponds to the way an 'infinite-four-arm event' is typically defined in the percolation literature: as two disjoint paths joining a point to infinity which are separated by two 'dual' paths joining the same point to infinity. For other choices of M ∞ (such as M ∞ = C) the level set at the level of an infinite-four-arm saddle may look quite different, so that this terminology is less intuitive.
Let us explain the importance of four-arm saddles. If f is aperiodic, then its lower connected saddle points are defined so as to correspond to an increase in the number of excursion set components in R 2 . However if such a saddle point is four-arm in B(R), then this increase cannot be observed from inside B(R) since the excursion sets that are created when passing through the saddle intersect ∂ B(R) (see Fig. 9). The case when f is singly periodic is similar (in the case that f is doubly periodic we do not have to worry about such saddles). Infinite-four-arm saddle points will be relevant when we prove Lemma 3.2. Fortunately we can control the number of four-arm saddles which occur, in terms of the boundary behaviour of f .  Intuitively this bound follows because as we raise the level past a saddle which is four-arm in B(R) or an infinite-four-arm saddle we separate two components of { f ≥ } ∩ ∂ B(R) (i.e. the boundary components are no longer in the same excursion set component). The total number of such separations which may occur is bounded above by the number of boundary excursion components at different levels, which in turn is bounded by the number of tangent points. This is formalised in the proof below.
We proceed by induction. Let A − denote the component of { f ≥ f (x 1 ) − } containing x 1 . Since x is lower connected and four-arm in M, A − will contain some S (1) j (it may contain two such elements: see Fig. 10), and so we make a preliminary association between x 1 and S (1) j . We now assume inductively that for each x 1 , . . . , x n we have either a preliminary association with some S (n) j (which differs across points) or a final association to some tangent point with value less than f (x n ).
We consider x n+1 and let j which has a preliminary association with some x i (where i ≤ n). Then A − must be simple, as it is a subset of S . We then make a preliminary association between x n+1 and S (n+1) p and make a new preliminary association between x i and S (n+1) q (so we remove the old association between x i and S (n) j ). If A − is contained in some S (n) j which does not have a preliminary association or if A − is not contained in any S (n) j then we choose some S (n+1) p contained in A − (which must exist since x n+1 is four-arm in M) and make a preliminary association between this and x n+1 . Now suppose that x i , where i ≤ n, has a preliminary association with S (n) j . If the maximum level of any point in S (n) j is less than f (x n+1 ) then we make a final association between x i and the highest such point in ∂ M ∩ S (n) j , which is by definition a tangent point. We also remove the preliminary association between x i and S for some k, and we then make a new preliminary association between x i and S (n+1) k . We perform this process for all preliminary associations which remain after the step described in the previous paragraph, which completes the inductive step. This algorithm will cease after the m-th step, at which stage if a point x i has a preliminary association with some S  Next we suppose that x is a saddle point, so that by Lemma 4.2, A c− is homotopy equivalent to A c+ with a 1-cell attached. So in particular, A c+ consists of either one or two components. First we note that if A c− is simple and does not intersect ∂ M, then any components of A c+ have both these properties, therefore N simple (c − ) − N simple (c+ ) ≤ 0. If A c− is simple but intersects ∂ M, then by repeating the sequential argument in the first paragraph of this proof we see that A c+ must intersect ∂ M and so cannot consist of two simple components disjoint from the boundary. Furthermore, if A c− is not simple, then by the above homotopy, A c+ cannot consist of two simple components. Combining these two observations show that N simple Finally let x be a tangent point, so that A c− is homotopy equivalent to A c+ with a k-cell attached for some k ∈ {0, 1, 2}. In each case, A c+ has at most two simple components disjoint from ∂ M, so |N simple (c − ) − N simple (c + )| ≤ 2. Now suppose that there are no critical or tangent points at level . Since  4N tang (∂ B(R)). Applying the first part of this lemma to each of the N E S,R terms here then completes the proof.
For the periodic cases, the argument is a little more technical since we cannot apply Lemma 4.6 to B(R) directly. Instead, we tile B(R) with translated parallelograms or rectangles, apply Lemma 4.6 to f T or f C(n) on each translated domain and aggregate the results.
Proof (Lemma 2.5 in the case of periodic functions) Let f be a doubly or singly periodic function satisfying the assumptions of Lemma 2.5. It suffices to prove the excursion set relationship, since the level set relationship will follow by the same argument as in the aperiodic case.
Suppose f is doubly periodic with periodic vectors y, z ∈ R 2 and associated parallelogram P. By assumption, f T is a Morse function almost surely. Suppose that { f T ≥ } has N simple simple components.
Let A be a component of { f ≥ } and A be the corresponding component of { f T ≥ }. If A is compact, then clearly A cannot contain a non-T-contractible loop, and so it is simple. Since f | ∂ P has a finite number of local extrema, { f ≥ } ∩ ∂ P has a finite number of components. Therefore if A intersects the boundary of a translated parallelogram P + n 1 y + n 2 z where n 1 , n 2 ∈ Z then it must contain the translation of one of these boundary components. So if A is unbounded, it must contain a path joining two translated versions of the same boundary component, which implies that A is not simple.
We can tile R 2 with the translated parallelograms P + n 1 y + n 2 z where n 1 , n 2 ∈ Z. We associate each bounded component A of { f ≥ } to a particular translation P +n 1 y+n 2 z of P such that (P +n 1 y+n 2 z)∩ A = ∅, N simple components are mapped to each translated parallelogram and if A is associated with P, then A + n 1 y + n 2 z is associated with P + n 1 y + n 2 z. If A is a component of { f ≥ } contained in B(R) then it must be associated with a parallelogram which intersects B(R), therefore where d = diam(P). Similarly, each compact component of { f ≥ } which is associated with a parallelogram inside B(R) must be in B(R) unless it intersects ∂ B(R), but the number of such intersections is bounded above by the number of tangent points of f to B(R). Therefore The number of translated parallelograms contained in the ball of a given radius can be approximated by a generalisation of Gauss' circle problem. It is shown in [11] that as R → ∞. So combining this with (20) and (21) as R → ∞.
Now suppose that { f T ≥ } contains m local maxima. Reasoning as above, it is clear that with an analogous result holding for lower connected saddles. Applying Lemma 4.6 to f T (i.e. with the setting M ∞ = M = T) then shows that as R → ∞ where the constant in the O(·) notation is independent of , as required (although it may depend on d).
Now suppose that f is singly periodic with periodic vector y = (y 1 , 0) and define R(n) = [ny 1 , (n + 1)y 1 ) × R. By assumption f almost surely has a finite number of tangent points in [0, y 1 ] × {n 1 } or {n 2 y 1 } × [−R, R] for any n 1 , n 2 ∈ Z or R > 0. Repeating the argument from the doubly periodic case then shows that the compact components of { f ≥ } correspond precisely to the simple components of { f C ≥ }.
As before, we can associate each compact component of { f ≥ } with a translated rectangle R(n) for n ∈ Z, in such a way that if A is a compact component of { f ≥ } which is mapped to R(n) then A intersects R(n), and if A = A + y then A maps to R(n + 1).
We recall our earlier definitions: if S(n 1 , n 2 ) = [n 1 y 1 , (n 1 + 1)y 1 ] × [n 2 , n 2 + 1] where n 1 , n 2 ∈ Z, then we define B int (x, R) to be the union over n 1 , n 2 ∈ Z of all S(n 1 , n 2 ) contained in B(x, R). Similarly we define B ext (x, R) to be the union over n 1 , n 2 ∈ Z of S(n 1 , n 2 ) which intersect B(x, R). We now fix R > 0 and define m n for n ∈ Z to be the largest integer such that and note that under the standard quotient relationship, this rectangle can be identified with C(m n ).
Let N simple ( f C(m n ) , ) denote the number of simple components of { f | C(m n ) ≥ } which do not intersect ∂C(m n ). By the mapping described above, there will be at least N simple ( f C(m n ) , ) compact components of { f ≥ } which are associated with R(n) and intersect [ny 1 , (n + 1)y 1 ]×[−m n , m n ]. If any of these components intersect ∂ B int (0, R), then they will do so at distinct components of { f | ∂ B int (0,R) ≥ } and so in particular the number of such components intersecting the boundary is at most the number of tangent points on ∂ B int (0, R). Therefore We can then apply Lemma 4.6 to f C(m n ) (i.e. with the setting M ∞ = C and M = C(m n )) for each n to get this inequality in terms of local maxima and lower connected saddles points. Let r = max{ √ 2, y 1 }. Then since B int (0, R) covers B(R − r ) we see that This gives the required lower bound for N E S (R, ). The upper bound follows from a very similar argument. First we define M n to be the largest integer such that Then and by applying Lemma 4.6 to each cylinder we see that as required.
We next complete the proof of Lemma 3.2. Again we shall separate the argument into the aperiodic case and the periodic cases, so that the reader interested only in the aperiodic case can access the simplest version of the argument. Since in this case we work only with the manifolds R 2 and B(R), we replace the condition that an excursion set component A is simple with the equivalent condition that it is bounded.
Proof (Lemma 3.2 in the case of aperiodic fields) Let f be an aperiodic stationary field satisfying Conditions 2.3. Let g : R 2 → R be a Morse realisation of f and assume that g| B(n) is Morse for all n ∈ N. Let x be a lower connected saddle of g at level .
The first step is to show that, with probability one, x is not also an upper connected saddle of g. For > 0 sufficiently small, let A − be the component of {g ≥ − } containing x and A + = A − ∩{g ≥ + }. From the definition of a lower connected saddle point, we have three cases to consider. In the first two cases we use the fact that g is Morse to (deterministically) rule out the possibility of x being upper connected; we show that the third case occurs with zero probability and can therefore be neglected. 2) A − is bounded and A + has two bounded components. The arguments in the previous case are also valid in this case where S is chosen to be either of the two components of A + .

1)A − is unbounded and
3)A − is unbounded and A + is bounded. We will show that when f is aperiodic, it almost surely has no saddle points of this type. We fix R > 0 and let g R = g| B(R) which we assume is Morse. For y ∈ R 2 we define A y,− to be the component of {g ≥ g(y) − } containing y and A y, = A y,− ∩ {g ≥ g(y) + }. Let x 1 , . . . , x n be the saddle points of g R for which A x i ,− is unbounded but A x i , is bounded for all > 0 sufficiently small. We then choose a fixed sufficiently small that this condition holds for each i and such that g R has no other critical points and no tangent points in are non-overlapping for i = 1, . . . , n. Finally we assume that the x i are ordered so that g(x i ) < g(x j ) for i < j.
Since A x i ,− is unbounded for each i, we see that is non-empty for each i. Since there are no tangent points with level in [g R (x i ) − , g R (x i ) + ], we can repeat the arguments in the proof of Lemma 4.6 to show that is non-empty and has the same number of components as B − i for each i = 1, . . . , n (see Fig. 11). Suppose there exists a point y ∈ B − i ∩ B − j for some i < j. Since A x j ,− is unbounded and path connected, there exists an unbounded path started at y and contained in {g ≥ g(x j ) − )}. Since y ∈ A x i , which is path connected, this path must also be contained in A x i , which is bounded. This contradiction implies that B − 1 , . . . , B − n are disjoint. Since each of these sets must have a local maximum, we see that n ≤ N tang (∂ B(R)). Since f is stationary, we know that the expected number of lower connected saddle points x contained in B(R) for which A x,− is unbounded and A x, is bounded for > 0 sufficiently small is c 0 R 2 for some c 0 ≥ 0. However by the argument just given, this number is bounded above by the expected number of tangent points of f to ∂ B(R) which is O(R) by Conditions 2.3. Hence c 0 = 0 so f almost surely has no saddle points of this type.  We have shown that the set of lower connected saddles of f is almost surely disjoint from the set of upper connected saddles. Using the same notation as above, we will now show that any saddle point of g must be lower connected or upper connected, completing the proof of the lemma. We now show that in the first case x will be lower connected. By symmetry, this will imply that in the second case, x is upper connected and so complete the proof for aperiodic fields. We fix a path γ in {g < } connecting the components of {g < }∩ N . From the existence of γ it is clear that A \{x} has two components B 1 , B 2 and that without loss of generality, B 1 is bounded (see Fig. 12).
For > 0 sufficiently small, A + therefore has at least one bounded component B 1 ∩{g ≥ + }. If A − is compact for some > 0 then clearly B 2 is also bounded and so A + has two bounded components so x is lower connected. If A − is unbounded for arbitrarily small we can assume that A + has an unbounded component, by the argument given above, and so x is again lower connected. We note that this also proves the equivalence of Definition 3.1 in the aperiodic case and Definition 1.7.  C depending on whether f is doubly periodic or singly periodic. We also assume that g| C(n) is Morse for all n ∈ N in the singly periodic case. Let x be a lower connected saddle of g at level .
Again the first step is to show that, with probability one, x is not an upper connected saddle of g. By reducing we can ensure that A + has two components, one of which is simple, and so both end-points of γ must be contained in the same component of M c ∩ {g ≤ − }. This will then be contained in a component of {g ≤ − } which we denote by B. If B is not simple, then clearly B ∪ γ is not simple. Suppose that B is simple but B ∪ γ is not simple, then B ∪ γ must contain a loop, denoted η, which is not M-contractible and in particular η must intersect γ . However B is path connected and B ∪ γ 'surrounds' a region which is homotopy equivalent to S and so must be simple (see Fig. 13 for an example of these sets in the doubly periodic case).
Therefore η is homotopy equivalent to a loop which is contained in B and so B is not simple, which is a contradiction. We have shown that B ∪ γ is simple if and only if B is, therefore M c ∩ {g ≤ − } and M c ∩ {g ≤ + } have the same number of simple components. This holds for any M c sufficiently large (possibly after reducing ) but if x were upper connected we could find M c such that M c ∩ {g ≤ − } has one more simple component than M c ∩ {g ≤ + } (again, we note that this argument is stated more formally in the proof of Lemma 4.6). We conclude that x is not an upper connected saddle point.

2)
A − is simple and A + has two simple components. Once again, the arguments in the previous case are also valid in this case where S is chosen to be either of the two components of A + .  A + is simple (in particular, connected). We can repeat the argument given in this section of the proof for aperiodic fields to show that if f is singly periodic then A − must be bounded. Clearly in the doubly periodic case, A − ⊂ T must also be bounded. We therefore choose a compact domain M c ⊂ M of the form T or C(n) for some n which contains a neighbourhood of A − . Let B 1 , . . . , B n denote the components of {g ≤ + } ∩ M c which intersect A + . Since A + is simple, we know that at most one of these sets is not simple and the remainder are. (One of the sets will 'surround' A + which in turn will 'surround' the remaining sets.) Without loss of generality, we assume B 1 is the 'surrounding' set, which may not be simple. By Theorem 4.2, A − is homotopy equivalent to A + with a 1-cell, denoted γ attached. Clearly γ is contained in B 1 and B 1 is not simple, otherwise A + ∪ γ would be simple. If M = C then B 1 \γ has two components and each of these components contains a loop which is not C-contractible (since A + ∪ γ is compact and so separated from ∞ and −∞ by B 1 \γ ). If M = T then B 1 \γ may have one or two components, in either case each component will contain a non-Tcontractible loop (to see this, consider a path on either side of γ which then traverses the boundary of A + ). This means that passing through the saddle point x can only create non-simple components of {g ≤ − }, so x is not upper connected.

3)A − is not simple and
This completes the proof that the sets of upper and lower connected saddle points of f are almost surely disjoint. We continue to use the notation defined above, and we now show that any saddle point of g must be lower connected or upper connected. As in the aperiodic case, Lemma 4.5 and Conditions 2.3 allow us to conclude that f almost surely has no infinite four-arm saddles.
Suppose that f is doubly periodic so that g is defined on the torus T. We fix a small neighbourhood N of x such that {g > } ∩ N has precisely two components denoted N 1 , N 2 . If γ : [0, 1] → T is a path contained in {g > } ∪ {x}, we say that γ cuts through x if x ∈ γ [0, 1] and for every t ∈ (0, 1) such that γ (t) = x, for > 0 sufficiently small γ ((t − , t)) ⊂ N 1 and γ ((t, t + )) ⊂ N 2 (intuitively, this means that the image of γ just before hitting x is always on the same side of the saddle point; see Fig. 14). We make an analogous definition for paths contained in {g < } ∪ {x}.
First we suppose that there exists a T-contractible loop γ in {g < }∪{x} which cuts through x. Since x is a saddle point, it is clear that A + must have two components, one of which is surrounded by γ and so is simple. If A − is simple, then both components of A + must be simple. If A − is not simple, we know from Lemma 4.6 that A + can contain at most one simple component. In either case, we see that x is lower connected. By symmetry, x is upper connected if there exists a T-contractible loop in {g > } ∪ {x} cutting through x. Now we suppose that all loops cutting through x are non-T-contractible. We define C + to be the union of {x} and the components of {g > } which x is in the closure of. We define C − to be the analogous set for {g < }. At least one of C + or C − must contain a non-contractible loop cutting through x (in order to stop any paths in C + cutting through x from joining to form a loop, there must be a loop in T\C + cutting through x which blocks them, and since x is the only critical point at level , such a 'blocking' loop can be found in C − ). We also note that both C + and C − must contain non-T-contractible loops (if C + was simple, then we could find a T-contractible loop in C − cutting through x, and this argument is symmetric).
Suppose that γ 1 is a loop contained in C + which cuts through x. If every non-Tcontractible loop in C + intersects x, then C + \{x} is simple, so A + is simple and hence x is lower connected. Therefore we may assume that C + contains a non-Tcontractible loop γ 2 which does not intersect x. Let η 1 and η 2 be two loops which generate the fundamental group of T, we denote the concatenation of n copies of η 1 and m copies of η 2 by nη 1 + mη 2 . We suppose γ 1 n 1 η 1 + m 1 η 2 and γ 2 n 2 η 1 + m 2 η 2 where denotes being T-homotopic. It is known (see [7,Section 1.2.3]) that any loop homotopic to n 1 η 1 + m 1 η 2 must intersect any loop homotopic to n 2 η 1 + m 2 η 2 unless n 1 m 2 = n 2 m 1 .
If n 1 m 2 = n 2 m 1 then we know that any non-contractible loop in C − must intersect either γ 1 or γ 2 . Since C − ∩ C + = {x} clearly such a loop must intersect γ 1 at x. By the previous paragraph this means that C − \{x} is simple and so x is upper connected.
If n 1 m 2 = n 2 m 1 then by choosing a path ξ in C + \{x} which joins γ 1 to γ 2 we consider the concatenated loop n 2 γ 1 + ξ + n 1 (−γ 2 ) + (−ξ) where − denotes inverting the direction of the path. This path is contained in C + and cuts through x (since γ 1 does but γ 2 and ξ do not hit x) and since n 1 m 2 = n 2 m 1 this loop is T-contractible. However this contradicts the above supposition, so this case is not possible. This completes the proof in the doubly periodic case. The proof in the singly periodic case is simply a repetition of parts of the proof for aperiodic and doubly periodic fields, so we omit it.
All that remains is to complete the calculations for the special class of degenerate fields.
Proof (Proposition 1.20) We recall that a random variable Y is Rayleigh distributed with parameter σ > 0 if for all x ≥ 0, and we denote this as Y ∼ Ray(σ ). If Y ∼ Ray(σ ) and θ is an independent uniform-[0, 2π ] random variable then it is well known that Y cos(θ ) ∼  N (0, σ 2 ).
Then f has the representation where X 0 ∼ N (0, α), Y 1 ∼ Ray( √ β), Y 2 ∼ Ray( √ γ ), the θ i are uniformly distributed on [0, 2π ] and all of these random variables are independent. (By well known properties of the Rayleigh distribution, the field defined by the right hand side of (22) will have Gaussian finite dimensional distributions and simple calculations show that the covariance function of this field is κ.) Let and P = {tu + sv : t, s ∈ [0, 1]}. Then P is the parallelogram associated with f and we consider f T to be the restriction of f to P when P is identified with the two dimensional torus. By rotating the axes, we may assume that K 1 > 0 (since this does not affect the definition of c L S ). Some basic calculations show that, on the event f T has four critical points which occur at different levels and are all non-degenerate. Therefore f T is almost surely a Morse function. Moreover, the critical points of f T −X 0 occur at levels −Y 1 − Y 2 , −|Y 1 − Y 2 |, |Y 1 − Y 2 |, Y 1 + Y 2 . Clearly the critical points at levels Y 1 + Y 2 and −Y 1 − Y 2 are a local maximum and local minimum respectively. We can use Lemma 4.6 to characterise the other two critical points and the number of simple components of { f T − X 0 ≥ }, denoted N simple ( ), at different levels.
Specifically, since { f T − X 0 ≥ Y 1 + Y 2 } = ∅ and { f T − X 0 ≥ −Y 1 − Y 2 } = T, we see that N simple ( ) = 0 whenever | | > Y 1 + Y 2 . Then, since the critical points at level Y 1 + Y 2 and −Y 1 − Y 2 must be a local maximum and local minimum respectively, applying Lemma 4.6 shows that N simple ( ) does not change as passes through −Y 1 − Y 2 and decreases by one as passes through Y 1 + Y 2 . Therefore Now assume that θ 1 = θ 2 = 0. Then by traversing the parallelogram P across one of its edges (depending on which of Y 1 , Y 2 is bigger) we can find a closed path which is not T-contractible, on which f T − X 0 is bounded below by |Y 1 − Y 2 |, so in particular is positive. This shows that { f T − X 0 ≥ 0} has a non-simple component (for general values of θ 1 , θ 2 , such a path will also exist, but it will be translated). Since { f T − X 0 ≥ |Y 1 − Y 2 |} consists of a single simple component, this implies that the critical point at level |Y 1 − Y 2 | must be a lower connected saddle point. Then by Lemma 4.6 we see that It is then clear that Arguments given in the proof of Lemma 2.5 for periodic functions show that Arguing as in the proof of Lemma 2.4 it can be shown that E(N tang (∂ B(R)) = O(R) as R → ∞. Then for ∈ R E N E S,R ( ) π R 2 − 1 Area(P) so that N E S,R ( )/(π R 2 ) converges in L 1 to the random variable 1 Area(P) 1 −X 0 ∈(|Y 1 −Y 2 |,Y 1 +Y 2 ] which is non-deterministic provided α > 0 or > 0.
In particular, this convergence shows that We can repeat these arguments for the sum of upper and lower excursion sets to derive the corresponding expression for c L S ( ).