Skip to main content
Help us understand how you use our websites. Take part in our 30 minute study now.

Youth smoking and anti-smoking policies in North Dakota: a system dynamics simulation study



The current study utilizes system dynamics to model the determinants of youth smoking and simulate effects of anti-smoking policies in the context of North Dakota, a state with one of the lowest cigarette tax rates in the USA.


An explanatory model was built to replicate historical trends in the youth smoking rate. Three different policies were simulated: 1) an increase in cigarette excise taxes; 2) increased funding for CDC-recommended comprehensive tobacco control programs; and 3) enforcement of increased retailer compliance with age restrictions on cigarette sales.


The explanatory model successfully replicated historical trends in adolescent smoking behavior in North Dakota from 1992 to 2014. The policy model showed that increasing taxes to $2.20 per pack starting in 2015 was the most effective of the three policies, producing a 32.6% reduction in youth smoking rate by 2032. Other policies reduced smoking by a much lesser degree (7.0 and 3.2% for comprehensive tobacco control program funding and retailer compliance, respectively). The effects of each policy were additive.


System dynamics modeling suggests that increasing cigarette excise taxes are particularly effective at reducing adolescent smoking rates. More generally, system dynamics offers an important complement to conventional analysis of observational data.


Various anti-smoking policies and programs have been effective at reducing smoking rates [1,2,3]. However, there is still progress to be made, particularly among adolescents since the majority of smokers initiate in their adolescent years [4, 5]. This is a critical period during which to intervene to reduce the smoking-related burden in the general adult population [2, 6, 7]. Successful strategies for reducing smoking rates include excise taxes on cigarettes [2, 8,9,10,11] and comprehensive tobacco control programs (which include media campaigns, cessation resources, and interventions) [1, 3]. Another common policy is legislating age restrictions on purchasing cigarettes, though the evidence is mixed as to whether age restrictions in isolation are actually effective in decreasing smoking prevalence [12,13,14]. A related issue is retailer compliance [15, 16]; enforcement efforts aimed at retailers are successful at reducing sales to minors, but this does not necessarily translate into reduced smoking [17].

A prominent challenge in evaluating the effectiveness of these policies is the uncertainty in their true causal effect on smoking rates, which stems from the limitations of observational data. Specifically, it is rare that an appropriate control group is available for a given policy; and even in this case, the true effects of the policy itself are difficult to tease apart from other influences of the outcome, such as larger economic trends. Thus, confounding is a pervasive challenge in accurately evaluating anti-smoking policies. A further challenge in existing policy research is the limited ability to account for feedback in complex systems. For example, there is a positive feedback loop between smoking behavior and nicotine dependence: smoking behavior leads to the development of nicotine dependence symptoms [18,19,20] and nicotine dependence in turn increases and perpetuates later smoking behavior [18, 21, 22]. However, conventional regression analyses artificially and inappropriately impose a uni-directional relationship, selecting one as the outcome and the other as a predictor, with only a few notable exceptions that examined the mutual directionality of this relationship [18, 19, 23].

Simulation methods can overcome many of the limitations found in observational data, and thus represent an important complementary tool to study anti-smoking policies. Previous simulation work has been done using discrete Markov models to estimate the reduction in smoking prevalence and smoking-attributable deaths that could result from tax increases and comprehensive tobacco control policies in both Louisiana [8], the Netherlands [24], and Brazil [25]. Additionally, system dynamics modeling has been utilized to model the effects of excise tax policy and smoking cessation interventions in New Zealand [26, 27] as well as anti-tobacco education in schools [28].

System dynamics has a wide range of applications. It decomposes complex systems and can make projections on various effects and consequences of a given decision. Thus, system dynamics is especially useful for studying health policy, agriculture, economics, and climate. Examining the effects of smoking policies among adolescents can be readily applied to system dynamics methodology given that this problem is complex, dynamic, and contains multiple feedback loops. A system dynamics model, which is essentially a formalized set of hypotheses (e.g. about the causes of smoking behavior), can be used to test causal hypotheses and effects within the simulation framework. Thus, system dynamics can overcome many of the methodological limitations of existing studies on smoking policy, provided the model reasonably approximates reality, after calibrating it to the population of interest and performing a series of validation tests (Model Validation). In particular, system dynamics can manipulate a single variable at a time and examine its effect on the system, whereas this is difficult or impossible with statistical analysis of observational data due to the presence of confounders that covary with other variables. In addition, system dynamics is able to account for mutual relationships between variables (i.e. feedback loops), whereas conventional statistical analyses must impose a unidirectional relationship between the variables. Since smoking behavior is driven by a variety of influences and feedback loops, this makes system dynamics a promising tool for studying population-level trends in smoking behavior as well as policies aimed at reducing smoking rates. Very little system dynamics work has been done examining the effects of smoking policies among adolescent populations in particular.

The current study presents a system dynamics model of adolescent smoking behavior and anti-smoking policies. As commonly done in system dynamics, an “explanatory model” was built first with the goal of replicating trends in real data (youth smoking rates from 1992 to 2014); Fig. 1 shows this historical data as well as the simulation data from the explanatory model (described in detail below). Second, a “policy model” was built to evaluate and compare the effects of three separate anti-smoking “policies:” 1) an increase in excise taxes for cigarettes; 2) increased funding for comprehensive tobacco control programs as recommended by the Center for Disease Control (CDC) [3]; and 3) stricter enforcement of retailers’ compliance with age restrictions on cigarette sales. The model was tailored to the case study of North Dakota, USA, which has one of the lowest cigarette tax rates in the USA [29] at $0.44/pack—a rate that has remained unchanged since 1993—and at the time of this research was considering legislation that would have substantially increased that tax rate (North Dakota Tobacco Tax Increase, Initiated Statutory Measure 4). Each of these policies were evaluated with respect to their impact on the reduction in youth smoking rates approximately 20 years into the future.

Fig. 1

Actual vs. simulated data on past-month smoking prevalence. Values are averaged across middle and high school populations. Simulated data are from the “base” stock-and-flow model shown in Fig. 3. Actual data are from the Monitoring the Future (MTF) Study, 1992–2014


Causal loop diagram

A causal loop diagram (CLD) is a conceptual model that represents the formalized hypothesized causal relationships in a system, with an emphasis on relationships that form complete feedback loops. A CLD was constructed to represent well-known risk factors and causes of smoking among youth, namely, parental and peer smoking, secondhand smoke exposure, nicotine dependence, risk perception, and exposure to secondhand smoke. Though there are a multitude of risk factors for smoking, not all of them produce dynamic (as opposed to static or constant) effects; consistent with principles of system dynamics [30]. In other words, a CLD is intended to display the minimal essential set of feedback relationships in the more complex quantitative stock-and-flow simulation model.

Figure 2 shows the CLD that was developed. Five major feedback loops were identified that affect smoking behavior, with 3 being reinforcing/positive feedback loops (in which any change becomes amplified and self-perpetuating) and 2 being balancing/negative feedback loops (in which any change is counteracted and self-limited). Specifically, R1 shows the reinforcing effect of social pressure (from both peers and parents) on youth smoking: the more common smoking is among adolescents, the more susceptible never-smokers become to initiate smoking due to social norms [31], peer pressure [32, 33], behavioral modeling of smoking [34, 35], perceived availability of cigarettes [36, 37], etc.; and the loop is closed when these new smokers are added to the pool of peers who smoke and (eventually) parents who smoke. Similarly, R2 shows the effect of secondhand smoke on nicotine dependence among non-smoking youth, which makes them more likely to initiate smoking [38, 39]; the loop again is closed when these youth contribute to further secondhand smoke in the environment. Finally, R3 shows the positive loop introduced by nicotine dependence: smoking increases nicotine dependence [38, 39], and nicotine dependence perpetuates smoking [18, 19, 23].

Fig. 2

Causal loop diagram of smoking behavior. Curved arrows represent (hypothesized) causal relationships, and the polarity of each relationship is labeled as +/−. Feedback loops are labeled as B (balancing or negative feedback loop) and R (reinforcing or positive feedback loop) and numbered. ND: nicotine dependence

Balancing loops also exist that attenuate the reinforcing effects in the system (Fig. 2). Risk perception counteracts increases in youth smoking (loop B1): such an increase in smoking would result in higher prevalence of tobacco-related death and disease, which increases risk perception and lowers the likelihood that future generations of smokers will initiate. Risk perception also motivates the development and implementation of anti-smoking policies (B2) such as tax increases, which reduces susceptibility of initiation among non-smoking youth.

The contribution of these feedback loops to the dynamics of the simulation model are explained in Results.

Stock-and-flow simulation model

Stocks and flows

A detailed stock-and-flow diagram was constructed based on the CLD, and equations and parameters were entered for all variables. The stock-and-flow diagram thus represents the system of hypothesized causal relationships in the system (i.e. the causal factors that determine smoking behavior). System Dynamics models are broken up into a series of “stocks” (a terminology used in system dynamics that refers to some quantity that can accumulate, e.g. people in different stages of addiction, such as experimental smokers) and “flows” (transition rates between stocks, e.g. smoking initiation rates that describe transitioning from a never-smoker to an experimenter), and is especially appropriate for systems that contain positive or negative feedback relationships and show dynamic (i.e. changing over time) behavior (e.g. declining rates of cigarette use over time) [30]. The stocks represent mutually-exclusive categories of smokers: never-smokers, experimenters, former experimenters, current smokers, and ex-smokers. The distinction between experimenters and smokers was the commonly used criterion of 100 cigarettes/lifetime [40] for characterizing lifetime smoking status; thus, once an individual becomes a smoker, he/she cannot again become an experimenter. The distinction between current and former smoking behavior (for both experimenters and current smokers) was any smoking in the past 30 days.

These stocks were changed through their corresponding inflows and outflows: smoking initiation, smoking progression, cessation, and relapse. The modeling task, therefore, was to hypothesize how the risk factors influence these flows that in turn, change the levels of the stocks. Figure 3 shows a simplified version of the stock-and-flow model; the full model with all structure and equations is available online [41].

Fig. 3

Simplified stock-and-flow diagram of smoking behavior. Experimenters: smoked < 100 cigarettes/lifetime. Smokers: smoked ≥100 cigarettes/lifetime; current smokers/experimenters: smoked within past 30 days; ex-smokers: did not smoke within past 30 days

Model structure

Using variables and feedback relationships from the CLD, a stock-and-flow simulation was then developed to model the initiation and progression of smoking behavior among youth. The crux of the model is an “aging chain” representing different categories of smokers as “stocks” and transitions between these categories as “flows.” Different age groups representing middle school or younger (ages 11–14) and high school (ages 15–18) were modeled in parallel through the use of array variables, to reflect this common age distinction in national surveys such as Monitoring the Future [42].

The model was built and run using Stella (Systems Thinking, Experimental Learning Laboratory with Animation) Architect, version 1.3.1. Stella Architect is produced by Isee Systems, one of the leading developers of system dynamics software. This software is frequently used in educational settings [43] as well as research settings [44]. Stella reaches its results through numeric approximation of a series of differential equations. Specifically, values of stocks are calculated as a function of their initial value (a parameter set by the user) and the integration of net inflows and net outflows, which are in turn determined by their user-specified equations and other variable inputs. Stella uses numeric approximation at each time step (a user-determined fraction of the underlying time unit) and calculates throughout the chosen range of time.

Calibration to real data

The time horizon of the model was set to 1992–2032 (a 40-year range), so that the explanatory power of the model (simulation vs. real data) could be evaluated over the first approximately 20-year period, and projections could be made 20 years into the future.

Model parameters were selected to replicate actual trends in youth smoking behavior in North Dakota, USA over the observed period. North Dakota was selected as a case study because it has one of the lowest tax rates on cigarettes in the US at $0.44. At the time of this research, legislation was being considered to increase the excise tax by $1.00 or by $1.76 (North Dakota Tobacco Tax Increase, Initiated Statutory Measure 4). Parameters were selected using several methods, including using publicly available population-level statistics, using coefficient estimates from regression models in scientific literature, and (in cases where quantitative data were not available) estimating based on qualitative information, according to standard procedures for incorporating “soft variables” in system dynamics [30]. Parameters were validated based on the tests described below. All equations and parameters are available in the model which is available online [41].

Model validation

A series of model validation tests were conducted identify and eliminate errors in the model and increase the validity of the results. Model testing is conducted throughout the entire model building process. This process exposes limitations or errors in the model and allows for continuous improvements. Using a wide variety of model validation tests, rather than a single validity test, helps to ensure the validity of the model’s output. The current study primarily utilized three separate types of validation tests: boundary adequacy tests, structure tests, and stock-and-flow assessment tests. Some validation tests are based on judgment of the modeler (e.g. boundary of the model), while others are based on software output (e.g. match between simulation data and real data; unit errors).

First, boundary adequacy tests were performed to make decisions about inclusion/exclusion of many risk factors. Only variables that were directly relevant to youth smoking were considered. Variables were generally included if they were hypothesized to 1) causally affect other variables in the model, 2) be endogenous to the model (i.e. are affected by other variables contained within feedback loops in the model) and 2) contribute dynamically to the model. These inclusion criteria, for example, excluded variables such as demographic characteristics because 1) they may represent associations rather than causal relationships, 2) they cannot be affected by other variables in the model, and 3) despite being associated with an individual’s risk for smoking, they do not produce dynamic effects over time at the population level.

Next, a series of structure tests were performed. Equations in the model were motivated by, and compared with, logic and findings from scientific literature. For example, monthly smoking quantity was defined as smoking frequency (number of days smoked in 1 month) times cigarettes per day. Parameters were also assessed using sensitivity testing, which is especially important for qualitative or “soft” variables. This was done by re-running many simulations, each time with a different value of a given parameter; thus, the full range of plausible parameter values are tested, and the overall sensitivity of the model (i.e. outcomes of youth smoking rates) to changes in specific parameters was evaluated. All models will display some numeric sensitivity (i.e. at least a small change in the overall level of the outcome) but a more serious issue is sensitivity of the behavior mode itself (i.e. the shape of the outcome’s behavior over time) [30]. Therefore, this test determines whether the model can produce a consistent mode of behavior, regardless of the stochastic variability which results from variation in parameters.

Finally, the stock-and-flow model was assessed with respect to whether it could replicate trends in actual data drawn from national and state survey data, namely the prevalence of past-month smoking from 1992 to 2015 according to the Monitoring the Future (MTF) survey [42]. In SD, a good match between simulated and observed data is evidence for the appropriateness of the model; for this purpose, replicating the overall dynamics (i.e. the shape of the curve of smoking prevalence over time) is more important than achieving an exact numeric match. Our explanatory model was able to replicate the observed data well, as indicated by the consistency in periods of increase, periods of decrease, and timing the peak of past-month smoking prevalence (Fig. 1). The key feedback loops contributing to this behavior are explained in Results.

Anti-smoking policy tests

Once the main (explanatory) model was completed, three policies were tested: 1) an increase in the excise tax (from $0.44 to either $1.44 or $2.20, representing the actual and two proposed excise tax rates in North Dakota at the time of this research); 2) an increase in funding for comprehensive tobacco control programs (from $9.80 per capita—the CDC recommended level for North Dakota [45]—to $20 per capita, chosen somewhat arbitrarily as roughly double the 2014 value), which was modeled as affecting anti-tobacco media and promotions from the tobacco industry; and 3) improving compliance of retailers regarding sales to minors (from 78% compliance to high school adolescents and 73% compliance to middle school adolescents [46], both to 95% compliance).

Policies were added to the model as follows. Increases in the excise tax were modeled as reducing the average smoking quantity which in turn reduces initiation and progression rates, and separately as increasing cessation rates. Increased funding for comprehensive tobacco control programs was modeled as lowering tobacco marketing and as increasing perception of health risks. Improved retailer compliance was modeled as reducing perceived cigarette availability, which in turn reduces initiation and progression rates. Policies were implemented in simulation year 2015, and results were simulated through 2032 (allowing approximately 20 years of actual, historical data and 20 years of simulated projections). Each policy was tested in isolation in different simulation runs, and a final simulation run tested all policies together. Policies were evaluated with respect to the projected smoking prevalence resulting from each simulation run, in terms of the percent change in smoking prevalence among adolescents (averaged across middle and high school students). As mentioned in Model Validation, parameter uncertainty produces some uncertainty in the projected smoking prevalence, so the relative ranking of effectiveness across policies is of more use than the exact numerical projections, provided that the same ranking is preserved across sensitivity testing of different plausible parameter values.


Causal loop diagram

The CLD presented in Fig. 2 contains three major positive feedback loops (in which a change in smoking prevalence perpetuates itself via social pressure to smoke, exposure to secondhand smoke, and nicotine dependence) and two major negative feedback loops (in which a change smoking prevalence limits or counteracts itself through awareness of health risks and through enactment of anti-tobacco policies). In isolation, each of the positive feedback loops would produce exponential growth in smoking prevalence, and each of the negative feedback loops would produce exponential decay or goal-seeking behavior. In the full model, these positive and negative feedback loops interact in important ways to produce the simulated behavior that replicates the actual historical data (Fig. 1). The initial rise in smoking prevalence from 1992 to approximately 1995 is produced by the positive feedback in the system – namely the increase in parental and peer smoking and exposure to secondhand smoke, and their reinforcing effects on adolescent smoking rates via increased susceptibility to smoking and smoking initiation rates. However, negative feedback relationships dominate the system starting in the mid-1990’s, as smoking rates start to decelerate, then peak, then decline until the present. Specifically, the Tobacco Master Settlement Agreement in 1998 which banned cigarette advertisement to minors, alleviated the social pressure loop, which resulted in lower initiation rates and reduced susceptibility to smoking.

Simulation model

The explanatory simulation model (Fig. 3) contains the essential feedback loops presented in the CLD and described above, along with additional detail such as modules for tobacco marketing (which affects social pressure and price) and other tobacco products (whose effects are limited in this model to nicotine dependence). By and large, the feedback loops produce the behavior observed in Fig. 1 as explained in the CLD results section; and the stock-and-flow simulation model offers additional insight into how this system produces the observed behavior. That is, modules in the explanatory simulation model were turned “on” and “off” as part of the process of calibrating the simulated data with the real data. Specifically, the negative feedback loops due to heightened perception of health risks of smoking and due to the marketing restrictions from the Master Settlement Agreement were necessary for the model to better match the observed decline in smoking prevalence after the year 2000. Additionally, the model behavior was affected to a large degree by price effects, i.e. how cigarette prices affect smoking rates; adjusting cigarette price to the average income in North Dakota for each year of the simulation resulted in a closer match of the simulated data to the observed data. This highlights the importance of price effects to the system as a whole.

The full simulation model was subjected to the validity tests described in Methods. Sensitivity tests showed that the model was most sensitive to changes in the effect of social pressure on the initiation rate (producing at most a 25% change in the overall level of the youth smoking rate), and less sensitive to the exact value for other parameters. However, the shape of the behavior over simulated time did not change in any case, indicating that the model is fairly robust to changes in uncertain parameters.

Policy tests

Figure 4 shows the results of the policy simulations if they had been implemented in 2015, and projects the current (past-30 day) smoking rate until year 2032 (nearly 20 years after implementation). Though even the base case (no policy changes) shows a decline in the fraction of adolescents who currently smoke cigarettes, each policy further decreases current smoking prevalence. Increasing the excise tax from $0.44 (base case) to $1.44 and $2.20 decreased the simulated prevalence of past-month smoking by 22.3 and 32.6%, respectively, in 2032, relative to the 2032 base case projections (no policy). In contrast, increasing funding for comprehensive tobacco control programs from $9.8 per capita (base case) to $20 per capita decreased the simulated smoking prevalence by 7.0% of the 2032 base case prevalence, and increasing retailer compliance with sales restrictions to 95% (from the base case of 78% for middle school students and 73% for high school students) decreased the simulated smoking prevalence by 3.2% of the base case prevalence in 2032. Finally, all three policies implemented together (with the higher $2.20/pack excise tax) decreased the simulated smoking prevalence by 42.7% in 2032, consistent with an additive effect of each policy.

Fig. 4

Policy simulations for policy tests for: increases in the excise tax per pack of cigarettes (a); increasing per-capita funding for comprehensive tobacco control programs (b); increasing retailer compliance with sales laws (c); policies A – C implemented together (d). Line type indicates the policy, with 1 (solid line) representing the base case simulation in all panels. Y-axis shows the current (past-30 day) smoking rate, averaged across middle school and high school students


The current study used system dynamics simulation modeling to examine the possible impacts of three different anti-smoking policies on the rate of current smoking among youth in North Dakota. While all three policies were effective in reducing smoking rates nearly 20 years later according to this model, increasing excise taxes on cigarettes produced the greatest reduction in rates of current youth smoking.

These findings are consistent with previous findings from conventional policy analysis. In particular, there is strong evidence that increasing excise taxes is one of the most effective ways to reduce the youth smoking rate [2, 8, 9, 24]. In comparison, our findings suggest that comprehensive tobacco control programs—which include media campaigns, interventions, and cessation resources—also lower the youth smoking rate, but to a much lesser degree. This aligns with previous research showing that anti-tobacco media has a comparatively weaker effect on reducing smoking [1, 8]; an effect which has not consistently been found in all settings [1, 47]. The current finding that policies have additive effects also corroborate previous findings that comprehensive tobacco control is more effective when implemented along with tax increases [1, 8]. Finally, our findings that increasing retailer compliance produced a minimal decrease in smoking rates is consistent with previous research that increasing compliance has an uncertain effect on actual youth smoking behavior, despite being successful at decreasing sales to youth [16, 17].

The current study advances these policy findings in important ways, namely through the use of system dynamics modeling. In particular, these findings indicate that the anti-smoking policies considered—increasing the excise tax, increasing funding for comprehensive tobacco prevention programs, and increasing retailer compliance with sales laws—can have plausible causal effects on reducing the youth smoking rate. These findings are consistent with previous system dynamics and other simulation studies in that tax policy and other comprehensive tobacco policies are successful at lowering the smoking rate among the general population [8, 24, 26]. Here, we extend these findings to adolescents in particular, among whom tax policy is particularly effective at lowering smoking rates.

System dynamics methodology offers a number of advantages over existing research. True tests of causality are possible in simulations, and conclusions can only be drawn about reality to the extent that the model is an appropriate simplification of the real-world system. Continual model testing during the building process of the model helps build confidence in the model’s ability to simulate the essential components of a real-world system. In this way, simulation studies offer an important complement to conventional statistical analyses that are limited by observational data, making causality difficult or impossible to establish. Moreover, though the current model was calibrated to the case study of North Dakota, a system dynamics model can be tailored to other situations merely by changing the relevant parameters in the model. Thus, a single model, once constructed, can be easily adapted and re-used to study a wide variety of other hypothetical and actual situations, increasing its utility and impact. Additionally, simulation provides results much more quickly and inexpensively than in real life. Not only can this increase the pace of research, but this also could facilitate learning in a rapid, interactive context. For example, legislators could interact with the model in ways that allow them to learn the short-term and long-term effects of a given anti-smoking policy within the system; such simulation-based instructional tools have been very effective in other applications such as management [30, 48].

Limitations of this study must be acknowledged. The major limitation of system dynamics modeling is that the generalization of findings depends on the model appropriately capturing the real-world system. Since all models are by definition a simplification of the real world, it is difficult to evaluate the ways and extent to which the model must reflect reality (e.g. with respect to the variables included, or the desired match between real and simulated data). However, this model was subject to a variety of validation tests that are standard in system dynamics [30], and was found to be fairly robust. Thus, though some small deviations between simulated and real data are to be expected, the overall shape and policy conclusions are likely to be accurate. Related to this main limitation, the model exhibited some sensitivity to certain parameters for which finding precise objective data was difficult (e.g. the effect of social pressure on smoking initiation altered the overall smoking rates by up to 25%). Though the overall behavior of the model was not sensitive to these parameters, future versions of this model would be improved if more precise estimates for these parameters could be obtained. Additionally, the cost-effectiveness of each policy is also a critical consideration given the different costs that would be required for implementation, but that is outside the scope of this study. Finally, it is unclear how well the current model’s findings would generalize to other environments; testing these policies in a different situation or setting would require re-calibrating the parameters of the model appropriately.


The results of this system dynamics study suggest that increasing taxes on cigarettes are an especially effective way to reduce the youth smoking rate, consistent with previous research [8, 9, 24]. These simulation results offer an important complement to conventional statistical approaches, in that they can conclusively establish (hypothesized) causal effects.

Availability of data and materials

The system dynamics model created during this study is available in full online at: The datasets used to calibrate the model are publicly available from the Monitoring the Future Study, waves 1992–2013,



Center for Disease Control


Causal loop diagram


Monitoring the Future


  1. 1.

    Wakefield M, Chaloupka F. Effectiveness of comprehensive tobacco control programmes in reducing teenage smoking in the USA. Tob Control. 2000;9(2):177–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  2. 2.

    Aldrich MC, Hidalgo B, Widome R, Briss P, Brownson RC, Teutsch SM. The role of epidemiology in evidence-based policy making: a case study of tobacco use in youth. Ann Epidemiol. 2015;25(5):360–5.

    PubMed  Article  Google Scholar 

  3. 3.

    Centers for Disease Control and Prevention. Best Practices for Comprehensive Tobacco Control Programs—2014. Atlanta: U.S. Department of Health and Human Services, Centers for Disease Control and Prevention, National Center for chronic disease and prevention, Office on smoking and health; 2014.

  4. 4.

    Anderson CM, Burns DM, Dodd KW, Feuer EJ. Chapter 2: birth-cohort-specific estimates of smoking behaviors for the U.S population. Risk Anal. 2012;32(Suppl 1):S14–24.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Rath JM, Villanti AC, Abrams DB, Vallone DM. Patterns of tobacco use and dual use in US young adults: the missing link between youth prevention and adult cessation. J Environ Public Health. 2012;2012:679134.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    U.S. Department of health and human services. Preventing Tobacco Use among Youth and Young Adults: A Report of the Surgeon General. Atlanta, GA; 2012.

  7. 7.

    National Institute on Drug Abuse, National Institutes of Health, U.S. Department of Health and Human Services. Preventing drug abuse among children and adolescents: a research-based guide. Washington, DC; 2003.

  8. 8.

    Levy D, Fergus C, Rudov L, McCormick-Ricket I, Carton T. Tobacco policies in Louisiana: recommendations for future tobacco control investment from SimSmoke, a policy simulation model. Prev Sci. 2016;17(2):199–207.

    PubMed  Article  Google Scholar 

  9. 9.

    Bader P, Boisclair D, Ferrence R. Effects of tobacco taxation and pricing on smoking behavior in high risk populations: a knowledge synthesis. Int J Environ Res Public Health. 2011;8(11):4118.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Wilson N, Thomson G, Tobias M, Blakely T. How much downside? Quantifying the relative harm from tobacco taxation. J Epidemiol Community Health. 2004;58(6):451–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Bader P, Boisclair D, Ferrence R. Effects of tobacco taxation and pricing on smoking behavior in high risk populations: a knowledge synthesis. Int J Environ Res Public Health. 2011;8(11):4118–39.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Kuipers MA, Brandhof SD, Monshouwer K, Stronks K, Kunst AE. Impact of laws restricting the sale of tobacco to minors on adolescent smoking and perceived obtainability of cigarettes: an intervention-control pre-post study of 19 European Union countries. Addiction. 2017;112(2):320–9.

    PubMed  Article  Google Scholar 

  13. 13.

    Kessel Schneider S, Buka SL, Dash K, Winickoff JP, O'Donnell L. Community reductions in youth smoking after raising the minimum tobacco sales age to 21. Tob Control. 2016;25(3):355–9.

    PubMed  Article  Google Scholar 

  14. 14.

    Friend KB, Lipperman-Kreda S, Grube JW. The impact of local U.S. tobacco policies on youth tobacco use: a critical review. Open J Prev Med. 2011;1(2):34–43.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    DiFranza JR. State and federal compliance with the Synar amendment: federal fiscal year 1998. Arch Pediatr Adolesc Med. 2001;155(5):572–8.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Forster JL, Wolfson M. Youth access to tobacco: policies and politics. Annu Rev Public Health. 1998;19(1):203–35.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Rigotti NA, DiFranza JR, Chang Y, Tisdale T, Kemp B, Singer DE. The effect of enforcing tobacco-sales laws on adolescents' access to tobacco and smoking behavior. N Engl J Med. 1997;337(15):1044–51.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Doubeni CA, Reed G, Difranza JR. Early course of nicotine dependence in adolescent smokers. Pediatr. 2010;125(6):1127–33.

    Article  Google Scholar 

  19. 19.

    Hu MC, Griesler PC, Wall MM, Kandel DB. Reciprocal associations between cigarette consumption and DSM-IV nicotine dependence criteria in adolescent smokers. Addiction. 2014;109(9):1518–28.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Kandel DB, Hu MC, Griesler PC, Schaffran C. On the development of nicotine dependence in adolescence. Drug Alcohol Depend. 2007;91(1):26–39.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Dierker L, Hedeker D, Rose J, Selya A, Mermelstein R. Early emerging nicotine dependence symptoms in adolescence predict daily smoking in young adulthood. Drug Alcohol Depend. 2015;151:267–71.

    PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    DiFranza JR, Savageau JA, Rigotti NA, Fletcher K, Ockene JK, McNeill AD, et al. Development of symptoms of tobacco dependence in youths: 30 month follow up data from the DANDY study. Tob Control. 2002;11(3):228–35.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Selya AS, Rose JS, Dierker L, Hedeker D, Mermelstein RJ. Evaluating the mutual pathways among electronic cigarette use, conventional smoking and nicotine dependence. Addiction. 2018;113(2):325–33.

    PubMed  Article  Google Scholar 

  24. 24.

    Nagelhout GE, Levy DT, Blackman K, Currie L, Clancy L, Willemsen MC. The effect of tobacco control policies on smoking prevalence and smoking-attributable deaths. Findings from the Netherlands SimSmoke tobacco control policy simulation model. Addiction. 2012;107(2):407–16.

    PubMed  Article  Google Scholar 

  25. 25.

    Levy D, de Almeida LM, Szklo A. The Brazil SimSmoke policy simulation model: the effect of strong tobacco control policies on smoking prevalence and smoking-attributable deaths in a middle income nation. PLoS Med. 2012;9(11):e1001336.

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Cavana RY, Clifford LV. Demonstrating the utility of system dynamics for public policy analysis in New Zealand: the case of excise tax policy on tobacco. Syst Dyn Rev. 2006;22(4):321–48.

    Article  Google Scholar 

  27. 27.

    Tobias MI, Cavana RY, Bloomfield A. Application of a system dynamics model to inform Investment in Smoking Cessation Services in New Zealand. Am J Public Health. 2010;100(7):1274–81.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Tengs TO, Osgood ND, Chen LL. The cost-effectiveness of intensive national school-based anti-tobacco education: results from the tobacco policy model. Prev Med. 2001;33(6):558–70.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Rauschenberger R. State and local taxes: an overview and comparative guide. Bismarck, ND: North Dakota Office of State Tax Commissioner; 2016.

  30. 30.

    Sterman J. Business dynamics: systems thinking and modeling for a complex world. New York: The McGraw-Hill Companies, Inc.; 2000.

    Google Scholar 

  31. 31.

    Eisenberg ME, Forster JL. Adolescent smoking behavior: measures of social norms. Am J Prev Med. 2003;25(2):122–8.

    PubMed  Article  Google Scholar 

  32. 32.

    Avenevoli S, Merikangas KR. Familial influences on adolescent smoking. Addiction. 2003;98(Suppl 1):1–20.

    PubMed  Article  Google Scholar 

  33. 33.

    Vink JM, Willemsen G, Boomsma DI. The association of current smoking behavior with the smoking behavior of parents, siblings, friends and spouses. Addiction. 2003;98(7):923–31.

    PubMed  Article  Google Scholar 

  34. 34.

    Kalesan B, Stine J, Alberg AJ. The joint influence of parental modeling and positive parental concern on cigarette smoking in middle and high school students. J Sch Health. 2006;76(8):402–7 quiz 38-9.

    PubMed  Article  Google Scholar 

  35. 35.

    Hoving C, Reubsaet A, de Vries H. Predictors of smoking stage transitions for adolescent boys and girls. Prev Med. 2007;44(6):485–9.

    PubMed  Article  Google Scholar 

  36. 36.

    Doubeni CA, Li W, Fouayzi H, DiFranza JR. Perceived accessibility of cigarettes among youth: a prospective cohort study. Am J Prev Med. 2009;36(3):239–42.

    PubMed  Article  Google Scholar 

  37. 37.

    Woodruff SI, Candelaria JI, Laniado-Laborin R, Sallis JF, Villasenor A. Availability of cigarettes as a risk factor for trial smoking in adolescents. Am J Health Behav. 2003;27(1):84–8.

    PubMed  Article  Google Scholar 

  38. 38.

    Belanger M, O'Loughlin J, Okoli CT, McGrath JJ, Setia M, Guyon L, et al. Nicotine dependence symptoms among young never-smokers exposed to secondhand tobacco smoke. Addict Behav. 2008;33(12):1557–63.

    PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Okoli CT, Browning S, Rayens MK, Hahn EJ. Secondhand tobacco smoke exposure, nicotine dependence, and smoking cessation. Public health Nurs (Boston, Mass). 2008;25(1):46–56.

    Article  Google Scholar 

  40. 40.

    Centers for Disease Control. Cigarette smoking among adults- United States, 2000. 2002.

  41. 41.

    Ivanov O. , Bachman A., Selya A., Wheat D. Youth -smoking and anti-smoking in North Dakota. 2019. Accessed 22 July 2019.

  42. 42.

    Miech R, Johnston L, O'Malley P, Bachman J, Schulenberg J. Monitoring the future national survey results on drug use, 1975–2015: volume I, secondary school students. Ann Arbor: Institute for Social Research, The University of Michigan; 2016.

    Google Scholar 

  43. 43.

    Doerr HM. Stella ten years later, a review of the literature. Int J Comput Math Learn. 1996;1(2):201–24.

    Article  Google Scholar 

  44. 44.

    Gambardella PK, Polk DE, Lounsbury DW, Levine RL. A co-flow structure for goal directed internal change. Syst Dyn Rev. 2017;33(1):34–8.

    PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Centers for Disease Control and Prevention (CDC). CDC Best Practices for Comprehensive Tobacco Control Programs. 2014.

  46. 46.

    Centers for Disease Control and Prevention. Preventing tobacco use among young people: A report of the surgeon general. 1994.

  47. 47.

    Niederdeppe J, Farrelly MC, Wenter D. Media advocacy, tobacco control policy change and teen smoking in Florida. Tob Control. 2007;16(1):47–52.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Bakken B, Gould J, Kim D. Experimentation in learning organizations: a management flight simulator approach. Eur J Oper Res. 1992;59(1):167–82.

    Article  Google Scholar 

Download references


Not applicable.

Consent for publications

Not applicable.


This work was supported by the National Institute on Drug Abuse (grant numbers L40 DA042431, 1F32DA033742-01A1), the National Institute for General Medical Sciences (grant number 1P20GM121341–01), and the Norwegian Centre for International Cooperation in Education (SIU; grant numbers PNA-2014/10088 and NNA-2016/10023). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH, NIDA, NIGMS, or SIU. These funding bodies had no role in the study design, data analysis, interpretation of data, or the writing of the manuscript.

Author information




AS conceptualized the study, interpreted findings, and drafted the manuscript. OI constructed and conceptualized the initial system dynamics model, interpreted findings, and drafted parts of the manuscript. AB modified and ran the system dynamics model, interpreted findings, and drafted parts of the manuscript. DW interpreted findings and drafted parts of the manuscript. All authors revised the manuscript and approved the final version for publication.

Corresponding author

Correspondence to Arielle S. Selya.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Selya, A.S., Ivanov, O., Bachman, A. et al. Youth smoking and anti-smoking policies in North Dakota: a system dynamics simulation study. Subst Abuse Treat Prev Policy 14, 34 (2019).

Download citation


  • Adolescent
  • Smoking
  • Simulation methods
  • System dynamics
  • Tobacco policy