Issue cover
Research Article | 29 May 2026

Early-life thermal stress and mid-growth light intensity as key environmental predictors of lameness in commercial Cobb 500 broilers under tropical closed-house production systems

Parita Noisamran ORCID , Wiriya Loongyai ORCID , Choawit Rakangthong ORCID , and Chaiyapoom Bunchasak ORCID Show more
VETERINARY WORLD | Article No. 29 | pg no. 2233-2251 | Vol. 19, Issue 5 | DOI: 10.14202/vetworld.2026.2233-2251
Cited by 0

Cite this Article

  • APA
  • MLA
  • Chicago
  • Vancouver
  • Harvard

                            
                        

ABSTRACT

Background and Aim: Lameness is a major welfare and economic concern in fast-growing broiler chickens, particularly under tropical intensive production systems where environmental stressors are pronounced. Despite its multifactorial etiology, large-scale field evidence quantifying the combined effects of housing microclimate on lameness in Southeast Asian closed-house systems remains limited. This study aimed to identify age-specific environmental predictors of lameness in commercial Cobb 500 broilers raised under tropical closed-house conditions.

Materials and Methods: An observational study was conducted using data from 296 broiler flocks (11.45 million birds) across five commercial farms in Eastern Thailand during 2023. Lameness incidence was derived from culling records. Environmental parameters included temperature, relative humidity, temperature–humidity index (THI), light intensity, and ammonia concentration. Statistical analyses comprised descriptive statistics, General Linear Model, Pearson’s correlation, and multivariable regression, with assessment of multicollinearity.

Results: The median lameness incidence was 1.00% (interquartile range: 0.76–1.38%), with significant variation among farms (p < 0.0001). Early-life thermal conditions were strongly associated with lameness, with day 1 maximum temperature showing a positive correlation (r = 0.54, p < 0.0001). Diurnal temperature variation on day 7 exhibited a negative association (r = −0.36, p < 0.0001), indicating the importance of thermal stability. Light intensity demonstrated the strongest association during mid-growth, particularly on day 14 (r = 0.56, p < 0.0001). Ammonia concentration became a significant factor during late production stages (day 28 onward; r = 0.31–0.41, p < 0.0001). THI values consistently exceeded recommended thresholds, indicating sustained thermal stress. The final multivariable model explained 42.9% of the variation in lameness incidence (adjusted R² = 0.416), identifying day 7 diurnal temperature variation and day 14 light intensity as independent predictors.

Conclusion: Lameness in tropical broiler systems is strongly influenced by age-specific environmental conditions. Early thermal stress and excessive mid-growth light intensity emerged as critical predictors, while late-stage ammonia accumulation further increased risk. These findings highlight the importance of precise microclimate management, particularly during early and mid-growth phases, to improve broiler welfare and reduce economic losses in tropical closed-house production systems.

Keywords: ammonia concentration, broiler chickens, environmental management, lameness, light intensity, temperature fluctuation, Thailand, welfare.

INTRODUCTION

In the past five decades, intensive genetic selection has substantially increased broiler growth rates, with modern commercial strains achieving daily weight gains of 65–100 g, more than three times those reported in the 1970s [1, 2]. Cobb 500 broilers typically reach approximately 2 kg of body weight at around 31 days of age, reflecting a high growth rate [3]. While this genetic progress has markedly improved production efficiency, it has concurrently increased susceptibility to metabolic disturbances, musculoskeletal disorders, and compromised bone integrity. Consequently, skeletal and locomotor disorders, particularly lameness, remain among the most critical welfare and economic challenges in the global broiler industry. Recent estimates indicate that leg disorders contribute to 15–30% of production losses due to impaired growth, increased culling, carcass condemnation, and reduced welfare, depending on the production system and genetic line [4, 5]. Birds affected by lameness experience reduced mobility, behavioral restriction, and chronic pain, raising major welfare concerns and negatively influencing consumer perception and compliance with international animal welfare standards [6, 7].

Lameness in broilers is commonly evaluated using gait scoring systems, most notably the Bristol gait scoring system (scores 0–5) [8] and the U.S. gait scoring system (scores 0–2) [9]. In commercial practice, a lameness incidence below 2% is generally considered acceptable, whereas higher values often indicate compromised management, environmental stress, or genetic predisposition [10]. Among the underlying pathological conditions, bacterial chondronecrosis with osteomyelitis (BCO) is one of the most prevalent causes of severe lameness in fast-growing broilers, particularly under high mechanical loads and microbial challenges [11]. Lameness development is multifactorial, involving both intrinsic factors (e.g., genetics and growth rate) and extrinsic factors, among which housing microclimate plays a central role in determining broiler leg health and welfare [5].

Environmental conditions within broiler houses, particularly temperature, light intensity, and air quality, are increasingly recognized as critical extrinsic drivers of locomotor health. Recent studies conducted in tropical and subtropical production systems have demonstrated that microclimate variables, such as ambient temperature, humidity, and thermal indices, strongly influence broiler welfare and performance [1215]. However, most existing research has evaluated these environmental factors independently, despite their frequent interaction under commercial conditions.

Despite the growing body of evidence on environmental influences on broiler welfare, several critical gaps remain. First, existing studies have largely focused on isolated environmental factors rather than evaluating their combined and interactive effects under commercial production conditions. Second, there is a lack of large-scale, multi-farm field studies conducted under tropical closed-house systems, where environmental stressors are more pronounced and management practices differ from temperate or open-sided systems. Third, limited information is available regarding age-specific susceptibility to environmental stressors, particularly during key developmental phases when skeletal growth and metabolic demands are highest. Furthermore, no published study has simultaneously quantified the integrated effects of temperature patterns, light intensity, and ammonia concentration on lameness incidence in commercial broiler production systems in Southeast Asia. This lack of comprehensive, real-world data limits the ability to develop precise, stage-specific management strategies to mitigate lameness.

Therefore, the aim of this study was to evaluate the combined influence of key environmental variables, including temperature, temperature-related indices, light intensity, and ammonia concentration, on lameness incidence in commercial Cobb 500 broilers raised under tropical closed-house conditions in Eastern Thailand. Specifically, this study sought to identify age-specific environmental predictors and critical risk windows associated with lameness development, using a large-scale, multi-farm dataset under standardized commercial management. The findings are intended to provide practical, evidence-based insights to improve microclimate management, enhance broiler welfare, and reduce economic losses in tropical poultry production systems.

MATERIALS AND METHODS

Ethical approval

This study used retrospective and prospective field data collected from commercial broiler farms operating under standard industry management practices. No experimental manipulation, invasive procedures, or handling of live animals was performed specifically for research purposes. Therefore, formal approval from an Institutional Animal Care and Use Committee (IACUC) or equivalent ethical review board was not required.

All procedures complied with the national livestock and animal welfare regulations of Thailand and adhered to internationally recognized animal welfare standards, including the principles of the Five Freedoms. The participating farms were certified under the LRQA Farm First assurance scheme, which ensures compliance with global animal welfare, biosecurity, and management guidelines.

Routine flock management, including health monitoring, culling decisions, and environmental control, was performed by trained farm personnel under the supervision of licensed veterinarians employed by the integrator company. Birds identified with moderate to severe lameness were humanely removed from the flock in accordance with established welfare protocols to prevent unnecessary suffering.

Data used in this study were obtained exclusively from routine farm records and automated environmental monitoring systems. No additional interventions were introduced for research purposes. Farm identities were anonymized (Farms 1–5) to ensure confidentiality and to comply with data protection policies agreed upon with the participating commercial entities.

The study design and data usage were reviewed and approved internally by the Department of Animal Science, Faculty of Agriculture, Kasetsart University, Bangkok, Thailand, to ensure compliance with ethical standards for the use of production data in research. The study adhered to accepted ethical principles for observational research involving animals and did not compromise animal welfare at any stage.

Study period and location

This observational study was conducted between January and December 2023 in the eastern region of Thailand.

Study design

Data on environmental conditions and poultry performance were gathered using Aviapp® (Huvepharma, Bangkok, Thailand), a cloud-based digital tool designed for monitoring poultry production. While Aviapp® can also track flock health indicators, this study focused solely on analyzing environmental factors (such as temperature, light intensity, and ammonia levels) and performance-related data. Although Aviapp allows direct mobile data entry, this study used a structured protocol. Farm workers recorded information manually, which house supervisors then entered into Aviapp®. The farm veterinarian and manager reviewed the data before final submission to the cloud system, ensuring multi-level quality control and minimizing errors.

At first, 296 commercial broiler flocks were evaluated. To maintain consistent data for statistical analysis, any flock with missing lameness records or essential environmental details was removed from the study. However, no flocks were excluded due to performance issues, disease status, or environmental conditions, which helped reduce selection bias.

The investigation covered five commercial broiler farms operating under standardized intensive production systems. All farms were managed by the same integrator and followed uniform feeding, health, and biosecurity protocols to ensure consistency in management. Each farm utilized closed-house systems equipped with double flooring and evaporative cooling mechanisms to maintain optimal thermal comfort. Houses were tunnel-ventilated under negative pressure, with wet-pad panels automatically controlled for temperature and humidity. House dimensions were 120 × 22.5 m for Farm 1, 125 × 36 m for Farm 2, 125 × 22.5 m for Farm 3, and 125 × 45 m for Farms 4 and 5. Final stocking density at slaughter was similar across farms, ranging from 25.5 to 26.54 kg/m², minimizing potential confounding effects on lameness outcomes. All houses used fresh rice hulls as litter, fully replaced each cycle, with mean litter allocation ranging from 5.17 to 5.98 kg/bird across farms.

A total of 296 flocks (representing 11,449,800 Cobb 500 broiler chickens) were included. The flock distribution by farm was as follows: Farm 1: 1,662,400 broilers (60 flocks); Farm 2: 2,772,500 broilers (60 flocks); Farm 3: 2,552,600 broilers (92 flocks); Farm 4: 1,665,500 broilers (36 flocks); Farm 5: 2,796,800 broilers (48 flocks). Broiler flocks on five farms included male, female, and mixed-sex groups. Males had higher lameness rates than females (1.18% vs. 0.96%, p = 0.0002). Because sex was closely linked to growth rate and body weight and was unevenly distributed, it was excluded from the final multivariable regression to avoid multicollinearity and was kept only for descriptive analysis.

All broiler chicks were sourced from multiple commercial hatcheries within the same integrator system using parent stock from its breeder farms. Mean breeder age was similar across farms (approximately 39.22–41.16 weeks). Because hatchery sources and parent stock origins overlapped, hatchery-specific analysis was not feasible in this study.

All diets were centrally manufactured at the same commercial feed mill and formulated (corn–soybean meal-based) to meet or exceed Cobb 500 nutritional recommendations. A four-phase feeding program was applied uniformly across farms: Starter 1 (crumble) from 1 to 15 days of age, Starter 2 (crumble) from 14 to 22 days, Grower (pellet) from 23 days to 7 days before slaughter, and Finisher (pellet) during the last 7 days before slaughter. Standardized levels of calcium, phosphorus, and vitamin D were applied across all farms, and no additional mineral or vitamin supplementation was used at the farm level. This centralized feed production ensured nutritional uniformity among flocks. Vaccination schedules and medication protocols were standardized across farms in accordance with the integrator’s veterinary guidelines.

Lighting programs were standardized across farms but implemented with age-dependent adjustments. All farms provided near-continuous light (23–24 h/day) during the first week after placement, followed by a gradual reduction in photoperiod during the growing period. From approximately 2 to 3 weeks of age, the daily light duration was reduced to 18–20 h with a consolidated dark period to promote rest. Thereafter, birds were maintained under 16–18 h of light per day until slaughter. Although the overall structure of the lighting program was uniform, the exact timing of light-on and light-off cycles varied slightly among farms (e.g., lights-on at 21:00–22:00 and dark periods ranging from 2 to 7 h/day depending on farm and age).

Ventilation and temperature management were controlled automatically in all houses using a Temp 1200 Plus controller. Target temperatures were set at 32–33°C at placement and gradually decreased with age to 26–27°C by 3–4 weeks and to 21–22°C by slaughter age. Evaporative cooling pads were activated when the house temperature exceeded the target by approximately 0.3–1.0°C. Minimum and maximum ventilation were achieved through staged operation of tunnel fans, ranging from 1–3 fans at minimum to 17–35 fans at maximum, depending on farm and house size. Heater set points were typically maintained 0.1–0.5°C below the target temperature during brooding.

While all farms adhered to the same integrator guidelines and used the same controller model, specific environmental settings varied slightly among farms. Differences were observed in cooling-pad differential settings, heater offsets, and fan staging patterns, primarily due to variations in house dimensions. Larger houses (Farms 4 and 5) operated a greater number of tunnel fans compared with smaller houses (Farms 1–3), resulting in farm-specific ventilation capacity despite a standardized overall microclimate management strategy.

Assessment of lameness

Daily flock inspections were conducted by trained farm personnel as part of standard management routines. Birds exhibiting visible lameness, such as limping, reluctance to move, or an abnormal gait, were manually captured and removed from the flock. The number of lame birds culled per day was recorded at the flock level and used to calculate the percentage of lameness using the following equation:

Lameness (%) = (Number of lame birds culled / Total birds in flock) × 100

Lameness was assessed using a standardized commercial visual scoring system adapted from established gait scoring methods. Birds showing moderate to severe gait impairment (approximately Bristol gait score ≥3) were culled for welfare reasons. Farm personnel received annual training from company veterinarians using reference images and videos, with periodic supervisory checks to reduce inter-observer variability.

Only birds with clear locomotor impairments leading to culling were classified as lame to ensure accuracy and reduce observer bias. Subclinical or borderline cases were excluded unless removal was required for welfare reasons.

Environmental temperature measurement

Ambient temperature and relative humidity (RH) were recorded using integrated digital sensors of the Temp 1200 Plus controller system installed at the central zone of each house at bird level. One sensor per house was used as part of the automated climate control system and was routinely calibrated by technical staff and positioned at bird level (approximately 30 cm above the floor). Although multi-point monitoring was not available, the central location was considered representative of average flock exposure under uniform tunnel-ventilated conditions.

Temperature was recorded in °C on days 1, 7, 14, 21, 28, 35, and the day before slaughter. For each sampling day, the daily maximum (Tmax) and daily minimum (Tmin) temperatures were recorded, and the diurnal temperature variation (ΔT) was calculated as:

ΔT = Tmax − Tmin

Relative humidity (RH, %) was recorded from the same house controller sensors at the corresponding sampling days. Descriptive statistics of RH at the time of the highest and lowest daily temperatures are provided in Supplementary Table 2. These temperature–humidity data were used to compute the temperature–humidity index (THI).

THI was calculated according to the following equation:

THI = (0.8 × T) + [RH/100 × (T − 14.4)] + 46.4

where T represents ambient temperature (°C) and RH represents relative humidity (%). For THI at high temperature, T = Tmax; for THI at low temperature, T = Tmin, using the corresponding RH recorded at the same time point.

Light intensity measurement

Light intensity (lux) was measured with a calibrated handheld digital lux meter between 09:00 and 10:00 h during the active photoperiod, when the artificial lighting systems were fully operational. Measurements were performed at bird level (approximately 30 cm above the floor) at the central location of the house to avoid edge or reflection bias.

Readings were obtained on days 1, 7, 14, 21, 28, 35, and the day before slaughter to represent changes in illumination throughout the production cycle. Light intensity was measured at two standardized locations in each house (central and rear zones). Because rear measurements were influenced by external light entering through ventilation openings and fan inlets, only central-zone readings were used for statistical analysis to minimize measurement bias.

Ammonia concentration measurement

Ammonia concentrations (NH3, ppm) were determined with a portable gas detector (ToxiRAE Pro, RAE Systems, Honeywell, San Jose, USA). The device was calibrated at the factory and then recalibrated annually by qualified technicians. Per the manufacturer’s guidelines, the instrument’s accuracy is within approximately ±2 ppm or ±10% of the reading, whichever value is higher.

Measurements were taken once per day at a central bird-level location (≈30 cm above the floor) between 09:00 and 10:00 h, reflecting typical exposure under tunnel ventilation. These values were used as the daily flock-level ammonia concentration. Instruments were calibrated before use. Multi-point sampling was not routinely conducted, which is a limitation. Ammonia levels above 10–15 ppm exceeded recommended welfare thresholds for broilers [3].

Statistical analysis

All statistical analyses were conducted using raw lameness percentages (0–100%) and performed using SAS University Edition (SAS Institute Inc., Cary, NC, USA) following the SAS/STAT® guidelines [16]. Prior to analysis of variance (ANOVA), lameness percentage was tested for normality using the Shapiro–Wilk test and for homogeneity of variance using Levene’s test. Analyses were performed using raw percentages (0–100%), and raw means are reported for interpretation.

A General Linear Model (GLM) was used to evaluate differences in lameness incidence among farms. In this model, the dependent variable was lameness percentage, and the fixed independent factor was Farm (five levels). The GLM model was expressed as:

Yij = μ + τi + εij

where Yij is the observation for treatment i, μ is the overall mean, τi is the treatment effect (farm), and εij is the random error term.

When significant effects were detected (p < 0.05), mean comparisons were performed using Duncan’s New Multiple Range Test (DMRT), applied as an exploratory post hoc comparison for descriptive purposes in a commercial field dataset.

Pearson’s correlation analysis was applied to assess the relationships between lameness incidence (dependent variable) and environmental parameters (independent variables: maximum temperature, minimum temperature, diurnal temperature variation, light intensity, and ammonia concentration). The Pearson correlation coefficient (r) was calculated as:

The significance of each correlation coefficient was tested using a t-test:

A p < 0.05 was considered statistically significant, and < 0.01 was considered highly significant. The interpretation of correlation strength followed the criteria of Best method [17] as shown in Table 1.

Table 1. Rule of thumb for interpreting the size of a Pearson’s correlation coefficient [17].

Size of correlationInterpretation
0.81 to 1.00 (−0.81 to −1.00)Very high positive (negative) correlation
0.51 to 0.80 (−0.51 to −0.80)Moderate positive (negative) correlation
0.21 to 0.50 (−0.21 to −0.50)Low positive (negative) correlation
0.01 to 0.20 (−0.01 to −0.20)Very low positive (negative) correlation
0.00Negligible correlation

To explore potential environmental predictors of lameness, a stepwise multiple linear regression was first performed using PROC REG in SAS. The dependent variable was flock-level lameness percentage (0–100%). Normality and homogeneity assumptions were evaluated prior to analysis, and the candidate independent variables included the highest and lowest temperatures, diurnal temperature difference, light intensity, and ammonia concentration on each measurement day. Variables were entered into and removed from the model based on statistical criteria (entry and stay p-value < 0.15), partial R², Mallows’ C(p), and biological plausibility.

A stepwise regression procedure was used to screen variables for candidate predictors. Due to multicollinearity and biological considerations, the stepwise output was not used as the final model but served only for variable preselection (Supplementary Table 3).

Based on a priori biological relevance and the correlation results, a final multivariable regression model was then fitted with five predictors: the highest temperature on day 1, the diurnal temperature difference on day 7, the light intensity on day 14, and the ammonia concentrations on days 28 and 35. Model fit was assessed using the coefficient of determination (R² and adjusted R²), F-test for overall significance, root mean square error (RMSE), and examination of residual plots to verify linear model assumptions.

For the multiple regression model, multicollinearity was assessed using variance inflation factors (VIF), with values <5 indicating acceptable levels. Residual independence was evaluated using the Durbin–Watson statistic (0.98), suggesting mild positive autocorrelation. Residual normality was examined using Q–Q plots, which largely followed the theoretical normal distribution with minor upper-tail deviations that were not considered to meaningfully affect model validity.

Power of analysis

A post hoc power analysis was conducted for the primary outcome (% lameness) comparing farms. Because farm differences were analyzed using a General Linear Model (GLM) within a one-way ANOVA framework, statistical power was calculated using G*Power 3.1 (F tests, fixed effects, omnibus, one-way ANOVA). Based on five groups (farms), a total sample size of 296 flocks, and α = 0.05, the effect size (Cohen’s f = 0.49) indicated a large effect. The achieved statistical power was >0.99, demonstrating that the study had sufficient power to detect differences in % lameness among farms.

Additionally, a retrospective power analysis for correlation analysis was performed using G*Power (two-tailed, α = 0.05). With a total sample size of 296 flocks, the study had >0.99 power to detect a moderate correlation (r = 0.30), indicating sufficient sensitivity to identify meaningful associations between environmental factors and % lameness.

RESULTS AND DISCUSSION

Lameness distribution and farm-level variation

Lameness information from five farms, including N, mean, SD, SEM, min–max, and 95% CI, is presented in Table 2, while the distribution across 296 flocks is illustrated in Figure 1. The median lameness incidence was 1.00%, with the 25th and 75th percentiles at 0.76% and 1.38%, respectively (Supplementary Table 1 and Figure 1). Most flocks exhibited lameness levels below 1.5%, although a small proportion exceeded 2.0%.

Table 2. Descriptive statistics of lameness incidence in broiler farms.

FarmnMean (%)SDSEMMin–max95% CI
1600.680.410.050.14–1.500.58–0.79
2601.480.470.060.68–2.511.36–1.60
3921.020.980.040.37–2.520.94–1.09
4361.560.530.090.57–3.161.38–1.74
5480.860.200.030.42–1.250.80–0.92

CI = Confidence interval, n = Number of flocks, SD = Standard deviation, SEM = Standard error of the mean.

Figure 1

Figure 1. Distribution of lameness incidence (%) across 296 broiler flocks. Figure 1 illustrates the frequency distribution of flock-level lameness incidence across 296 commercial broiler flocks. The histogram shows that most flocks had lameness levels between 0.75% and 1.50%, with the highest frequency around 1.00%. A small proportion of flocks exceeded 2.00%, indicating outliers with elevated lameness incidence. Overall, the distribution shows a slight right skew, suggesting that while most flocks maintained relatively low lameness levels, a subset experienced substantially higher lameness prevalence. This distribution supports the use of parametric statistical analyses after confirmation of approximate normality.

Compared with previous commercial broiler studies, in which gait abnormalities vary widely due to genetics, management, and environment, this study found a relatively low incidence. While intensively reared flocks often show clinically detectable lameness ranging from a few percent up to over 20%, the levels observed here are low for commercial settings [18].

At the farm level (Table 3), Farm 4 and Farm 2 showed the highest mean lameness rates (1.56 ± 0.09% and 1.48 ± 0.06%, respectively) and did not differ significantly from each other (p > 0.05), but were significantly higher than all other farms (p < 0.05). Farm 3 (1.02 ± 0.04%), Farm 5 (0.86 ± 0.03%), and Farm 1 (0.68 ± 0.05%) had lower mean values, each significantly different per Duncan’s test. Environmental or management factors at Farms 4 and 2 may increase locomotor disorder risk compared to the other farms

Table 3. Analysis of variance and Duncan’s multiple range test for lameness percentage by farm.

FarmnMean lameness (%) ± SEM95% CI
1600.68 ± 0.05ᵃ0.58–0.79
2601.48 ± 0.06ᵈ1.36–1.60
3921.02 ± 0.04ᶜ0.94–1.09
4361.56 ± 0.09ᵈ1.38–1.74
5480.86 ± 0.03ᵇ0.80–0.92

ANOVA: F(4, 291) = 45.88, p < 0.0001. Values with different superscript letters within the same column differ significantly (p < 0.05) according to Duncan’s multiple range test. ANOVA = Analysis of variance, CI = Confidence interval, n = Number of flocks, SEM = Standard error of the mean.

Even slight increases in lameness can lower profits by reducing growth, flock uniformity, and increasing mortality, culling, and carcass condemnation. Broiler welfare studies link leg weakness to economic loss and poorer welfare outcomes [18].

Furthermore, the use of a single centralized feed mill with standardized diet formulation minimized nutritional variability among farms, making dietary factors unlikely to confound the associations between environmental conditions and lameness.

Environmental temperature and lameness relationship

Environmental temperature greatly affects broiler performance and welfare, especially in tropical regions where heat often exceeds comfortable levels. This study assessed daily high, low, and ΔT (Table 4) to quantify thermal load in commercial settings, and correlation analysis revealed significant age-dependent links between temperature and lameness (Table 5). RH, which impacts heat load, was summarized (Supplementary Table 2) and factored into THI calculations (Table 6) to reflect combined temperature and humidity exposure.

Table 4. Descriptive statistics of environmental temperature by farm, including the highest temperature, the lowest temperature, and the diurnal temperature.

DayFarmnMean ± SEM (°C)SDMin–max (°C)95% CI (°C)
Highest temperature
115733.52 ± 0.070.5232.5–35.133.38–33.66
25235.47 ± 0.030.2434.5–35.935.41–35.54
38833.98 ± 0.020.2133.6–34.933.93–34.02
43534.83 ± 0.130.8032.9–36.734.56–35.10
54833.97 ± 0.060.3933.0–34.533.86–34.08
715832.37 ± 0.141.0330.9–36.232.10–32.65
26033.11 ± 0.040.3031.1–33.533.04–33.19
39233.23 ± 0.050.4432.0–33.833.14–33.32
43633.38 ± 0.271.5930.4–39.232.84–33.92
54832.86 ± 0.100.6630.6–33.532.67–33.06
1416031.45 ± 0.100.7429.7–33.231.26–31.64
26032.69 ± 0.131.0130.6–34.232.43–32.95
39231.93 ± 0.090.8529.7–33.631.76–32.10
43631.95 ± 0.160.9530.2–34.431.63–32.28
54831.96 ± 0.080.5330.8–32.831.81–32.12
2116030.79 ± 0.120.9228.8–32.830.55–31.03
26031.74 ± 0.090.7329.4–33.131.55–31.93
39230.56 ± 0.151.4426.0–33.330.26–30.86
43630.66 ± 0.241.4328.5–33.830.18–31.15
54831.11 ± 0.070.5029.9–31.730.96–31.25
2816030.23 ± 0.181.4024.2–32.829.87–30.59
26030.57 ± 0.070.5829.7–33.230.42–30.72
39229.62 ± 0.151.4125.0–33.429.33–29.91
43629.69 ± 0.181.1027.7–32.829.32–30.06
54830.32 ± 0.070.4829.1–31.830.18–30.46
3516029.95 ± 0.080.6328.6–31.229.79–30.11
26029.75 ± 0.141.0526.1–32.129.47–30.02
39228.79 ± 0.141.3325.1–31.228.52–29.07
43628.88 ± 0.150.8927.0–30.928.57–29.18
54829.01 ± 0.090.6427.1–30.828.82–29.19
Slaughter16028.72 ± 0.050.4228.1–30.128.61–28.83
26029.16 ± 0.141.0626.5–30.828.88–29.43
39227.93 ± 0.080.7626.8–30.327.78–28.09
43528.49 ± 0.201.1826.2–31.728.09–28.90
54828.40 ± 0.110.7826.0–30.628.18–28.63
Lowest temperature
115831.67 ± 0.090.6929.9–33.231.49–31.85
25233.49 ± 0.130.1833.0–33.933.44–33.54
38733.43 ± 0.040.3931.7–33.833.34–33.51
43533.63 ± 0.160.9531.2–35.433.31–33.96
54832.92 ± 0.080.5431.0–33.532.76–33.07
715829.58 ± 0.171.3026.5–32.529.24–29.92
25932.23 ± 0.020.1932.0–32.732.18–32.28
39232.28 ± 0.090.8530.2–33.332.10–32.46
43631.74 ± 0.211.2429.2–34.531.31–32.17
54831.56 ± 0.151.0429.0–32.631.25–31.86
1415828.80 ± 0.120.8927.2–31.128.57–29.03
26028.32 ± 0.110.8926.2–29.828.09–28.55
39228.97 ± 0.121.1326.0–31.228.73–29.21
43330.44 ± 0.191.0828.0–32.530.06–30.83
54829.84 ± 0.080.5328.7–30.729.68–29.99
2115728.11 ± 0.151.1225.1–30.527.82–28.41
26026.76 ± 0.120.9523.4–28.026.51–27.00
39227.88 ± 0.131.2822.5–29.627.61–28.14
43628.62 ± 0.140.8626.9–30.828.33–28.91
54829.02 ± 0.060.4127.8–30.028.90–29.14
2815727.30 ± 0.191.4223.3–32.526.92–27.67
26026.27 ± 0.080.5825.1–27.926.12–26.42
39226.82 ± 0.151.4623.0–28.726.52–27.12
43627.97 ± 0.231.3826.0–32.127.51–28.43
54828.07 ± 0.120.8225.0–29.727.83–28.31
3515826.41 ± 0.171.3124.3–30.526.06–26.75
26025.47 ± 0.292.2520.7–28.324.88–26.06
39225.88 ± 0.151.4221.5–27.525.58–26.17
43627.38 ± 0.191.1425.2–30.126.99–27.76
54827.31 ± 0.191.3223.9–28.726.93–27.69
Slaughter15825.76 ± 0.141.1024.2–29.025.48–26.05
26024.99 ± 0.211.6221.8–28.124.57–25.41
39225.96 ± 0.111.1021.1–28.425.73–26.18
43625.81 ± 0.211.2320.3–27.525.39–26.23
54827.10 ± 0.120.8524.2–28.626.85–27.34
Diurnal temperature
11571.85 ± 0.110.860.0–4.21.62–2.08
2521.98 ± 0.040.310.7–2.61.90–2.07
3870.51 ± 0.050.480.0–2.70.41–0.61
4341.24 ± 0.120.720.2–3.10.99–1.48
5481.05 ± 0.050.350.4–2.20.95–1.16
71582.73 ± 0.141.070.9–5.52.51–3.07
2590.92 ± 0.030.200.3–1.30.87–0.97
3920.95 ± 0.070.640.0–2.50.82–1.09
4351.72 ± 0.221.270.5–7.61.28–2.15
5481.31 ± 0.090.650.5–3.41.12–1.50
141582.68 ± 0.130.950.4–4.82.41–2.96
2604.38 ± 0.171.282.0–7.44.04–4.71
3922.98 ± 0.141.380.3–6.62.67–3.25
4341.49 ± 0.241.380.4–6.40.99–1.98
5482.13 ± 0.060.411.0–3.02.01–2.25
211572.69 ± 0.201.500.0–7.22.27–3.07
2604.99 ± 0.110.883.1–7.44.76–5.21
3922.68 ± 0.141.390.5–5.42.40–2.97
4362.04 ± 0.211.250.8–5.01.62–2.47
5482.09 ± 0.050.371.1–3.11.98–2.20
281572.83 ± 0.161.220.0–5.52.51–3.15
2604.29 ± 0.110.862.0–6.74.07–4.51
3922.80 ± 0.141.310.7–6.32.53–3.07
4361.72 ± 0.181.060.5–4.51.36–2.08
5482.25 ± 0.100.671.2–4.42.05–2.44
351583.53 ± 0.181.360.0–6.13.17–3.89
2604.28 ± 0.221.732.0–9.23.83–4.73
3922.92 ± 0.131.260.3–6.12.65–3.18
4361.50 ± 0.130.770.4–4.21.24–1.76
5481.70 ± 0.151.060.5–5.01.39–2.01
Slaughter1582.97 ± 0.131.010.0–4.92.70–3.23
2604.16 ± 0.151.142.2–8.63.87–4.46
3921.98 ± 0.111.010.3–6.31.77–2.19
4352.70 ± 0.291.730.8–7.82.10–3.29
5481.31 ± 0.090.620.4–4.11.13–1.49

CI = Confidence interval, n = Number of flocks, SD = Standard deviation, SEM = Standard error of the mean.

Table 5. Pearson’s correlation coefficients between the highest temperature, the lowest temperature, the diurnal temperature, and the total lameness incidence across all farms and rearing days.

VariableDaynr95% CIp-value
Highest temperature12800.540.46 to 0.62< 0.0001
72920.130.01 to 0.240.0271
142930.380.28 to 0.47< 0.0001
212940.220.10 to 0.320.0002
28294−0.03−0.14 to 0.090.6251
352940.06−0.05 to 0.170.3405
Slaughter2930.09−0.02 to 0.210.1062
Lowest temperature12800.500.41 to 0.59< 0.0001
72900.380.28 to 0.48< 0.0001
14289−0.08−0.18 to 0.040.1723
21291−0.12−0.24 to −0.010.0347
28291−0.16−0.27 to −0.040.0085
35292−0.05−0.16 to 0.060.3883
Slaughter292−0.10−0.21 to 0.020.1022
Diurnal temperature1278−0.01−0.12 to 0.110.9155
7292−0.36−0.46 to −0.26< 0.0001
142900.320.21 to 0.41< 0.0001
212930.260.15 to 0.36< 0.0001
282930.170.05 to 0.280.0042
352940.10−0.01 to 0.210.0807
Slaughter2930.330.22 to 0.42< 0.0001

CI = Confidence interval, n = Number of flocks, r = Pearson’s correlation coefficient.

Table 6. Descriptive statistics of temperature–humidity index at the highest and lowest temperatures by farm.

DayFarmnMean ± SEM (%)SDMin–max (%)95% CI (%)
Temperature–humidity index at highest temperature
115785.28 ± 0.272.0679.24–89.7384.73–85.83
25285.33 ± 0.050.3284.09–85.9885.24–85.42
39383.82 ± 0.413.9873.52–91.1483.00–84.64
43584.03 ± 0.251.4680.85–86.8283.52–84.53
54883.43 ± 0.231.5981.00–86.8882.97–83.89
715884.64 ± 0.312.3879.45–93.4884.02–85.27
26082.55 ± 0.050.4279.86–83.3282.44–82.65
39684.76 ± 0.302.8978.37–89.5284.17–85.35
43583.22 ± 0.422.5078.03–89.3282.36–84.08
54882.85 ± 0.352.4277.52–86.9282.15–83.55
1416083.70 ± 0.272.1179.89–88.8283.16–84.25
26082.43 ± 0.181.3779.56–84.8182.07–82.78
39683.59 ± 0.313.0279.32–90.5682.98–84.21
43581.70 ± 0.372.2076.60–85.9380.94–82.45
54882.37 ± 0.352.4577.40–86.5581.66–83.08
2116083.47 ± 0.241.8380.38–87.2283.00–83.94
26081.81 ± 0.181.4377.62–85.5081.44–82.17
39681.66 ± 0.323.1275.55–89.0481.03–82.29
43681.45 ± 0.533.1677.31–90.4380.38–82.51
54881.29 ± 0.251.7676.98–84.6880.78–81.80
2816082.15 ± 0.332.5372.41–87.2681.49–82.80
26081.09 ± 0.141.0679.02–85.2680.81–81.36
39680.94 ± 0.313.0872.62–89.5680.32–81.57
43680.11 ± 0.352.0974.58–86.5979.41–80.82
54880.88 ± 0.362.5076.01–84.7780.15–81.61
3516081.88 ± 0.221.7274.91–84.2781.43–82.32
26080.61 ± 0.322.5072.99–85.3279.96–81.25
39679.72 ± 0.323.1874.06–86.6579.07–80.36
43679.87 ± 0.301.8076.51–83.0079.26–80.48
54880.13 ± 0.332.3074.91–83.5579.47–80.80
Slaughter16080.45 ± 0.151.1678.18–82.8580.16–80.75
26081.19 ± 0.261.9972.85–84.0880.68–81.71
39678.86 ± 0.201.9774.84–85.0578.46–79.26
43678.86 ± 0.563.3570.23–85.7977.73–79.99
54879.80 ± 0.322.2074.04–82.3779.16–80.43
Temperature–humidity index at lowest temperature
115884.25 ± 0.221.6780.15–87.8583.81–84.69
25282.84 ± 0.030.2582.30–83.4082.77–82.91
38684.52 ± 0.312.9279.42–93.0883.90–85.15
43582.39 ± 0.211.2680.22–85.4481.95–82.82
54882.22 ± 0.211.4479.89–84.8381.80–82.63
715881.48 ± 0.332.5477.41–88.2080.81–82.15
26081.58 ± 0.040.3481.08–82.6881.49–81.66
39684.02 ± 0.312.9977.86–88.6783.41–84.63
43680.75 ± 0.331.9777.07–87.6180.08–81.41
54881.31 ± 0.362.4673.14–85.3080.59–82.02
1415880.49 ± 0.322.4274.78–84.3679.85–81.13
26077.03 ± 0.181.3973.82–79.5176.67–77.38
39679.69 ± 0.212.0774.56–83.0479.27–80.11
43479.69 ± 0.472.7269.42–83.5278.74–80.64
54879.67 ± 0.302.1075.00–82.6479.07–80.28
2115880.07 ± 0.322.4374.48–85.2879.43–80.71
26075.59 ± 0.241.8269.70–78.5975.12–76.06
39678.43 ± 0.222.1169.61–82.0578.01–78.86
43678.47 ± 0.331.9575.16–83.8677.81–79.13
54878.98 ± 0.231.5675.05–81.6178.52–79.43
2815878.55 ± 0.433.2471.69–88.2977.70–79.40
26075.62 ± 0.151.2072.80–78.4275.31–75.93
39476.99 ± 0.272.6270.88–81.9376.46–77.53
43677.91 ± 0.392.3174.17–85.3777.13–78.69
54878.15 ± 0.362.4670.82–81.6677.44–78.87
3515877.00 ± 0.352.6771.99–81.6276.30–77.70
26074.97 ± 0.513.9666.31–79.4773.95–75.99
39675.68 ± 0.252.4968.24–79.8875.18–76.19
43677.60 ± 0.271.6275.50–81.3477.05–78.14
54877.94 ± 0.402.7571.51–81.9977.14–78.74
Slaughter15876.13 ± 0.403.0568.63–82.5475.32–76.93
26075.00 ± 0.342.6569.80–80.0074.31–75.68
39675.85 ± 0.212.0568.69–81.3475.43–76.26
43575.27 ± 0.362.1167.12–78.1974.54–76.00
54878.02 ± 0.332.3071.55–80.9477.35–78.68

THI = Temperature–humidity index, CI = Confidence interval, n = Number of flocks, SD = Standard deviation, SEM = Standard error of the mean.

Across farms, mean maximum temperatures ranged from 33.52–35.47°C on day 1 and 27.93–29.16°C at slaughter (Table 4), with Farm 2 recording the highest values at both time points (35.47°C and 29.16°C). Mean minimum temperatures ranged from 31.67–33.63°C on day 1 and 24.99–27.10°C at slaughter. Mean ΔT ranged from 0.51 to 4.99°C, indicating substantial differences in thermal stability among farms. Farms 2 and 4 consistently showed higher temperatures, particularly during the early rearing period (days 1–21), with Farm 2 recording the highest day-1 temperature (35.47°C), followed by Farm 4 (34.83°C). The Cobb 500 Broiler Management Guide recommends air temperatures of 32–33°C at placement, decreasing to 27–29°C by days 8–14 and 25–27°C by days 15–21 [3]. In this study, all farms recorded temperatures 3–5°C above these targets, especially in Farms 2 and 4. High temperatures persisted beyond placement throughout brooding and early starter phases, resulting in prolonged exposure during critical skeletal development stages.

Correlation analysis revealed significant age-dependent links between temperature and lameness (Table 5). As shown in Figure 2, the highest daily temperature was strongly correlated with lameness on day 1 (r = 0.54, p < 0.0001), with weaker but still significant associations on days 7, 14, and 21, which became non-significant later. These results suggest that higher temperatures during brooding and early starter phases increase risk for locomotor disorders in tropical commercial settings.

Figure 2

Figure 2. Relationship between the highest temperature on day 1 and lameness incidence. Figure 2 illustrates the association between the highest ambient temperature on day 1 and flock-level lameness incidence. Each point represents one flock, and the dotted line indicates the fitted linear regression. A clear positive trend is observed, indicating that higher brooding temperatures were associated with increased lameness prevalence. Most observations clustered between 33.0°C and 35.5°C, while flocks exposed to temperatures above 35.5°C tended to exhibit higher lameness levels. This pattern supports the significant positive correlation observed in the statistical analysis.

Lowest temperature showed a dynamic pattern. Positive correlations with lameness were observed on days 1 and 7 (p < 0.001) and remained significant after FDR adjustment (Supplementary Table 4). After FDR correction, the day 21 association was not significant, whereas the negative correlation at day 28 remained significant (q = 0.03040), indicating a stronger late-stage effect. ΔT showed stage-specific patterns: a negative correlation on day 7 (r = −0.36; q = 0.00025) shifted to positive from day 14 onward, remaining significant at day 28 (q = 0.00871), suggesting that large day–night temperature swings may increase lameness risk.

THI was calculated to reflect the combined thermal load. THI remained relatively high throughout the production cycle, particularly during early and mid-growth (Table 6). Based on Cobb-Vantress guidelines, reference THI values are approximately 74–81 during brooding and 67–72 during finishing. Observed THI values were generally higher than these references, indicating sustained exposure to elevated heat–humidity load that may contribute to heat stress and impaired musculoskeletal development. Similar associations between high THI, reduced performance, and increased leg disorders have been reported in tropical production systems, supporting the broader relevance of these findings.

Table 7 shows no significant association between high-temperature THI and lameness (p > 0.05), whereas low-temperature THI showed a notable negative correlation during days 14–28 (r = −0.22 to −0.26, p < 0.01), emphasizing the importance of humidity and temperature effects at night. Early skeletal and thermoregulatory immaturity, combined with rapid muscle growth after day 14, increases broilers’ sensitivity to environmental conditions. Temperatures often exceeded Cobb-Vantress recommendations, particularly on Farms 2 and 4 during days 8–14, potentially impairing bone development and increasing the risk of leg disorders. Early-life heat exposure, unstable diurnal temperatures, and unfavorable nighttime THI were linked to lameness, with day 1 temperature and day 14 light intensity serving as practical early-warning indicators.

Table 7. Pearson’s correlation coefficients between THI at the highest and lowest temperatures and total lameness incidence across all farms and rearing days.

VariableDaynr95% CIp-value
THI at highest temperature1283−0.00−0.12 to 0.110.98
7293−0.11−0.23 to 0.010.05
14295−0.01−0.21 to 0.110.91
212960.03−0.08 to 0.150.57
28296−0.04−0.15 to 0.070.50
352960.00−0.11 to 0.120.94
Slaughter2960.07−0.05 to 0.180.26
THI at lowest temperature1275−0.04−0.15 to 0.080.55
72940.04−0.08 to 0.150.52
14292−0.26−0.38 to −0.15< 0.0001
21294−0.23−0.34 to −0.12< 0.0001
28292−0.22−0.32 to −0.100.0002
35294−0.05−0.16 to 0.070.41
Slaughter293−0.18−0.29 to −0.070.0019

THI = Temperature–humidity index, CI = Confidence interval, n = Number of flocks, r = Pearson’s correlation coefficient.

Light intensity and locomotor health

Light intensity varied substantially across the five farms and frequently exceeded the recommended management levels for broilers (Table 8). During the brooding period (days 1–7), Farms 2 and 4 recorded the highest mean light intensities (43.69 and 48.96 lux, respectively), whereas Farms 1, 3, and 5 maintained lower levels between 22.72 and 34.17 lux. From day 14 onwards, Farms 2 and 4 consistently sustained high light levels (>50 lux). In contrast, Farm 3 maintained the lowest levels throughout the production cycle, averaging 31–32 lux. By slaughter day, Farms 1, 2, and 4 still had light intensities exceeding 59 lux, while Farm 3 and Farm 5 remained around 32–33 lux. Notably, Farms 1, 3, and 5, which maintained comparatively lower light intensities, also exhibited lower lameness incidence, further supporting the link between excessive light exposure and leg health problems.

Table 8. Descriptive statistics of light intensity in broiler houses across farms.

DayFarmnMean ± SEM (lux)SDMin–max (lux)95% CI (lux)
116022.72 ± 0.141.0920–2622.43–23.00
26043.69 ± 0.735.6634.3–53.442.23–45.15
37631.24 ± 0.383.3225–3630.48–32.00
43648.96 ± 1.7410.4530.2–64.645.43–52.50
54834.17 ± 0.573.9820–3833.01–35.32
716025.45 ± 0.544.2221–4524.36–26.54
26058.32 ± 0.856.6150–78.256.61–60.03
37631.72 ± 0.373.2125–3831.00–32.46
43652.86 ± 1.046.2542.3–62.850.75–54.97
54834.08 ± 0.513.5420–3833.06–35.11
1416028.17 ± 0.866.7022–5026.44–29.90
26058.45 ± 0.715.5050–7157.03–59.87
37631.72 ± 0.363.1225–3631.01–32.44
43653.56 ± 1.016.0841.3–6851.50–55.61
54834.58 ± 0.584.0320–3933.41–35.75
2116034.58 ± 1.179.0322–5232.25–36.91
26059.43 ± 0.816.3151.5–7857.80–61.05
37231.61 ± 0.373.1426–3830.87–32.35
43655.99 ± 1.016.0544.2–66.453.94–58.04
54834.44 ± 0.614.2020–4033.22–35.66
2816043.27 ± 1.6012.4025–6540.06–46.47
26061.14 ± 0.584.5252.1–6559.97–62.31
36032.25 ± 0.372.8425–3831.52–32.98
43658.52 ± 0.915.4544.7–65.556.67–60.36
54834.02 ± 0.543.7121–3732.94–35.10
3516053.98 ± 1.5812.2625–7550.82–57.15
26062.00 ± 0.665.1050.3–67.260.69–63.31
35632.27 ± 0.423.1425–3831.43–33.11
43657.09 ± 0.845.0648.4–6555.38–58.80
54835.31 ± 0.573.9722–3834.16–36.46
Slaughter16066.82 ± 1.7213.2925–9563.38–70.25
26063.23 ± 0.675.1951.1–66.961.89–64.57
35632.68 ± 0.362.7325–3831.95–33.41
43659.81 ± 0.714.2949.5–67.658.36–61.26
54833.27 ± 0.513.5121–3932.25–34.29

CI = Confidence interval, n = Number of flocks, SD = Standard deviation, SEM = Standard error of the mean.

Correlation analysis confirmed that light intensity was strongly associated with lameness incidence across all rearing stages (Table 9). This relationship is illustrated in Figure 3, showing that higher light intensity on day 14 was associated with increased lameness incidence. The strongest correlations were observed at day 14 (r = 0.56, p < 0.0001) and day 21 (r = 0.52, p < 0.0001), coinciding with the period of rapid skeletal mineralization and subsequent accelerated muscle growth. Even at day 1, light intensity was positively correlated with lameness (r = 0.41, p < 0.0001), suggesting that inappropriate lighting conditions from the earliest stages can predispose broilers to locomotor disorders. This period corresponds to the transition phase in which skeletal mineralization may lag behind rapid muscle accretion, making birds particularly vulnerable to excessive mechanical loading induced by high activity levels under intense illumination.

Table 9. Pearson’s correlation coefficients between light intensity and total lameness incidence across all farms and rearing days.

VariableDaynr95% CIp-value
Light intensity12700.410.31–0.50< 0.0001
72740.520.43–0.60< 0.0001
142700.560.47–0.64< 0.0001
212720.520.43–0.60< 0.0001
282640.520.42–0.61< 0.0001
352720.330.22–0.43< 0.0001
Slaughter2700.150.03–0.270.012

CI = Confidence interval, n = Number of flocks, r = Pearson’s correlation coefficient.

Figure 3

Figure 3. Relationship between light intensity on day 14 and lameness incidence. Figure 3 shows the relationship between light intensity on day 14 and lameness incidence. The scatterplot demonstrates a strong positive association, with increasing light intensity corresponding to higher lameness prevalence. Two main clusters are evident, representing moderate lighting regimes of approximately 25–35 lux and high lighting regimes of >50 lux. Flocks exposed to higher light intensities consistently exhibited elevated lameness levels. The regression line confirms this trend, consistent with the strong correlation and multivariable regression results.

These findings deviate markedly from the Cobb Genetics guidelines [19], which recommend a minimum of 25 lux during the first 0–7 days to stimulate early feed intake, followed by a gradual reduction to 5–10 lux from day 7 onwards to balance welfare and performance. In addition, the variation in light intensity across the floor should not exceed 20%. The persistently high light levels in Farms 2 and 4, particularly beyond day 14, exceeded these recommendations and likely contributed to higher lameness incidence. During days 8–21, measured light intensities in Farms 2 and 4 frequently reached 50–59 lux, exceeding recommended levels by approximately two- to four-fold [3]. Such prolonged exposure to excessive illumination during periods of rapid skeletal growth may have amplified mechanical loading on immature leg structures.

Recent research provides biological support for the observed association between light environment and leg health in broilers. Higher light intensities (>30–50 lux) have been shown to increase bird activity, such as walking, foraging, and competitive interactions, thereby elevating mechanical loading on immature skeletal structures and increasing the likelihood of locomotor problems, particularly in fast-growing strains [20, 21]. Prolonged or continuous light exposure further reduces opportunities for rest and compromises cartilage turnover and bone remodeling, especially during periods of rapid weight gain [5, 20]. Conversely, suboptimal light levels during brooding (<20 lux) may limit feed- and water-locating behavior, delay early skeletal development, and increase long-term susceptibility to lameness [5, 21]. In this study, Farms 2 and 4 exhibited characteristics consistent with excessive light exposure, while Farm 1 occasionally approached suboptimal levels during brooding, aligning with these documented risks.

Overall, the results demonstrate that inappropriate light-intensity management, particularly excessive illumination, during both the early brooding period (days 0–7) and the critical skeletal development phase (days 14–21) was strongly associated with increased lameness incidence. Recent studies show that high light intensity elevates activity levels and mechanical loading on immature leg structures, while insufficient dark periods reduce resting opportunities and impair skeletal recovery [20, 21]. These findings emphasize the importance of adopting appropriate lighting programs, including reducing light intensity after the first week and ensuring structured dark periods to support rest, bone remodeling, and overall locomotor health. Aligning farm practices with breeder recommendations is therefore essential to minimizing leg strain and improving broiler welfare.

Ammonia concentration and late-stage risk

Ammonia concentrations gradually increased with bird age across all farms (Table 10). According to EFSA and international welfare guidelines, ammonia concentrations above 15.00 ppm are associated with respiratory irritation, reduced performance, and compromised welfare. During the first two weeks, mean ammonia levels remained low. At day 7, concentrations ranged between 3.15 and 4.35 ppm, and by day 14 between 4.57 and 6.00 ppm, all well below the recommended threshold of 10.00–15.00 ppm [3].

Table 10. Descriptive statistics of ammonia concentration in broiler houses across farms.

DayFarmnMean ± SEM (ppm)SDMin–max (ppm)95% CI (ppm)
71603.15 ± 0.110.842–52.93–3.37
2603.55 ± 0.070.533–53.41–3.69
3823.63 ± 0.141.241–63.36–3.91
4344.35 ± 0.160.922–74.03–4.67
5483.31 ± 0.100.722–43.10–3.52
141604.57 ± 0.151.182–84.26–4.87
2605.60 ± 0.080.623–65.44–5.76
3825.70 ± 0.191.704–95.32–6.07
4346.00 ± 0.130.794–85.73–6.27
5484.69 ± 0.090.663–54.50–4.88
211606.30 ± 0.161.273–95.97–6.63
2607.12 ± 0.080.586–86.97–7.27
3827.15 ± 0.252.225–116.66–7.63
4347.06 ± 0.150.895–96.75–7.37
5486.21 ± 0.120.824–75.97–6.45
281607.90 ± 0.161.264–107.57–8.23
26011.95 ± 0.131.0211–1611.69–12.21
3629.52 ± 0.282.186–138.96–10.07
43411.09 ± 0.120.7110–1210.84–11.34
54810.67 ± 0.281.975–1210.09–11.24
351609.73 ± 0.191.496–129.35–10.12
26013.95 ± 0.161.2712–1613.62–14.28
35011.96 ± 0.171.1810–1411.63–12.29
43412.35 ± 0.271.559–1511.81–12.90
54811.79 ± 0.251.767–1311.28–12.30
Slaughter16012.03 ± 0.151.138–1511.74–12.33
26016.15 ± 0.080.6315–1715.99–16.31
35015.34 ± 0.120.8713–1715.09–15.59
43414.97 ± 0.160.9012–1614.66–15.29
54815.67 ± 0.171.1910–1615.32–16.01

NH₃ = Ammonia, CI = Confidence interval, n = Number of flocks, SD = Standard deviation, SEM = Standard error of the mean.

From day 21 onwards, ammonia concentrations progressively increased, reaching approximately 6.21–7.15 ppm. By day 28, mean values ranged from 7.90 to 11.95 ppm, with several farms exceeding 10.00 ppm. At day 35, concentrations further increased to 9.73–13.95 ppm, approaching the upper welfare threshold. By slaughter day, mean ammonia levels reached 16.15 ppm in Farm 2 and 15.34 ppm in Farm 3, surpassing the 15.00 ppm welfare threshold, while Farm 4 approached this limit (14.97 ppm). These findings indicate that although ammonia remained within acceptable limits during early rearing, late-stage accumulation exposed birds to suboptimal air quality conditions.

The correlation analysis (Table 11) further supported these observations. Ammonia concentrations showed no significant relationship with lameness incidence during the early stages (day 7–21), with correlation coefficients close to zero (r = 0.06 to 0.03, p > 0.3). From day 28 onwards, however, the relationship became evident, as illustrated in Figure 4. Positive and significant correlations were observed (day 28: r = 0.31, p < 0.0001; day 35: r = 0.41, p < 0.0001; slaughter day: r = 0.39, p < 0.0001). This suggests that ammonia accumulation beyond 10 ppm during the late growing period was associated with increased lameness incidence.

Table 11. Pearson’s correlation coefficient between ammonia concentration and total lameness incidence across all farms and rearing days.

VariableDaynr95% CIp-value
Ammonia concentration72840.06−0.06 to 0.170.3555
14284−0.01−0.13 to 0.100.8163
212840.03−0.08 to 0.150.5674
282640.310.20–0.42< 0.0001
352520.410.31–0.51< 0.0001
Slaughter2520.390.28–0.49< 0.0001

CI = Confidence interval, n = Number of flocks, r = Pearson’s correlation coefficient.

Figure 4

Figure 4. Relationship between ammonia concentration on day 28 and lameness incidence. Figure 4 presents the association between ammonia concentration on day 28 and lameness incidence. Although a weak positive trend is visible, the data show substantial variability, particularly at moderate ammonia levels of 8–12 ppm. Several flocks with similar ammonia concentrations exhibited markedly different lameness outcomes. The regression line indicates only a modest association, which is consistent with the non-significant effect of ammonia in the final multivariable model.

The biological mechanisms underlying this relationship are well-documented. High ammonia exposure can cause respiratory irritation and lung damage, leading to hypoxia, which indirectly impairs skeletal development [3, 22]. Elevated ammonia levels also compromise litter quality, predisposing birds to footpad dermatitis and increased joint pressure, which subsequently contribute to locomotor impairment [3, 2325]. Additionally, poor air quality due to high ammonia levels can increase susceptibility to bacterial infections, such as bacterial chondronecrosis with osteomyelitis (BCO), thereby exacerbating leg health issues [22]. However, detailed litter quality parameters were not systematically evaluated. Although litter quantity per bird was similar across farms, variation in litter condition may have affected ammonia accumulation and leg health, potentially confounding the interpretation of the results. Taken together, these results indicate that ammonia did not pose a major risk for lameness during the brooding and early rearing stages, as concentrations remained below the critical threshold of 10 ppm. However, progressive ammonia accumulation in the later stages (>28 days) significantly increased the risk of lameness, particularly in Farms 2, 4, and 5, where ammonia levels exceeded or approached 15 ppm. These findings underscore the importance of maintaining effective ventilation and litter management throughout the production cycle to prevent late-stage ammonia buildup and its detrimental effects on leg health.

Prolonged exposure to suboptimal temperatures, excessive light intensity, and elevated ammonia levels may increase economic losses through reduced growth, increased culling, and greater carcass rejection due to locomotor disorders. After Benjamini–Hochberg FDR adjustment, key associations (day 1 Tmax, day 7 ΔT, day 14 light intensity, and late-stage ammonia) remained significant (Supplementary Table 4), supporting the robustness of these relationships.

Multivariable regression model

The stepwise selection procedure (Supplementary Table 3) identified early thermal conditions and mid-cycle light intensity as candidate predictors of lameness. However, to maintain biological interpretability and avoid multicollinearity among correlated environmental variables, a final multivariable model was constructed a priori using five predictors: highest temperature on day 1, diurnal temperature difference on day 7, light intensity on day 14, and ammonia concentrations on days 28 and 35 (Table 12). The estimated regression equation was:

Lameness (%) = −0.137 + 0.0122 T_maxD1 − 0.0829 ΔT_D7 + 0.0197 L_D14 + 0.0381 NH3_D28 − 0.0234 NH3_D35

where T_maxD1 is the highest temperature on day 1 (°C), ΔT_D7 is the diurnal temperature on day 7 (°C), L_D14 is the light intensity on day 14 (lux), and NH3_D28 and NH3_D35 are ammonia concentrations on days 28 and 35 (ppm).

Table 12. Multiple regression model evaluating environmental predictors of lameness (%) in commercial broilers.

Variableβ (Estimate)SEtp-valueVIF
Intercept−0.1370.679−0.200.840
Highest temperature (Day 1)0.01220.02040.600.5511.24
Diurnal temperature (Day 7)−0.08290.0264−3.140.00191.28
Light intensity (Day 14)0.01970.002448.07< 0.00011.58
Ammonia concentration (Day 28)0.03810.02661.430.1534.36
Ammonia concentration (Day 35)−0.02340.0276−0.850.3994.62

The final regression model incorporated five environmental predictors selected a priori based on biological relevance and previous correlation analysis. Light intensity on day 14 and diurnal temperature on day 7 were statistically significant predictors of lameness, whereas highest temperature on day 1 and ammonia concentration on days 28 and 35 were not statistically significant in the multivariable model.

R² = 0.429, adjusted R² = 0.416, F(5, 222) = 33.38, p < 0.0001, root MSE = 0.400, Durbin–Watson = 0.98, n = 228.

SE = Standard error, VIF = Variance inflation factor.

The final multivariable regression model (Table 12) explained 42.9% of the variation in flock-level lameness percentage (R² = 0.429; adjusted R² = 0.416; F(5,222) = 33.38, p < 0.0001; Root MSE = 0.400 % units). A greater diurnal temperature difference on day 7 (β = −0.0829, p = 0.0019) was linked to less lameness, while higher light intensity on day 14 (β = 0.0197, p < 0.0001) predicted more lameness. Maximum temperature on day 1 and ammonia levels on days 28 and 35 were not significant predictors. Variance inflation factors (VIFs) ranged from 1.24 to 4.62, indicating no significant multicollinearity. The Durbin–Watson statistic of 0.98 suggested mild positive autocorrelation. Q–Q plots in Figure 5 show near-normal residuals, with minor upper-tail deviations unlikely to impact model validity.

Figure 5

Figure 5. Quantile–quantile plot of flock-level lameness (%) showing approximate normal distribution of residuals. Figure 5 presents the normal quantile–quantile plot for flock-level lameness percentage. Each point represents one flock plotted against the theoretical normal distribution quantiles. Most observations closely follow the reference line, indicating that the lameness data approximate a normal distribution. Minor deviations are observed at the upper tail, where a few flocks exhibited relatively higher lameness incidence. However, these deviations were limited and did not substantially violate normality assumptions. This pattern supports the suitability of applying the General Linear Model and linear regression analyses to the dataset.

Although the final models explained a substantial proportion of the variation in lameness, a considerable amount of unexplained variability remained. Potential interaction effects between early thermal stress, light intensity, and late-stage air quality were not formally tested in the present model and warrant further investigation in future studies. Other unmeasured factors such as within-strain genetic differences, flock-level variation in stocking density, diet formulation (e.g., calcium, phosphorus, and vitamin D levels), litter management routines, and infectious causes of leg disorders may also have contributed to lameness but were not available in the present dataset. Therefore, the associations reported here should be interpreted as the environmental component of a broader multifactorial process rather than as definitive causal effects.

These results suggest that a practical farm-level risk framework using early-life maximum temperature (day 1), diurnal temperature (day 7), mid-cycle light intensity (day 14), and late-stage ammonia may support proactive management in tropical closed-house systems. However, the framework should be validated in additional flocks and tested for interaction effects before use as a predictive tool.

CONCLUSION

This study demonstrated that lameness in commercial Cobb 500 broilers raised under tropical closed-house conditions was associated with age-specific environmental stressors. The median lameness incidence was 1.00%, with most flocks remaining below 1.5%; however, significant farm-level variation was observed. Early-life thermal stress, particularly higher day 1 temperature, was positively associated with lameness, while day 7 diurnal temperature variation and day 14 light intensity emerged as key predictors in the final multivariable model. Late-stage NH3 accumulation after day 28 was also positively correlated with lameness in univariable analysis, indicating that air quality becomes increasingly important as birds approach market age.

The findings have direct practical implications for tropical broiler production. Maintaining brooding temperatures near breeder recommendations, avoiding excessive light intensity during the rapid skeletal development phase, and improving ventilation and litter management during late growth may reduce the risk of lameness. In particular, monitoring day 1 temperature, day 7 diurnal temperature variation, day 14 light intensity, and late-stage NH3 concentration may provide a practical early-warning framework for farm managers and veterinarians.

The major strength of this study is its large-scale commercial dataset, including 296 flocks and more than 11.45 million broilers across five closed-house farms managed under a standardized integrator system. This design reduced variation related to feeding, genetics, health management, and biosecurity, allowing clearer evaluation of environmental risk factors under real-world tropical production conditions.

However, the study also had limitations. The lameness diagnosis was based on clinical culling records without postmortem confirmation, and interobserver reliability was not formally assessed. Environmental measurements were mainly collected at central house locations, and within-house variation, litter condition, detailed infectious causes, and economic outcomes were not fully evaluated. The regression model explained approximately 43% of the variation in lameness, indicating that additional biological, management, and health-related factors likely contributed to the observed outcomes.

Future studies should validate these findings in independent commercial flocks and across different seasons, stocking densities, ventilation systems, and genetic lines. Further research should also include multi-point environmental sensing, litter quality assessment, postmortem confirmation of leg lesions, economic impact analysis, and mixed-effects or longitudinal modeling to test interactions among early thermal stress, lighting, and late-stage air quality.

In conclusion, lameness in tropical closed-house broiler systems is a multifactorial condition strongly influenced by time-specific environmental exposure. Optimizing early thermal management, controlling mid-growth light intensity, and preventing late-stage NH3 buildup may improve broiler welfare, reduce culling losses, and support more sustainable commercial poultry production.

DATA AVAILABILITY

The supplementary data can be made available from the corresponding author upon request.

AUTHORS’ CONTRIBUTIONS

PN and WL: Conceived and designed the study. PN: Data collection, preliminary data organization, and drafting of the original manuscript. PN and CR: Data analysis and statistical interpretation. CB: Supervision of the research, scientific guidance, and critical revision of the manuscript. PN, WL, CR, and CB: Review and editing of the manuscript. All authors have read and approved the final version of the manuscript.

COMPETING INTERESTS

The authors declare that they have no competing interests.

PUBLISHER’S NOTE

Veterinary World remains neutral with regard to jurisdictional claims in the published institutional affiliations.

ACKNOWLEDGMENTS

The authors would like to thank the commercial broiler farms in Eastern Thailand for granting access to farm data. The authors are also grateful to the Department of Animal Science, Faculty of Agriculture, Kasetsart University, Bangkok, Thailand, for academic support and technical assistance throughout the study. The authors acknowledge Huvepharma (Thailand) Ltd., Bangkok, Thailand, for supporting access to the Aviapp® program, which was used solely as a data recording and management tool.

REFERENCES

  1. Tallentire CW, Leinonen I, Kyriazakis I. Artificial selection for improved energy efficiency is reaching its limits in broiler chickens. Sci Rep 2018;8(1):1168. [Google Scholar]
  2. Zuidhof MJ, Schneider BL, Carney VL, Korver DR, Robinson FE. Growth, efficiency, and yield of commercial broilers from 1957, 1978, and 2005. Poult Sci 2014;93(12):2970-2982. [Google Scholar]
  3. Broiler performance and nutrition supplement. Siloam Springs (AR): Cobb-Vantress Inc; 2022. [Google Scholar]
  4. Wilcox CH, Sandilands V, Mayasari N, Asmara IY, Anang A. A literature review of broiler chicken welfare, husbandry, and assessment. Worlds Poult Sci J 2024;80(1):3-32. [Google Scholar]
  5. Liu KL, He YF, Xu BW, Lin LX, Chen P, Iqbal MK. Leg disorders in broiler chickens: A review of current knowledge. Anim Biotechnol 2023;34(9):5124-5138. [Google Scholar]
  6. Riddell C. Skeletal deformities in poultry. Adv Vet Sci Comp Med 1981;25:277-310. [Google Scholar]
  7. Scientific opinion on the influence of genetic parameters on the welfare and the resistance to stress of commercial broilers. EFSA J 2010;8(7):1666. [Google Scholar]
  8. Kestin SC, Knowles TG, Tinch AE, Gregory NG. Prevalence of leg weakness in broiler chickens and its relationship with genotype. Vet Rec 1992;131(9):190-194. [Google Scholar]
  9. Webster AB, Fairchild BD, Cummings TS, Stayer PA. Validation of a three-point gait-scoring system for field assessment of walking ability of commercial broilers. J Appl Poult Res 2008;17(4):529-539. [Google Scholar]
  10. Shim MY, Karnuah AB, Anthony NB, Pesti GM, Aggrey SE. The effects of broiler chicken growth rate on valgus, varus, and tibial dyschondroplasia. Poult Sci 2012;91(1):62-65. [Google Scholar]
  11. Anthney A, Do ADT, Alrubaye AA. Bacterial chondronecrosis with osteomyelitis lameness in broiler chickens and its implications for welfare, meat safety, and quality: A review. Front Physiol 2024;15:1452318. [Google Scholar]
  12. Quintana-Ospina GA, Alfaro-Wisaquillo MC, Oviedo-Rondon EO, Ruiz-Ramirez JR, Bernal-Arango LC, Martinez-Bernal GD. Effect of environmental and farm-associated factors on live performance parameters of broilers raised under commercial tropical conditions. Animals (Basel) 2023;13(21):3312. [Google Scholar]
  13. Kpomasse CC, Oke OE, Houndonougbo FM, Tona K. Broiler production challenges in the tropics: A review. Vet Med Sci 2021;7(3):831-842. [Google Scholar]
  14. Sánchez-Casanova RE, Sarmiento-Franco L, Phillips CJ. The effects of providing outdoor access to broilers in the tropics on their behaviour and stress responses. Animals (Basel) 2022;12(15):1917. [Google Scholar]
  15. Apalowo OO, Ekunseitan DA, Fasina YO. Impact of heat stress on broiler chicken production. Poultry 2024;3(2):107-128. [Google Scholar]
  16. SAS/STAT®15.3 user's guide. Cary (NC): SAS Institute Inc; 2023. [Google Scholar]
  17. Best JW. Research in education. Englewood Cliffs (NJ): Prentice-Hall; 1977. [Google Scholar]
  18. Granquist EG, Vasdal G, de Jong IC, Moe RO. Lameness and its relationship with health and production measures in broiler chickens. Animal 2019;13(10):2365-2372. [Google Scholar]
  19. Cobb broiler management guide. Siloam Springs (AR): Cobb-Vantress; 2021. [Google Scholar]
  20. Kang SW, Christensen KD, Kidd MT, Orlowski SK, Clark J. Effects of a variable light intensity lighting program on the welfare and performance of commercial broiler chickens. Front Physiol 2023;14:1059055. [Google Scholar]
  21. Zhou S, Thornton T, Gan H, Tabler T, Zhao Y. Impact of lighting intensity on welfare and performance in broiler chickens. Animals (Basel) 2025;15(22):3348. [Google Scholar]
  22. Zhou Y, Zhang M, Zhao X, Feng J. Ammonia exposure induced intestinal inflammation injury mediated by intestinal microbiota in broiler chickens via TLR4/TNF-αsignaling pathway. Ecotoxicol Environ Saf 2021;226:112832. [Google Scholar]
  23. Barbosa LV, Lima NDDS, Barros JDSG, Moura DJ, Estellés F, Ramón-Moragues A. Predicting risk of ammonia exposure in broiler housing: Correlation with incidence of health issues. Animals (Basel) 2024;14:615. [Google Scholar]
  24. Kwon BY, Park J, Kim DH, Lee KW. Assessment of welfare problems in broilers: Focus on musculoskeletal problems associated with their rapid growth. Animals (Basel) 2024;14(7):1116. [Google Scholar]
  25. Mocz F, Berthelot M, Michel V, Contreras-Jodar A, Guinebretière M. The assessment of litter quality in broiler chickens: Validity, inter-assessor reliability, and intra-assessor repeatability of three visual scoring systems. Br Poult Sci 2025;66(2):147-154. [Google Scholar]