Apply suggestions

src/field/gaussian/fft.jl

@@ -17,7 +17,7 @@ inverse Fourier transform.
 * `neighborhood` - Search neighborhood (default to `nothing`)
 * `distance`     - Distance used to find nearest neighbors (default to `Euclidean()`)
 
-### References
+## References
 
 * Gutjahr 1997. [General joint conditional simulations using a fast
   Fourier transform method](https://link.springer.com/article/10.1007/BF02769641)

src/field/gaussian/lu.jl

@@ -16,7 +16,7 @@ Gaussian are drawn via LU factorization.
 * `correlation`   - Correlation coefficient between two covariates (default to `0`)
 * `init`          - Data initialization method (default to `NearestInit()`)
 
-### References
+## References
 
 * Alabert 1987. [The practice of fast conditional simulations
   through the LU decomposition of the covariance matrix]

src/field/gaussian/seq.jl

@@ -24,7 +24,7 @@ For each location in the process `path`, a maximum number of
 neighbors `maxneighbors` is used to fit the conditional Gaussian
 distribution. The neighbors are searched according to a `neighborhood`.
 
-### References
+## References
 
 * Gomez-Hernandez & Journel 1993. [Joint Sequential Simulation of
   MultiGaussian Fields](https://link.springer.com/chapter/10.1007/978-94-011-1739-5_8)

src/field/lindgren.jl

@@ -11,7 +11,7 @@ and `sill` (total variance) as described in Lindgren 2011.
 The process relies relies on a discretization of the Laplace-Beltrami 
 operator on meshes and is adequate for highly curved domains (e.g. surfaces).
 
-### References
+## References
 
 * Lindgren et al. 2011. [An explicit link between Gaussian fields and
   Gaussian Markov random fields: the stochastic partial differential

src/field/strata.jl

@@ -16,10 +16,9 @@ as described in Hoffimann 2018.
 * `fillbase`    - Fill value for the bottom layer (default to `NaN`)
 * `filltop`     - Fill value for the top layer (default to `NaN`)
 
-### References
+## References
 
-Hoffimann 2018. *Morphodynamic analysis and statistical
-synthesis of geormorphic data.*
+* Hoffimann 2018. *Morphodynamic analysis and statistical synthesis of geormorphic data.*
 """
 struct StrataProcess{E,S,ST,N,FB,FT} <: FieldProcess
   environment::E

src/field/turing.jl

@@ -14,9 +14,9 @@ Turing process as described in Turing 1952.
 * `edge`   - edge condition (default to `nothing`)
 * `iter`   - number of iterations (default to `100`)
 
-### References
+## References
 
-Turing 1952. *The chemical basis of morphogenesis.*
+* Turing 1952. [The chemical basis of morphogenesis](https://royalsocietypublishing.org/doi/10.1098/rstb.1952.0012)
 """
 struct TuringProcess{P,B,E,I} <: FieldProcess
   params::P