The new analysis, published in the journal Subterranean Biology, compiled 291 measurements of oxygen consumption rates from 37 species spanning five taxonomic classes, from tiny copepods weighing less than a microgram to the olm, a blind amphibian that can outweigh a small mammal. The dataset drew on 23 studies published between 1971 and 2023, with research effort accelerating sharply after the year 2000, and covered groundwater fauna from both the Nearctic and Palearctic regions, including sites in Canada, the United States, Mexico, France, Portugal, Spain, Italy and Slovenia. The findings reveal that groundwater species obey the same universal metabolic scaling rules that govern life everywhere else, yet respond to temperature in a way that sets them apart from their surface-dwelling relatives.

Metabolic rate, the energy an organism burns over time, is central to virtually every aspect of animal ecology. It determines how much food an animal needs, how fast it grows, how quickly it reproduces and how much heat it produces as an inevitable byproduct of the biochemical machinery of life. For ectotherms, animals whose body temperature follows that of their surroundings, metabolism is also profoundly shaped by the environment. Within permissive thermal ranges, a rise of 10 degrees Celsius often doubles or even triples metabolic rates, an effect quantified by the familiar Q10 coefficient or by the Arrhenius activation energy, which typically falls between 0.5 and 0.8 electron volts for most ectotherms, equivalent to a 7 to 12 percent increase in metabolism per degree Celsius.

Metabolic rate also scales with body size according to the allometric law, in which standard metabolic rate equals a taxon-specific constant multiplied by body mass raised to a scaling exponent. When Kleiber first described this relationship in 1932, he found an exponent close to three-quarters, and many invertebrate taxa cluster between 0.66 and 0.80, a shared characteristic whose underlying reason remains one of biology’s enduring puzzles. Exponents below one indicate non-isometric scaling, meaning metabolic rates do not double when body mass doubles. Proposed explanations range from heat exchange with the environment, since larger animals have proportionally less surface area relative to their volume, to the nonlinear scaling of circulatory and respiratory systems and to cellular and mitochondrial constraints.

To test whether groundwater species follow this universal pattern, the researchers systematically searched the Web of Science and grey literature, screening more than 2,200 papers to isolate studies that quantified oxygen consumption in groundwater metazoans. For the scaling analysis, they assembled data on 10 stygobitic species across 64 measurements, including two copepods, three isopods, four amphipods of the genus Niphargus and the cave-dwelling amphibian Proteus anguinus. Only measurements from resting, fasted adults, unexposed to stressors and kept within temperatures matching their native habitats, were included, and one value per species per study was retained to prevent statistical bias. The result was striking: body mass explained 90.8 percent of the linear variation in oxygen consumption, with a slope of 0.80, squarely within the range seen in crustaceans and other surface taxa. Groundwater life, it seems, adheres to the universal allometric law without deviation.

Yet the picture changes when temperature enters the equation. The team compiled 288 measurements from 19 stygobitic species, including copepods, isopods, amphipods, crayfish, the olm and the cavefish Astyanax mexicanus, and compared them with data from 12 non-stygobitic species drawn from the same literature. Because aquatic subterranean organisms naturally cope with periodic oxygen deficiency, shifting from normoxia above 3 milligrams per liter to hypoxia below 0.3, the analysis included measurements taken under both oxygen regimes. Oxygen consumption rates were converted to a standard unit of microliters of oxygen per milligram of dry weight per hour, and the thermodynamic response across the range from minus 2 to 30 degrees Celsius was calculated from Arrhenius plots, with the slope estimated separately for low and high temperature ranges.

The results followed a classic bell-shaped thermal profile, with metabolic rates rising up to roughly 16 to 18 degrees Celsius and declining thereafter. In the low temperature range, the activation energy for stygobitic species was 94.7 kilojoules per mole, well below the 132.2 kilojoules per mole calculated for non-stygobitic species. This means that at low to moderate temperatures, the metabolisms of groundwater species are less responsive to warming than those of surface species, a pattern consistent with the stable, cold conditions of their habitats and possibly with reduced costs for mitochondrial biosynthesis and maintenance. In the energy-poor subterranean realm, where every joule counts, a metabolism that resists temperature-driven acceleration may be an adaptation for conserving scarce resources.

But above the thermal break, the relationship reverses dramatically. In the high temperature range from 19 to 30 degrees Celsius, the activation energy for stygobitic species plunged to minus 232.7 kilojoules per mole, compared with minus 83.6 for non-stygobitic species, indicating that metabolic rates in groundwater species decline far faster with warming than in surface taxa. The authors interpret this as evidence that hypogean species possess a less efficient enzyme system at elevated temperatures, with respiratory enzymes becoming inactivated outside their usual thermal range and leading to an irreversible loss of function. Surface species, which live in environments with large temperature fluctuations, may carry a more robust molecular defense against heat stress. Comparative genomics of terrestrial subterranean and surface species supports this idea, showing that surface taxa express a notably larger number of genes in response to thermal stress, suggesting a stronger heat shock response capability, although this has not yet been demonstrated for aquatic groundwater fauna.

The activation energies reported in the study, roughly 0.98 and 1.3 electron volts for the low range and minus 2.41 and minus 0.86 electron volts for the high range, exceed the average of 0.48 electron volts calculated across 314 aquatic and terrestrial ectothermic species within their permissive temperature range. This suggests that the temperature ranges investigated may not entirely fall within the permissive thermal limits of the species analyzed, and that even modest global warming scenarios could impose disproportionate metabolic impairments on groundwater species compared with surface organisms. The finding aligns with a recent meta-analysis on the effects of climate change on subterranean ecosystems and adds a quantitative physiological foundation to growing conservation concerns.

The study also uncovered remarkable metabolic variability among closely related groundwater groups. Stygobitic cyclopoid copepods showed oxygen consumption rates roughly six times higher than those of stygobitic harpacticoids, despite both belonging to the same class, indicating that factors beyond temperature and body mass shape metabolic rates. In some stygobitic species, including the copepod Diacyclops belgicus, the isopod Proasellus lusitanicus and the amphipod Gammarus acherondytes, oxygen consumption appears mass-independent altogether, a pattern regarded as an adaptive trait to energy-limited habitats where individuals may grow larger without increasing their metabolic rates. Weight gain in these animals is attributed to metabolically inactive reserves, and stygobitic species possess larger fat stores than their surface relatives, allowing them to endure prolonged fasting. However, low metabolic variability within a species may heighten population vulnerability to anthropogenic disturbances, adding a further layer of fragility.

The authors acknowledge limitations inherent in any synthesis of published data. Inconsistent reporting of body size forced calculations and estimations, including the assumption that dry weight equals 15 percent of wet mass for amphipods, and pooling species into broad ecological groups may mask species-specific variation. Certain species are overrepresented in particular temperature classes, and the dataset contains no hyporheic species, the groundwater fauna living in the saturated sediments beneath rivers, despite their integral connection to surface ecosystems. Future research, the authors argue, should expand into lesser-studied habitats, employ experimental approaches to directly assess the impact of environmental stressors, and integrate energy budgets, food web dynamics and ecosystem functioning into metabolic studies.

For conservation biologists, the implications are potentially far-reaching. Understanding the metabolic responses of stygobitic species to temperature change can inform strategies to protect these animals in regions predicted to experience significant warming, and the identification of a physiological thermal break around 18 degrees Celsius offers a concrete threshold for managers of cave and aquifer systems. Groundwater species, with poor dispersal abilities and high habitat specialization, will depend almost entirely on their physiological capacity to cope with environmental change. As the planet warms, the slow, stable metabolisms that have sustained life in the dark for millions of years may prove to be the very trait that determines whether it can persist.

Subject of Research: groundwater ecology

Article Title: Metabolic rates of groundwater species as a function of body mass and temperature

Article References: Di Lorenzo, T., Mori, N., & Simčič, T. (2024). Metabolic rates of groundwater species as a function of body mass and temperature. Subterranean Biology, 49, 53-74. https://doi.org/10.3897/subtbiol.49.125131

Image Credits: AI Generated

DOI: 10.3897/subtbiol.49.125131

Keywords: stygobitic species, groundwater, metabolic rate, oxygen consumption, allometric scaling, Arrhenius activation energy, temperature, climate change, Subterranean Biology, ectotherms, cave fauna, thermal tolerance

The research team, led by Sofia Eirini Chatoutsidou of the Technical University of Crete together with colleagues from the N.C.S.R. Demokritos research center and the National and Kapodistrian University of Athens, built their simulation around a realistic exposure pattern. They assumed eight hours per day of breathing diluted diesel exhaust, chosen to approximate an occupational setting such as a depot, loading dock, or roadside workplace, followed by sixteen hours of ordinary urban air pollution. That split matters because the two aerosols differ sharply in character. Diesel exhaust particles are soot agglomerates roughly 70 to 80 nanometers across, studded with polycyclic aromatic hydrocarbons, semi-volatile organics, and metals, while urban particles are fewer in number and spread across a broader size range.

The particle data came from measurements in Toronto, where urban air was sampled at a central intersection of two high-traffic roads, and from a heavy-duty diesel engine run at low load and intermediate speed, conditions that mimic urban cruising. Using a Scanning Mobility Particle Sizer, the diesel exhaust showed a number concentration of about 5.37 × 10⁴ particles per cubic centimeter, peaking near 70 to 80 nanometers, while urban air carried roughly 1.76 × 10⁴ particles per cubic centimeter. The authors checked these values against the wider literature and found them plausible: vehicles without diesel particulate filters can emit particles in the range of 10⁶ to 10⁷ per cubic centimeter, whereas filter-equipped vehicles drop to around 10⁴, though active filter regeneration events can spike concentrations back up by orders of magnitude.

To track where those particles go, the team used ExDoM2, a dosimetry model built on the semi-empirical respiratory tract framework of the International Commission on Radiological Protection. The model divides the human respiratory tract into five regions: the anterior nose, the posterior nasal passages with pharynx and larynx, the bronchial region, the bronchiolar region, and the deep alveolar-interstitial zone where gas exchange happens. Each region acts as a series of filters, and deposition in each depends on particle size, breathing rate, and the balance between aerodynamic mechanisms like impaction and settling and thermodynamic mechanisms like diffusion. For nanoparticles, diffusion dominates, which is precisely why particles smaller than 100 nanometers penetrate so effectively into the lung’s deepest recesses.

The headline result of the deposition simulations is stark. For an adult male of normal weight engaged in light exercise such as walking, 69 percent of inhaled nanoparticles from either diesel exhaust or urban air deposited in the alveolar-interstitial region, and combining that with the tracheobronchial zone means 86 percent of the inhaled nanoparticles ended up in the lungs after a single day. In absolute terms, one hour of diesel exposure delivered about 2.7 × 10¹⁰ particles to the lungs, corresponding to a surface dose of 4.15 × 10¹⁴ square nanometers and a mass dose of 9.9 micrograms. Surface area, rather than mass, turns out to be the critical metric, because rodent toxicology studies have shown that the tumor response to poorly soluble particles tracks the accumulated particle surface in lung tissue.

Clearance told an equally important story. During an eight-hour diesel exposure, 19.3 micrograms deposited in the extrathoracic region and 79.2 micrograms in the lungs, but after clearance only 9.4 and 71.6 micrograms respectively remained. That means the upper airways shed just over half of their deposited load, while the lungs cleared a mere 9.5 percent. Meanwhile, small fractions of the dose migrated elsewhere: 10.5 micrograms traveled to the esophagus, 5.9 micrograms were absorbed into the blood, and a trace amount reached the lymph nodes. The deep lung, with its slow macrophage-mediated clearance, effectively becomes a reservoir for diesel soot, and it is exactly this retained burden that the model projects forward over years.

When the simulations were extended to long-term exposure, the retained lung mass followed an exponential curve that reached a plateau after roughly 1,000 days, or about three years of constant exposure. Projected out to 10,000 days, approximately 27 years, the retained dose stabilized at around 7,786 micrograms for combined diesel and urban exposure, compared with 7,761 micrograms for diesel alone, a difference so small it underscores how overwhelmingly the eight-hour diesel window dominates the total burden. The saturated dose was roughly 77 times the initial one-day dose. The authors caution that this plateau reflects the assumption of unchanging exposure conditions; in reality, shifts in concentration or episodes of very high exposure can overwhelm alveolar macrophages, the lung’s first-line scavenger cells, and change the deposition picture entirely.

The most innovative step came next. Using the allometric scaling approach recommended by the Environmental Protection Agency, in which toxicity doses between species scale with body mass to the power of three-quarters, the team converted the long-term retained human lung dose into an animal equivalent dose. Because the target tissue is the lung itself, the scaling used lung masses rather than body masses: about 1,000 grams for a human lung versus 1.5 grams for a rat lung, a ratio of roughly 667. Smaller animals have higher metabolic rates and correspondingly higher dose accumulation, so the equivalent rat dose came out higher than the human dose. The converted values fell in the lower range of the experimental dose-response curve compiled from classic rat inhalation studies by Oberdörster, Miller, and others, where the logistic fit predicts a lung tumor proportion of 2.42 percent, with a 95 percent confidence interval of 0.073 to 6.12 percent.

Notably, this tumor estimate was essentially identical across all the receptor types the model examined, including normal-weight and overweight men and women at different activity levels. Variations in inhalation rate and body composition did shift the deposited dose, with males under light exercise receiving the highest doses among normal-weight individuals, while overweight females sitting at the same activity level received higher doses than their normal-weight counterparts. But once those doses were scaled to the rat equivalent, they all landed in the flat lower portion of the logistic curve, where the predicted tumor proportion barely changes. The implication is that, within the physiological variation explored here, exposure conditions, not personal characteristics, are the dominant driver of the modeled carcinogenic risk.

The study has clear limitations that the authors acknowledge. It relied on a single particle size distribution for all simulations, whereas real-world exposure varies constantly; it assumed fixed receptor characteristics and constant exposure conditions; and the allometric and clearance formulations carry their own uncertainties, as does fitting a logistic curve to a limited set of in vivo observations. Still, the work represents a meaningful bridge between two worlds that rarely speak directly: computational dosimetry of human inhalation and experimental rodent carcinogenicity. By quantifying exactly how much diesel soot accumulates in human lungs over a working lifetime and showing that the scaled equivalent sits within the tumor-producing range observed in rats, the study gives risk assessors a concrete, mechanistically grounded link between everyday traffic pollution and the cancer endpoint that regulators care about most. For the millions of workers who spend their shifts beside idling diesel engines, the message is hard to ignore.

Subject of Research: Respiratory tract dosimetry of diesel exhaust particles and its allometric link to lung tumor responses in rats

Article Title: Linking human respiratory tract dose from diesel exhaust particles to lung tumor responses observed in rats

Article References: Chatoutsidou, S. E., Chalvatzaki, E., Mammi-Galani, E., Diapouli, L., Gini, M. I., Manousakas, M. I., Kasdagli, M.-I., Samoli, E., Katsouyanni, K., Eleftheriadis, K., & Lazaridis, M. (2026). Linking human respiratory tract dose from diesel exhaust particles to lung tumor responses observed in rats. Air Quality, Atmosphere & Health, 19(10), Article 226. https://doi.org/10.1007/s11869-026-02115-0

Image Credits: AI Generated

DOI: 10.1007/s11869-026-02115-0

Keywords: diesel exhaust particles, lung dosimetry, nanoparticles, allometric scaling, lung cancer, air pollution, respiratory tract, clearance, occupational exposure, risk assessment, ExDoM2 model, rat toxicology

Allometric growth is a concept borrowed from biology, where it describes how, for example, the surface area of an animal scales with its body volume. In its classical mathematical form, the law states that the relative growth rate of one measure is proportional to the relative growth rate of another, with a constant coefficient. The solution to this differential equation is a power function, and the exponent of that function is the allometric scaling exponent. When John Q. Stewart and later Naroll and von Bertalanffy imported these ideas into the study of cities and urbanization, scientists assumed the exponents should be explainable by Euclidean dimensions. Urban population was treated as a three-dimensional measure and urban area as a two-dimensional one, predicting a scaling exponent of two-thirds. Where two measures shared the same dimension, the expected exponent was exactly one.

Observations refused to cooperate. Decades of empirical calculations on real cities produced scaling exponents that were neither integers nor simple ratios of integers, and this mismatch created what the literature calls the dimensional dilemma. Scientists faced an unpalatable choice: either abandon the law of allometric growth, which clearly described real patterns in the data, or find a more sophisticated geometry to explain it. The escape route, Chen argues, was opened by the arrival of fractal geometry in the 1970s and 1980s, through the work of Benoit Mandelbrot and, in urban contexts, Michael Batty and Pierre Frankhauser. Once cities were recognized as fractal objects, sprawling, self-similar structures whose detailed shapes repeat across scales, the exotic values of scaling exponents suddenly became interpretable.

The theoretical core of the new study rests on the principle of dimensional consistency, an idea stretching back to ancient Greek mathematics. Measures of different dimensions, such as length, area, and volume, cannot form simple proportional relationships; to construct such relationships, the dimensions must be brought into consistency through exponents. When two measures of a complex system are proportional, the scaling exponent connecting them must equal the ratio of their dimensions. If those dimensions are Euclidean, the exponent should be a recognizable fraction like two-thirds or three-halves. But if the exponent is not such a ratio, then at least one of the underlying measures must be fractal, possessing a dimension greater than its topological dimension and typically non-integer. Chen turns this logic into a diagnostic table: integer or simple-fraction exponents consistent with Euclidean geometry point to ordinary measures, while values such as three-quarters or 0.85 betray the presence of fractal structure.

To make the argument concrete, Chen assembled empirical evidence from Chinese cities, drawing on census data from 2000 and 2010 covering the thirty-one regions of the Chinese mainland. Two allometric relationships were examined: the scaling between the total urban population of each region and the number of cities it contains, and the scaling between total urban population and the population of each region’s central city. In every case, the estimated exponents fell well below one, clustering around values from roughly 0.64 to 0.72 depending on the year, the relationship, and the estimation method. None of these numbers can be produced by any ratio of Euclidean dimensions. Under the dimensional consistency principle, the only coherent conclusion is that the measures involved are fractal, and the exponents are ratios of fractal dimensions.

The analysis is careful about statistics in ways that matter for reproducibility. Chen distinguishes between the three worlds of science, the real world of cities, the mathematical world of deductive reasoning, and the computational world of data and algorithms, and shows that parameter estimates depend on which method is used. The least squares regression of population on city number yields a different exponent from the regression of city number on population, and the product of the two exponents equals the goodness of fit, a relationship confirmed in the data. To handle this asymmetry, the study employs the reduced major axis method, which averages the two directions of regression and yields modified exponents. Municipalities directly under the central government, which contain only a single city and therefore behave unlike multi-city provinces, were identified as outliers using standardized residuals, scatter plots, and K-means clustering, and were excluded from the fitting in a documented and principled way.

Beyond the headline relationships, the study connects allometric scaling to other cornerstones of urban science. Zipf’s law, the celebrated rank-size rule stating that city populations follow a power-law distribution, turns out to be mathematically intertwined with cross-sectional allometry: a pair of correlated Zipf distributions for city population and city area derives the transversal allometric model, and the Pareto exponent of the size distribution can itself be read as a fractal dimension. In this framework, cross-sectional allometric scaling is a secondary law flowing from rank-order scaling, and the allometric exponent between two size measures equals the ratio of their two Pareto-derived fractal dimensions. Spatial allometry, measured through concentric circles drawn around a city center, similarly yields exponents that are ratios of radial fractal dimensions describing how population density and built-up land fall away from the core.

The work does not stop at geometry; it also ventures into why scaling laws exist at all. Chen reviews candidate mechanisms, including self-organized criticality, proportional random growth, and preferential attachment, and then advocates a dual explanation grounded in two mathematically equivalent principles: entropy maximization at the macro level and utility maximization at the micro level. A single entropy-maximizing process generates an exponential distribution, but a pair of coupled entropy-maximizing processes, one governing the growing number of cities and one governing the growth of each city, combines two exponential functions into a power law. The same structure emerges from microeconomic reasoning, where individuals and organizations maximizing their satisfaction produce the same aggregate pattern. Intriguingly, this equivalence suggests that the socio-economic character of a city is partly legible through its physical, fractal properties.

The research even carries policy implications. In the multivariate Cobb-Douglas analysis, the contribution of city number to total urban population grew stronger between 2000 and 2010, while the contribution of central-city population weakened, indicating that many small and medium-sized cities collectively drive urbanization more than a few giant metropolises. This is the long tail effect of power-law distributions: when the fractal dimension of the city-size distribution exceeds one, the numerous small settlements at the tail of the hierarchy can absorb enormous urban populations. For Chen, the practical lesson is that developing small and medium-sized cities first, and then tuning the largest cities, is the mathematically informed path to raising urbanization, a conclusion that elevates a seemingly abstract fractal ratio into a tool for planning the urban future. The study’s remaining open question, the possible connection between allometric growth and the third defining property of fractals, entropy conservation, is flagged honestly as unsolved, leaving an inviting frontier for the next chapter of urban science.

Subject of Research: Fractal dimension as the mathematical basis of allometric scaling exponents in urban systems

Article Title: Fractal dimension accounts for allometric scaling exponents of cities

Article References: Chen, Y. (2026). Fractal dimension accounts for allometric scaling exponents of cities. Heliyon, 12(15), Article e45409. https://doi.org/10.1016/j.heliyon.2026.e45409

Image Credits: AI Generated

DOI: 10.1016/j.heliyon.2026.e45409

Keywords: fractal dimension, allometric scaling, urban science, power laws, Zipf's law, Chinese cities, urbanization, dimensional consistency, entropy maximization, city size distribution, Heliyon, complex systems

The study drew on data from a six-year, school-based physical activity initiative conducted in south-eastern Norway. Of 351 children whose parents provided consent, 330 completed at least one valid maximal treadmill test, an unusually rigorous approach in pediatric research where field estimates of fitness are far more common than direct measurement. Because valid peak oxygen uptake data from 2016 were available from only three of the nine participating schools following a technical data loss, the researchers based their inferential analyses on measurements collected in 2015, 2017, 2018, 2019 and 2020, providing a longitudinal window that spans most of primary school.

At the heart of the research lies a methodological question that has divided pediatric exercise scientists for decades: how should oxygen uptake be expressed when comparing children of different sizes? Raw absolute values in liters per minute naturally favor older, heavier children. Dividing by total body mass is the most common correction, but it penalizes children with higher body fat, who contribute mass but little metabolic machinery for running. Allometric scaling, which uses exponentials of body mass, and scaling to fat-free mass, which isolates metabolically active tissue, offer alternatives. The HOPP team, led by Asgeir Mamen of Kristiania University College, together with Julianna Buer of the Norwegian School of Sport Sciences and Per Morten Fredriksen of the University of Inland Norway, designated fat-free-mass-scaled peak oxygen uptake as their prespecified primary outcome and tracked all four expressions simultaneously.

The statistical architecture reflected the complexity of the design. Linear mixed models included test year, school allocation group, sex, and the interaction between test year and group as fixed effects, with participant identity as a random intercept and test year as the repeated factor. The primary model incorporated 904 observations from 314 children, a sample that gives the analysis substantial power to detect developmental trends. School-level models were run as sensitivity analyses to probe whether results were being driven by the clustering of children within individual schools, a known hazard of school-based interventions.

The headline result was unambiguous in direction, if nuanced in interpretation. Fat-free-mass-scaled peak oxygen uptake was significantly associated with test year, allocation group, sex, and the test year by group interaction, with all p-values at or below 0.030. From fifth to sixth grade, the adjusted mean fell by 5.73 milliliters per kilogram of fat-free mass per minute in the intervention group and by 5.11 in the control group. Critically, the difference between those group-specific changes was just 0.62 milliliters, with a 95 percent confidence interval spanning from minus 4.25 to 3.00 and a p-value of 0.736. In other words, the school-based physical activity program appeared to make no measurable difference to the final-year decline.

That equivalence between intervention and control schools is perhaps the study’s most sobering finding. The HOPP project was designed as a pedagogical experiment in which schools adopted enhanced physical activity programming, and the reasonable hope was that extra activity would buffer children against fitness losses as they approach adolescence. Instead, the data suggest the late-primary-school decline was resistant to the intervention as delivered. Body-mass-relative and allometrically scaled oxygen uptake also fell significantly in both groups during the final year, while absolute peak oxygen uptake, unadjusted for growth, held roughly steady, hinting that increases in raw aerobic capacity were simply failing to keep pace with children’s rapid physical development.

Sex differences emerged across every expression of fitness, with boys showing higher adjusted values than girls. Yet the magnitude of that gap depended heavily on how the data were scaled. When oxygen uptake was normalized to fat-free mass, which strips away the influence of differences in body composition between the sexes, the sex difference was substantially attenuated. This observation reinforces a growing consensus in pediatric physiology: apparent fitness gaps between boys and girls are partly artifacts of scaling choices and body composition rather than true differences in the oxidative capacity of muscle tissue.

The school-level sensitivity analyses added another layer of nuance. Changes in the three body-size-adjusted expressions of fitness were negative at all nine schools, suggesting the final-year decline was a pervasive phenomenon rather than a quirk of any particular classroom, neighborhood or teaching staff. Such consistency strengthens the biological plausibility of the trend, while also underscoring that whatever combination of growth, maturation and behavior drives the decline, it operates broadly across the study population and overwhelms school-level variation in programming.

The authors are candid about the limitations that temper firm causal conclusions. The study design could not separate growth-related developmental changes from shifts in behavior associated with the COVID-19 pandemic, which overlapped with the later measurement years and is known to have disrupted children’s physical activity worldwide, nor from the possibility that the composition of participants changed over time. A retrospective trial registration, dated 20 June 2015 on ClinicalTrials.gov as NCT02495714, and the 2016 data gap further complicate the inferential picture. Still, the convergence of negative findings across multiple scaling methods and all nine schools makes the late-primary decline difficult to dismiss as statistical noise.

For researchers, the message is methodological: conclusions about childhood aerobic fitness hinge materially on scaling method, body composition and school-level heterogeneity, and studies that report only a single expression of peak oxygen uptake may tell an incomplete or even misleading story. For parents, educators and policymakers, the findings add to mounting evidence that the transition out of primary school is a vulnerable window for children’s cardiovascular health, one that conventional school-based activity interventions, at least as implemented in HOPP, did not shield against. Identifying what does work during those years, and distinguishing genuine developmental biology from modifiable behavior, now stands as an urgent task for pediatric exercise science.

Subject of Research: Longitudinal development of peak oxygen uptake in 6- to 12-year-old children

Article Title: Longitudinal development of oxygen uptake in 6- to 12-year-old children: the Health Oriented Pedagogical Project (HOPP)

Article References: Mamen, A., Buer, J., & Fredriksen, P. M. (2026). Longitudinal development of oxygen uptake in 6- to 12-year-old children: the Health Oriented Pedagogical Project (HOPP). BMC Pediatrics. https://doi.org/10.1186/s12887-026-07631-7

Image Credits: AI Generated

DOI: 10.1186/s12887-026-07631-7

Keywords: peak oxygen uptake, aerobic fitness, children, fat-free mass, allometric scaling, treadmill running, primary school, physical activity, longitudinal study, pediatric health, HOPP, Norway