  Precipitation is a large-scale, spatio-temporally heterogeneous phenomenon,
with frequent anomalies exhibiting unusually high or low values. We use Markov
Random Fields (MRFs) to detect spatio-temporally coherent anomalies in gridded
annual rainfall data across India from 1901-2005. MRFs are undirected graphical
models where each node is associated with a \{location,year\} pair, with edges
connecting nodes representing adjacent locations or years. Some nodes represent
observations of precipitation, while the rest represent unobserved
(\emph{latent}) states that can take one of three values: high/low/normal. The
MRF represents a probability distribution over the variables, using \emph{node
potential} and \emph{edge potential} functions defined on nodes and edges of
the graph. Optimal values of latent state variables are estimated by maximizing
the posterior probability of the observations, using Gibbs sampling. Edge
potentials enforce spatial and temporal coherence, and node potentials
influence threshold for anomalies by affecting the prior probabilities of the
states. The model can be tuned to recover anomalies detected by threshold-based
methods. The competing influences of spatial and temporal coherence can be
adjusted through edge potentials.
  We study spatio-temporal properties of rainfall anomalies discovered by this
method, using suitable measures. We identify nonstationarities in occurrence of
positive and negative anomalies between the first and second halves of the 20th
century. We find that between these periods, there has been decrease in
rainfall during June-September (JJAS) and an increase during other months.
These effects are highlighted prominently in the statistics of anomalies.
Properties of anomalies learnt from this approach could present tests of
regional-scale rainfall simulations by climate models and statistical
simulators.
