Original article | Open | Published:
Does waterfall aerosol influence mucosal immunity and chronic stress? A randomized controlled clinical trial
Journal of Physiological Anthropologyvolume 36, Article number: 10 (2017)
The specific microclimate of alpine waterfalls with high levels of ionized water aerosols has been suggested to trigger beneficial immunological and psychological effects. In the present three-armed randomized controlled clinical study, we focused on effects on (i) immunological reagibility, on (ii) physiological stress responses, and on (iii) stress-related psychological parameters.
People with moderate to high stress levels (n = 65) spent an active sojourn with daily hiking tours in the National Park Hohe Tauern (Großkirchheim, Austria). Half of the group was exposed to water aerosol of an alpine waterfall for 1 h/day (first arm, n = 33), whereas the other half spent the same time at a distant site (second arm, n = 32). A third arm (control, n = 26) had no intervention (except vaccination) and stayed at home, maintaining their usual lifestyle. The effect of the interventions on the immune system was tested by oral vaccination with an approved cholera vaccine and measuring specific salivary IgA antibody titers. Lung function was determined by peak expiratory flow measurement. Electric skin conductance, heart rate, and adaption of respiration rate were assessed as physiological stress parameters. Psychological stress-related parameters were analyzed by questionnaires and scales.
Compared to the control group, both intervention groups showed improvement of the lung function and of most physiological stress test parameters. Analysis of the mucosal immune response revealed a waterfall-specific beneficial effect with elevated IgA titers in the waterfall group. In line with these results, exposure to waterfall revealed an additional benefit concerning psychological parameters such as subjective stress perception (measured via visual analog scale), the Global Severity Index (GSI), and the Positive Symptom Total (PST).
Our study provides new data, which strongly support an “added value” of exposure to waterfall microclimate when combined with a therapeutic sojourn at high altitude including regular physical activity.
Therapies at high altitude and green exercise are known to affect a variety of physiological and immunological parameters. These include neurovegetative, cardiovascular or thermoregulation mechanisms [1, 2], and also the reduction of inflammatory diseases, and the induction of balancing immunomodulatory effects [3–5].
High-altitude climate therapy is also characterized as a successful alternative medical treatment for respiratory and allergic diseases such as bronchial asthma, atopic dermatitis, psoriasis, or eczema [6–8].
In addition to high altitude, European alpine regions host numerous waterfalls, which produce inhalable, negatively charged nano-water particles known as “Lenard ions” or ballo-electric ions. Negative air ions nearby waterfalls (hereinafter we call them “ionosols”) are generated by aerosolization of water droplets at an obstacle, aqueous surface, or by aerodynamic break-up during free fall, undergoing charge redistribution forming “water bags” with negatively charged surfaces and fracturing into micro-bubbles due to water shearing. After breaking up, smaller fragments carry negative charge and remain in the air for some time carried by air stream. Their diameter is between 1.5 and10 nm, whereby 2-nm-sized negative ions were the most abundant. Lifetime of ionosols is long enough so that they can be inhaled. The remaining bigger fragments are positive and precipitate to the ground [9, 10].
“Negative air ions” are positively correlated with relative humidity and their concentration is higher in natural environments compared in urban environment [11–14]. Waterfalls not only produce high levels of negative air ions, they also create a specific microbiological atmosphere by spreading of microbes in the impact zone of a waterfall, which may influence immunological and physiological parameters. Negative air ions also interact with phytoncides released from trees which in turn can interact with microbes, altogether influencing the environmental atmosphere . Phytoncides may also act directly on the innate immune system, as they have been shown to increase the level of natural killer cells .
Positive immunomodulating effects of the microclimate in close proximity to the impact zone of alpine waterfalls have already proven beneficial effects for the treatment of allergic asthma . In this previous study, half of a group of asthmatic allergic children (n = 54) was exposed to water aerosol of an alpine waterfall for 1 h/day over 3 weeks, whereas the other half spent this hour at a water aerosol free control site. The exposure to water aerosol had a long-lasting beneficial effect on asthma symptoms and lung function, accompanied by a decreased inflammatory immune status as indicated by elevated levels of IL-10 and regulatory T cells.
Mechanisms underlying the immunomodulatory effect of high altitude include the normalization of eosinophil levels, balancing of the TH2/Treg cell ratio, increase of regulatory cytokines (e.g., IL-10), and even alteration of immunoglobulin class switching [3, 4, 16, 17]. Furthermore, alterations of adrenocorticotropic hormone (ACTH) and glucocorticoid secretion indicate that high altitude affects both neuronal as well as immune circuits . Negative air ions are also associated with psychological well-being, e.g., due to an increased serotonin level . Considering, it is conceivable that high-altitude climate may also affect stress and its pathologic immunological, physiological, and psychical consequences. Stress is known to induce endocrine changes via mediators such as glucocorticoid hormones and catecholamines and thus also to influence humoral and cellular immunity . In a well-balanced but fragile molecular and cellular network, the immune, endocrine, and nervous system affect cognitive performance, well-being, and behavior, and, in general, also maintain health [20, 21]. For example, people with high occupational stress, like caregivers, show impaired specific and unspecific immune responses, both concerning humoral and cellular immunity [22, 23], and chronic stress is associated with a shift from a TH1 to a more TH2-type immune response cytokine pattern . Similar T cell response types are also observed in the development of asthma, which is not only influenced by stress, anxiety, and depression but also affected by the microbiome . Also, mucosal immune responses show a clear correlation with stress. As a consequence of chronic stress, the abundance of salivary IgA is decreased [26–29], and even more subtle psychological parameters (like disgruntlement) seem to influence salivary IgA levels [30–32]. As secretory IgA is present in all mucosal surfaces, it is of crucial importance for the first line of defense of the immune system [33, 34]. Therefore, monitoring specific antibodies in a workload situation, e.g., during a vaccine-induced immune response, represents a suitable model to estimate the immune status of an individual. Several studies demonstrate the strong impact of psychological stress on a developing humoral response, e.g., chronic burden (like nursing care of a sick member of the family) interferes with the humoral immune response to an influenza vaccination [33, 35–39].
Typical and useful physiological parameters indicative for stress and emotional strain are peripheral physiological signals such as electric skin conductance (SC), respiration rate (RSP), and heart rate (HR), all being measurable in response to an artificial stressor [40–42]. The conductance of the skin is dependent of the quantity of sweat, which is controlled by the sympathetic nervous system and positively correlates with stress. Also, respiration reflects the psychological and emotional state, e.g., in relaxed situations, respiration is slow and regular, whereas fast and irregular rhythms of breath are ascribed to stress emotions . Similarly, the heart rate and its variability correlate with comforting impacts and defensive reactions. Stress and increased sympathetic activity normally lead to an increase of heart rate, whereas parasympathetic activity has the opposite effect [44, 45]; thus, effective stress adaption is known to increase general well-being .
Based on the data showing an immunomodulatory effect of the waterfall aerosols, we hypothesized that the waterfall microclimate may exert a positive effect on different stress parameters.
Therefore, the present randomized controlled clinical study investigated whether in combination with a high-altitude climate therapy, waterfall-generated ions would provide an effective and cost-efficient therapy to mitigate various stress-related symptoms and strengthening immunity.
As a trigger of many secondary diseases such as neuropsychiatric, cardiovascular, or metabolic disorders, stress induces enormous economic costs [47, 48]. Interestingly, nature itself may preserve a sustainable health resource to counteract this trigger.
Ninety-one persons working in care professions (19–61 years old; 44 men, 47 women) were enrolled in the study. People working in care professions are known as high-risk occupational groups concerning psychological and immunological consequences of stress [19, 24]. Requirement for study inclusion was moderate to high stress levels as defined by the Trier Inventory for the Assessment of Chronic Stress (TICS) questionnaire, assessed 7 days before beginning of the study (Fig. 1), and no previous cholera vaccination. Exclusion criteria were fever, chronic or severe diseases, including immune dysfunctions, and treatment with antibiotics or immunosuppressants. The study was approved by the ethical committee of Salzburg (E1066/29042009).
The study was set up as a randomized controlled clinical trial with three arms.
Except for the non-intervention (control) group (CO, n = 26), all participants spent 1 week at the village of Großkirchheim (Carinthia, Austria) located 1024 m above sea level, hosted in hotels and receiving the same meals. For 1 h each day, the groups were separated for intervention into a waterfall- (WF+, n = 33) and a “non-waterfall” group (WF−, n = 32); group separation was kept identical throughout the study. Individuals of the control group stayed at home (without any intervention, except vaccination), following their usual daily life.
The probands of the WF+ group were exposed for 1 h to the Gartl Waterfall at 1000 m a.s.l., whereas the WF− group spent this hour at a site about 1000 m away from the waterfall, but at the same altitude. The average concentration of ionosols at this “control” site was 840 ions/cm3 (mean air humidity 72.3% ±23.1, mean temperature 15.3 °C ±2.5), compared to the average 42660 ions/cm3 with maxima up to 57510 ions/cm3 at the waterfall exposed site (mean air humidity 84.6% ±9.9, mean temperature 16.9 °C ±4.1) . In order to exclude effects based on different trainings, this location was selected at the same altitude as the waterfall, with identical cumulative elevation gain and distance from the starting point.
Except for this 1-h intervention, participants of the WF+ and WF− group had the same daily routine, and both groups went through a daily hiking program at an altitude between 1400 and 2900 m, with structured hikes of at least 600 meters of cumulative elevation gain.
All participants received an oral cholera vaccine on day 0 and day 6 of the study (Fig. 1).
The study took place from July to August 2009; follow-up examinations were performed after 2 months (Fig. 1). Physiological and psychological data from day 0 and day 6 were evaluated on site in a medical field laboratory in Großkirchheim, Carinthia. Follow-up examinations on day 16 and day 66 were held in the Paracelsus Medical University Salzburg. Randomization was performed with “Random allocation software” program .
Environmental parameters of the Gartl Waterfall
The waterfall is located in the heart of the National Park Hohe Tauern in Carinthia (Austria), situated in the community Großkirchheim in the upper Möll Valley. It is easily accessible via a gently rising 1.5-km walk from the village center and well established for waterfall exposition. Its approximate drop height is 50 m in two cascades and the average water flow is 0.6 m3/s (max 0.9 m3/s). Despite the moderate drop height, the waterfall features very high concentrations of ionosols (mean 42660 ions/cm3, max 57510 ions/cm3) floating the whole valley beneath the impact zone. Air humidity near the impact zone is 60–90%, consisting of respirable small water drops, whereas ozone levels are virtually zero.
Mucosal immunity was tested by oral vaccination with DUKORAL® from SBL Vaccines, which was given on days 0 and 6 (Fig. 1). Dosage and mode of application were according to the package information leaflet [50, 51].
Collecting and analysis of saliva
Saliva (2–3 ml) was collected from all subjects in sterile 15-ml plastic tubes on days 1, 6, 16, and 66 in the morning (Fig. 1), immediately after waking up (thus containing the matutinal IgA concentration peak) . The tubes were stored at –80 °C until preparation and analysis. To determine DUKORAL® vaccine specific salivary IgA concentrations, 96-well high-bind, flat-bottomed immunoplates (FluoroNunc, Roskilde, Denmark) were coated overnight at 4 °C with DUKORAL® vaccine antigens (1 mg recombinant cholera toxin B subunit (rCTB)) diluted 1:100 with PBS (1×), corresponding to 0.01 mg per well. After blocking for 1 h at room temperature with 200 μl/well PTB (PBS/0.1% Tween 20/0.5% bovine serum albumin (BSA)), the plates were washed with PBS/0.1% Tween 20, and saliva was added at a dilution of 1:100 in 100 μl/well PTB for specific IgA and incubated for 1 h at RT. The plates were washed with PBS/0.1% Tween 20 and incubated with 100 μl/well horseradish peroxidase (HRPO)-conjugated goat anti-human IgA (BioRad, Austria), diluted 1:1000 in PTB for 1 h at RT. Again, the plates were washed (with PBS) and the reaction was detected by adding luminol (Boehringer Mannheim, Germany) (5 ml solution A (100 mM glycine pH 10 (NAOH)/0.4 mM luminol/8 mM iodophenol), 5 mL solution B (0.12% H2O2 in dH2O), and 10 ml dH2O). Chemiluminescence was measured by a Promega ELISA-plate luminometer (Promega, Mannheim, Germany). Luminescence was calculated as photon counts per second.
Lung function and physiological stress test
The respiratory parameter analyzed was peak expiratory flow (PEF) by means of a forced ex-in maneuver. All tests were performed with EasyOneTM Plus Diagnostic from ndd Medical Technologies by aid of two trained scientists.
A computer-guided physiological stress test was performed to measure heart rate (three-point ECG), respiration rate (respiration sensor), and skin conductance (finger clip) during a 3-min baseline phase, a 1-min stressor phase (optical and acoustical stressors on the computer screen), and a 4-min post-stress recovery phase using a NeXus 10 multimodal device from Mind Media (Mind Media BV, Netherlands). The arithmetic mean of the respective physiological signal during the stressor phase was set to 100%. The arithmetic mean of the post-stress recovery phase was compared to the 100% line and used to measure the stress adaption of the peripheral signals heart rate, respiration rate, and skin conductance.
Questionnaires and scales
Trier Inventory for the Assessment of Chronic Stress
This questionnaire was used to differentiate between various facets of chronic stress, measured by the retrospective rate of stress events in six areas (worries, work overload, work discontent, lack of social recognition, social stress, and intrusive memories) on a five-point rating scale .
Visual analog scale
VAS was performed as a measurement tool for subjective stress perception. On basis of a linear scale, participants can state their current stress level by indicating a position along a continuous line between the two end-points “no stress at all” and “unbearable stress”. The analog aspect of linear scales creates preferably metrical characteristics instead of discrete scales .
Symptom Check List–SCL 90
The SCL 90 questionnaire enables to analyze treatment or progress of a variety of symptom dimensions, such as somatization, obsessive-compulsive, interpersonal sensitivity, depression, anxiety, hostility, phobic anxiety, aggression, paranoid ideation, and psychoticism. The test provides an overview and displays intensity of symptoms at the moment and in progress, including three global indices, i.e., Global Severity Index (GSI) (describing overall psychological distress), Positive Symptom Distress Index (PSDI) (describing intensity of symptoms) and Positive Symptom Total (PST) (describing the number of self-reported symptoms) .
Maslach Burnout Inventory (MBI-D©)
The original Maslach Burnout Inventory was developed to measure intensity and incidence of perceived burnout in care professions. The authorized German version MBI-D, closely aligned with the original, surveys three components of burnout syndrome based on several specific items: emotional exhaustion, depersonalization, and personal accomplishment .
List of complaints (KSb–S BL)
This questionnaire is an independent part of the test-battery KSb–S (the german abbreviation KSb–S stands for “clinical self-assessment-scales”, BL for “list of complaints”, PSYCHIS Munic), which captures the degree of subjective impairment due to physical, e.g., joint pains and general complaints (e.g., faintness). Combination of quantity and levels of manifestations (heavy–moderate–little–not at all) result in a total interference of complaints, reflecting a constitutional change over time .
Mental state scale (KSb–S Bf–S)
Like the BL, the Bf–S is also an autonomous part of the test-battery KSb–S which is arranged to reflect the change of the present, psychical, and subjective mental state in an aggregate value, figuring the whole range of normal and pathologic changes of well-being .
This questionnaire figures an individual’s current recovery-stress state using retrospective (last 3 days) information of stressful situations, reactions, and recovery activities in a self-assessment test. Answers are queried in a seven-point scale, ranging from never to always .
All analyses were performed using the IBM SPSS statistics version 22 (IBM, NY, USA, http://www.spss.com/). Different linear mixed models, three-armed (LMM1, WF+, and WF− group compared to CO group) and two-armed (LMM2, WF+ group compared to Wf− group) were used for statistical analysis of all variables with more than two time points and high individual variations. Unlike a linear regression analysis, a linear mixed model is able to detect changes despite high individual variations in the IgA levels and questionnaires. Time, treatment, and interaction of treatment and time were set as fixed factors. To take individual differences into account, the patient IDs were included in the models as random effect. A third linear mixed model (LMM3) was used to detect changes over time for the main outcome variable of secretory IgA (sIgA) levels, excluding the interaction of time and treatment, with only time and treatment as fixed factors of each. Linear regressions were performed to detect intervention effects in all physiologic parameters and the visual analog scale, lacking of high-individual variations and inset at two time points. To identify significant differences of the baselines between the three study groups, comparison of means was done with ANOVA including Bonferroni correction for parametric data; Wilcoxon signed-rank test or Mann–Whitney U test was applied to non-parametric data.
Figures are shown as means ± standard deviation (SD). Statistical significance was expressed as p ≤ 0.05 (*) or as p ≤ 0.01 (**).
Furthermore, all study data were subjected to Kolmogorov–Smirnov (with Lilliefors correction of significance) and Shapiro–Wilk normality tests—70% of the data generated passed the tests.
The randomized, controlled clinical study comprised 91 participants working in care professions consisting of 47 females (mean BMI 22.50 ± 3.51 SD) and 44 males (mean BMI 24.33 ± 3.16 SD), ranging from 19–61 (30.93 mean ± 10.39 SD) years of age. No significant differences were found between baseline values of both intervention groups (WF+ and WF– group). The control group (CO) had a bias in age and was slightly younger compared to both other groups.
Response to mucosal vaccination
The vaccination protocol led to a detectable antigen-specific mucosal immune response in approximately one half of the immunized subjects. Responders and non-responders were defined by measuring vaccine-specific salivary sIgA titers. A value of three standard deviations above the IgA titer of the serum taken prior to immunization (pre-serum) served as baseline. Participants with values above baseline were determined as responders (∑n = 48 (n (WF+) = 16; n (WF–) = 18; n (CO) = 14)). Data of non-responders (∑n = 43, equates ~47%) got excluded from statistical analysis of specific IgA levels.
A comparison of IgA titers of both intervention groups (WF+ and WF–) with the non-intervention control group with a linear mixed model (LMM1) shows no significant difference but a statistical trend (p = 0.055) on day 16 of the WF+ group, compared to that of the CO group (Fig. 2). The strongest percentage increase at all time points can be found in the WF+ group (Table 1).
The specific impact of the waterfall ionosols was evaluated with a linear mixed model (LMM2) of the two intervention groups WF+ and WF–. The result indicates a significant increase of salivary IgA titers in the WF+ group on day 6 and day 66 (Fig. 3).
A linear mixed model excluding interaction of time and treatment (LMM3) showed significant changes over time in both intervention groups at all time points (day 6 p = 0.000, day 16 p = 0.000, day 66 p = 0.001).
Subjective stress perception
The visual analog scale assessing the subjective stress perception was applied on day 1 and day 6. A significantly lower stress level could be measured in the WF+ group after the intervention week on day 6, compared in the control group (lin. regression, p = 0.000, R 2 = 0.615), and a clear trend to a waterfall-specific beneficial influence when compared to the WF– group (lin. regression, p = 0.054, R 2 = 0.179) (Fig. 4).
Psychological stress parameters and burnout symptoms
Linear mixed model (LMM1) analysis of the SCL-90 questionnaire revealed improvement in six of the ten symptom dimensions in both intervention groups after 6 days (compared in the control group): aggression, day 6: WF+ p = 0.000; WF– p = 0.006, day 66: WF+ p = 0.058; obsessive-compulsive, day 6: WF+ p = 0.005; WF– p = 0.023; depression, day 6: WF+ p = 0.006; WF− p = 0.036; paranoid ideation, day 6: WF+ p = 0.006, day 66: WF+ p = 0.023; phobic anxiety, day 6: WF+ p = 0.005; WF− p = 0.042, day 66: WF+ p = 0.002; interpersonal sensitivity, day 6: WF+ p = 0.034 (Fig. 5). No differences could be detected with respect to the parameters somatization, anxiety, hostility, and psychoticism (data not shown) . Significant long-lasting effects (day 66) could be measured in the waterfall group for two parameters (phobic anxiety and paranoid ideation). Again, the data indicate a clear trend of a beneficial effect of the waterfall environment.
Furthermore, on day 6, both the Global Severity Index assessing overall psychological distress, as well as the Positive Symptom Total, which displays the number of self-reported symptoms, were significantly decreased in the waterfall group only (p = 0.009 for GSI and p = 0.03 for PST), again indicating a waterfall-specific effect (Fig. 6). The Positive Symptom Distress Index, which describes the intensity of symptoms, elicited no differences over time between the three study groups (data not shown).
Measuring incidence and severity of burnout by means of the Maslach Burnout Inventory (MBI-D) elicited a significant improvement in the WF+ group in the category depersonalization, compared in the control group on day 66 (day 6 p = 0.064, day 66 p = 0.002). No differences could be measured concerning the other two components, i.e., emotional exhaustion and personal accomplishment over time in all groups. With respect to the list of complaints (BL), reflecting constitutional changes during the trial, the WF+ group showed a significant positive long-term effect on day 66 (p = 0.011). Both intervention groups show a decrease of complaints on day 16 (WF+ p = 0.086, WF− p = 0.057), and the mental state scale (Bf–S) clearly indicates enhanced well-being short-term effects on day 6 (WF+ p = 0.051, WF− p = 0.06). The recovery and stress questionnaire (EBF) revealed less stress (day 6, p = 0.068) and improved recovery only in the waterfall group (day 6 p = 0.07, day 66 p = 0.036). All statistical analyses were done with LMM1 (data not shown).
Lung function and stress-associated peripheral physiological parameters
No significant differences between the groups regarding all determined physiological parameters and peripheral signals could be detected at the beginning of the study. Exposure to the waterfall had a significant positive effect on lung function as measured via peak expiratory flow (PEF) on day 6 (WF+ to CO p = 0.023, R 2 = 0.346). No significant effect could be detected in the WF− group (WF− to CO p = 0.136, R 2 = 0.047) or between the two intervention groups after 1 week of intervention (WF+ to WF− p = 0.359, R 2 = 0.079) (Table 2, Fig. 7).
All peripheral signals of the physiological stress test decreased in the W+ and W− group after 1 week. The results show the arithmetic mean in percentage of the 4-min post-stress recovery phase, compared to a 100% reference line, which reflects the arithmetic mean of the 1-min stressor phase in percentage. Results below 100% indicate the ability to reconstitute a parasympathetic physiological state, to calm down, and to relax after a stressor within the 4-min post-stress recovery phase. Compared to the control group, the respiration rate was decreased significantly in both intervention groups (WF+ to CO p = 0.016, R 2 = 0.164; WF− to CO p = 0.004, R 2 = 0.169). Similarly, skin conductance indicating sympathetic activity and sweat production was reduced after 1 week of waterfall exposition (WF+ to CO p = 0.017, R 2 = 0.151), with a trend to a waterfall-specific additional effect (WF+ to WF− p = 0.056, R 2 = 0.187). The heart rate was significantly lowered in both intervention groups (WF+ to CO p = 0.011, R 2 = 0.421; WF− to CO p = 0.000, R 2 = 0.492) (Table 2, Fig. 7).
The present controlled and randomized clinical study addressed the question whether the specific environment of a waterfall provides additional beneficial effects for prophylactic and therapeutic stress management when combined with high-altitude climate therapy and physical activity (mountain hiking in the protected area National Park Hohe Tauern in Austria). As stress represents a complex phenomenon including psychological, physiological, and immunological effects, we investigated selected parameters of all of these three aspects. For this purpose, two study groups with moderate to high stress levels spent an active 1-week sojourn in the National Park Hohe Tauern with identical parameters concerning daily mountain hiking, accommodation and food supply, except for 1 h/day. One group spent this time at the waterfall, the other at the same sea level but free from the waterfall-specific environment. A third study arm was recruited as a control group, keeping their normal daily life without any intervention.
The results of our study are in line with publications indicating positive health effects of green exercise and high-altitude climate therapy, which have been already proven to contribute to the improvement of respiratory and/or allergic diseases [3–5, 16]. Our data suggest to expand the applicability of this simple and cost-effective health provision for treatment of stress-related symptoms. Interestingly, several of the investigated parameters were significantly changed with a 1-h/day exposure to the environment of a waterfall. Most strikingly, beyond psychological and physiological stress-related symptoms, the specific humoral sIgA immune response after oral vaccination was significantly improved by the additional waterfall exposure. Secretory IgA is present in all mucosal surfaces; it is of crucial importance for the first line of defense of the immune system at mucosal barriers [33, 34].
Currently, we know that the immense power of the falling water cracks small water droplets in a specific way, resulting in small nano particles carrying a negative charge, whereas the created bigger fragments are charged positively. The specific environment around waterfalls is characterized by these mainly negatively charged “Lenard’s ions” hovering in the air, while the positive-charged fragments sink to the ground [10, 60–62]. The airborne nano-aeorosol is assumed to trigger a variety of biological effects, e.g., mild activation of the immune system, inducing a balance between TH1/TH2 immunity, stabilizing the autonomous nervous system, and improving blood flow [16, 63–65]. In a mouse model, water-generated negative ions have been shown to enhance cytotoxic activity of natural killer cells . It is also conceivable that the high concentration of negative air ions near waterfalls could influence the human microbiota. Furthermore, the waterfall environment itself provides a specific microbiologic atmosphere, which may directly effect the microbiota of skin and mucosal surfaces . This waterfall-altered microbiome could act as an immunological adjuvant and thus be responsible for the observed effects on mucosal immunity.
Negative air ions have also been shown to influence psychologic well-being by an increased serotonin level  and to enhance positive emotional processing in seasonal affective disorder  and have also been associated with lower depression scores . An “added value” concerning the mood-raising effect may simply originate from the extraordinary beauty of the waterfall itself .
Numerous curing and healing effects are ascribed to waterfalls in ancient traditions and folk wisdom in many regions of the world. The present randomized, controlled clinical study provides evidence for an added value of a daily 1-h stay for 1 week in a waterfall environment in combination with green exercise and high-altitude climate therapy. Furthermore, the data point to an influence of the waterfall ionosols on complex “psychoneuroimmunological” regulatory circuits.
The results of this study also provides a rational basis for further research with an aim to develop schedules for new and effective prophylactic therapies for high-risk patients suffering from psychological and physiological stress symptoms.
Cepilova Z, Porubska M, Slamova J. Climate therapy of uveitis in a mountain environment. Cesk Oftalmol. 1993;49:90–4.
Butykova L, Kolesar J, Turecekova A, Najvarova S. Changes in the elimination of histamine, noradrenaline, adrenaline and vanilmandelic acid during high-mountain climatic therapy at Strbske Pleso. Fysiatr Revmatol Vestn. 1973;51:103–7.
Karagiannidis C, Hense G, Rueckert B, Mantel PY, Ichters B, Blaser K, Menz G, Schmidt-Weber CB. High-altitude climate therapy reduces local airway inflammation and modulates lymphocyte activation. Scand J Immunol. 2006;63:304–10.
Simon HU, Grotzer M, Nikolaizik WH, Blaser K, Schoni MH. High altitude climate therapy reduces peripheral blood T lymphocyte activation, eosinophilia, and bronchial obstruction in children with house-dust mite allergic asthma. Pediatr Pulmonol. 1994;17:304–11.
van Leeuwen WS. Bronchial asthma in relation to climate. Proc R Soc Med. 1924;17:19–26.
a Porta B, Barandun J, Wuthrich B. Atopic neurodermatitis—therapy in high altitude climate. Praxis (Bern 1994). 2000;89:1147–53.
Engst R, Vocks E. High-mountain climate therapy for skin diseases and allergies—mode of action, therapeutic results, and immunologic effects. Rehabilitation (Stuttg). 2000;39:215–22.
Rijssenbeek-Nouwens LH, Bel EH. High-altitude treatment: a therapeutic option for patients with severe, refractory asthma? Clin Exp Allergy. 2011;41:775–82.
Kolarž P, Gaisberger M, Madl P, Hofmann W, Ritter M, Hartl A. Characterization of ions at Alpine waterfalls. Atmos Chem Phys. 2012;12:3687–97.
Laakso L, Hirsikko A, Gronholm T, Kulmala M, Luts A, Parts TE. Waterfalls as sources of small charged aerosol particles. Atmos Chem Phys. 2007;7:2271–5.
Wang XL, Li CR, Xu JW, Hu DM, Zhao ZL, Zhang LD. Air negative ion concentration in different modes of courtyard forests in southern mountainous areas of Jinan, Shandong Province of East China. Ying Yong Sheng Tai Xue Bao. 2013;24:373–8.
Skromulis A, Noviks G. Atmospheric light air ion concentrations and related meteorologic factors in Rezekne city, Latvia. J Environ Biol. 2012;33:455–62.
Mao GX, Cao YB, Lan XG, He ZH, Chen ZM, Wang YZ, Hu XL, Lv YD, Wang GF, Yan J. Therapeutic effect of forest bathing on human hypertension in the elderly. J Cardiol. 2012;60:495–502.
Craig JM, Logan AC, Prescott SL. Natural environments, nature relatedness and the ecological theater: connecting satellites and sequencing to shinrin-yoku. J Physiol Anthropol. 2016;35:1.
Li Q, Kawada T. Effect of forest environments on human natural killer (NK) activity. Int J Immunopathol Pharmacol. 2011;24:39S–44.
Gaisberger M, Sanovic R, Dobias H, Kolarž P, Moder A, Thalhamer J, Selimovic A, Huttegger I, Ritter M, Hartl A. Effects of ionized waterfall aerosol on pediatric allergic asthma. J Asthma. 2012;49:830–8.
Hummelshoj L, Ryder LP, Poulsen LK. The role of the interleukin-10 subfamily members in immunoglobulin production by human B cells. Scand J Immunol. 2006;64:40–7.
Wu CF, Lai CH, Chu HJ, Lin WH. Evaluating and mapping of spatial air ion quality patterns in a residential garden using a geostatistic method. Int J Environ Res Public Health. 2011;8:2304–19.
Glaser R, Kiecolt-Glaser J. How stress damages immune system and health. Discov Med. 2005;5:165–9.
Harbuz MS, Chover-Gonzalez AJ, Jessop DS. Hypothalamo-pituitary-adrenal axis and chronic immune activation. Ann N Y Acad Sci. 2003;992:99–106.
Calcagni E, Elenkov I. Stress system activity, innate and T helper cytokines, and susceptibility to immune-related diseases. Ann N Y Acad Sci. 2006;1069:62–76.
Levy MM. Caring for the caregiver. Crit Care Clin. 2004;20:541–7. xi.
Lambert VA, Lambert CE, Ito M. Workplace stressors, ways of coping and demographic characteristics as predictors of physical and mental health of Japanese hospital nurses. Int J Nurs Stud. 2004;41:85–97.
Glaser R, MacCallum RC, Laskowski BF, Malarkey WB, Sheridan JF, Kiecolt-Glaser JK. Evidence for a shift in the Th-1 to Th-2 cytokine response associated with chronic stress and aging. J Gerontol A Biol Sci Med Sci. 2001;56:M477–82.
Trueba AF, Ritz T, Trueba G. The role of the microbiome in the relationship of asthma and affective disorders. Adv Exp Med Biol. 2016;874:263–88.
Deinzer R, Kleineidam C, Stiller-Winkler R, Idel H, Bachg D. Prolonged reduction of salivary immunoglobulin A (sIgA) after a major academic exam. Int J Psychophysiol. 2000;37:219–32.
Deinzer R, Schuller N. Dynamics of stress-related decrease of salivary immunoglobulin A (sIgA): relationship to symptoms of the common cold and studying behavior. Behav Med. 1998;23:161–9.
Jemmott 3rd JB, Magloire K. Academic stress, social support, and secretory immunoglobulin A. J Pers Soc Psychol. 1988;55:803–10.
Mouton C, Fillion L, Tawadros E, Tessier R. Salivary IgA is a weak stress marker. Behav Med. 1989;15:179–85.
Stone AA, Cox DS, Valdimarsdottir H, Jandorf L, Neale JM. Evidence that secretory IgA antibody is associated with daily mood. J Pers Soc Psychol. 1987;52:988–93.
Stone AA, Cox DS, Valdimarsdottir H, Neale JM. Secretory IgA as a measure of immunocompetence. J Human Stress. 1987;13:136–40.
Stone AA, Marco CA, Cruise CE, Cox DS, Neale JM. Are stress-induced immunological changes mediated by mood? A closer look at how both desirable and undesirable daily events influence sIgA antibody. Int J Behav Med. 1996;3:1–13.
Burns VE, Drayson M, Ring C, Carroll D. Perceived stress and psychological well-being are associated with antibody status after meningitis C conjugate vaccination. Psychosom Med. 2002;64:963–70.
Mestecky J, Russell MW, Jackson S, Brown TA. The human IgA system: a reassessment. Clin Immunol Immunopathol. 1986;40:105–14.
Kiecolt-Glaser JK, Glaser R, Gravenstein S, Malarkey WB, Sheridan J. Chronic stress alters the immune response to influenza virus vaccine in older adults. Proc Natl Acad Sci U S A. 1996;93:3043–7.
Miller GE, Cohen S, Pressman S, Barkin A, Rabin BS, Treanor JJ. Psychological stress and antibody response to influenza vaccination: when is the critical period for stress, and how does it get inside the body? Psychosom Med. 2004;66:215–23.
Vedhara K, Cox NK, Wilcock GK, Perks P, Hunt M, Anderson S, Lightman SL, Shanks NM. Chronic stress in elderly carers of dementia patients and antibody response to influenza vaccination. Lancet. 1999;353:627–31.
Glaser R, Kiecolt-Glaser JK, Bonneau RH, Malarkey W, Kennedy S, Hughes J. Stress-induced modulation of the immune response to recombinant hepatitis B vaccine. Psychosom Med. 1992;54:22–9.
Morag M, Morag A, Reichenberg A, Lerer B, Yirmiya R. Psychological variables as predictors of rubella antibody titers and fatigue—a prospective, double blind study. J Psychiatr Res. 1999;33:389–95.
Yumatov EA, Pertsov SS, Dudnik EN, Kramm MN, Strelkov NO. Monitoring of human physiological functions during the daily night sleep. Patol Fiziol Eksp Ter. 2015;3:47–54.
Azam MA, Katz J, Fashler SR, Changoor T, Azargive S, Ritvo P. Heart rate variability is enhanced in controls but not maladaptive perfectionists during brief mindfulness meditation following stress-induction: a stratified-randomized trial. Int J Psychophysiol. 2015;98:27–34.
Dieleman GC, Huizink AC, Tulen JH, Utens EM, Tiemeier H. Stress reactivity predicts symptom improvement in children with anxiety disorders. J Affect Disord. 2016;196:190–9.
Ritz T, Dahme B, Dubois AB, Folgering H, Fritz GK, Harver A, Kotses H, Lehrer PM, Ring C, Steptoe A, Van de Woestijne KP. Guidelines for mechanical lung function measurements in psychophysiology. Psychophysiology. 2002;39:546–67.
Hosseini SA, Khalilzadeh MA, Naghibi-Sistani MB, Homam SM. Emotional stress recognition using a new fusion link between electroencephalogram and peripheral signals. Iran J Neurol. 2015;14:142–51.
Del Rosso S, Abreu L, Webb HE, Zouhal H, Boullosa DA. Stress markers during a rally car competition. J Strength Cond Res. 2016;30:605–14.
Pusenjak N, Grad A, Tusak M, Leskovsek M, Schwarzlin R. Can biofeedback training of psychophysiological responses enhance athletes’ sport performance? A practitioner’s perspective. Phys Sportsmed. 2015;43:287–99.
Bejean S, Sultan-Taieb H. Modeling the economic burden of diseases imputable to stress at work. Eur J Health Econ. 2005;6:16–23.
Kalia M. Assessing the economic impact of stress—the modern day hidden epidemic. Metabolism. 2002;51:49–53.
Saghaei M. Random allocation software for parallel group randomized trials. BMC Med Res Methodol. 2004;4:26.
Cholera vaccines: WHO position paper. Wkly Epidemiol Rec. 2010; 85:117-128.
Rahman A, Rashu R, Bhuiyan TR, Chowdhury F, Khan AI, Islam K, LaRocque RC, Ryan ET, Wrammert J, Calderwood SB, et al. Antibody-secreting cell responses after Vibrio cholerae O1 infection and oral cholera vaccination in adults in Bangladesh. Clin Vaccine Immunol. 2013;20:1592–8.
Hucklebridge F, Clow A, Evans P. The relationship between salivary secretory immunoglobulin A and cortisol: neuroendocrine response to awakening and the diurnal cycle. Int J Psychophysiol. 1998;31:69–76.
Schulz PSW. The Trier Inventory for the Assessment of Chronic Stress (TICS): Scale construction, statistical testing, and validation of the scale work overload. In: Diagnostica, vol. 45. pp. 8-19. Goettingen: Hogrefe Verlag; 1999. p. 8-19.
Reips UD, Funke F. Interval-level measurement with visual analogue scales in Internet-based research: VAS Generator. Behav Res Methods. 2008;40:699–704.
Franke GH. SCL-90-R. Die Symptom-Checkliste von Derogatis, Deutsche Version (Manual). Göttingen: Beltz Test GmbH; 2002.
Büssing A, Perrar KM. Die Messung von Burnout. Untersuchung einer deutschen Fassung des Maslach Burnout Inventory (MBI-D). Diagnostica. 1992;38:328–53.
Zerssen DDB-LaT, 25, 1908 (1971). Die Beschwerden-Liste als Test. Therapiewoche 21. 1971; 25:1908-1914.
Dv Z, Koeller D-M, Rey E-R. Die Befindlichkeitsskala (B-S)—ein einfaches Instrument zur Objektivierung von Befindlichkeits-Störungen, insbesondere im Rahmen von Längsschnitt-Untersuchungen. Arzneim Forsch. 1970;20:915–8.
Kallus K. Der Erholungs-Belastungs-Fragebogen. Manual. Frankfurt/M: Swets Test Services; 1995.
Luts A, Parts TE, Laakso L, Hirsikko A, Gronholm T, Kulmala M. Some air electricity phenomena caused by waterfalls: correlative study of the spectra. Atmos Res. 2009;91:229–37.
Parts TE, Luts A, Laakso L, Hirsikko A, Gronholm T, Kulmala M. Chemical composition of waterfall-induced air ions: spectrometry vs. simulations. Boreal Environ Res. 2007;12:409–20.
Vostrikov AA, Drozdov SV, Rudnev VS, Kurkina LI. Molecular dynamics study of neutral and charged water clusters. Comput Mater Sci. 2006;35:254–60.
Krueger AP, Reed EJ. Biological impact of small air ions. Science. 1976;193:1209–13.
Palti Y, De Nour E, Abrahamov A. The effect of atmospheric ions on the respiratory system of infants. Pediatrics. 1966;38:405–11.
Takahashi K, Otsuki T, Mase A, Kawado T, Kotani M, Ami K, Matsushima H, Nishimura Y, Miura Y, Murakami S, et al. Negatively-charged air conditions and responses of the human psycho-neuro-endocrino-immune network. Environ Int. 2008;34:765–72.
Yamada R, Yanoma S, Akaike M, Tsuburaya A, Sugimasa Y, Takemiya S, Motohashi H, Rino Y, Takanashi Y, Imada T. Water-generated negative air ions activate NK cell and inhibit carcinogenesis in mice. Cancer Lett. 2006;239:190–7.
Harmer CJ, Charles M, McTavish S, Favaron E, Cowen PJ. Negative ion treatment increases positive emotional processing in seasonal affective disorder. Psychol Med. 2012;42:1605–12.
Perez V, Alexander DD, Bailey WH. Air ions and mood outcomes: a review and meta-analysis. BMC Psychiatry. 2013;13:29.
Joye Y, Bolderdijk JW. An exploratory study into the effects of extraordinary nature on emotions, mood, and prosociality. Front Psychol. 2014;5:1577.
We greatly appreciate all our study participants which enabled the implementation of the study. We want to thank Peter Suntinger, deputy mayor of Grosskirchheim, Austria, and Franziska Fellner, community of Grosskircheim, for their important support. Many thanks to General Director Dr. Christian Heu, Grossglockner Hochalpenstrassen AG for his important input and discussion. We also want to thank Mag. Peter Rupitsch, director of the Hohe Tauern National Park, and the National Park Rangers Konrad Mariacher and Johanna Dullnig, for their help and guidance. We also want to thank Antonia Suntinger for her medical support, Veronika Granig for her input, and Mag. Brigitta Bernhard for the discussion of the psychological data. Also, many thanks to the Hotel Post, Hotel Schlosswirt, Hotel Marx, and Hotel Dorfwirt for their great hospitality.
This work was funded by the Oesterreichische Forschungsfoerderungsgesellschaft FFG, KNet Water Project 223.
CG carried out the statistical analysis, data analysis, visualization and interpretation, laboratory analysis, and preparation of the manuscript. MG participated in the study design and made operative involvement in conducting the study on site, laboratory analysis, data warehousing, and data analysis. JP participated in the statistical analysis. MR carried out the medical advising and proofreading of the manuscript. PK executed the physical characterization of the waterfall environment. CP made operative involvement in conducting the study on site and carried out the proofreading. JT participated in the preparation and proofreading of the manuscript. AH conceived of the study idea and the study design, made operative involvement in conducting the study on site, and participated in the laboratory analysis, data warehousing, data analysis, and interpretation of the data. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
About this article
- Ionized water aerosol
- High-altitude climate therapy
- Chronic stress
- Green exercise
- Burnout prevention
- Mucosal immune response
- Mountain hiking
- Alpine space