Variation in salinity in soil in the Yangtze delta was analysed geostatistically.
Strong lateral and vertical trends and covariances of the residuals were estimated by REML.
The anisotropy could be represented by a non-affine transformation.
Space–time covariance functions led to satisfactory fitting.
Kriged predictions were mapped as three-dimensional charts.