Participants (n = 5017) contributed an average of 3.9 measurements per person (n = 19 705 total) over the course of 20 years. A lowess smoother was used to calculate the smoothed average. FEV1 indicates forced expiratory volume in first second of expiration; FVC, forced vital capacity.
Associations between continuous current and lifetime exposure measurements and pulmonary function were modeled via cubic splines (see “Methods”). All 4 exposure measurements were included in each model (one model each for forced expiratory volume in the first second of expiration [FEV1] and forced vital capacity [FVC]). Mixed models with a random intercept and a random 3-knot age spline were used to adjust for repeated measures, and fixed effects were included for year, center and center-year (their interaction), race-sex category, education, and asthma; cubic splines for age, height, waist circumference, secondhand smoke exposure, and exposure to airborne particulate matter less than 10 microns and less than 2.5 microns in size; and interactions between the age spline variables and race-sex, asthma, waist spline variables, and height spline variables. Point estimates and confidence intervals for slopes and net associations at different exposure levels are provided in Table 3.
Pletcher MJ, Vittinghoff E, Kalhan R, Richman J, Safford M, Sidney S, Lin F, Kertesz S. Association Between Marijuana Exposure and Pulmonary Function Over 20 Years. JAMA. 2012;307(2):173–181. doi:10.1001/jama.2011.1961
Author Affiliations: Department of Epidemiology and Biostatistics (Drs Pletcher and Vittinghoff and Mr Lin) and Division of General Internal Medicine, Department of Medicine (Dr Pletcher), University of California, San Francisco; Asthma-COPD Program, Division of Pulmonary and Critical Care Medicine, Northwestern University Feinberg School of Medicine, Chicago, Illinois (Dr Kalhan); Department of Surgery (Dr Richman) and Division of Preventive Medicine (Drs Safford and Kertesz), University of Alabama at Birmingham; Center for Surgical, Medical and Acute Care Research and Transitions, Veterans Affairs Medical Center, Birmingham (Drs Richman and Kertesz); and Division of Research, Kaiser Permanente of Northern California, Oakland (Dr Sidney).
Context Marijuana smoke contains many of the same constituents as tobacco smoke, but whether it has similar adverse effects on pulmonary function is unclear.
Objective To analyze associations between marijuana (both current and lifetime exposure) and pulmonary function.
Design, Setting, and Participants The Coronary Artery Risk Development in Young Adults (CARDIA) study, a longitudinal study collecting repeated measurements of pulmonary function and smoking over 20 years (March 26, 1985-August 19, 2006) in a cohort of 5115 men and women in 4 US cities. Mixed linear modeling was used to account for individual age-based trajectories of pulmonary function and other covariates including tobacco use, which was analyzed in parallel as a positive control. Lifetime exposure to marijuana joints was expressed in joint-years, with 1 joint-year of exposure equivalent to smoking 365 joints or filled pipe bowls.
Main Outcome Measures Forced expiratory volume in the first second of expiration (FEV1) and forced vital capacity (FVC).
Results Marijuana exposure was nearly as common as tobacco exposure but was mostly light (median, 2-3 episodes per month). Tobacco exposure, both current and lifetime, was linearly associated with lower FEV1 and FVC. In contrast, the association between marijuana exposure and pulmonary function was nonlinear (P < .001): at low levels of exposure, FEV1 increased by 13 mL/joint-year (95% CI, 6.4 to 20; P < .001) and FVC by 20 mL/joint-year (95% CI, 12 to 27; P < .001), but at higher levels of exposure, these associations leveled or even reversed. The slope for FEV1 was −2.2 mL/joint-year (95% CI, −4.6 to 0.3; P = .08) at more than 10 joint-years and −3.2 mL per marijuana smoking episode/mo (95% CI, −5.8 to −0.6; P = .02) at more than 20 episodes/mo. With very heavy marijuana use, the net association with FEV1 was not significantly different from baseline, and the net association with FVC remained significantly greater than baseline (eg, at 20 joint-years, 76 mL [95% CI, 34 to 117]; P < .001).
Conclusion Occasional and low cumulative marijuana use was not associated with adverse effects on pulmonary function.
Exposure to tobacco smoke causes lung damage with clinical consequences that include respiratory symptoms, chronic obstructive pulmonary disease, and lung cancer.1,2 Chronic obstructive pulmonary disease and lung cancer are leading causes of death,2,3 and smoking tobacco cigarettes is the most important preventable cause of death in the United States.4,5
Marijuana smoke contains many of the same constituents as tobacco smoke,6 but it is unclear whether smoking marijuana causes pulmonary damage similar to that caused by tobacco. Prior studies of marijuana smokers have demonstrated consistent evidence of airway mucosal injury and inflammation7- 9 as well as increased respiratory symptoms such as cough, phlegm production, and wheeze, similar to that seen in tobacco smokers.10- 12 However, analyses of pulmonary function and lung disease have failed to detect clear adverse effects of marijuana use on pulmonary function.10- 13 It is possible that cumulative damage to the lungs from years of marijuana use could be masked by short-term effects; prior analyses have not attempted to disentangle these factors. Smoking marijuana is increasingly common in the United States,14 and understanding whether it causes lasting damage to lung function has important implications for public health messaging and medical use of marijuana.15,16
The Coronary Artery Risk Development in Young Adults (CARDIA) study collected repeated measures of tobacco and marijuana smoking as well as pulmonary function over the course of 20 years (March 26, 1985-August 19, 2006) in more than 5000 study participants. We estimated both current intensity and lifetime cumulative exposure to tobacco and marijuana smoking and analyzed their associations with spirometric measures of pulmonary function over the 20 years of follow-up.
CARDIA is a longitudinal study designed to measure risk factors for coronary artery disease in a cohort of black and white women and men (n = 5115) aged 18 through 30 years and healthy at enrollment in 1985.17,18 Participants were sampled from 4 US communities without selection for smoking behaviors and comprise a broad cross-section of typical tobacco and marijuana use patterns.
With the written informed consent of participants and the approval of institutional review boards at each study center (Oakland, Chicago, Minneapolis, and Birmingham), participants underwent a baseline examination and 6 follow-up examinations, with 69% retention at year 20. Pulmonary function testing was performed at years 0, 2, 5, 10, and 20. For this investigation, we included all visits for which pulmonary function, smoking behavior, secondhand smoke exposure, height, and waist circumference were available.
Current intensity of tobacco use (cigarettes smoked per day) was assessed at each examination. These data, along with baseline examination data on past years of smoking, were used to estimate cumulative lifetime exposure to cigarettes in terms of pack-years, with 1 pack-year of exposure equivalent to 7300 cigarettes (1 year × 365 days/y × 1 pack/d × 20 cigarettes/pack). Misclassification of smoking exposure by self-report, measured by comparisons with serum cotinine levels, is uncommon.19
Current intensity of marijuana use (episodes in the last 30 days) was also assessed at each examination. Using baseline examination data on past lifetime exposure to marijuana, current intensity of marijuana use, and another question designed to assess number of joints or filled pipe bowls smoked per episode (eMethods), we calculated total lifetime exposure to marijuana joints in joint-years, with 1 joint-year of exposure equivalent to 365 joints or filled pipe bowls smoked (1 year × 365 days/y × 1 joint/d), as described previously.20
Study outcomes were forced expiratory volume in the first second of expiration (FEV1) and forced vital capacity (FVC) measured by forced spirometry. These were collected using a Collins Survey 8-L water-sealed spirometer and an Eagle II microprocessor (years 0, 2, 5, and 10) and then an OMI rolling seal spirometer (year 20). A comparability study performed among 25 participants demonstrated an average difference of less than 1% for both measurements. Standard quality control and testing procedures were maintained according to established guidelines.21,22
CARDIA was designed to recruit approximately equal numbers of self-identified “black, not Hispanic” and “white, not Hispanic” men and women to ensure an adequate sample of the largest minority group in the United States at that time. Height and waist circumference were measured at each examination. As a proxy for socioeconomic status, we used the maximum educational grade attained for each participant. Secondhand smoke exposure in hours per week (sum of exposure in the home, small enclosed spaces, and large spaces) was assessed at each examination, with linear interpolation for missing data. Asthma was self-reported at each examination; we used the baseline assessment. We obtained average annual city-specific levels of airborne particulate matter less than 10 microns and less than 2.5 microns in size23 around the 4 CARDIA study centers from the Environmental Protection Agency24 (eMethods).
Participants were categorized by whether they ever reported current use of tobacco, marijuana, or both at a CARDIA examination and compared across these categories using descriptive statistics. We then categorized participants according to degree of current and lifetime tobacco and marijuana exposure at each examination and described pulmonary function (FEV1 and FVC) across categories before and after adjustment. Tests of trend and interaction were performed in fully adjusted models.
The categorized exposure models described above represent a standard approach to multivariable-adjusted association testing. Categorization models, however, use necessarily arbitrary category thresholds and do not take full advantage of the continuous exposure measurements for estimation or adjustment purposes. To fully explore and test potential nonlinear associations, we modeled tobacco and marijuana exposure variables as flexible cubic splines (eMethods) in adjusted models to allow associations with pulmonary function to take different shapes at lower vs higher levels of exposure.25
For each adjusted analysis described above, we used mixed models accounting for repeated measures of pulmonary function within participants, with a random intercept and a random 3-knot age spline within each individual and an unstructured variance-covariance structure. Fully adjusted models included fixed effects for year, center, and center-year (their interaction), race-sex category, education, and asthma; cubic splines for age, height, waist circumference, secondhand smoke exposure, and exposure to airborne particulate matter less than 10 microns and less than 2.5 microns in size; and interactions between the age-spline variables and race-sex, asthma, waist-spline variables, and height-spline variables to allow for differing flexible age-based trajectories of pulmonary function for participants with differing characteristics. Models were queried to produce adjusted estimates of slope (reflecting the incremental difference in pulmonary function observed with additional tobacco or marijuana smoking) and net association (reflecting the net observed difference between persons with a particular level of consumption and persons with none) at various points along the association curve. All analyses were performed using Stata version 11 and used 2-sided tests for significance at the .05 level, with 95% CIs.
The 5115 CARDIA participants recruited in 1985-1986 contributed 20 777 total visits that included pulmonary function testing. Of these, 959 visits were excluded for lack of complete information on smoking behavior, 114 for lack of height or waist measurements, and 1 for an unknown visit date, leaving 19 703 visits (95%) with complete data from 5016 participants (98%). Participants contributed 3.9 visits/participant on average; attrition was more common in tobacco smokers but not associated with marijuana use. FEV1 and FVC varied across participants, increased slightly with age through the late 20s, and declined slowly thereafter (Figure 1).
More than half of participants (54%; mean age at baseline, 25 years) reported current marijuana smoking, tobacco smoking, or both at 1 or more examinations (Table 1). Smoking patterns differed by race and sex, with black women most likely to smoke tobacco only, white men most likely to smoke marijuana only, and black men most likely to smoke both. Tobacco smokers tended to have lower education and income and to be slightly shorter and less active, whereas marijuana smokers tended to be taller and more active. The median intensity of tobacco use in tobacco smokers was substantially higher (8-9 cigarettes/d) than the median intensity of marijuana use in marijuana smokers (2-3 episodes in the last 30 days). Although marijuana and tobacco exposures were strongly correlated, our sample included 91 participants with no tobacco exposure and more than 10 joint-years of marijuana exposure (contributing 153 observations of pulmonary function), 40 (56 observations) of whom had more than 20 joint-years of exposure.
In fully adjusted models that considered 4-level categorizations of current and lifetime exposure to tobacco and marijuana, tobacco smoking (both current and lifetime) was associated with a lower FEV1 and current smoking with a lower FVC (Table 2). For example, compared with zero exposure, FEV1 was 63 mL lower (95% CI, −89 to −36; P < .001 for trend) and FVC was 69 mL lower (95% CI, −97 to −41; P < .001 for trend) with current tobacco exposure of more than 20 cigarettes per day and 101 mL lower (95% CI, −136 to −65; P < .001 for trend) with lifetime tobacco exposure of more than 20 pack-years.
In contrast, exposure to marijuana (both current and lifetime) was associated with higher FVC and lifetime exposure with higher FEV1. For example, compared with zero exposure, FVC increased with greater lifetime exposure in joint-years (P = .01 for trend) and FEV1 increased with greater lifetime exposure of up to 10 joint-years and then declined to 36 mL (95% CI, −6.5 to 79) greater than the zero exposure level (P = .049 for trend). FVC increased with smoking intensity up to 20 marijuana smoking episodes in the past 30 days and then declined to 20 mL greater than the zero exposure level (P = .03 for trend). We found no statistically significant interactions between tobacco and marijuana exposure for either FEV1 or FVC.
When we modeled current and lifetime tobacco and marijuana exposure as continuous exposures and permitted flexible nonlinear associations (via splines), we again found strong, dose-related associations (P < .001) between increasing exposure to tobacco and lower FEV1 and FVC (Figure 2), with no evidence of nonlinearity (Table 3). Declining slopes ranged as steep as −2.8 mL (95% CI, −4.8 to −0.7; P = .007) per additional cigarette smoked per day and −7.0 mL (95% CI, −10 to −3.7; P < .001) per additional pack-year for FEV1 and were of similar magnitude for FVC (Table 3). At 50 pack-years of exposure, FEV1 was on average 332 mL lower (95% CI, −401 to −263; P < .001) and FVC was 229 mL lower (95% CI, −310 to −147; P < .001), compared with no exposure.
For marijuana, we found strong statistical evidence that associations between marijuana use and pulmonary function were nonlinear (Figure 2, Table 3). At low lifetime exposure levels, increasing marijuana use was associated with a steep increase in both FEV1 (13 mL/joint-year higher [95% CI, 6.4 to 20], P < .001) and FVC (20 mL/joint-year higher [95% CI, 12 to 27], P < .001), but at higher levels of exposure (>7 joint-years), the slope leveled or even turned downward. At more than 10 joint-years of lifetime exposure, we found a nonsignificant decline in FEV1 (−2.2 mL/joint-year [95% CI, −4.8 to 0.3], P = .08) but a significant decline in FEV1 at more than 20 episodes of marijuana use per month (−3.2 mL/episode [95% CI, −5.8 to −0.6], P = .02). Although net associations with FEV1 became negative at very high exposure levels (>40 joint-years or >25 episodes/mo), these negative deflections were not statistically significant (Table 3). FVC remained significantly elevated in even heavy users (eg, 76 mL [95% CI, 34 to 117; P < .001] at 20 joint-years).
In this 20-year study of marijuana and pulmonary function, we confirmed the expected reductions in FEV1 and FVC from tobacco use. In contrast, marijuana use was associated with higher FEV1 and FVC at the low levels of exposure typical for most marijuana users. With up to 7 joint-years of lifetime exposure (eg, 1 joint/d for 7 years or 1 joint/wk for 49 years), we found no evidence that increasing exposure to marijuana adversely affects pulmonary function. This association, however, was nonlinear: at higher exposure levels, we found a leveling off or even a reversal in this association, especially for FEV1. Although our sample contained insufficient numbers of heavy users to confirm a detrimental effect of very heavy marijuana use on pulmonary function, our findings suggest this possibility.
The associations we found between tobacco and pulmonary function are consistent with a large body of prior research on the adverse pulmonary consequences of tobacco smoking. The high prevalence of tobacco smoking, the wide range of exposure intensity among smokers, and the legality of tobacco have made tobacco smoking an easy target for observational epidemiology. Exposure predicts reduced expiratory flow and air trapping, gas-exchange abnormalities, and emphysema,1 and smoking cessation interventions reduce the rate of FEV1 decline in smokers26 (ie, these associations are likely causal). Our findings of a linear dose-response relationship showing lower FEV1 and FVC with increasing tobacco exposure, consistent with prior findings, represent a positive control for our study of the association between marijuana smoking and pulmonary function.
Prior studies of marijuana smoking and pulmonary function have yielded apparently conflicting results.10- 13 Many studies have focused on FEV1:FVC ratio, lower values of which suggest the presence of airway obstruction, and have found either no association10,20,27 or lower FEV1:FVC ratios with marijuana use.28- 32 Lower FEV1:FVC ratios in marijuana smokers, however, can be explained at least partly by a tendency toward higher FVC or total lung capacity.28,29,32 A recent longitudinal study, which demonstrated significantly higher FVC and total lung capacity with marijuana exposure, strongly supports this notion,13,20 as does our study.
The potential association of marijuana smoking with FEV1 has been even less clear. Tobacco smoking reduces FEV1, but despite the similarities in the constituents of marijuana smoke and tobacco smoke and our a priori expectations that marijuana smoking might have similar effects, prior research has not demonstrated this. In studies that report FEV1 in association with marijuana use, findings have mostly been null,20,28,32- 35 although one study reported the apparently paradoxical finding of a lower FEV1 with past marijuana use but a nonsignificantly higher FEV1 with current use.29
Our study suggests a way to reconcile these findings. Because of the many thousands of measurements obtained over 20 years among more than 5000 participants with a wide range of smoking habits, we could simultaneously account for levels of current and past lifetime use of both marijuana and tobacco and test for nonlinearity in their associations with pulmonary function to disentangle short-term and long-term effects. We found highly significant nonlinearity, with a positive association for both FEV1 and FVC at low levels of exposure that reversed in direction toward a possibly negative association for FEV1 at higher levels of exposure (Figure 2 and slopes in Table 3). These findings could explain the paradox previously noted regarding past and current use29 and are also consistent with the average null association reported in studies20,28,32- 35 that either dichotomized marijuana exposure (user/nonuser)28- 31,33,36 or constrained the association to be linear across all levels of exposure.10,20,32,35 When we looked at “marijuana only” smokers (Table 2), we also found a null association with FEV1 and FVC. Only after parsing the association at different levels of exposure, with careful control for confounding, did the suggestion emerge of a negative association for FEV1 at high levels of exposure.
These findings suggest that marijuana smoking could influence pulmonary function via multiple mechanisms. To explain the higher FVC previously observed in marijuana smokers,20,32 some investigators have proposed that the deep inspiratory maneuvers practiced by marijuana smokers could stretch the lungs,13,20 resulting in larger lung volumes.20,32 Another speculative possibility is strengthening of chest wall musculature or another “training” effect that allows marijuana users to inspire more fully (closer to total lung capacity) on spirometry testing. A nondestructive stretch or training effect is consistent with previously reported findings in marijuana smokers of lower lung density32 and a lack of emphysematous change32 or diminished diffusion capacity.20,27,32,36 This mechanism would explain our FVC results and could explain the positive deflection of FEV1. The functional effects of this association on lung health or respiratory function in daily life are unclear.13 An alternate explanation is the acute bronchodilatory effect of marijuana use that has been directly observed in some studies.11 This effect, however, is transient (lasting approximately 60 minutes11) and seems unlikely to explain higher lung volumes measured during the CARDIA examination unless many marijuana users smoked immediately before the examination.
The suggestion of a negative association with FEV1 at higher exposure levels could reflect mixing of this putative stretch/training effect with a second mechanism operating on a different time-exposure scale. A negative association with heavy exposure to marijuana smoke aligns with our a priori hypothesis that marijuana smoking should produce damage to the airways and accelerated loss of lung function similar to that caused by tobacco smoking. Hypothetically speaking, a positive effect from marijuana in the short term (the stretch/training effect) and a negative effect in the long term (damage from smoke exposure) should result in a nonlinear association such as the one we observed. According to this explanation, the predominant effect for FEV1 at very high exposure (more than 40 joint-years) reflects cumulative damage; the predominant effect for FVC at all levels of exposure is from the stretch/training mechanism.
Our study has limitations. Although CARDIA offers longitudinal spirometry measurements, it lacked body plethysmographic measurements of static lung volumes (total lung capacity and residual volume) and measures of diffusing capacity and radiographic emphysema. A minority of our participants reported very high levels of marijuana exposure (and a smaller minority of these were nonsmokers of tobacco), so our estimates at high marijuana exposure levels are imprecise. The self-reported measures of marijuana and tobacco smoking are certain to include recall error, both random and systematic, and do not include any indication of smoking method (joint, pipe, “bong”, etc). It is unlikely, however, that such error would differentially occur in association with pulmonary function, and nondifferential error would most likely bias results toward the null. Our mixed modeling approach is ideal for filtering out random error and taking advantage of individual-level correlations in the data.
As with any observational analysis, unmeasured or inadequately modeled confounding effects could be mixed with our estimates, but the extensive covariate measurements and large sample in our study permitted more extensive efforts to control confounding than were possible in previous studies. This study addressed respiratory exposure to marijuana and not exposure by ingestion. Recent increases in the potency of marijuana are unlikely to have influenced our estimates, because we did not detect an interaction of marijuana and pulmonary function by calendar time.
Marijuana may have beneficial effects on pain control, appetite, mood, and management of other chronic symptoms.15,16 Our findings suggest that occasional use of marijuana for these or other purposes may not be associated with adverse consequences on pulmonary function. It is more difficult to estimate the potential effects of regular heavy use, because this pattern of use is relatively rare in our study sample; however, our findings do suggest an accelerated decline in pulmonary function with heavy use and a resulting need for caution and moderation when marijuana use is considered.
Corresponding Author: Mark J. Pletcher, MD, MPH, Department of Epidemiology and Biostatistics, University of California, San Francisco, 185 Berry St, Ste 5700, San Francisco, CA 94107 (email@example.com).
Author Contributions: Dr Pletcher had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Pletcher, Richman, Safford.
Acquisition of data: Sidney.
Analysis and interpretation of data: Pletcher, Vittinghoff, Kalhan, Richman, Safford, Lin, Kertesz.
Drafting of the manuscript: Pletcher, Safford.
Critical revision of the manuscript for important intellectual content: Pletcher, Vittinghoff, Kalhan, Richman, Safford, Sidney, Lin, Kertesz.
Statistical analysis: Pletcher, Vittinghoff, Richman, Lin.
Obtained funding: Pletcher, Sidney, Kertesz.
Administrative, technical, or material support: Kertesz.
Study supervision: Kertesz.
Conflict of Interest Disclosures: The authors have completed and submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Dr Kalhan reported serving as a consultant for Boehringer-Ingelheim, Forest Laboratories, and AstraZeneca; receiving honoraria for lectures from GlaxoSmithKline and AstraZeneca; receiving honoraria for development of educational materials from Quantia Communications and Medscape Education; and receiving industry-sponsored grants from GlaxoSmithKline and Boehringer-Ingelheim. Dr Kertesz reported chairing a committee that advised the Drug Treatment Task Force for the Chief Justice of the State of Alabama and that he is an employee of the Department of Veterans Affairs. No other authors reported disclosures.
Funding/Support: This study was supported by the National Institute on Drug Abuse (R01-DA-025067) and the National Heart, Lung, and Blood Institute (N01-HC-95095 and N01-HC-48047).
Role of the Sponsors: The National Institute on Drug Abuse funded this analysis but did not participate in CARDIA or review the manuscript. The National Heart, Lung, and Blood Institute helped design CARDIA and funds data collection, supports a Publications and Presentations Committee that reviews and approves all publications, and provides a representative who sits on that committee, but does not otherwise control publication.
Disclaimer: The views expressed in this article do not reflect positions of the Department of Veterans Affairs or of any other entity of the federal government.
Online-Only Material: The Author Video Interview is available here.