Skip to main content


Int J Public Health, 01 April 2021

What Regional Living Conditions Affect Individual Smoking of Adults in Russia

Sergey A. Maksimov
Sergey A. Maksimov*Svetlana A. ShalnovaSvetlana A. ShalnovaYulia A. BalanovaYulia A. BalanovaVladimir A. KutsenkoVladimir A. KutsenkoSvetlana E. EvstifeevaSvetlana E. EvstifeevaAsiia E. ImaevaAsiia E. ImaevaOksana M. DrapkinaOksana M. Drapkina
  • Department of Epidemiology of Chronic Non-Communicable Diseases, National Medical Research Center for Therapy and Preventive Medicine of the Ministry of Healthcare of the Russian Federation, Moscow, Russia

Objectives: Our study evaluated the impact of a wide range of characteristics of large administrative regions on the individual level of cigarette smoking in the Russian adult population.

Methods: The pool of participants included 20,303 individuals aged 25–64 years. We applied 64 characteristics of the 12 Russian regions under study for 2010–2014. Using principal component analysis, we deduced five evidence-based composite indices of the regions. We applied the generalized estimating equation to determine associations between the regional indices and the individual level of smoking.

Results: The increased Industrial index in the region is associated with the probability of smoking (odds ratio = 1.15; 95% confidence interval = 1.06–1.24). The other indices show associations with smoking only in separate gender and educational groups. Surprisingly, it was found that the Economic index has no associations with the probability of smoking.

Conclusion: We evaluated the key associations of the territorial indices with the individual probability of smoking, as well as the mutual influence between the territorial indices and individual factors.


Smoking causes the death of millions of people every year. Specifically, in 2017, about eight million people worldwide died of tobacco-related diseases (1). Numerous international prospective studies and meta-analyses attest to the significant role of smoking in the development of a wide range of serious chronic non-communicable diseases: Coronary heart disease and brain hemorrhage (2), diabetes (3), and oncological diseases (4). According to the World Health Organization, the global trend consists of a decreasing degree of smoking incidence throughout the world: In 2000, 33.3% of the world population was smoking; by 2015, however, this figure dropped to 24.9%. At the same time, the dynamic of such a decrease varies significantly depending on demographic and regional characteristics (5). The degree of smoking incidence among men in the majority of high-income countries began to decline in the mid-1990s, while the degree of smoking incidence in many low and middle-income countries remained unchanged or even increased. The degree of smoking incidence among women has also decreased in most countries, although the decline began later and was slower compared to that of men (6).

At the individual level, smoking addiction significantly varies depending on gender, age, and a number of socio-economic characteristics, such as education, income level, marital status, and profession (5, 7, 8). However, even with adjustments for individual characteristics, there are territorial differences in the degree of smoking incidence depending on the scale: International (7, 8), regional (9), administrative, and district (10).

In the mid-1990s, several leading epidemiologists started emphasizing the need to switch from the phenomenology of the prevalence of non-communicable diseases and their risk factors, to the development of procedures for identifying associated cause-and-effect relationships (1113). This gave a boost to the study of new epidemiological theories, including the socio-ecological model, which provides for biological and social production of diseases as equally significant components. This resulted in a large number of epidemiological surveys that studied the multi-level impact of territorial characteristics on individual health. For the past two decades, there have also been numerous studies of how the territory of residence affects the smoking status irrespective of individual characteristics (14, 15). It is significant that the impact of territorial living conditions on individual health indicators can vary considerably depending on the scale of the selected territory (1619). Most of the “smoking geography” studies focus on small territorial entities: Postal index zones, city districts, municipalities, and areas with a population of under 50,000–100,000 people (14, 15, 2027). There are considerably fewer works devoted to the research of territorial smoking predictors at the international level (2830). Moreover, studies of “smoking geography” covering territorial entities between countries and small areas, that is, large regions, provinces, and states within a country, are even more sporadic (3133). Apparently, such studies are likely to be of interest mostly for countries with a large territory.

Such surveys of “smoking geography” essentially focus on studying one or several social or economic characteristics selected as territorial predictors: The Gini index (23, 28, 31, 33), crime rate (20), Gross Domestic Product (28, 30), education level (21), per capita income (22, 23), unemployment rate (22), and poverty headcount (34). However, in reality, there are numerous complex relationships between territorial characteristics (17, 35). Therefore, associations of a specific territorial characteristic with health indicators may not indicate a cause-and-effect relationship, but a relationship mediated by other territorial characteristics. Thus, many researchers of “smoking geography” apply composite socio-economic indices calculated based on several social and economic territorial characteristics at once (14, 15, 24, 25, 28, 32). These studies allow to evaluate socio-economic territorial predictors of health status. However, despite the fact that the socio-economic environment is viewed as one of the most significant predictors of health status (and smoking in particular), and rightfully so, one, nevertheless, cannot but take into account the possible impact of other territorial characteristics. Once this line of reasoning is applied, it may be of interest to focus on those few studies that have employed an empirical approach to evaluate territorial characteristics, with subsequent analysis of their impact on individual health indicators (36), including smoking status (26, 27).

Finally, as part of the rationale for this study, it should not go unspoken that in Russia, there has been no analysis of the impact of territorial characteristics on individual health indicators, including smoking. A certain number of studies show geographical differences in the degree of smoking incidence (37, 38); moreover, they show a geographical trend of the degree of smoking incidence in Russia increasing from south to north and from west to east (39, 40). However, the reasons for such patterns have not been analyzed.

Thus, although it has been found that the socio-economic environment affects smoking habits, there is, however, no convincing evidence with regard to large territorial entities and other, non-socio-economic, characteristics. With this study, we aimed to evaluate the impact of the characteristics of large administrative regions on individual probability of smoking in the Russian adult population. Herein, we intended to find answers to the following questions: 1) What main groups (indices) of characteristics that describe Russian regions from different perspectives—geographic, demographic, social, economic, industrial, environmental, etc.—can be singled out? 2) What and how do regional living conditions affect the individual probability of smoking in the Russian population? 3) Do regional living conditions have any particular impact on smoking status depending on individual characteristics of the Russian people; in other words, are there any interactions between regional and individual levels?


Sample Description

For the purpose of analysis, we used data from the cross-sectional phase of the epidemiological study “Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation” (ESSE-RF) conducted in 2013–2014. A total of 21,923 individuals aged 25–64 years were examined. More detailed information on the ESSE-RF sampling and study protocol was given previously (41). In brief, the study was conducted in 13 regions of the Russian Federation. The sample was drawn based on the Kish method, which provides for systematic, multi-step, random community-based sampling on the premises of medical and preventive treatment facilities. The study was carried out in accordance with the standards of Good Clinical Practice and the principles of the Declaration of Helsinki. The study protocols were approved by the Ethics Committee of the National Medical Research Center for Therapy and Preventive Medicine (Moscow), the National Medical Research Center of Cardiology (Moscow), and the Almazov National Medical Research Centre (St. Petersburg), as well as by collaborating centers in the regions where this study was conducted. All participants gave their written informed consent prior to being included in the study. The response rate was approximately 80%, with some variations across the study regions.

St. Petersburg (1,588 people) was excluded from the final sample, since its regional characteristics are substantially different from those of the other 12 regions. The city of St. Petersburg is classified as a separate administrative territory in the Russian Federation, while the other 12 regions are large territories that include both urban and rural areas. Figure 1 shows the geographical location of the 12 regions participating in the study.


FIGURE 1. Regions of Russia participating in the study (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014). 1, Vologda; 2, Ivanovo; 3, Voronezh; 4, Volgograd; 5, Vladikavkaz; 6, Samara; 7, Orenburg; 8, Tyumen; 9, Tomsk; 10, Kemerovo; 11, Krasnoyarsk; 12, Vladivostok.

Another 32 individuals (0.2%) without a recorded smoking status were excluded from the remaining sample. As regards other individual factors, 265 people (1.3%) did not provide data on income, 152 individuals (0.7%) on marital status, and 15 individuals (0.1%) on the level of education. The missing data for these three factors were reconstructed using the k-nearest neighbor algorithm. Data imputation was carried out based on the input parameters, i.e., region, location, gender, and age. Thus, the final sample size with complete data (i.e., initial and reconstructed) amounted to 20,303 individuals. The general characteristics of the sample are given in Table 1.


TABLE 1. Baseline characteristics of the study participants (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).

Individual Variables and Responses

From individual variables, we selected the socio-economic and demographic characteristics as having the highest evidenced impact on smoking. These included gender, age, level of education (at non-university/university level), marital status (family/no family), and income level. Location was represented by urban or rural area. The income level was evaluated indirectly based on three questions characterizing the share of income spent on food, the respondents’ opinion on the family budget, and prosperity as compared to other families. Each question offered five possible answers that ranged in points from 1 (the “poorest” answer) to 5 (the “richest” answer). The point total was used to calculate terciles, which in turn were used to break the income level down into three categories: “Low,” from 3.0 to 7.2 points; “Middle,” from 8.0 to 10.3 points; and “High,” from 11.0 to 15.0 points. The current smoking status was construed as a response: Smokes (smoking of one or more cigarettes per day) or does not smoke (including quitted smoking). All individual variables were obtained by means of face-to-face interviews.

Regional Variables

We used figures provided on the official website of the Federal State Statistics Service of Russia for the 12 regions of the Russian Federation under study that give a measure of the following aspects: Climatic and geographic (4 factors), demographic (10 factors), social (13 factors), economic (9 factors), environmental (2 factors), industrial (6 factors), medical and infrastructure (3 factors), health status of the population (10 factors), and alcohol consumption (6 factors). The majority of the figures provided were for the period 2010–2014, except for gross regional product and per capita household consumption, which were for the period 2010–2013. Regional characteristics were averaged over the given period; a total of 64 regional characteristics were analyzed. A complete definition, average, minimum and maximum values, and standard deviations of all of the regional characteristics used are presented in Supplementary Table S1.

Regional characteristics are substantially interconnected, both contextually and statistically. In order to reduce the dimensionality of the data, we used principal component analysis (PCA), which is widely used to determine the composite indices of territories of residences based on several reference characteristics (14, 35, 4244). Our aim was to derive a parsimonious set of factors that captures the shared variance of a broad spectrum of region structural characteristics. We used varimax orthogonal rotation to simplify the structure of factors and to improve their interpretability. We identified factors with a share of explained variance of over 5%, with subsequent assessment of the gap according to the Cattell scree plot. Since there is a strong correlation between regional characteristics, we recognized a factor loading of | ≥0.65 | as a substantial contribution to the structure of factors. All of the resulting factor scores were standardized to have a mean of zero and a standard deviation of one.

We identified a total of five integrated factors (Table 2) that explain a cumulative 77.6% of the total variance. The first factor comprises ten characteristics: Average per capita consumption of vodka, wine, low-alcohol beverages and brandy, average annual air temperature (negative load on the factor), timberland area, per capita crime rate, geographical latitude of the regional center, proportion of decrepit and dilapidated housing, and proportion of students in the second and third shifts. The second factor comprises five characteristics: Natural population growth (negative load), birth (negative load) and mortality rates, proportion of population over unemployable age, and mortality from respiratory diseases. The third factor comprises eight characteristics: Rates of production of minerals and electric power, mortality from tuberculosis, infectious diseases and external causes, proportion of people in the region working under harmful working conditions, population size of the region, and emissions into the atmosphere. The fourth factor comprises five characteristics: Number of people employed in fisheries, average per capita volume of paid services, average per capita number of cars, male/female ratio (negative load), and geographical longitude of the regional center. The fifth factor comprises five characteristics: Per capita retail turnover, per capita household consumption, Gini index, per capita income of the population, and level of manufacturing in the region. These factors are easily interpreted, except for factor 4. Given the prevailing loads and for the purposes of this study, we named the identified factors as follows: Factor 1—Socio-geographic index; factor 2—Demographic index; factor 3—Industrial index; factor 4—Mixed index; and factor 5—Economic index.


TABLE 2. Factor loadings of the principal regional indices identified (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).

Statistical Analysis

We used bivariate statistics (proportions and chi-squared test) to summarize the characteristics of the samples and the distribution of the socio-economic status indicators. We used Pearson correlation to assess the interconnections between regional characteristics. The survey data were represented by a complex two-level sample with individual and regional characteristics, which requires the application of appropriate methods of statistical analysis. Studies based on complex cluster samples commonly use mixed regression models (for example, a generalized linear mixed model); however, a number of surveys have proven that application of marginal approaches that provide for a more robust and valid inference can also be successful (45, 46). Therefore, we used the generalized estimating equation (GEE) with constant standard errors to determine associations between regional indices and individual smoking levels, with due regard to the nested data structure (i.e., individuals in the regions). We completed several sets of logistic models of smoking probability that included a calculation of the odds ratio (OR) and the Wald statistic. The “zero” model included individual variables only. Model 1 included individual variables and all regional indices. Next, we evaluated interactions of the most important individual variables (i.e., gender, age, and level of education) and all of the regional indices. Since we detected substantial interactions between the individual variables of “gender” and “level of education” and the regional indices, we performed a separate analysis of Model 1 for men/women and various levels of education for the purpose of a better interpretation of such interactions. For the purposes of descriptive statistics, two-way analysis of categorical variables, correlation analysis, and PCA, we used Statistica Version 10.0 (Statsoft Inc., United States), and for the GEE, we used SPSS Version 22 (IBM Corp., United States).


The “zero” model showed a impact of gender, level of education, and age on the probability of smoking (Table 3). Model 1, which also includes the regional indices, showed a impact of individual income level. With regard to the regional indices, an increased probability of smoking is associated with an increased Socio-Geographic index (OR = 1.12: 1.03–1.22) and an increased Industrial index (OR = 1.15: 1.06–1.24). Chi-squared estimates of the likelihood Type III test are as follows, in descending order: age, 143.9; level of education, 104.2; gender, 90.1; Industrial index, 12.9; income level, 8.6; and Socio-geographic index, 7.0.


TABLE 3. Multivariate association of individual and regional characteristics with smoking (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).

We detected interactions between the regional indices and individual factors in terms of their impact on the probability of smoking (Table 4). There is interactions between gender and all of the regional indices, except for the Economic index, between age and the Socio-geographic index, and between education and the Socio-geographic and Mixed indices. Judging by the OR values, the impact of the indices on individual smoking is weaker or negative for men, people with university-level education, and people of older age.


TABLE 4. Interaction effects on smoking individual and regional characteristics (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).

Separate modeling shows that gender is a substantial factor in terms of the associations between individual smoking probability and characteristics of the region of residence (Table 5). Thus, the direct impact of the Socio-geographic and Demographic indices is typical for women only (respectively, OR = 1.47: 1.35–1.60 and OR = 1.18: 1.01–1.36), but not for men. The Industrial index directly affects the probability of smoking in both gender groups; however, the impact is more substantial for women (OR = 1.21: 1.10–1.33) than for men (OR = 1.08: 1.01–1.16). The Mixed index shows opposite trends based on gender: A direct impact on women (OR = 1.16: 1.10–1.22) and the reverse on men (OR = 0.90: 0.86–0.94).


TABLE 5. Multivariate association of individual and regional characteristics with smoking (gender stratification) (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).

In terms of the level of education (Table 6), the direct impact of the Socio-geographic index is typical for non-university education (OR = 1.19: 1.09–1.29), but not for university-level education. An inverse association between the Mixed index and the probability of smoking, being close to statistically significant, is found in individuals with university-level education (OR = 0.97: 0.94–1.00), but not with non-university education. It should be noted that the associations between the Industrial index and the probability of smoking are equal in their direction and intensity for both gradations of education.


TABLE 6. Multivariate association of individual and regional characteristics with smoking (education stratification) (“Epidemiology of Cardiovascular Diseases in the Regions of the Russian Federation”, Russia, 2013–2014).


The results we obtained allowed us to answer all of the questions posed, and provided an overall measure of the impact of region-specific characteristics on the degree of smoking incidence. Based on the 64 territorial characteristics, we determined five latent factors that form the main part of the variance: The Socio-geographic, Demographic, Industrial, Mixed, and Economic indices. Of these, the Economic index, surprisingly, showed no associations with the probability of smoking, in either the total sample size or the stratification analyses.

The Industrial index shows the most consistent associations, both in terms of the primary effects and in interactions with individual factors. Living in regions with a high level of mineral extraction, electric power production, poor labor conditions for the majority of workers (which may be related to the high mortality rate from external causes), and high levels of emissions into the atmosphere from stationary sources increases the individual probability of smoking.

The Socio-geographic, Demographic, and Mixed indices show substantial interactions with gender (most intense), age, and level of education. There are associations between the high levels of alcohol consumption in the region, accompanied by a high crime rate and deterioration of certain social conditions (i.e., quality of housing and the educational environment for children), as well as the climatic and geographic location of the regions (i.e., further to the north and colder). Living in such regions is associated with an increased individual probability of smoking among women, people with low-level education or qualifications, as well as older people.

Living in demographically depressed regions with a low birth rate, a high mortality rate, and, as a result, a negative rate of natural increase, as well as a large proportion of older people in the general population, is associated with an increased individual probability of smoking among women.

Finally, the Mixed index is the most incomprehensible, both in terms of its interpretation and its correlation with the probability of smoking. We found correlations between the development of fisheries and fish farms, the high volume of paid services in the region, the large number of private cars in the region, as well as the increased proportion of women in the general population, as well as the geographic location of the region (i.e., eastern longitude). Living in such regions increases the probability of smoking among women and, on the contrary, reduces it among men and people with university-level education.

The contribution of the selected territorial indices to the individual probability of smoking is rather small, especially in comparison to the substantial contribution of individual factors (i.e., age, level of education, and gender). This is consistent with the results of other multi-level studies of the degree of smoking incidence (22, 26). Nevertheless, this contribution is statistically significant, and may be of interest for the practical development of preventive measures and for monitoring of the degree of smoking incidence in administrative territories.

Comparison to Similar Multi-Level Studies

Conducting a comparative analysis of the obtained results and the results of other similar multi-level studies is not straightforward, since there are substantial differences in terms of the analyzed territorial characteristics and/or the indices determined on such a basis. First, it should be noted that there is a rather high level of differentiation for four of the five regional indices that we determined. Thus, it is not only the widely analyzed social and economic distinguishing characteristics, but also other regional characteristics that are not usually studied in such surveys that make a significant contribution to the overall variance. These are geographic characteristics, in terms of the industrial, environmental, and certain demographic factors. Most commonly, multi-level studies are based on either specific socio-economic characteristics (for example, the Gini index or the crime rate) or on a priori calculated composite socio-economic indices (for example, the Swedish Care Need Index (47) or the Japanese Areal Deprivation Index (25). Even those few studies that use an empirical assessment of territorial characteristics include mostly socio-economic characteristics in their fundamental lists (26, 27).

Studies covering districts, boroughs, and municipalities often show no associations between smoking and socio-economic characteristics, for example, the crime rate (20), the unemployment rate (22), the average per capita income (22, 23), or the Gini index (23). Furthermore, a number of studies, surprisingly, show positive associations between socio-economic characteristics and smoking, which means that an improvement of the socio-economic situation caused an increase in the degree of smoking incidence (14, 26). However, a rather large number of studies show the traditionally expected increase in the degree of smoking incidence amid the deteriorating socio-economic environment in terms of specific characteristics (22, 34) and, in particular, in terms of the composite socio-economic indices (15, 25, 27, 47).

At the same time, all of the few studies covering the territorial level of states and provinces show no correlation between the probability of smoking and socio-economic characteristics, for example, the Gini index for adolescents (33) and elderly people (31) and the socio-economic index for pregnant women (24). We also found that the Economic index did not show a correlation with the probability of smoking, either in terms of the total sample size or in terms of the specific characteristics of gender, age, and level of education. At the same time, the contribution of a number of social characteristics is evident, especially for women.

Review of the Results Obtained in Other Russian Studies

As noted above, earlier multi-center epidemiological studies in Russia (in 1993 and 2003–2004) showed a shift in the degree of smoking incidence toward the northern and eastern regions of the country (39, 40). The subsequent study conducted in 2013–2014 (we used the same individual data for our analysis) showed a slight change in the epidemiological situation in terms of descriptive statistics (38, 39). It showed an increased degree of smoking incidence from north to south and from west to east among men, and from south to north and from west to east among women. Certainly, the geographical location itself cannot be a predictor of the degree of smoking incidence, and the observed trends are also dependent on other regional characteristics. Our results are a “step forward” toward understanding these trends. First, we applied a statistical method that allows for reliable evaluation of complex data with a hierarchical (nested) structure, which increases the reliability of the “smoking geography” evaluation. Second, based on the latent factors identified with the help of PCA, the geographic characteristics correlate with other regional characteristics that presumably determine territorial associations with the degree of smoking incidence. The results show that the increased probability of smoking among women, as well as among older people in the northern regions of Russia, is associated with the respective deterioration of the social environment and, in particular, with the high levels of alcohol consumption, crime rate, and poor quality of housing, as well as the educational environment for children. At the same time, from west to east, the probability of smoking increases among women, but decreases among men and people with university-level education. This trend is also associated with the other territorial characteristics included in the Mixed index, which, unfortunately, is difficult to interpret.

Advantages and Limitations

This study is the first Russian multi-level analysis of hierarchical data on the “geography” of the degree of smoking incidence. For that matter, this study is the first Russia multi-level analysis of the “geography” of any health indicator whatsoever. We performed this analysis based on a large amount of data. In order to obtain scientific evidence, we used modern suitable methods of statistical analysis, for example, the GEE, for multi-level analysis purposes.

It should be noted that our empirical approach to the identification of territorial indices has not yet been used widely enough in similar studies. From this perspective, our study provides new scientific data, including for the purposes of the possible application of new methodological approaches in this particular scientific field.

Finally, it should be noted that we used an evaluation of the interactions between the impact of the territorial indices and individual factors on the probability of smoking for the purposes of our work. This allowed us to clarify the identified patterns in terms of the total sample size (that is, the key effects). Such evaluation of the interactions between the territorial and individual factors is rarely found in similar studies.

It should be noted that the shortcomings of the analysis include the relatively small number of regions (12 in total), that is, territorial entities under study. Clearly, this weakens the analysis in terms of using poor quantitative scales of territorial characteristics, which in turn may affect the resulting territorial indices, as well as the analysis of the key effects and the interactions between the territorial indices and the probability of smoking. On the other hand, such surveys are not unusual, especially when studying large territorial entities (33).

Another shortcoming is that the analysis also resulted in the applicability of a Mixed index that is difficult to interpret. However, obtaining such latent factors that are difficult to interpret is often the “flip side of the coin” of using statistical methods to reduce the dimensionality of data in the context of empirical approaches. It should be noted that the “socio-economic indices,” the “deprivation indices,” and the “district welfare indices” widely applied in similar studies also often use PCA to obtain the first key factor speculatively designated as socio-economic.


The analysis results allowed us to evaluate the impact of regional characteristics on the individual probability of smoking in a cross-sectional study of the Russian population. Based on the data collected by the official state statistics authorities of the Russian Federation, we determined the main groups (indices) of characteristics describing the Russian regions from different perspectives. The empirical approach we applied to determine the territorial indices is still quite new and poorly known in the context of multi-level studies, which means that this study provides new scientific data. We evaluated the key associations of the territorial indices with the probability of individual smoking, as well as interactions between the territorial indices and individual factors (i.e., gender, age, and level of education). The results we obtained provide, for a first-time, multi-level evaluation of the health status in Russia from the perspective of environmental epidemiology. Furthermore, this allowed us to provide a well-founded description of the “smoking geography” in Russia and to add the Russian results to the pool of similar global data. Since the degree and the dynamic of smoking incidence in Russia make this country one of the most adverse in the world (along with other countries of Eastern Europe), the results we obtained will be of interest for the purposes of healthcare management and preventive medicine in Russia.

Data Availability Statement

The datasets presented in this article are not readily available because of prohibition of transferring data to third parties. Requests to access the datasets should be directed to Svetlana Shalnova,

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the National Medical Research Center for Therapy and Preventive Medicine (Moscow). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

SAM contributed to conception and design of the study. YAB, SEE, and AEI organized the database. SAM and VAK performed the statistical analysis. SAM, SAS and OMD wrote the first draft of the manuscript. SAM, VAK, and SAS wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.


The study was conducted as part of a state assignment in National Medical Research Center for Therapy and Preventive Medicine of the Ministry of Healthcare of the Russian Federation. This research received no external funding.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The authors are sincerely grateful to those who assisted in the study and data collection: S. Bojcov, R. Oganov, A. Deev, A. Kapustina, G. Muromceva (Moscow); O. Rotar, E. Shljahto, E. Baranova, A. Konradi (St. Petersburg); N. Kulakova, V. Nevzorova, N. Shestakova, M. Mokshina, L. Rodionova (Vladivostok); G. Tolparov (Vladikavkaz); V. Il’in, A. Shabunova, K. Kalashnikov, O. Kalachikova, A. Popov (Vologda); S. Nedogoda, E. Chumachek, A. Ledjaeva (Volgograd); G. Furmenko, T. Chernyh, V. Ovsjannikova, L. Bondarcov (Voronezh); O. Belova, S. Romanchuk, O. Nazarova, O. Shutemova (Ivanovo); O. Barbarash, G. Artamonova, E. Indukaeva, T. Mulerova, A. Skripchenko, N. Cherkass, M. Tabakaev, J. Danil’chenko (Kemerovo); J. Grinshtejn, M. Petrova, L. Danilova, A. Evsjukov, V. Shabalin, R. Ruf, A. Kosinova, I. Filonenko, O. Bajkova (Krasnoyarsk); R. Libis, E. Lopina, I. Basyrova (Orenburg); D. Dupljakov, S. Gudkova, N. Cherepanova (Samara); I. Trubacheva, V. Kaveshnikov, R. Karpov, V. Serebrjakova (Tomsk); A. Efanov, I. Medvedeva, M. Storozhok, S. Shalaev (Tyumen).

Supplementary Material

The Supplementary Material for this article can be found online at:


1.Institute for Health Metrics and Evaluation (IHME). Findings from the global burden of disease study 2017. Seattle, United States: Institute for Health Metrics and Evaluation (IHME). (2018). Available at: (Accessed July 09, 2020).

2. Hackshaw, A, Morris, JK, Boniface, S, Tang, JL, and Milenković, D. Low cigarette consumption and risk of coronary heart disease and stroke: meta-analysis of 141 cohort studies in 55 study reports. BMJ (2018). 360:j5855. doi:10.1136/bmj.j5855

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Willi, C, Bodenmann, P, Ghali, WA, Faris, PD, and Cornuz, J. Active smoking and the risk of type 2 diabetes. JAMA (2007). 298(22):2654–64. doi:10.1001/jama.298.22.2654

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Pesch, B, Kendzia, B, Gustavsson, P, Jöckel, KH, Johnen, G, Pohlabeln, H, et al. Cigarette smoking and lung cancer-relative risk estimates for the major histological types from a pooled analysis of case-control studies. Int J Cancer (2012). 131(5):1210–9. doi:10.1002/ijc.27339

PubMed Abstract | CrossRef Full Text | Google Scholar

5.World Health Organisation. WHO global report on trends in prevalence of tobacco use 2000–2025. 3rd Edn. Geneva: World Health Organization. (2019). 121 p. Available at: (Accessed July 09, 2020).

6. Islami, F, Stoklosa, M, Drope, J, and Jemal, A. Global and regional patterns of tobacco smoking and tobacco control policies. Eur Urol Focus (2015). 1(1):3–16. doi:10.1016/j.euf.2014.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Degenhardt, L, Chiu, WT, Sampson, N, Kessler, RC, Anthony, JC, Angermeyer, M, et al. Toward a global view of alcohol, tobacco, cannabis, and cocaine use: findings from the world mental health surveys. Plos Med (2008). 5(7):e141. doi:10.1371/journal.pmed.0050141

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sinha, DN, Suliankatchi, RA, Amarchand, R, and Krishnan, A. Prevalence and sociodemographic determinants of any tobacco use and dual use in six countries of the WHO South-East Asia region: findings from the demographic and health surveys. Nicotine Tob Res (2016). 18(5):750–6. doi:10.1093/ntr/ntv286

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hu, SS, Homa, DM, Wang, T, Gomez, Y, Walton, K, Lu, H, et al. State-specific patterns of cigarette smoking, smokeless tobacco use, and e-cigarette use among adults-United States, 2016. Prev Chronic Dis (2019). 16:E17. doi:10.5888/pcd16.180362

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Snyder, RE, Rajan, JV, Costa, F, Lima, H, Calcagno, JI, Couto, RD, et al. Differences in the prevalence of non-communicable disease between slum dwellers and the general population in a large urban area in Brazil. Trop Med Infect Dis (2017). 2(3):E47. doi:10.3390/tropicalmed2030047

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Krieger, N. Epidemiology and the web of causation: has anyone seen the spider? Soc Sci Med (1994). 39(7):887–903. doi:10.1016/0277-9536(94)90202-x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Davey Smith, G. Post-modern epidemiology: when methods meet matter. Am J Epidemiol (2019). 188(8):1410–9. doi:10.1093/aje/kwz064

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Susser, M, and Susser, E. Choosing a future for epidemiology: II. From black box to Chinese boxes and eco-epidemiology. Am J Public Health (1996). 86(5):674–7. doi:10.2105/ajph.86.5.674

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Fleischer, NL, Thrasher, JF, Sáenz de Miera Juárez, B, Myriam Reynales-Shigematsu, L, Arillo-Santillán, E, Osman, A, et al. Neighbourhood deprivation and smoking and quit behaviour among smokers in Mexico: findings from the ITC Mexico survey. Tob Control (2015). 24(Suppl. 3):56–63. doi:10.1136/tobaccocontrol-2013-051495

CrossRef Full Text | Google Scholar

15. Algren, MH, Bak, CK, Berg-Beckhoff, G, and Andersen, PT. Health-risk behaviour in deprived neighbourhoods compared with non-deprived neighbourhoods: a systematic literature review of quantitative observational studies. PLoS One (2015). 10(10):e0139297. doi:10.1371/journal.pone.0139297

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ghenadenik, AE, Frohlich, KL, and Gauvin, L. Beyond smoking prevalence: exploring the variability of associations between neighborhood exposures across two nested spatial units and two-year smoking trajectory among young adults. Int J Environ Res Public Health (2016). 13(1):E106. doi:10.3390/ijerph13010106

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Diez-Roux, AV, Kiefe, CI, Jacobs, DR, Haan, M, Jackson, SA, Nieto, FJ, et al. Area characteristics and individual-level socioeconomic position indicators in three population-based epidemiologic studies. Ann Epidemiol (2001). 11(6):395–405. doi:10.1016/s1047-2797(01)00221-6

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bao, H, Zhang, L, Wang, L, Zhang, M, Zhao, Z, Fang, L, et al. Significant variations in the cervical cancer screening rate in China by individual-level and geographical measures of socioeconomic status: a multilevel model analysis of a nationally representative survey dataset. Cancer Med (2018). 7(5):2089–100. doi:10.1002/cam4.1321

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Schuurman, N, Bell, N, Dunn, JR, and Oliver, L. Deprivation indices, population health and geography: an evaluation of the spatial effectiveness of indices at multiple scales. J Urban Health (2007). 84(4):591–603. doi:10.1007/s11524-007-9193-3

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Felker-Kantor, EA, Cunningham-Myrie, C, Greene, LG, Lyew-Ayee, P, Atkinson, U, Abel, W, et al. Neighborhood crime, disorder and substance use in the caribbean context: Jamaica national drug use prevalence survey 2016. PLoS One (2019). 14(11):e0224516. doi:10.1371/journal.pone.0224516

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Räisänen, S, Kramer, MR, Gissler, M, Saari, J, Hakulinen-Viitanen, T, and Heinonen, S. Smoking during pregnancy was up to 70% more common in the most deprived municipalities - a multilevel analysis of all singleton births during 2005–2010 in Finland. Prev Med (2014). 67:6–11. doi:10.1016/j.ypmed.2014.06.026

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hamad, R, Brown, DM, and Basu, S. The association of county-level socioeconomic factors with individual tobacco and alcohol use: a longitudinal study of United States adults. BMC Public Health (2019). 19(1):390. doi:10.1186/s12889-019-6700-x

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Galea, S, Ahern, J, Tracy, M, and Vlahov, D. Neighborhood income and income distribution and the use of cigarettes, alcohol, and marijuana. Am J Prev Med (2007). 32(Suppl. 6):S195–202. doi:10.1016/j.amepre.2007.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lian, M, Madden, PA, Lynskey, MT, Colditz, GA, Lessov-Schlaggar, CN, and Schootmanet, M. Geographic variation in maternal smoking during pregnancy in the Missouri adolescent female twin study (MOAFTS). PLoS One (2016). 11(4):e0153930. doi:10.1371/journal.pone.0153930

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hanibuchi, T, and Nakaya, T. Associations of neighborhood socioeconomic conditions with self-rated health, mental distress, and health behaviors: a nationwide cross-sectional study in Japan. Prev Med Rep (2020). 18:101075. doi:10.1016/j.pmedr.2020.101075

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Pasman, JA, Verweij, KJH, Abdellaoui, A, Hottenga, JJ, Fedko, IO, Willemsen, G, et al. Substance use: interplay between polygenic risk and neighborhood environment. Drug Alcohol Depend (2020). 209:107948. doi:10.1016/j.drugalcdep.2020.107948

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Cambron, C, Kosterman, R, Catalano, RF, Guttmannova, K, and Hawkins, JD. Neighborhood, family, and peer factors associated with early adolescent smoking and alcohol use. J Youth Adolesc (2018). 47(2):369–82. doi:10.1007/s10964-017-0728-y

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mutumba, M, and Schulenberg, JE. Tobacco and alcohol use among youth in low and middle income countries: a multi-country analysis on the influence of structural and micro-level factors. Subst Use Misuse (2019). 54(3):396–411. doi:10.1080/10826084.2018.1497063

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Pförtner, TK, Rathmann, K, Moor, I, Kunst, AE, and Richter, M. Social inequalities in adolescent smoking: a cross-national perspective of the role of individual and macro-structural factors. Bundesgesundheitsblatt Gesundheitsforschung Gesundheitsschutz (2016). 59(2):206–16. doi:10.1007/s00103-015-2280-9

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Pampel, FC, and Denney, JT. Cross-national sources of health inequality: education and tobacco use in the world health survey. Demography (2011). 48(2):653–74. doi:10.1007/s13524-011-0027-2

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Cohen, SA, Greaney, ML, and Klassen, AC. A “Swiss paradox” in the United States? Level of spatial aggregation changes the association between income inequality and morbidity for older Americans. Int J Health Geogr (2019). 18(1):28. doi:10.1186/s12942-019-0192-x

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Shmulewitz, D, Stohl, M, Keyes, KM, Brown, Q, Saha, TD, and Hasin, D. Effects of state-level tobacco environment on cigarette smoking are stronger among those with individual-level risk factors. Nicotine Tob Res (2016). 18(10):2020–30. doi:10.1093/ntr/ntw114

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Quon, EC, and McGrath, JJ. Province-level income inequality and health outcomes in Canadian adolescents. J Pediatr Psychol (2015). 40(2):251–61. doi:10.1093/jpepsy/jsu089

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kravitz-Wirtz, N. A discrete-time analysis of the effects of more prolonged exposure to neighborhood poverty on the risk of smoking initiation by age 25. Soc Sci Med (2016). 148:79–92. doi:10.1016/j.socscimed.2015.11.027

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Lian, M, Struthers, J, and Liu, Y. Statistical assessment of neighborhood socioeconomic deprivation environment in spatial epidemiologic studies. Open J Stat (2016). 6(3):436–42. doi:10.4236/ojs.2016.63039

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Morenoff, JD, House, JS, Hansen, BB, Williams, DR, Kaplan, GA, and Hunte, HE. Understanding social disparities in hypertension prevalence, awareness, treatment, and control: the role of neighborhood context. Soc Sci Med (2007). 65:1853–66. doi:10.1016/j.socscimed.2007.05.038

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Andreeva, EA, Pokhaznikova, MA, and Kuznetsova, OY. The prevalence of smoking among residents of two cities in the North-West region of Russia according to the RESPECT international study. Profilakticheskaya Meditsina (2020). 23(1):92–9. doi:10.17116/profmed20202301199

CrossRef Full Text | Google Scholar

38. Balanova, YA, Kontsevaya, AV, Shalnova, SA, Deev, AD, Artamonova, GV, Gatagonova, DV, et al. Prevalence of behavioral risk factors for cardiovascular disease in the Russian population: results of the ESSE-RF epidemiological study. Profilakticheskaya Meditsina (2014). 17(5):42–52.

Google Scholar

39. Balanova, YA, Shalnova, SA, Deev, AD, Kapustina, AV, Konstantinov, VV, and Boitsov, SA. Smoking prevalence in Russia. What has changed over 20 years? Profilakticheskaya Meditsina (2015). 18(6):47–52. doi:10.17116/profmed201518647-52

CrossRef Full Text | Google Scholar

40. Shalnova, SA, Deev, AD, and Oganov, RG. Prevalence of smoking in Russia. Survey results of a national representative population sample. Profilakticheskaya Meditsina (1998). 1(3):9–12.

Google Scholar

41. Boitsov, SA, Chazov, EI, Shlyakhto, EV, Shalnova, SA, Konradi, AO, Karpov, YA, et al. Scientific organizing committee of the ESSE-RF. Epidemiology of cardiovascular diseases in different regions of Russia (ESSE-RF). The rationale for and design of the study. Profilakticheskaya Meditsina (2013). 16(6):25–34.

Google Scholar

42. Messer, LC, Laraia, BA, Kaufman, JS, Eyster, J, Holzman, C, Culhane, J, et al. The development of a standardized neighborhood deprivation index. J Urban Health (2006). 83(6):1041–62. doi:10.1007/s11524-006-9094-x

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Schederecker, F, Kurz, C, Fairburn, J, and Maier, W. Do alternative weighting approaches for an Index of multiple deprivation change the association with mortality? A sensitivity analysis from Germany. BMJ Open (2019). 9(8):e028553. doi:10.1136/bmjopen-2018-028553

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lamnisos, D, Lambrianidou, G, and Middleton, N. Small-area socioeconomic deprivation indices in Cyprus: development and association with premature mortality. BMC Public Health (2019). 19(1):627. doi:10.1186/s12889-019-6973-0

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Hubbard, AE, Ahern, J, Fleischer, NL, Van der Laan, M, Lippman, SA, Jewell, N, et al. To GEE or not to GEE: comparing population average and mixed models for estimating the associations between neighborhood risk factors and health. Epidemiol (2010). 21(4):467–74. doi:10.1097/EDE.0b013e3181caeb90

CrossRef Full Text | Google Scholar

46. Hanley, JA, Negassa, A, Edwardes, MD, and Forrester, JE. Statistical analysis of correlated data using generalized estimating equations: an orientation. Am J Epidemiol (2003). 157(4):364–75. doi:10.1093/aje/kwf215

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Cubbin, C, Sundquist, K, Ahlén, H, Johansson, SE, Winkleby, MA, and Sundquist, J. Neighborhood deprivation and cardiovascular disease risk factors: protective and harmful effects. Scand J Public Health (2006). 34(3):228–37. doi:10.1080/14034940500327935

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: smoking, regional environment, health geographic, multi-level study, Russia

Citation: Maksimov SA, Shalnova SA, Balanova YA, Kutsenko VA, Evstifeeva SE, Imaeva AE and Drapkina OM (2021) What Regional Living Conditions Affect Individual Smoking of Adults in Russia. Int J Public Health 66:599570. doi: 10.3389/ijph.2021.599570

Received: 27 August 2020; Accepted: 01 March 2021;
Published: 01 April 2021.

Edited by:

Mariachiara Di Cesare, Middlesex University, United Kingdom

Copyright © 2021 Maksimov, Shalnova, Balanova, Kutsenko, Evstifeeva, Imaeva and Drapkina This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sergey A. Maksimov,

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.