In this paragraph we will concentrate on measures for which
Note that because of Eq. (1) p(t,x) can be regarded as the (defective) probability distribution of patches over local population density x.
Following Metz & Diekmann [9] we can represent the metapopulation dynamics by the forward Kolmogorov functional partial differential equation
for which we must define a boundary condition at x=0.
When total disasters may occur, i.e. , empty patches are formed at a rate
When a total disaster does not result in a loss of habitat quality, the patch will be re-colonized by dispersers instantly. The resulting boundary condition is
which simply states that the flux into habitable patches of zero population size is the same as the rate of total disasters. A result of habitat quality being unaffected by total disasters is that there can be no empty patches. This is in contrast to what was assumed by Gyllenberg & Hanski [4]. In their model, patches hit by a total disaster are set aside in a pool of empty patches ( ), instead of contributing to the pool of suitable patches. Colonization of empty patches, at the arbitrarily level x=1, happens due to the law of mass action involving dispersers and empty patches. The rate of change in the fraction of empty patches is then
resulting in the boundary condition
which is dependent on the density of dispersers.
Suppose now that there is a probability u that after a total disaster a patch becomes unsuitable for habitation, and that the time required to regain habitat quality is taken from an exponential distribution with parameter , we then get
and
independent of D.
This is just a simple model describing habitat recovery. More complex models might require the incorporation of an age distribution of empty patches. The resulting model will certainly be different from the model of Gyllenberg & Hanski, except for the case where the size of the instream of dispersers affects the recovery of habitat, i.e. is a function of D.
The problem becomes well posed if we specify the initial conditions. For ease of reference we repeat the full model