Introduction

In the context of global aging, the number of people aged 65 and above is rapidly increasing, and the prevalence of chronic diseases in the elderly population is rising sharply1. What is more concerning is that the majority of elderly patients with chronic diseases often suffer from multiple chronic diseases concurrently2, as, chronic diseases are often characterized by complex etiologies, insidious onsets, long duration and difficulties in treatment. The World Health Organization (WHO) defines the presence of two or more chronic diseases as chronic disease comorbidity3. Compared to having a single chronic disease, comorbidity not only seriously reduces the quality of life of the middle-aged and the elderly, but also increases the economic burden on families and society4. In China, studies have shown that by 2018,the rate of chronic disease comorbidity among middle-aged and the elderly had reached 55.77%5. making it a pressing public health challenge.

The high incidence of multiple chronic diseases in China’s rural areas echoes the challenge of global aging. The 'healthy ageing framework advocated by the World Health Organization emphasizes that integrating lifestyle interventions, such as physical activity, with primary care is a core strategy to address the burden of comorbidities. Emerging evidence suggests that higher-intensity physical activity may impose stronger physiological adaptations that are critical for mitigating comorbidity risk. For instance, high-intensity activity (e.g., aerobic exercise at 70–85% of maximal heart rate) enhances insulin sensitivity more effectively than moderate-intensity exercise in patients with metabolic disorders6,7, which is particularly relevant for hypertension-diabetes comorbidity prevalent in rural elderly populations. In response, many countries have adopted to promote healthy lifestyles, including physical activity interventions, as a novel approach to combat the rising prevalence of comorbidity8. While prior research has established the benefits of exercise for single chronic diseases (e.g., hypertension or diabetes), evidence on how physical activity intensity specifically impacts comorbidity patterns—such as the co-occurrence of hypertension and diabetes—remains sparse, particularly in rural settings. Despite significant investment in biomedical technology, these efforts have had limited success in rural areas, where awareness, treatment, and control rates of comorbidity remain substantially lower than urban counterparts. Recent national health expenditure reports indicate that rural chronic disease screening coverage remains below 40% despite a 12% annual increase in biomedical technology investments9.This calls for a shift in approach—one that emphasizes healthy aging, maintaining basic health, and improving the quality of life for the elderly through multidisciplinary, multimodal strategies that integrate healthcare, physical activity, and the establishment of healthy lifestyle systems10.

With the popularization of the concept of “Exercise is Medicine”, intervening in the intensity of physical activity which can influence the prevalence of comorbidity has become one of the basic public health strategies in a majority of countries, because Compared to medication or surgical interventions, physical activity interventions do not require additional medical costs (such as drug purchases or hospitalizations) and may indirectly reduce long-term health care costs by reducing chronic disease complications11. Numerous studies have shown that appropriate physical activity contributes to improved health outcomes in comorbid populations. However, these studies have mainly focused on the effects of physical activity in the elderly on a particular chronic disease, with limited evidence exploring the relationship between physical activity intensity and comorbidity12. Previous studies focused on the impact of physical activity on a single chronic disease, such as hypertension or diabetes, but there is still a lack of large-scale empirical evidence on the association between different intensity of physical activity and multiple chronic diseases, such as hypertension combined with diabetes. its adoption and feasibility in elderly rural populations—often constrained by limited infrastructure or cultural preferences for low-intensity activities like walking—remain underexplored. This study is the first to systematically compare the difference between high intensity and low intensity activity in rural elderly population to fill this gap.

Given this gap, public health policymakers tended to take physical activity interventions as an auxiliary means of controlling comorbidity in the elderly, rather than recognizing its pivotal role in the prevention and treatment. Surveys of Chinese rural health administrators suggest that 68% prioritized medication adherence over structured PA programs13. This highlights a potential underutilization of PA in current policies. This study aims to address this gap by examining the relationship between physical activity intensity and comorbidity prevalence employing a propensity score matching methodology. By exploring the effect of physical activity intensity on comorbidity, while minimizing the influence of potential disturbance variables, this study provides a more accurate assessment of the net benefits of different physical activity levels. This study will provide valuable theoretical insights for health promotion among the elderly in rural areas and contribute to the development of a new public health service model.

Methods

Procedures and study participants

From July 1 to August 31, 2023, we conducted a cross-sectional survey in Jia County, Henan Province (one of the national demonstration areas for comprehensive prevention and control of chronic diseases). We conducted the survey simultaneously in batches, centered on the townships, with villages or communities as the units, covering 377 villages or communities across 15 townships in Jia County. Our research team consists of 70 highly trained students with a background in public health and 4–5 professional public health persons from each township. During the survey, each of the 377 villages or communities dispatched 1 to 2 general practitioners to assist throughout the investigation. Quality control personnel regularly conduct quality control check on the questionnaire to ensure the reliability and validity of the collected data. These examinations included a random examination of questionnaire records and a feedback session with the assessment team. Utilizing whole cluster sampling, we assessed hypertensive individuals aged 65 and above registered under the National Basic Public Health Hypertension Management program. Participants gave written informed consent, answered an interview questionnaire, and underwent a physical examination by trained health personnel. A total of 21,379 residents with diabetes and hypertension were included. Inclusion criteria: (1) age ≥ 65 years, clear consciousness, no communication barriers; (2) no motor disabilities, reaching the village clinics independently; (3) cooperating with the completion of questionnaires and physical examination; (4) permanent residents of the study site. Exclusion criteria: (1) those who suffered from other serious diseases and could not cooperate in completing the study N = 479 (2) those who were missing basic information N = 3507 (3) those who could not express themselves N = 598. 16,795 participants were finally included in the study. All participants provided informed consent before data collection and all data is kept strictly confidential.

Assessment of covariates

A structured questionnaire was used to collect information through face-to-face interviews, namely (1) basic demographic characteristics: sex, age, household registration, education level, marital status, occupation, annual household income, and self-reported health status; (2) disease status: Seven chronic diseases were surveyed (hypertension, diabetes, bone and joint diseases, COPD, digestive diseases, heart disease, and stroke). Comorbidity in this study was operationally defined as the concurrent diagnosis of hypertension and diabetes. In addition, we used the Age-adjusted Charlson Comorbidity Index to assess the severity of comorbidities by stratification (see Table TS9 for score details). The higher the score, the higher the severity of comorbidities14. In this study, 1 represents no comorbidity: Hypertension or Diabetes alone; 2 represents comorbidity: Hypertension and diabetes. Other chronic diseases (e.g., COPD, heart disease, etc.) were exclusively used to calculate the Age-adjusted Charlson Comorbidity Index (ACCI). The ACCI score reflects the cumulative burden of non-focal diseases but was independent of the primary comorbidity outcome. (3) physical examination indicators: height, weight, BMI was calculated as weight in kilograms divided by squared height in meters (kg/m2) (4) Lifestyle indicators: Physical activity intensity: referring to the International Physical Activity Questionnaire (IPAQ)15, which is a more widely used questionnaire for measuring the physical activity level of adults worldwide, and has good validity and reliability in our empirical studies. We used rigorous quality control procedures (e.g., trained interviewers, repeated verifications) aimed to minimize reporting errors. The situation of exercises was that an individual’s total weekly energy expenditure was calculated based on the number of days per week that the individual was engaged in low/moderate/high levels of physical activity, the length of the day, and the metabolic equivalent (MET) of each level of activity (low MET was assigned 3.3, moderate MET was assigned 4.0, and high MET was assigned 8.0). The formula is 1 week of exercise energy expenditure = MET corresponding to exercise level × duration of daily activity × number of activity days per week. And according to the International Physical Activity Questionnaire (IPAQ), the physical activity level of middle-aged and the elderly were categorized into high physical activity (assignment: ≥ 1500 MET-min/w = 3) medium physical activity (assignment: ≥ 5 days per week, each time ≥ 30 min = 2) and low physical activity (assignment: dissatisfy the medium–high physical activity = 1).

Statistical analysis

Data were entered using Epi Data 3.0. Counting data were expressed as percentages. Measurement data that conformed to normal distribution were described by mean ± standard deviation, and when they did not conform to normal distribution, they were expressed by quartiles (P25,P75). The χ2 test was used, and multivariate analyses were performed using logistic regression models. Using these to examine factors associated with the prevalence of comorbidity among the elderly while adjusting for potential confounders. In addition, we conducted a multiple linear regression analysis using SPSS to examine the association between the Age-adjusted Charlson Comorbidity Index (ACCI) and physical activity intensity. PSM was performed to minimize potential confounding bias by Logistic to select the sets of covariates whenever possible (p < 0.05). A one-to-one nearest neighbor matching algorithm was applied with a caliper value of 0.03 in our study. The Probit regression model was performed to estimate the propensity scores of participants with moderate- to low-intensity physical activity and high-intensity physical activity. Finally, 3387 were selected from 16,795 participants, and the differences in comorbidity between high-intensity physical activity and other levels of physical activity intensity were compared by using 1:1 propensity score matching, after controlling for confounders. A bilateral test was utilized, and p < 0.05 was considered statistically significant.

Results

General characteristics of participants

A total of 16,795 participants were included in this study, aged from 65 to 104, with a mean age of 73.1 ± 5.58. Among them, 7534 (44.9%) were male and 9261 (55.1%) were female. The elderly with low-intensity physical activity were 2944 (17.5%), with moderate-intensity physical activity were 10,464 (62.3%), and with high-intensity physical activity were 3387 (20.2%). 7434 (44.3%) suffered from one chronic disease, and 9361 (55.7%) suffered from comorbidity. Demographic characteristics are shown in Table 1.

Table 1 Sociodemographic characteristics.

Multivariate logistic regression analysis on the prevalence of comorbidity

Multivariate logistic regression analysis was performed with whether the participants had comorbidity as the dependent variables (assignment: 0 = no 1 = yes), and intensity of physical activity as the categorical independent variables. The result showed that Participants engaging in high-intensity physical activity had a 14.3% reduced odds (OR 0.857, 95% CI [0.812–0.905], p < 0.001) of comorbidity compared to those in the low-intensity group. The absolute risk reduction (ARR) was 13.2%, indicating that one comorbidity case could be prevented for every 7.6 elderly individuals transitioning to high-intensity activity. In addition, gender, age, income, occupation, and educational background was also influencing factors for comorbidity among the elderly in rural areas (p < 0.05), as shown in Table 2, Fig. 1, TS10.

Table 2 The associations between Physical activity and comorbidity prevalence.
Fig. 1
figure 1

Forest plot of regression results for physical activity.

Differences in the prevalence of comorbidity among the elderly after PSM

Moderate-intensity physical activity was not associated with the prevalence of comorbidity in the elderly. Participants who were in moderate-intensity physical activity were combined with those who were in low-intensity physical activity for better propensity score matching. In total, 6774 participants from 16,795 participants were matched using PSM. The intensity of physical activity was significantly associated with the prevalence of comorbidity after PSM for gender, age, income, occupation, marital status, education, and BMI (p < 0.001). Based on the balanced samples, the differences in the prevalence of comorbidity between low- and moderate-intensity physical activity participants and high-intensity physical activity participants were shown in Fig. 1, which suggested that the prevalence of comorbidity in high-intensity physical activity participants (45.4%, 95% CI [0.531–0.65]) was still lower than that of the low- and moderate-intensity physical activity population (58.6%, 95% CI [0.872–1.03]) by 13.2% (p < 0.05). After adjusting for potential confounders in the multiple linear regression analysis, physical activity intensity showed a coefficient of − 0.82 (p < 0.001) in relation to the Age-adjusted Charlson Comorbidity Index (ACCI). In other words, as physical activity intensity increases, the ACCI score decreases, indicating a negative correlation between physical activity intensity and comorbidity burden (as shown in TS8) Notably, this finding is consistent with the results obtained from the PSM analysis. Figure 3 illustrates sex- and age-stratified differences in comorbidity prevalence after propensity score matching. While high-intensity activity males aged 65–70 showed higher comorbidity rates than females (54.2% vs. 48.9%, p < 0.05), this pattern reversed in older cohorts: females exhibited progressively greater comorbidity risks than males from ages 71–75 (58.8% vs. 52.1%), 76–80 (61.3% vs. 55.7%), and ≥ 81 years (63.1% vs. 57.4%, all p < 0.01).As shown in Figs. 2 and 3. (Balance test in TS2 and TS4).

Fig. 2
figure 2

Prevalence of comorbidity between high-intensity physical activity and low-intensity physical activity participants before and after PSM.

Fig. 3
figure 3

Prevalence of comorbidity by age group and sex in rural areas after PSM.

Discussion

Physical activity status of the elderly in rural areas

The results of this study showed that the physical activity level of the elderly in rural areas was low. Only 20.2% of the elderly performed high-intensity physical activity every week, while 17.5% of the elderly did not perform physical activity every week or only performed low-intensity physical activity. 62.3% of the elderly engaged in medium-intensity physical activity. This phenomenon may be attributed to a multitude of factors, including aging, declining physical functions, and a sedentary lifestyle16, which makes the willingness of the elderly to engage in physical activity gradually decrease, and most of the elderly are only willing to do simple physical activity, such as walking, housework, and other activities. Only a minority of the elderly are willing to jogging, cycling, and other high-intensity physical activities. Walking is a relatively safer and easier way to exercise for the elderly, which may be related to the lack of knowledge and guidance on other forms of physical activity17.

Some studies have shown that the situation of chronic disease comorbidity among the elderly in China is rigorous. With the aggravation of population aging, the rate of chronic disease comorbidity among the middle-aged and elderly may continue to increase in the future18. Extensive studies have shown that there is an effective relation between the amount of physical activity and health status, with increasing physical activity will bring about cumulative health benefits19. Regular physical activity plays a positive role in preventing disease and promoting health by reducing mortality, lowering the risk of comorbidity, and improving the quality of life20. Concurrently, appropriate physical activity is beneficial to the physiology of the elderly, the functional activity of the brain and nervous system, and the quality of sleep21, and helps to improve physical and mental health. Therefore, village clinics could provide health education to the elderly about the benefits of increasing physical activity for their physical and mental health, as well as education and guidance on safety and precautions for physical activity. The elderly may benefit from active learning of simple motor skills, which could enhance their interest in exercise and facilitate the formation of a regular exercise daily routine22.

While our findings highlight the benefits of high-intensity physical activity in reducing comorbidity among ambulatory elderly populations, they may not generalize to individuals with severe mobility limitations. Policymakers should consider supplementary interventions tailored to non-ambulatory elderly populations, who face distinct barriers to physical activity and higher chronic disease risks. Such as accessible community exercise programs and home-based rehabilitation.

Comorbidity status of the elderly in rural areas

This study found that the prevalence of comorbidity among the elderly in rural areas was 55.7%. This result is slightly higher than the co-prevalence rate of chronic diseases in the elderly obtained by Xu Xiaobing et al. based on CHARLS data analysis (45.92%)23 and lower than the result obtained by Sun Ying et al. using the data of elderly inpatients in a hospital (83.74%)24. However, it is lower than the results of foreign studies, for instance, the U.S. National Health Interview Survey (NHIS) showed that the prevalence of comorbidity in the elderly over 65 years old was 63.7%25, and the comorbidity rate of the elderly of advanced age in the community of Kawasaki City, Japan was 94.7%26. Although the prevalence of comorbidity varies from study to study, it indicated that the prevalence of comorbidity was higher among the elderly in rural areas, and the situation is not optimistic. This phenomenon that the risk of comorbidity among the elderly is increasing may be attributed to the fact that the proportion of the elderly population is rising as the aging process is accelerating, and their physical functioning is declining as well as the body’s susceptibility to disease is increasing27. Simultaneously, the elderly in rural areas generally have low education levels, limited awareness of healthy lifestyles, insufficient physical exercise, unhealthy diet, with smoking and excessive drinking, and other bad lifestyles are more common28, increasing the risk of the higher prevalence of comorbidity of the elderly in rural areas. Moreover, economic reasons may also be one of the reasons for the increasing prevalence of comorbidity among the elderly in rural areas. The elderly in rural areas have a single source of income, mainly relying on renting out their cultivated land or selling their food crops as well as government subsidies to obtain income. Chronic comorbidity requires long-term medication and regular checkups, but most of the elderly can only choose cheap therapeutic drugs, and it is even harder for them to undergo regular checkups, which makes it difficult to ensure the effectiveness of drug treatment29. Moreover, due to the age-related decline in self-control and memory, it is easy for the elderly to forget to take their medication, which undoubtedly increases the burden of disease on the elderly30.

Factors influencing the decrease in the prevalence of comorbidity in participants with high-intensity physical activity

This study found that the prevalence of comorbidity decreased by 13.2% in participants who performed high-intensity physical activity after matching the variables of gender, age, occupation, education, income, marital status, and BMI. it means that with the increase of intensity of physical activity, the probability of comorbidities decreases in the elderly. This decrease may be related to the fact that the level of physical activity in the elderly has an impact on bone mineral density. Studies have found that the elderly with insufficient physical activity mostly have abnormal bone density compared with those with sufficient physical activity31. Previous studies have established a close relationship between physical activity levels and bone density, with higher activity levels being associated with a lower risk of low bone density32. The higher the physical activity levels are, the less prone to low bone density problems. Studies have shown that gradually increasing daily activity levels can prevent the decrease in bone density for postmenopausal elderly females; increasing physical activity and exercise can improve bone health for both elderly males and females33,34. Therefore, elderly individuals with insufficient physical activity should be encouraged to exercise more to maintain bone density and reduce the risk of fragility fractures. However, this does not mean that the elderly can carry out high-intensity physical activities without concern, and the exercise intensity and exercise time should be reasonably arranged according to their actual situation to avoid the possibility of injury caused by exercise.

Gender differences in comorbidity prevalence in rural areas

Our analysis revealed sex-specific differences in comorbidity prevalence across age groups (Fig. 3). High-intensity physical activity reduced comorbidity prevalence overall, but age- and sex-stratified analyses revealed nuanced patterns. While males aged 65–70 engaging in high-intensity activity paradoxically exhibited higher comorbidity rates than females (54.2% vs. 48.9%), this trend reversed sharply after age 70, with older females (≥ 71 years) experiencing disproportionately greater risks than males in the same age groups (p < 0.01).This divergence likely stems from sex-specific biological trajectories (e.g., postmenopausal metabolic shifts in women) and rural occupational exposures. Which is similar to the findings of Yuan Yupeng35. In rural areas, where physical labor is predominant, men are mainly engaged in higher-intensity agricultural activities, such as construction, and heavy lifting, while women often perform lower-intensity housework, resulting in lower overall activity for women. In addition, we also found that the prevalence of comorbidity among females who engaged in high-intensity physical activity increased with age. The reason for this may be related to female-specific physiological factors36. Domestic and foreign studies consistently agreed that this phenomenon may be related to female-specific physiological factors and longer average life expectancy. On the one hand, Females are at increased risk of chronic comorbidity due to pregnancy, menopause, and so forth, resulting in poorer health for females compared to males. On the other hand, female bear multiple responsibilities, including caring for the family, providing for elderly parents, and farming, while the medical services available to them are more limited37. As women age, it is no longer appropriate for them to engage in high-intensity physical activity. It is recommended that females actively participate in other types of physical activity, such as square dancing, traditional Chinese medicine health exercises, Baduanjin exercise, and other moderate-intensity exercise programs, to increase their level of physical activity38.

Occupational differences in the prevalence of comorbidity in rural areas

The results of the study showed that different occupations also had an impact on the prevalence of comorbidity among the elderly in rural areas. In the present study, 43.7% of the elderly people whose occupation was farming did not suffer from comorbidity, 45.1% of the elderly whose occupation was housework did not suffer from comorbidity, and 46% of the elderly with other occupations did not suffer from comorbidity. It indicated that farming in rural areas could be used as a way of high-intensity physical activity for the elderly, based on which they can appropriately participate in low-intensity physical activities to avoid the harm caused by overwork to the body. It is also necessary to pay more attention to the non-farmers39 and to increase the activity time and activity intensity of this group accordingly. For people engaged in non-agricultural labor, low-intensity physical activity should be supplemented with high-intensity physical activity with a weekly metabolic equivalent of ≥ 1500 MET-min/w. Surveys have found that non-farmers in rural areas generally spend a lot of time sitting still, and a large number of studies have shown that long sitting time is often significantly correlated with low levels of physical activity40. Health interventions may be necessary to encourage non-farmers to engage in more physical exercise, reduce sedentary behavior, and adopt healthier lifestyles41.

Strengths and weakness of this study

There are several strengths in our study. First, it is a big survey conducted among the elderly across an entire county, allowing for more comprehensive data on the aggregation of diseases in rural areas. Second, the large sample size and scientific methodology, particularly the use of PSM, we employed PSM analysis to exclude confounding variables when exploring the association between physical activity and comorbidity prevalence. This study goes beyond previous research in its scope and depth of analysis. Third, it specifically examined the effects of high-intensity physical activity on comorbidity, expanding the research on the relationship between physical activity and comorbidity. Last, The ACCI served as a supplementary tool to control for variations in overall disease burden arising from conditions other than hypertension and diabetes. For instance, participants with higher ACCI scores (indicating severe comorbidity beyond hypertension/diabetes) showed lower adherence to physical activity interventions. This stratification underscores that our findings on physical activity and hypertension-diabetes comorbidity remain statistically robust even in participants with additional chronic diseases. However, some limitations exist. First, while this study identified associations between gender, age, physical inactivity, and chronic diseases, it could not establish causality. Across-sectional nature of the study design may be another limitation that precluded the establishment of cause-and-effect associations. Second, the self-reported nature of physical activity data via the IPAQ short form introduces potential social desirability and recall biases. Participants may have overestimated their activity levels due to a tendency to report socially acceptable behaviours or inaccurately recalled the duration and intensity of their physical activities over the past 7-day period. These biases could lead to misclassification of high- versus low-intensity activity groups, thereby affecting the validity of the observed associations. The IPAQ’s 7-day recall period may also fail to capture long-term physical activity patterns, which could influence chronic disease outcomes. These limitations highlight the need for future studies to incorporate objective measures and longer monitoring periods to validate self-reported data and strengthen causal inferences. Third, while PSM was used to control for known confounders, unknown covariates may still have influenced the results. unmeasured variables such as dietary habits, medication use, and broader social determinants of health were not considered, which could influence comorbidity prevalence. Future research should prioritize the use of objective measures (e.g., accelerometers) to validate self-reported physical activity data and address these biases. Longitudinal designs are needed to clarify temporal relationships and disentangle the effects of activity intensity on comorbidity. Fourth, while we employed whole-cluster sampling to reflect real-world population structures in rural healthcare settings, our analysis did not formally assess or adjust for intra-cluster correlation (ICC). Although propensity score matching (PSM) helped control individual-level confounders, residual clustering effects may persist if health behaviors or environmental factors shared within villages influenced outcomes. This limitation suggests that estimates of physical activity benefits could be marginally overprecise. Future investigations should incorporate multilevel modeling to disentangle individual- and cluster-level effects on comorbidity risk among rural elderly populations.

Conclusions and implications

This study found that healthier individuals might engage in more high-intensity exercise in rural areas. Besides, the prevalence of comorbidity is influenced by a variety of factors, including gender, age, and occupation. It is recommended that the farming population maintain a balance between high-intensity physical activity (at least three days per week for 63 min per time each session) and low-intensity activity to optimize health benefits while avoiding overexertion. Moreover, we may pay more attention to the non-farming population, encouraging an increase in both the intensity and duration of physical activity to mitigate the harmful effects of a sedentary lifestyle. Additionally, Given the heightened comorbidity risk among women—particularly those aged ≥ 71—tailored interventions should address age- and sex-specific barriers. While older women may benefit from community-based exercise promoting more frequent engagement in physical activities and the development of healthy exercise habits. For younger men (65–70 years), appropriate increases in the intensity or duration of physical activity should be considered.