Solution (source code)

= Solution

The <Zero-boundary Gaussian free field> on the unit disc is the centered <Gaussian process> indexed by finite <Borel measures> of finite <Green energy>, with <covariance>
$$
\mathbb E[(h,\rho)(h,\nu)]=\iint_{\mathbb D\times\mathbb D}G_{\mathbb D}(x,y)\rho(dx)\nu(dy).
$$
This covariance determines all its <finite-dimensional distributions>.

The <Domain Markov property of the Gaussian free field> says that for every suitable open $U\subset\mathbb D$ one can write
$$
h=h_U+h^{\mathbb D\setminus U},
$$
where $h_U$ is a zero-boundary Gaussian free field on $U$, independent of $h^{\mathbb D\setminus U}$, while $h^{\mathbb D\setminus U}$ is harmonic on $U$ and carries the information from the field outside $U$.