Geospatial examination of lithium in drinking water and suicide mortality
© Helbich et al.; licensee BioMed Central Ltd. 2012
Received: 26 March 2012
Accepted: 24 May 2012
Published: 13 June 2012
Skip to main content
© Helbich et al.; licensee BioMed Central Ltd. 2012
Received: 26 March 2012
Accepted: 24 May 2012
Published: 13 June 2012
Lithium as a substance occurring naturally in food and drinking water may exert positive effects on mental health. In therapeutic doses, which are more than 100 times higher than natural daily intakes, lithium has been proven to be a mood-stabilizer and suicide preventive. This study examined whether natural lithium content in drinking water is regionally associated with lower suicide rates.
Previous statistical approaches were challenged by global and local spatial regression models taking spatial autocorrelation as well as non-stationarity into account. A Geographically Weighted Regression model was applied with significant independent variables as indicated by a spatial autoregressive model.
The association between lithium levels in drinking water and suicide mortality can be confirmed by the global spatial regression model. In addition, the local spatial regression model showed that the association was mainly driven by the eastern parts of Austria.
According to old anecdotic reports the results of this study support the hypothesis of positive effects of natural lithium intake on mental health. Both, the new methodological approach and the results relevant for health may open new avenues in the collaboration between Geographic Information Science, medicine, and even criminology, such as exploring the spatial association between violent or impulsive crime and lithium content in drinking water.
One of the first anecdotic evidence of effects of lithium containing water on mental health was mentioned by Cade . Since then, several studies have investigated the relationship between the naturally occurring lithium content in drinking water and suicide mortality [2–5], a relationship which is plausible since lithium is a well recognized mood-stabilizer and suicide preventive in psychiatric treatment nowadays. However, this relationship is also somewhat surprising, as the level of lithium for therapeutic use is multiple times higher than naturally occurring lithium. The first ecological study on the association of lithium in drinking water and suicide rates came from Texas, US. This study was based on aggregated data from 27 counties grouped into low, medium, and high lithium counties and statistically tested by t-tests to explore relationships between lithium and crime, suicide, etc. within each of the three groups . Twenty years later, utilizing more comprehensive statistics of weighted ordinary least-squares (OLS) regression, Ohgami et al.  confirmed the association using data on lithium in drinking water from 19 municipalities in Oita, a prefecture of Japan. This study sparked interest and was followed by two subsequent studies from Europe. Kabacs et al.  examined 47 subdivisions of counties in East England using a Pearson’s correlation coefficient to examine the association, but found results to be not significant. Finally, a nationwide study of all 99 districts of Austria based on the most comprehensive lithium measurement of 6,460 water samples, applied weighted univariate and multivariate OLS, again confirming the suicide preventive potential of natural lithium even after adjustment for confounding factors . However, the main limitation of all these studies is that they do not consider the unique properties of space, having serious consequences concerning the statistical validity of published results (see, e.g. [6–8]).
Essential special properties of geospatial data are spatial autocorrelation and spatial heterogeneity (non-stationarity). The first property implies a spatial association between an attribute value at a particular location and attribute values at other locations close by. The second property describes systematic spatial variation of attribute values across space [6, 7]. It is well established that these spatial effects must be taken into account when modeling the spatial relationship in a regression framework , such as the association between the amount of lithium in drinking water and mortality ratios or other explanatory variables across space. The consequence of ignoring both spatial effects makes statistical inference invalid and results in misleading conclusions (for a detailed discussion see, e.g. [6–8]).
Therefore, in contrast to all previous research [2–5] this study models both spatial autocorrelation and non-stationarity explicitly by means of global and local spatial statistical models and thus aims to challenge the hypothesis of lithium as a suicide preventive substance in natural distribution.
The same original data set as in  is used in this analysis. In total, 6,460 water samples from drinking water supplies from 99 Austrian districts were analysed for lithium. The average was 65.3 samples per district, with a range from 1 to 312 samples. The lowest measurable threshold lithium level by inductively coupled plasma optical emission spectrometry was 0.0033 mg/l. For the statistical calculations, lithium levels were averaged per district with a mean lithium level of 0.0113 mg/l (SD ± 0.027). Further, data on suicide rates adjusted as standardized mortality ratios (SMR; ), proportion of Roman Catholics, population density, and per capita income were obtained from Statistics Austria. Density of psychiatrists, general practitioners, and psychotherapists were obtained from the Austrian Medical Chamber and the Austrian Institute of Health (ÖBIG). Unemployment rates per district were obtained from the Austrian Public Employment Service (AMS) and were averaged for the available years 2005–2008. Further details were described previously in .
First, exploratory spatial data analysis is applied to demask spatial patterns which may contradict fundamental model assumptions of the subsequent regression. To start off, spatial autocorrelation and non-stationarity in the variables are explored using the global Moran’s I and the local G*-statistic. In this context “global” refers to testing for spatial autocorrelation for the entire study area at once and deriving a single value indicating whether spatial autocorrelation exists and what its strength is.
The Moran’s I is a widely used measure of global spatial autocorrelation, which tests whether there are some relationships between location and attribute values. A significant positive statistic indicates that nearby locations of similar attribute values are more spatially clustered than randomly distributed. In contrast, a significant negative statistic shows dissimilar values at nearby locations showing a more dispersed pattern.
This global measure is not capable to explore distinctive local features as well as the non-stationarity of a spatial process. Often, a non-stationary response variable is a first clue for spatially-varying relationships in a multivariate regression model. Therefore, Getis and Ord  have introduced the G*-statistic to detect clusters of high or low values by location. High positive values refer to “hot spots” and high negative values to “cold spots”. As such, “hot spots” can be described as areas, where districts with high levels are surrounded by other districts with high levels. In contrast, “cold spots” are clusters with low level districts surrounded by other low level districts. A crucial step in spatial modeling is the choice of an appropriate representation of space. Common choices for such representations are, for instance, contiguity, k-nearest neighbors, among others. For a summary see .
The classical ordinary least squares (OLS) model is widely used to model the global relationship between a response and one or more explanatory variables. OLS assumes, among other things that residuals are spatially independent. Residual autocorrelation captures unexplained similarities between neighboring districts, which can be a results of omitted variables or a misspecification of the regression model in general . Accounting for spatial effects reduces the magnitude of the prediction error, removes most of the systematic error, and leads to more reliable estimates . In this study the OLS model serves as the reference model.
Global models that account for spatial effects are spatial autoregressive models, introduced by Anselin . They comprise of two special cases, namely the spatial lag and the spatial error model. The spatial lag model extends the standard OLS regression model by including a spatially lagged dependent variable, which can be mostly interpreted as spill-over effects. Ignoring this lagged dependent variable leads the OLS model to be biased and inconsistent . This model will be referred to as the spatial autoregressive (SAR) model in the reminder of this article. The spatial error model addresses the presence of spatial autocorrelation by defining a spatial autoregressive process for the error term and, by doing so, captures unexplained similarities. Not considering this spatial process in the error term yields the OLS estimates to be inefficient, although unbiased . The Lagrange multiplier test statistic (LM) on the estimated OLS residuals helps to decide between these two alternative model specifications.
Global regression models assume a homogeneous behavior of the estimated parameters across space, which has often proven to be unrealistic. A way out is the application of the Geographically Weighted Regression (GWR), a local regression model that explores local spatial variations of the parameters . The GWR models spatial autocorrelation and spatial heterogeneity for subsets of the entire data set. Each subset is established around a regression point with near data points exhibiting a higher influence than more distant data points. The weighting is often based on a bisquare kernel function (e.g. ). The setting of an appropriate bandwidth length of the weighting function is crucial. The most common is the adaptive bandwidth, where its length is allowed to vary across space, depending on the density of the data points. In densely populated areas the kernel possesses a shorter bandwidth in contrast to regions with larger inter-point distances, where the bandwidth is longer . A kernel function with an adaptive bandwidth improves the goodness-of-fit, if data points are irregularly distributed across space. Despite the critique that the GWR is more suitable for exploratory analysis (see e.g. [17, 18] and references therein), it is a flexible model type to investigate spatially varying relationships, which is the focus of this research. GWR estimations are modeled for the same level of aggregation or spatial units as the original data and can be interpreted in the same way as OLS regression estimates.
In a step-by-step approach, an exploratory spatial data analysis, a global non-spatial regression model, a global spatial regression model, and finally a local spatial regression model were developed and applied to explore the association that lithium and other factors have on suicide mortality.
In the analysis to follow, the standardized mortality ratio (SMR) is the dependent variable. The eight independent variables include mean lithium level, population density, per capita income, proportion of Roman Catholics, psychiatrist density, psychotherapist density, general practitioner density, and unemployment rate.
Spearman correlation r s between standardized mortality ratios (SMRs) for suicide (2005–2009) and district characteristics
Lithium level, mean (mg/l)
Population density (per km2)
Per capita income (in 1,000 Euro)
Proportion of Roman Catholics, %
Unemployment rate, %
Psychiatrist density (per 10,000)
Psychotherapist density (per 10,000)
General practitioner density (per 10,000)
Global Moran’s I statistic for variables included in this study
Per capita income
Prop. of Roman Catholics
General practitioner density
Global regression estimates on standardized mortality ratios (SMRs) for suicide (2005–2009)
Per capita income
Prop. of Roman Catholics
Gen. practitioner density
Of the two different types of spatial autoregressive models, the (robust) LM test recommends the use of a spatial lag model (rob. LM: 8.477; p = 0.004). Again, variables are selected based on the full model and predictors that do not make a significant contribution on a 0.05 level are omitted stepwise. In the final global spatial lag model, both the lithium level as well as the proportion of Roman Catholics are significant predictors of the overall suicide SMR, the psychiatrist density is nearly significant at α = 0.1 level but its consideration in the model does not increase the AIC (Table 3, right column). Sensitivity analysis with different neigbourhood configurations confirm the robustnes of these results.
A Jarque Bera test (JB = 3.038; p = 0.219) shows that the residuals are normally distributed, homoscedastic (BP = 4.303; p = 0.231), and no longer significantly different from a spatially random distribution (i.e., they possess no spatial autocorrelation, Moran’s I: -0.060; p = 0.761). This means that the autoregressive component successfully captures spatial autocorrelation, which is indicated by a significant value of rho (ρ = 0.448; LR: 14.336; p < 0.001) and non-significant LM test for residual autocorrelation. In other words, residuals are now in compliance with the assumption of being spatially independent of each other. In addition, the resulting global spatial lag model showed much improvement over the final stepwise OLS model in terms of performance, based on the AIC values (−85 for the global spatial model versus −73 for the stepwise regression model). Overall, the SAR model explains more than 43% (Nagelkerke R 2 ) in the variation of the suicide SMR.
While the global spatial model is clearly an improvement over the aspatial OLS model, the relationships between the dependent and the independent variables remain stationary (i.e., constant) across the entire study area of Austria. However, multiple studies have indicated that such relationships are in fact non-stationary and thus vary across the study area (e.g., ). Local spatial regression models take such non-stationarity into account. Thus, the final step in this analysis involved estimating a GWR that used an adaptive, bisquare kernel function and the significant independent variables indicated by the SAR model.
Given the increasing possibilities of disease mapping due to a growing availability of digital data in medicine, the experiences of psychiatric epidemiology with geospatial analyses are growing . Up to now, it has been assumed that lithium in drinking water originated from natural sources such as lithium containing rocks . Natural lithium traces in water have already been considered by Cade  to have the potential to influence mental health. Since then, the mood-stabilizing effects of lithium became widely recognised by the psychiatric community  and in the meanwhile its suicide preventive properties are well documented . As a natural trace element, lithium is washed out by rain from rock and soil and dissolves in ground water and reaches the food-chain via the drinking water. In some geographic regions, its concentrations may reach up to 5.2 mg/l, reflecting a natural daily intake of lithium of up to 10 mg/l [24, 27]. Although, such daily doses of lithium are considerably lower than those used therapeutically, it is unknown to which extent intake of natural lithium levels may influence mental health or suicide mortality. Only one randomised placebo-controlled study showed favourable effects of low-dose lithium supplementation (0.4 mg daily) on mood, in a small sample of former drug users . Although, the effects of therapeutic doses of lithium (about 70 mg dissolved lithium per day) are well established, still only little is known about the health effects of natural or low lithium intake as examined in this study.
This study did not consider that lithium prescribed to patients may accumulate via waste water in ground water aquifers due to urinary excretion. There are several reports pointing to the public health problem of antibiotics leading to resistant bacteria in waste water [29, 30] as well as pharmaceuticals like carbamazepine, ibuprofen, and biphosphenol A in finished drinking water in Canada . Similarly, acetaminophen, caffeine, carbamazepine, codeine, and other substances are measureable in groundwater from the Los Angeles metropolitan area, when compared to other regions in California . In the United States, also Fluoxetine has been detected in 4.3% of all examined ground water sites . The effects of Fluoxetine and Carbamazepine concentrations found in sediment and water systems were strong enough to affect health and behavior in aquatic invertebrates [34, 35]. It is possible that lithium, which is used in psychiatry since more than 60 years, has accumulated in drinking water reservoirs and thus has contributed to the nonrandom distribution of lithium in drinking water across Austrian districts . With this in mind, the found association of lithium in drinking water and suicide mortality could be a function of regional lithium treatment rates rather than an effect of natural lithium occurrence. To answer this question, a further examination of this hypothesis is necessary and spatial regression models are most useful in this regards. Indeed, the authors are already in the process of preparing a project assessing lithium prescription rates at a geographical basis in Austria, which are available and can be purchased.
In particular, local spatial regressions allow to further challenge the hypotheses drawn from linear models. As shown in this study, the associations found between lithium in drinking water vary in strength by region. This underscores the necessity to examine large areas instead of selected regions. For example, the contradicting results in Ohgami at al.  and Kabacs et al. , which both examined only a random region of each country, a prefecture in Japan and East England, might have been biased by the selection of these regions. Such local spatial regression models are especially informative when multiple risk factors coalesce as causes of a disease. Then, as shown here, strongest risk factors may vary from region to region and such results may be used to inform policy makers and prevention initiatives, which than may be tailored to the local needs in addition to global spatial predictors. While the suicide mortality of the north-east of Austria is best described by the Proportion of Roman Catholics, the south-east is better described by lithium levels. The results suggest that in the south-eastern parts of Austria the density of psychiatrists coincidently exhibits a significant negative relationship with suicide mortality with lithium levels in water, which may be an argument for an influence of local lithium prescriptions by psychiatrists on lithium levels in drinking water.
The outcomes found here still need replications based on other countries data. Also, replications in other time periods would be informative to substantiate the possible effects of lithium on people’s health. The limitations of the results derived from such ecological models still include the problem of the ecological fallacy in particular, and the Modifiable Areal Unit Problem in general, and do not allow inferring direct causality from such associations but are of generative and explorative value for further basic research. The examined model of suicide mortality is an example for further interesting research questions. May regional patterns of violent crime incidence such as homicides and other violent or impulsive assaults be explained by lithium content in drinking water as previously hypothesized .
The research was funded by the Austrian Science Fund (FWF) through the Doctoral College GIScience (DK W 1237-N23). MH was funded by the Alexander von Humboldt Foundation.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.