Abstract

Household air pollution from solid cooking fuel use influences multiple health outcomes, but its association with body pain remains poorly understood. This was a longitudinal study of 8880 adults who participated in the China Health and Retirement Longitudinal Study (CHARLS) from 2011 to 2018. Household cooking fuels were extracted from the baseline household questionnaire. Transitions in cooking fuels from 2011 to 2018 were also identified. Body pain status was reported in the three waves of surveys conducted in 2011, 2015, and 2018. The associations between cooking fuel type, fuel transition, and pain site number were examined using generalized estimating equations. Among the 8880 participants, 41.4% () primarily used clean fuels for cooking, and 58.6% () used solid ones at baseline. Cooking with solid fuels was associated with more pain sites (incidence rate ratio (IRR): 1.14; 95% confidence interval (CI): 1.08 to 1.21), but a slower rate of pain sites increases from 2011 to 2018 (; 95% CI: 0.71 to 0.86, for fuels). Compared with those who persistently used clean fuels for cooking, the number of pain sites increased by 10% in participants who transiting from using solid to clean fuels (; 95% CI: 1.04 to 1.18), by 21% in those transiting from cooking with clean to solid fuels (: 95% CI: 1.08 to 1.35) and by 25% among those persistent using solid fuels for cooking (; 95% CI: 1.18 to 1.34). Our findings provided new evidence linking using solid fuels for cooking with more pain sites, but a slower rate of pain sites increases. Public health efforts should focus on fuel transition and take measures to help clean fuels spread.

1. Introduction

After adopting the United Nations Declaration on the Human Environment in 1972, air quality improvement has become a general trend. However, public health awareness of indoor air pollution has lagged behind outdoor air pollution. Billions of the world’s population are still exposed to cooking with solid fuels, a major source of indoor air pollution [1]. Compared with those using cleaner fuels (like gas or electricity), burning of solid fuels (such as coal, wood, or crop residues) in inefficiency simple stoves generates substantial air pollutants, including carbon monoxide, nitrogen dioxide, and particulate matter, and results in adverse consequences for air qualities [2], climate change [3], and human health [4, 5].

In 2019, household air pollution was ranked 9th in terms of attributable deaths among 69 risk factors considered potentially modifiable [6]. About 2.3 million deaths and 91.5 million disability-adjusted life years, most occurring in China and India, were due to household air pollution from polluting cooking fuel [7]. Epidemiological studies have linked cooking with solid fuels to a broad range of health problems, including respiratory illness [8], cardiovascular diseases [9, 10], arthritis [11], cognitive decline [12], and depression [13]. Existing studies showed that outdoor air pollution was associated with body pain, including headaches [14], chest [15], and angina pectoris [16], but little is known about the relationship between indoor air pollution and body pain.

Pain is defined by the International Association for the Study of Pain (IASP) as an unpleasant sensory and emotional experience associated with actual or potential tissue damage [17]. It is known that pain is associated with consequences that far exceed itself, including adverse biological (i.e., cognitive decline and disability) [18, 19], psychological (i.e., depression, anxiety, and stress) [20], and social conditions (i.e., economic costs, isolation, and social withdrawal) [21, 22]. According to the Global Burden of Disease Study 2019 [23], pain conditions such as low back pain, neck pain, and headache were the most important cause of years lived with disability. Indeed, the pain increases with age and is the main reason people seek healthcare services [24]. However, it is often neglected or shifted attention to diagnosed diseases after treatment. In this setting, additional attention to the pain is warranted.

2. Mechanisms, Materials, and Methods

2.1. Theoretical Mechanism

The classification of pain is various and depends on the perspectives. When considering the duration of pain, the cutoff of three months is used for defining acute and chronic pain. According to the cause of pain, the latter can be further divided into the following seven categories: chronic primary pain, chronic neuropathic pain, chronic secondary visceral pain, chronic posttraumatic and postsurgical pain, chronic secondary headache and orofacial pain, chronic cancer-related pain, and chronic secondary musculoskeletal pain [25]. In addition, pain can be classified into different types based on its location, such as headache, chest pain, and low back pain. As the site of pain is easily felt, screened, and located, many epidemiological studies used body parts to define the pain [2628]. The current study focused on the association between solid fuel usage and the number of pain sites.

2.1.1. Air Pollution and Musculoskeletal Pain

Musculoskeletal pain is generally defined as pain caused by joints, bones, muscles, and tendons. According to the Standardized Nordic Questionnaire [29], the affected body area(s) of musculoskeletal pain include neck, shoulder, upper back, elbows, lower back, wrists/hands, hips/thighs, knees, and ankles/feet. Previous studies have found that exposure to air pollution was associated with musculoskeletal pain [27], bone fractures [30], and osteoporosis [31]. Indeed, air pollution containing particulate matter leads to systemic inflammation [32], which is directly associated with immune cells and nociceptors and indirectly modulates pain intensity [26, 33]. In addition, the plausibility of an association between air pollution and musculoskeletal pain is supported by studies showing musculoskeletal pain associated with tobacco smoking [34, 35]. Carbon monoxide is not only a component of smoke but also an important component of ambient and household air pollution. It leads to reduced perfusion and malnutrition of tissues and causes these tissues to respond inefficiently to mechanical stress, which may impair the bone remodeling process [34].

2.1.2. Air Pollution and Pain in Other Parts of the Body

Evidence exists on the association between air pollution and headache and chest pain. In a randomized stove intervention trial [36], the headache was substantially reduced in the plancha group relative to the group using open fires for the follow-up period. A recent review [37] also has found that exposure to airborne pollutants, especially particulate matter and nitrates, is associated with headaches. As we breathe, the pollutants produced by traffic, industrial processes, and burning fossil fuels may enter the bloodstream via the lungs or move directly to the brain via the olfactory nerves. Reactive pollutant chemicals are then activated in the body and provoke inflammatory processes that may cause headaches or migraines [14, 38, 39]. Chest pain, where angina commonly occurs, is also related to air pollution [16, 40]. The pollutants may cause chest pain by affecting the cardiovascular system through systemic inflammation, systemic oxidative stress, thrombosis, and coagulation [41, 42]. Although, few studies [43] have found a relationship between air pollution and abdominal pain. As air pollutants can enter the body via the lungs and gut, growing evidence links its exposure to gastrointestinal and respiratory impairments and diseases [44]. The directly toxic effects starting from the lungs may influence the gastrointestinal system through the gut-lung axis, leading to gut microbiota effects [45], gut permeability effects, and even systemic inflammation.

2.1.3. The Development of Multisite Pain

Pain rarely occurs in a single anatomical location, and the more widespread the pain is, the higher the impact on functioning and health [46, 47]. The study of multisite pain and its underlying factors has gained scientific interest internationally. One possible mechanism for explaining the development of multisite pain is repeated stressors. Individual risk factors may act accumulatively through the concept of allostatic load, whereby repeated stressors increase the physiological stress response in individuals and further lead to widespread pain [48]. Another potential mechanism for explaining why pain becomes widespread is central sensitization [49], which is an increased synaptic response in nociceptive neurons in the central nervous system.

Based on the discussion above, ambient air pollution and smoking are important risk factors for pain localized in different regions, including the musculoskeletal system, chest, head, and abdomen. Few studies are exploring household air pollution and the trajectory of pain sites. Still, ambient air pollutants and tobacco smoking have a similar component with household air pollutants, like carbon monoxide and particulate matter, and this provided us with a theoretical basis for studying the relationship between household air pollution and pain. In addition, as the exposure may persist, single pain may develop into multisite pain. This study seeks to answer the following research questions: (1)Is exposure to solid fuels for cooking associated with more pain sites?(2)Does transition in solid fuels influence the trajectory of pain sites?(3)Is exposure to solid fuels for cooking associated with a specific pain site?

2.2. Study Design and Population

The China Health and Retirement Longitudinal Study (CHARLS) is an open-access dataset that enrolled a nationally representative sample of Chinese residents aged 45 and above [50]. Briefly, CHARLS enables the study of the older population’s health in China patterned after the US Health and Retirement Study (HRS) and related aging surveys worldwide. This national survey has collected rich information on the respondents’ social, economic, behavioral, and health. The first survey was conducted in 2011 with a response rate of 80.5%, and 17708 individuals were interviewed. Then, these participants are followed up once every two years, and data from the repeated surveys performed in 2011, 2013, 2015, and 2018 are available. The details of the design and methods of CHARLS have been described extensively. CHARLS has been approved by the Biomedical Ethics Review Committee of Peking University.

For the current study, we used data from three repeated surveys of the CHARLS (2011, 2015, and 2018) because information on body pain differed in the second survey (2013) and then limited participants to 45-75 years old at baseline. After excluding those missing data on cooking fuel in 2011 or 2018, without data on body pain across three repeated surveys, and missing baseline data on any other variables of interest, the final analytic sample included 8880 participants, as detailed in Figure 1.

2.3. Outcome Ascertainment

Body pain was defined using one item: are you often troubled with any body pains? If the answer was “yes,” then the next question would be asked: on what part of your body do you feel pain? Please list all body parts that are currently in pain (including head, shoulder, arm, wrist, fingers, chest, stomach, back, waist, buttocks, leg, knees, ankle, toes, and neck). The primary outcome of the study was the number of body pain sites.

Specific body pain site was the second outcome in this study. According to previous studies [51, 52], we categorized the above 15 pain sites into 7 regions: head, neck/shoulders, arms/hands, frontal torso/genitals, back, buttocks, and legs/feet.

2.4. Exposure Measurement

Household cooking fuels were extracted from the baseline household questionnaire, which asked, what is the main source of cooking fuel? The categories of this question include (1) solar, (2) coal, (3) natural gas, (4) liquefied petroleum gas, (5) electric, and (6) crop residue/wood burning. According to previous studies [53, 54], respondents in households that answered coal, crop residue, or wood burning were deemed primarily to use solid fuels for cooking. In contrast, those who used natural gas, marsh gas, liquefied petroleum gas, or electricity were considered to use clean fuels.

Transition in fuel use was proxied using data from 2011 to 2018, and we identified change in fuel use into four categories: (1) solid fuels in 2011/solid fuels in 2018; (2) solid fuels in 2011/clean fuels in 2018; (3) clean fuels in 2011/solid fuels in 2018; (4) clean fuels in 2011/clean fuels in 2018. The first condition was considered to be always using solid fuels, and the last was deemed to be always using clean fuels. Solid to clean fuels presented those who used solid fuels for cooking in 2011 but switched to clean fuels in 2018, whereas those who used clean fuels for cooking in 2011 but switched to solid fuels in 2018 were deemed clean to solid fuels.

2.5. Covariates

Based on previous studies [12, 53, 54], age, gender, education level, marital status, residence, currently working, smoking status, alcohol consumption, number of diseases, household consumption expenditure, number of people living in the household, and heating fuel types were chosen as covariates and acquired from the baseline questionnaire (survey in 2011). The categories of most covariates were according to the questionnaire, including gender (male or female), education level (less than lower secondary, upper secondary, and vocational training or tertiary), residence (rural village or urban community), currently working (yes or no), smoking status (never or current/ever), and alcohol consumption (never or current/ever). Currently working indicates the respondent engaged in any work in the past year. Smoking status indicates whether the respondent reports ever smoking, and alcohol consumption indicates whether the respondent has had any alcoholic beverages in the past. We reclassified five variables: marital status, number of diseases, household consumption expenditure, number of people living in the household, and heating fuel types. In the database, the categories of marital status include (1) married with spouse present, (2) married but not living with spouse temporarily for reasons such as work, (3) separated, (4) divorced, (5) widowed, or (6) never married. Marital status in the study was classified as married and unmarried. Married included the first two statuses, whereas unmarried included the last four statuses. Given that multimorbidity (having at least two diseases at the same time) is associated with a higher risk of mortality and adverse health outcomes [55], we categorized the number of conditions (coded as 0, 1, >1) in our analyses. Household consumption expenditure indicates the amount of total household consumption as aggregated from all consumption activities, including food and other activities, and is calculated as annual expenditures. We used household consumption expenditure to represent the economic level of the household, which was confirmed as a better measure than income in developing countries. According to other studies [5, 56], we classified it as a four-categorical variable by quartiles (first quartile, < CNY¥9201; second quartile, ¥9201 to <16700; third quartile, ¥16700 to <27760; fourth quartile, ≥ ¥27760). The number of people living in the household is the number of people coresiding with the respondent, including the respondent. Generally, the family structure of two persons is composed of the respondent and the spouse. In contrast, more than two people’s family structure includes the respondent, the spouse, and their children. Thus, the number of people living in the household was classified as one person, two people, or more than two people in the analyses. Similar to the categories of cooking fuels, we reclassified heating fuels as clean fuels, solid fuels, and other fuels.

2.6. Statistical Analysis

The baseline characteristics of the study sample are presented as mean (standard deviation (SD)) for continuous variables, mean (interquartile range (IQR)) for count data, and number (percentage) for categorical variables, respectively. We compared baseline characteristics of the study sample according to the type of cooking fuels using Pearson’s chi-square tests, analysis of variance (ANOVA), and Kruskal-Wallis rank-sum tests. In the primary analyses, the association between baseline polluting cooking fuels and the number of pain sites from 2011 to 2018 was analyzed using a repeated-measure negative binomial regression analysis by the generalized estimating equation (GEE) method with an unstructured correlation structure after the model selection (Table S1 & S2) [57, 58]. The results are presented as the incidence rate ratios (IRRs) with 95% confidence intervals (CIs). Three models adjusted for factors that might affect the association, including age, gender, marital, residence, education level, and currently working (model 1); further adding smoking status, alcohol consumption, the number of pain sites at baseline, and the number of diseases (model 2); and adding household consumption expenditure, the number of people living in the household, and heating fuels (model 3). We further included statistical interaction terms between survey year and cooking fuel type to investigate whether fuel type was related to the rate of change in pain site number over time. The association between cooking fuel’s transition (2011/2018) and the number of pain sites was also analyzed using a GEE model with a log link and negative binomial distribution and tested with the same covariates as our primary analyses.

A series of exploratory and sensitive analyses were performed to assess the robustness of the primary findings. First, we evaluated the influence of moderators on the association between fuel types and the number of pain sites by assessing interaction terms of fuel types and each covariate. Second, considering diverse effects on pain sites, we examined the associations between polluting cooking fuels and each pain site (separate models) by using a GEE model with a logit link and logistic distribution. Results are presented as odds ratio (OR) with 95% CI. Third, we compared the baseline characteristics of those excluded for missing and included data. We further performed multiple imputations by chained equations (MICE) and repeated the modeling analyses using the imputed datasets to address potential bias from missing data [59]. Fourth, we excluded those having chronic diseases at baseline and repeated the analyses. Fifth, we further controlled for physical activity in the model. Finally, we examined the associations of cooking fuel types and the number of pain sites using different working correlation structures to address potential bias from other working correlation structures. The detailed analyses are provided in Supplementary Methods.

Analyses were performed with Stata/SE 17.0, and a two-side value level of 0.05 was considered the threshold for statistical significance.

3. Results

3.1. Baseline Sample Characteristics

Of the 8880 eligible participants, 47.9% were male, and the mean baseline age was 57.1. Among these, 41.4% () primarily used clean fuel for cooking, and 58.6% () used solid fuels for cooking (Table 1). Participants who cooked with solid fuels were more likely to have lower levels of education, live in rural villages, have no employment, be in a lower quartile level of household consumption expenditure, use solid heating fuels, and have more diseases when compared to those using clean cooking fuels. Table S3 provides additional information about the distribution of cooking fuel types and fuel transitions from 2011 to 2018. In brief, from 2011 to 2018, 3324 (37.4%) participants persistently used clean fuels for cooking, 2829 (31.9%) participants kept using solid fuels for cooking, 2371 (26.7%) individuals transitioned from using solid to clean fuels, and 356 (4.0%) individuals transitioned from cooking with clean to solid fuels.

3.2. Association of Cooking with Solid Fuels and the Number of Pain Sites

Table 2 presents the association between cooking with solid fuels and the number of body pain sites. In model 1, individuals who used solid fuels for cooking experienced more pain sites than those who used clean fuels (; 95% CI: 1.39 to 1.56). When smoking status, alcohol consumption, the number of diseases, and the number of pain sites at baseline were added in model 2 and household consumption expenditure, the number of people living in the household and heating fuels were added in model 3, the parameter estimate for polluting cooking fuels weakened but was still related with the increased number of pain sites (model 2: , 95% CI: 1.13 to 1.25; model 3: , 95% CI: 1.08 to 1.21).

When considering transitions in cooking fuel, in model 1, the number of pain sites increased in participants who were transiting from using solid to clean fuels (; 95% CI: 1.27 to 1.47), transiting from cooking with clean to solid fuels (; 95% CI: 1.21 to 1.59), and those persistent using solid fuels for cooking (; 95% CI: 1.60 to 1.84) compared with those who persistently used clean fuels for cooking. Adjustment for more covariates in model 2 and model 3 resulted in weaker parameter estimates.

3.3. Association of Cooking with Solid Fuels and Short-Term Trajectories of Body Pain Sites

Results from the interaction terms suggest that the effect of using solid cooking fuels was modified by time (Table 3 and Figure 2). In model 1, individuals who used solid cooking fuels had a slower rate of pain site increase (; 95% CI: 0.80 to 0.98, for fuels; ; 95% CI: 0.64 to 0.77, for fuels). In the fully adjusted model, the interaction effect was attenuated in 2015 (; 95% CI: 0.90 to 1.14, f fuels) but remained statistically significant in 2018 (; 95% CI: 0.71 to 0.86, for fuels).

Similarly, a negative interaction was observed between cooking fuel transition and time on the body pain sites in the fully adjusted model (; 95% CI: 0.69 to 0.87, for to clean; ; 95% CI: 0.68 to 0.85, for solid); and the slower rate of body pain site increase was observed in those who persistently used solid fuels for cooking as well as those transited from cooking with solid to clean fuels compared with those persistently used clean fuels (Table 3 and Figure 3).

3.4. Subgroup Analyses and Sensitivity Analyses

In the subgroup analyses, statistically significant interactions ( for interaction <0.05) on the number of pain sites were detected for fuel types across subgroups of gender, currently working, smoking status, number of diseases, and household consumption expenditure, as detailed in Table S4. The effect of using solid cooking fuels on the number of pain sites was more pronounced among males, smokers, those not having jobs, those with the highest level of household consumption expenditure, or those with one chronic disease (Figure S1).

The sensitivity analyses yielded similar findings to our primary analyses (Tables 2 and 3). First, we assessed the impact of cooking with solid fuels on the specific pain, and additional information is presented in Table S5 and Figure S2. Compared to clean fuel use for cooking, using solid fuels was associated with higher odds of leg/feet pain (, 95% CI: 1.10 to 1.28), buttock pain (, 95% CI: 1.07 to 1.39), back pain (, 95% CI: 1.05 to 1.28), frontal torso/genital pain (, 95% CI: 1.04 to 1.21), arm/hand pain (, 95% CI: 1.17 to 1.39), neck/shoulder pain (, 95% CI: 1.03 to 1.21), and headache (, 95% CI: 1.05 to 1.25). Second, the baseline characteristics varied between the excluded and included participants (Table S6). We performed multiple imputations by chained equations (MICE) to address potential bias from missing data and repeated the modeling analyses using the imputed datasets. The associations were not materially altered (Table S7 & S8). Third, we excluded those with chronic diseases at baseline to reduce bias caused by diseases. The results were similar to those observed in our primary analysis (Table S9). Still, the interaction terms between exposure variables and time were not statistically significant (Table S10), suggesting the effect of chronic diseases on the rate at which the number of pains increased. Fourth, gathering and using solid fuels increase the risk of musculoskeletal injuries, which also cause body pain; thus, we further controlled for physical activity, and there were no changes in the statistical significance of associations (Table S11 & S12). Finally, to address potential bias from different working correlation structures, we repeated the modeling analyses using the other two types of working correlation structures, and the main results remained robust (Table S13).

4. Discussion

Our main findings were that cooking with solid fuels in this Chinese population-based study was associated with more pain sites, but a slower rate of pain sites increases. In addition, males, smokers, those not having jobs, those with the highest household consumption expenditure, and those with one chronic disease responded more strongly to the harmful effects of indoor air pollution. When fuel switching was further considered, our study observed a nonconventional energy transition from clean to solid fuels and highlighted its impact on pain.

The evidence on the association between air pollution and pain is based chiefly on outdoor air pollution or limited to a specific pain site. A stove intervention study [60] in Peru observed reduced headache in the intervention group, those with new chimney-equipped stoves, than those with mostly open-fire stoves in the control group. Another study in multiple countries showed that higher odds of angina pectoris were observed mainly with household solid fuel use [61]. Other researchers [27, 30] gave evidence of an association between air pollution and musculoskeletal disorders. By comparison, our study is the first to simultaneously investigate the relationship between using solid fuels and the number of pain sites and specific pain sites. It extends previous knowledge by demonstrating that solid cooking fuel use was associated with more pain sites. In addition, we examined the impact of fuel transition on pain and found that transiting from clean to solid fuels and from solid to clean fuels was associated with more pain sites than those who persist in using clean fuels. This finding suggested that despite later switching to clean fuels, it was impossible to quickly offset the hazard from earlier exposure to polluting fuels, which underscored the importance of clean energy.

The potential biological mechanisms by which cooking with solid fuels influences the type and number of pain sites remain unclear; however, clues can be taken from the impact of outdoor air pollution and tobacco smoking on pain. Ambient air pollutants and tobacco smoking have components similar to household air pollutants, like carbon monoxide and particulate matter, and this provided us with a theoretical basis for studying the relationship between household air pollution and pain. The effect of pollution on health is commonly attributable to substantial amounts of pollutants released during combustion, including respirable particulates, carbon monoxide, nitrogen oxides, benzene, and polyaromatic compounds [62]. Sustained exposure to these air pollutants, especially particulate matter, may lead to elevated concentrations of proinflammatory biomarkers in the systemic circulation [63, 64], sensitize nociceptors, and amplify nociceptive input to the central nervous system [65, 66]. Systemic inflammation thus contributes to pain, a cardinal feature of inflammation [67, 68]. While the pain results from inflammation are adaptive, it still needs to be reduced; if not, ongoing inflammation may cause central sensitization, modify pain perception thresholds, and transit acute localized pain into chronic widespread pain [69, 70]. Another possible mechanism for the relationship between using solid fuels and pain involves oxidative dysfunction. Some pollutants, such as nitrogen oxides and carbon monoxide, contribute to the intracellular formation of free radicals and cause oxidative damage to DNA, lipids, proteins, and cells [71, 72]. Previous studies [39, 73] supported the potential role of oxidative stress and related structures in the pathogenesis of migraine. Mounting evidence [73, 74] has also revealed that oxidative dysfunction is critically involved in the induction and maintenance of neuropathic pain. In addition, pain comorbidities may explain the link between cooking with solid fuels and pain. Exposure to solid fuels has been associated with a series of health consequences, including arthritis [11], cognitive decline [12], and depression [13]. Much evidence from epidemiological and animal models shows that patients with chronic pain often coexist with the above diseases [20, 51, 75], and body pain frequently presents as a result of a disease or an injury. Our study further explored a more precise and individual reflection of pain and observed that pain localized in the musculoskeletal system, including arm/hands, buttock, and legs/feet, exhibited the most prominent associations. Several potential mechanisms may explain the association, including systemic inflammation, changes in pain intensity, reduced perfusion, and malnutrition of tissues [26, 34]. We also observed a positive association between solid fuel use and pain in other body parts like the head, chest, and hips, and the mechanisms may differ from location to location. As we breathe, the pollutants produced by burning fossil fuels may enter different organs and provoke inflammatory processes that may cause pain in specific areas, such as headache via the nerves [37], chest pain via the cardiovascular system [40], or abdominal pain via the gut-lung axis [44]. Significantly, various pain syndromes may point towards similar underlying pathology and have a cumulative effect that subsequently contributes to widespread pain [47, 76]. Considering that the causal mechanisms are still hypothetical and pain is highly complex, more research is warranted to confirm.

Prior literature reported that individual and household characteristics could affect susceptibility to the health effects of using polluting fuels [53, 77, 78]. The current study also observed inconsistent associations between cooking with solid fuels and the number of pain sites across gender, currently working, smoking status, number of diseases, and household consumption expenditure. Compared with females, males were more susceptible to the effect of cooking with solid fuels on pain. One possible explanation is that females are more often exposed to solid fuels for cooking [79] and more frequently display pain symptoms [80]; this may mix the findings and contribute to a statistical difference. We also found that those not having jobs were more likely to be affected by the effect of cooking with solid fuels on pain than those having jobs, potentially due to having meals at home more frequently. Moreover, our study indicated that ever or current smokers were more likely to be affected by the fuel-pain association than nonsmokers. This may also partly explain the observed sex difference, with smokers almost being males. A smoking difference may be due to the additive effect of harmful pollutants. Studies found that smoking increased the odds of chronic pain onset and progression [81], and recurrent nicotine exposure resulted in physiological tolerance, leading to reduced pain inhibition and more excellent pain perception [82]. Strangely, those with the highest household consumption expenditure tended to be more susceptible to the hazardous effects of using solid fuels. As household consumption expenditure reflects household income to some extent, this finding may be attributed to the higher impact of income on pain than using solid fuels [83], with the group with the lowest level of household consumption expenditure having the most pain sites among those using clean fuels. Similarly, the more prominent association in those with one disease rather than those with multimorbidity or none may also occur due to the more significant effect of pain-comorbidity on pain [84]. A study showed that the number of pain sites tends to be stable in the general population [47]; thus, for those with multimorbidity that had the most pain sites when using clean fuels, it is hard to develop more sites.

Pain has often been deemed a state that will be stable, improved, or worsening; thus, repeated data measurements were used to investigate the development and maintenance of pain. Contrary to our expectations, using solid fuels was associated with a slower rate of pain site increase over seven years. Indeed, chronic pain tends to amplify over time, but the number of pain sites that most subjects can develop is stable [47]. This may explain the above findings. In the present study, those using clean fuels for cooking had fewer pain sites than those using solid fuels, which meant that the former was faster to reach the peak. On the other hand, those exposed to solid fuels were close to the top, which may limit the speed of pain development. The current study observed that the changes in pain in different groups from 2015 to 2018 tend to be similar, suggesting that other risk factors that did not occur may exist. Future studies can explore the trajectory of pain from a lifecycle perspective and a deeper understanding of the mechanisms of pain. Furthermore, when fuel switching was considered, pain trajectories became more complicated. The group transiting from clean to solid fuels had a faster rate of pain site increase than those from solid to clean fuels over seven years, further confirming the harmful effect of using solid fuels. Generally, household energy choices should change from solid to cleaner fuels as economic conditions improve. However, the combined use of multiple energy is still common in rural homes in China [85]. This phenomenon is attributed to poor economic conditions and is associated with specific cultural contexts [86]. China is in an earlier stage of epidemiological transition than high-income countries, which may make people’s behaviors in various patterns. Regardless, our data showed a nonconventional energy transition and highlighted its impact on pain; thus, it is important to explore further why such a strange way of switching fuels and take measures to help clean fuel spread.

The study’s strengths are the availability of data regarding fuel transition, repeated pain measurements, and systematic investigation of pain (number of pain sites and locations). However, there are several limitations to this study. First, the major limitation of this study was based on the self-reported information on body pain. To confirm the robustness of our findings, we used pain data from three repeated measurements (in 2011, 2015, and 2018) and added the analyses on the relationship between household use of polluting cooking fuels and different regional pain sites to get more information about the underlying mechanism of pathogenesis. The etiology and prognosis of pain are complex and believed to involve interactions between biological, psychological, and social. Thus, future research should define the pain via doctors’ diagnoses or using a body map for participants to report the anatomical localization of their pain to understand the underlying pathogenic mechanism. Second, the pain assessment in the second survey (in 2013), which asked participants about yesterday’s pain, differs from other survey years. Thus, the information on pain in the first (in 2011), third (in 2015), and fourth (in 2018) surveys were included in the study, leading to a large time interval between each pain assessment. Collecting continuous annual data on pain, which can help describe the trajectory of pain development, is better. Third, this study relied on one simple item to measure pain and did not provide information on the causes of pain, limiting us from exploring more pain traits and getting clinical implications, such as pain level, duration, and frequency. Fourth, the follow-up time in our study was relatively short, so it was hard to depict the life-course pain trajectories. For example, our study indicated that the number of pain sites might increase with time and stabilize at a level, but we failed to observe the subsequent development of pain. More longitudinal studies with longer observation times are needed to confirm our findings. Fifth, a few important risk factors were not assessed in this study, including ergonomic postures, stove types, duration of use, and outdoor air quality due to limited data. These risk factors need further verification in future research.

5. Conclusions

In summary, this study provided new evidence linking using solid fuels for cooking with more pain sites, but a slower rate of pain sites increases. Males, smokers, those not having jobs, those with the highest household consumption expenditure, and those with one chronic disease responded more strongly to the harmful effects of indoor air pollution. When fuel switching was further considered, our study observed a nonconventional energy transition from clean to solid fuels and highlighted its impact on pain. Compared with those who persistently used clean fuels for cooking, the number of pain sites increased by 10% in participants who were transiting from using solid to clean fuels, by 21% in those transiting from cooking with clean to solid fuels, and by 25% among those persistent using solid fuels for cooking. Deeper explores why such a strange way of switching fuels is important and takes measures to help clean fuels spread.

Data Availability

All data in CHARLS are maintained at Peking University and will be accessible publicly. The datasets generated during the current study were from the Harmonized CHARLS dataset and Codebook (http://charls.pku.edu.cn/).

Ethical Approval

The data collection in CHARLS was ethically approved by the Biomedical Ethics Review Committee of Peking University (approval number: IRB00001052-11015).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Authors’ Contributions

Yi Zhu was responsible for conceptualization, methodology, formal analysis, writing the original draft, reviewing, and editing. Lijin Chen was responsible for formal analysis, writing the original draft, reviewing, and editing. Honghong Feng was responsible for formal analysis, writing the original draft, reviewing, and editing. Esthefany Xu Zheng was responsible for writing, reviewing, and editing. Yixiang Huang was responsible for conceptualization, supervision, writing, reviewing, and editing. Yi Zhu and Lijin Chen contributed equally as co-first authors.

Supplementary Materials

A series of exploratory and sensitive analyses were performed to assess the robustness of the primary findings. The details of these results were provided in the supplementary part. We also briefly introduced the GEE model. The title of the additional content is as follows. (1) Supplementary Methods. (2) Table S1: model selection. (3) Table S2: correlation matrix. (4) Table S3: distribution of cooking fuel types and fuel transitions from 2011 to 2018. (5) Table S4a: interaction effects of fuel types and potential moderators on body pain sites. (6) Table S4b: interaction effects of fuel types and potential moderators on body pain sites (continued). (7) Table S5: distribution of the number of pain sites from 2011 to 2018 by cooking fuel types. (8) Table S6: baseline characteristics across the study sample and excluded individuals. (9) Table S7: association of fuel types and fuel transitions with body pain sites in imputed dataset. (10) Table S8: association of fuel types and fuel transitions with short-term trajectories of body pain in the imputed dataset. (11) Table S9: association of fuel types and fuel transitions with body pain sites excluding those having chronic diseases at baseline (). (12) Table S10: association of fuel types and fuel transitions with short-term trajectories of body pain excluding those having chronic diseases at baseline (). (13) Table S11: association of fuel types and fuel transitions with body pain sites with additional adjustment for physical activity at baseline (). (14) Table S12: association of fuel types and fuel transitions with short-term trajectories of body pain with additional adjustment for physical activity at baseline (). (15) Table S13: association between fuel types and body pain sites through different working correlation matrix. (16) Figure S1: subgroup analyses of the association between cooking fuels and the number of body pain sites. (17) Figure S2: association of cooking fuels with specific pain sites. (Supplementary Materials)