1. Introduction
Drinking safe water is a critical public health element worldwide. Improved water supply and sanitation, and better management of water resources, can boost countries’ economic growth and can contribute greatly to poverty reduction [
1]. Anyway, in the next future, the availability of safe water could be a worldwide issue. Environmental challenges such as climate change and extreme events, water scarcity, and pollution are increasingly becoming a global concern [
2]. Many studies have shown how the production and consumption of bottled water (BW) are contributors to these problem [
3,
4] s[
3,
4]. Others have shown how tap water (TW) would be preferable due to the high environmental impact of bottle manufacturing [
5]. Moreover, increasing evidence has shown a higher level of exposure to microplastics from BW versus TW [
6] Nowadays, in many countries, there is the possibility to have access to cheap and safe potable TW, but BW consumption has been increasing globally[
7]. Recently there has been a growing interest in investigating possible reasons for preferring TW or BW. Some of them have focused attention on health and safety concerns about TW [
8,
9], while others on organoleptic properties [
10,
11]. The European Region has diverse TW consumption patterns between and within countries, with the highest share of TW intake observed in northern Europe. Moreover, different sociodemographic characteristics associated with TW or BW consumption have been investigated, delineating a complex scenario[
12]; in particular, gender differences have been highlighted, with women drinking more BW [
13,
14,
15] since women seem to perceive TW riskier but at the same time, they have more general environmental awareness. Moreover, ethnicity and culture, as well as self-perceived health status and diet choices are among the key drivers of TW and city dwellers tend to drink less TW, with the lowest consumption observed at restaurants and the highest in people’s homes [
16]. Additionally, research shows that free access to safe and acceptable drinking water may be beneficial in reducing sugary drinks[
17].
The aim of this study is to evaluate the consumption of TW and its association with selected sociodemographic- characteristics, self-perceived health status, and body Mass index of consumers at a nationwide level in Italy.
2. Materials and Methods
Information on tap and BW consumption has been obtained from the multi-purpose survey on families “Aspects of daily life”, carried out by the Italian National Statistics Institute (ISTAT). The survey is carried out on a yearly basis in a representative sample of the Italian population and is part of the integrated system of multipurpose surveys on families aimed at detecting a plurality of behavioral dimensions and segments of daily life. The questionnaire is standardized, and it is administered both in face-to-face mode by trained interviewers and in self-compilation and it has been used in several studies[
18,
19]. We analyzed data coming from the 2021 edition of the survey, including 45,597 subjects and 20.000 families. In the specific case, the populations for the current surveys, have been identified by the Italian National Statistics Institute within the set of municipalities which has been divided into two subsets: Municipalities with larger demographic size constitute a separate stratum and are defined as Self-Representative (SR); the remaining municipalities are defined as Non Self-Representative (NSR) and are divided, based on demographic size, into strata of equal breadth. From these strata, the sample municipalities (two for each stratum) were selected with probabilities proportional to their size. For each of the municipalities involved in the survey (SR and NSR), a cluster sampling is carried out: the clusters—the families—are randomly selected from the registry list, and all the members belonging to the actual family are surveyed. The minimum number of sample families for each municipality has been set at 24. The families are selected for each sample municipality from the theoretical sample selected for the Master Sample for each family included in the sample, the characteristics under investigation of all actual members belonging to the same family are recorded. The size of the theoretical sample in terms of families, set at the national level primarily based on cost and operational criteria, is approximately 24,000 families. The number of involved sample municipalities should not exceed 900, to allow for effective control and supervision. The allocation of the sample of families and municipalities, among the various regions, has therefore been calculated by adopting a compromise criterion to ensure both the reliability of estimates at the national level and that of estimates within each of the territorial domains[
20] The following variables have been included in the analyses: educational level (graduate/postgraduate degree, high school diploma, middle school diploma, primary school diploma/none), age (<20 years old, 20-44 years old, 45-59 years old, 60-74 years old, ≥75 years old), sex (males/females), economical resources in last 12 months (1 = inadequate, 2 = scarse, 3 = adequate, and 4 = optimal), geographical area of residency (North-West, North-East, Center, South, The Islands), concern about waste production (yes/no), concern about climate change (yes/no), consumption of carbonated drinks excluding water (1=more than 1 L/day, 2=between ½ / 1 L day, 3=less than ½ L, 4=less than occasionally, 5=occasionally, 6= no use), alcohol consumption (1=yes, 0=no), consumption of vegetables (1=more than once a day, 2=once a day, 3=sometimes a week, 4=less than once a week, 5=never), consumption of snacks (1=more than once a day, 2=once a day, 3=sometimes a week, 4=less than once a week, 5=never), body mass index (BMI ; 1=underweight, 2=normal, 3 overweight, 4= obese), self-perceived health status(1=satisfactory, 0=not satisfactory). Bivariate analyses were performed to study the association of TW consumption with relevant variables using Chi-square tests. Logistic regression models have been built to adjust for confounders and to evaluate the factors independently associated with TW consumption (1, if TW consumption is present; 0 if not). Regarding health outcomes, we considered self-perceived health status satisfaction level (dichotomized in two levels 1=satisfactory, 0=not satisfactory) and BMI (dichotomized in two 1=obese, 0= not obese). The level of significance was set at 0.05. Analyses have been performed with STATA, version 15 (Stata Corp.).
3. Results
Among the selected sample, people who usually drink TW were 19,674 representing the 43.15 (95% confidence interval, CI 42.69-43.60). In terms of prevalence, bivariate analysis (
Table 1) has highlighted a higher prevalence of TW consumption in people aged 45 to 59 years old (N=4942, 25.12%, p-value<0.05) with no statistically significant differences among males and females. People with a graduate/post-graduate degree diploma showed the highest prevalence of consumption (N=3111, 48.56%, p-value<0.05) as well as people with optimal economic resources (N=356, 50%, p-value<0.05). Considering concern for environmental issues our sample showed a higher prevalence of people consuming TW among those that were concerned about waste production (N=7963, 44.21%, p-value<0.05) and among those that were concerned about climate change (N=9586, 44.66%, p-value<0.05). Regarding the Geographical Area of Residency, there was a higher prevalence of TW consumption in the North-East regions of Italy (N=5966, 64.98%, p-value<0.05). Underweight people showed higher prevalence of TW assumption (N=509, 46.27% p-value<0.05) as well as those who drink less than occasionally carbonated drinks (N=1932, p-value<0.05), drink alcohol (N=12519, 44.66% p-value<0.05), consume vegetables more than once a day (N=3584, 46.83% p-value<0.05) and snacks less than once a week (N=7265, 45.10%, p-value=0.05). Regarding perceived health status highest prevalence of TW consumption has been shown from those with a satisfying level (N=14323, 44.03% p-value<0.05)
Multilevel regression analysis (
Table 2) confirmed that there are no statistically significant differences between males and females in terms of the likelihood of drinking TW (females OR=0.98, C.I. 0.94-1.03). All age classes from 20 to 74 years old are less likely to drink TW than people younger than 20 years old but only one class reached statistical significance (20-44 OR=0.75, C.I. 0.64-0.88). Regarding economic resources, there is an increasing trend of TW consumption but none of the several different categories reached statistical significance. On the other hand, our model highlights also how the level of schooling may influence TW assumption. Compared to the highest level of education (Graduate/Post-graduate) other educational classes show a low likelihood of consuming TW with a statical significance level (High School diploma OR= 0.8, C.I. 0.75-0.85, Middle School Diploma OR=0.72, C.I. 0.67-0.78, Elementary school diploma/none OR=0.7, C.I. 0.64-0.77). About the Geographical Area of residence, we did observe statical significant differences among different regions with increased consumption in Northwest (baseline) and North-East (OR=2.17, C.I.= 2.04-2.32) and less consumption in Center (OR=0.85, C.I.=0.79-0.9), South (OR=0.69, C.I.=0.65-0.74) and The Islands (OR 0.25, C.I.=0.23-0.28). In our model, we considered also some variables associated with attention to environmental issues. A concern about waste production was associated with an increased likelihood of assuming TW (OR=1.09, C.I.=1.04-1.14). A similar trend of likelihood has been shown for concern for climate change but is not statistically significant (OR=1.03, C.I.=0.99-1.08).
We then performed a multilevel regression analysis with perceived health status as the main outcome (
Table 3) and possible variables that were independently associated: consumption of tap water showed no statistically significant association (OR=0.93, C.I. 0.99-1.12) with satisfying perceived health status. This perception showed a decreasing trend of likelihood with age increase, and it reached statistical significance all among different age classes (20-44, OR=0.65, C.I. 0.49-0.87; 45-59, OR=0.34, C.I. 0.25-0.45; 60-74, OR=0.21, C.I. 0.16-0.28; ≥75 OR=0.12, C.I. 0.09-0.16). The same decreasing trend has been found regarding lower levels of education but only people with middle school diploma (OR=0.89, C.I. 0.81-0.98) and with elementary school diploma/none diploma (OR=0.66, C.I. 0.59-0.74) reached statistical significance. A low likelihood of a satisfying perceived health status has been found also with a decreasing availability of economic resources in the last 12 months, particularly for people who declared few resources (OR=0.51, C.I. 0.39-0.66) and insufficient (OR=0.33, C.I. =0.25-0.44). Again, self-perceived health status was worse in Central, Southern and island regions (OR=0.87, C.I. 0.8-0.95; OR=0.8, C.I. 0.74-0.87; OR=0.74, C.I: =0.67-0.82), as well as in underweight people (OR=0.69, C.I. 0.58-0.82), overweight (OR=0.91, C.I. 0.86-0.97) and obese people (OR=0.61 C.I. 0.56-0.66) with respect to normal weight subjects. Regarding consumption of carbonated drinks (excluding water) compared to those who declared to consume more than 1 liter per day, all other subjects showed a positive association trend with perceived health status that reached statical significance for people who declared a consumption from 1/2L to 1L/day (OR=1.41, C.I. 1.01-1.96) and an occasional consumption (OR=1.41, C.I. 1.06-1.87); as well, the presence of alcohol consumption was found to be positively associated with perceived health status (OR=1.33, C.I. 1.25-1.42). Considering the assumption of vegetables and consumption of snacks the likelihood of association with a satisfying perceived health status showed opposite trends since the lower the assumption of vegetables the lowest odds of association and this resulted to be significant for “less than once a week” (OR=0.77, C.I. 0.68-0.88) and “never” (OR=0.71, C.I. 0.6-0.86), while the lower the consumption of snacks, compared with those who declare consumption “more than once a day”, the highest odds of associations with a statistical significance all among classes of assumption (“Once a day” OR=1.,46 C.I. 1.04-2.07; “Sometimes a week” OR=1.,91 C.I. 1.4-2.59; Less than once a week OR=1.,93 C.I. 1.43-2.63; Never OR=1.,84 C.I. 1.35-2.5).
We then performed a multilevel regression analysis with obesity as the main outcome (
Table 4) and possible variables that were independently associated: consumption of TW showed a statistically significant low likelihood association (OR=1.05, C.I. 0.86-0.99) of being obese. The obese status also showed an increasing trend of likelihood with aging, and it reached statistical significance among all different age classes (20-44 years, OR=1.71, C.I. 1.15-2.54; 45-59, OR=2.31, C.I. 1.56-3.42; 60-74, OR=2.64, C.I. 1.77-3.91; ≥75 OR=1.93, C.I. 1.28-2.9). The same trend has been found regarding lower levels of education statistically significant for high school diplomas (OR=1.39, C.I. 1.23-1.57), middle school diplomas (OR=1.5, C.I. 1.32-1.7), and elementary school diploma/none diploma (OR=1.69, C.I. 1.46-1.95). A low likelihood of satisfaction with self-perceived health status has been found also with a decreasing availability of economic resources in the last 12 months, particularly for people who declared few resources (OR=0.51, C.I. 0.39-0.66) and insufficient (OR=0.33, C.I. =0.25-0.44). An increasing trend of association for being obese in southern parts of Italy, compared to northern ones, has been found but only the south geographical area reached a statically significant association (OR=1.14, C.I. 1.04-1.26). A satisfying perceived health status was negatively associated with the likelihood of being obese (OR=0.68, C.I. 0.62-0.73) as well as alcohol consumption (OR=0.89, C.I. 0.82-0.96). Regarding consumption of carbonated drinks (excluding water) compared to those who declared to consume more than 1 liter per day, all other subjects showed a negative association trend with being obese but none of the classes reached statical significance. Considering the assumption of vegetables and consumption of snacks the likelihood of association with being obese showed opposite trends since the lower the assumption of vegetables the highest odds of association, while the lower the consumption of snacks, compared with those who declare consumption “more than once a day”, the lowest odds of associations. However, the association for vegetables and snacks did not reach statistical significance for any of the classes of consumption.
4. Discussion
Safe drinking water represents a worldwide major concern because of several factors such as pollution and climate change that are heavily impacting human health all around the globe, in particular, among residents of developing countries. Impacts on surface water and groundwater resources and water-related illnesses are increasing, especially under changing climate scenarios such as diversity in rainfall patterns, increasing temperature, flash floods, severe droughts, heat waves, and heavy precipitation [
21]. Recently, it has been shown a link between health and environmental impacts and drinking water choices, since it has been estimated that the environmental impact of BW is 1400-3500 higher than TW[
22]. On the other hand, BW consumption has sharply increased in the last years worldwide [
23] and interestingly the recent increase in BW use globally has been driven by an increase in demand in low- and middle-income countries (LMICs), despite parallel increases in access to piped water in some countries [
24]. To better understand these trends, it’s extremely important to study possible variables associated with increased trap water consumption that can allow policy decision-makers to take adequate actions. Our study fits in this contest since we investigated several demographic and socioeconomic variables that can be associated with increased consumption of TW. Unlike other studies[
13,
14] our regression model does not observe a statistically significant association of TW consumption with sex; on the other hand, regarding age, the younger the age classes were less likely TW consumers, in particular people aged 20-44 years old. This result agrees with other studies [
12]. Interestingly, lower levels of schooling compared to graduate/post-graduate level, showed a reduced frequency of TW consumption and this is partially in disagreement with other findings [
25]. It is possible that in our country, people with higher levels of schooling have an increased awareness about the safety and importance of consuming TW also for environmental issues, as testified by our result in terms of the presence of concern about waste production since we found a statistically significant association of increased likelihood of TW consumption for this variable. So, generally speaking, an increased level of knowledge may be associated with increased awareness and environmentally friendly behaviors. It must be highlighted that concern for climate change showed a similar trend, but it did not reach statistical significance in terms of association of TW assumption, indicating a complex scenario in terms of knowledge of environmental issues by the general population. An important consideration regards the Geographical area of residency: in the regression model we observed a statistically significant difference among different geographical areas and we found a decreasing trend of TW consumption from North to South, particularly evident on the Islands where there is the lowest likelihood of TW consumption. This result may be explained by the lower presence of natural water sources between the Northern and the Southern parts of Italy, especially for islands like Sicily and Sardinia where water is particularly scarce(Fresh Water Resources - Italy - Climatechangepost.Com).
Looking at health outcomes we focused on the perceived health status declared by subjects and we studied possible association with a positive attitude and consumption of TW together with other variables included in the logistic model. Even if TW assumption showed an increased likelihood of association with a satisfying perceived health status, the association did not reach statistical significance. On the other hand, other variables showed a coherent association in agreement also with studies conducted by others since the likelihood of perceiving a satisfying health status decreases with age[
27], with underweight or overweight/obese BMI classes[
28], with lower availability of economic resources[
29] and lower levels of education[
30,
31]. Moreover, we found a worsening perception of own health status from North to south as found by others[
32,
33]. Also, consumption of snacks and consumption of vegetables showed interesting trends with a low likelihood of a satisfying perceived health status associated with lowering the assumption of vegetables that reached statical significance for the lowest ones identified as “less than once a week” and “never”. Regarding the consumption of snacks, we found that the lowest classes of assumption of this type of food were associated with an increased likelihood of perceived satisfying health status that resulted to be statistically significant all among the several classes compared to the baseline of subjects who declared a consumption of more than “once a day”. Both results are in accordance with other studies[
34,
35]. Regarding others liquid intake besides TW, a decreasing consumption of carbonated drinks, compared to those who declared an assumption of more than 1L/day, showed a positively associated trend to an increased likelihood of a satisfying perceived health status, resulting particularly significant for the “occasionally” class of consumption in agreement with other works that found an association with depression and low perceived health status [
36,
37]. Taken together, these considerations may allow us to speculate that the model proposed describes quite accurately the reality of variables independently associated with a satisfying perceived health level including consumption of TW that showed a positive likelihood of association even if not statically significant.
Looking at BMI we found that consumption of TW showed a possible protective effect for being obese. This could be explained by the fact that individuals who prefer to consume TW are more likely to have both a healthier lifestyle and lower consumption of refined, carbonated, and sugary beverages[
17,
38] which are responsible for weight gain[
39], cardiovascular issues[
40], diabetes[
41] and overall poorer glycemic control[
42]. Also, this model showed coherent associations with other independent variables that have been confirmed by other studies since we found an increased likelihood of being obese with aging[
43,
44], with lower educational level[
44,
45], and for those subjects with a geographical area of residency in the southern parts of Italy[
46].
To our knowledge, this is the first study in Italy, to take into consideration the consumption of TW and possible sociodemographic variables associated and health outcomes such as perceived health status and BMI. Moreover, our data cover a nationwide sample of thousands of subjects, and interviews and data collection are conducted under rigorous methodological methods since they have been carried out by ISTAT-trained personnel.
Nevertheless, some limitations should be acknowledged. Given the nature of this study, as a prevalence study rather than a longitudinal one, certain details regarding individuals’ dietary habits might lack precision. This is testified, for example, by the absence of the exact amount of TW consumption that could have been more informative about water attitudes. Also, we did not have specific information on possible specific pathological conditions, that could have been associated with the assumption of TW, but only aggregated variables like the perceived health status that we used in our models. Moreover, all data must be considered declarative generating a possible declarative or recall bias in our sample. Lastly, it’s important to note that the official language of ISTAT, the government institution conducting the survey, is Italian. The survey administration was not declared in other languages, which may have resulted in a selection bias. This bias could have excluded individuals, even residents in Italy, who did not fully understand the Italian language.