Biogeochemical Feedbacks Research Articles

Organic ligands in whale excrement support iron availability and reduce copper toxicity to the surface ocean

Nutrient recycling by marine megafauna is a key ecosystem service that has been disturbed by anthropogenic activity. While some hypotheses attribute Southern Ocean ecosystem restructuring to disruptions in micronutrient cycling after the elimination of two million great whales, there is little knowledge of trace metal lability in whale excrement. Here we measured high concentrations of dissolved iron and copper in five baleen whale fecal samples and characterized micromolar levels of organic metal-binding ligands as a proxy for their availability. The iron-ligand pool consisted of weakly-binding ligands and intermediate-binding ligands which enhanced iron stability and potential bioavailability. In comparison, 47 novel strongly-binding metallophores dominated copper-binding, curtailing its potential toxicity. These results illustrate how marine megafauna transform prey biomass into highly-labile micronutrients that they inject directly into the surface ocean, a mechanism whaling reduced by over 90%. Thus, the rapid restructuring of pelagic ecosystems through overharvesting may cause large biogeochemical feedbacks, altering primary productivity and carbon sequestration processes in the ocean.

North‐East Peri‐Tethyan Water Column Deoxygenation and Euxinia at the Paleocene Eocene Thermal Maximum

AbstractThe Paleocene–Eocene Thermal Maximum (PETM) is associated with climatic change and biological turnover. It shares features with the Oceanic Anoxic Events (OAEs) of the Mesozoic, such as transient global warming and biogeochemical perturbations. However, the PETM experienced a more muted expansion of marine anoxia compared to the Mesozoic OAEs (especially OAE2), with geographically limited evidence for photic zone euxinia (PZE). We explore the extent and drivers of marine deoxygenation during the PETM using biomarkers for water column euxinia and anoxia as well as an intermediate complexity Earth system model (cGEnIE). These reveal that the water column in the North‐East Peri‐Tethys became anoxic, with euxinic conditions reaching the photic zone (PZE) during the PETM. Our model shows that euxinia developed due to a global increase in the ocean nutrient inventory with concomitant oxygen consumption, similar to findings for OAE2. The particularly strong regional response in the NE Peri‐Tethys appears to arise from a combination of global CO2‐weathering forcing, regionally restricted circulation and upwelling of sulphidic thermocline waters. Unlike OAE2, anoxia and PZE do not become widespread in our PETM simulations, consistent with new and existing geochemical and biological proxy data. This globally muted response could result from reduced biogeochemical feedbacks to climate forcing relative to the mid‐Cretaceous climate. Our observations suggest that similar mechanisms operated in response to disparate Cenozoic (PETM) and Mesozoic (OAEs) transient global warming events, while also highlighting that background conditions are crucial in modulating the sensitivity of Earth's system to them.

Riverine nutrient impact on global ocean nitrogen cycle feedbacks and marine primary production in an Earth system model

Abstract. Riverine nutrient export is an important process in marine coastal biogeochemistry and also impacts global marine biology. The nitrogen cycle is a key player here. Internal feedbacks are shown to regulate not only nitrogen distribution, but also primary production and thereby oxygen concentrations. Phosphorus is another essential nutrient and interacts with the nitrogen cycle via different feedback mechanisms. After a previous study of the marine nitrogen cycle response to riverine nitrogen supply, here we include phosphorus from river export with different phosphorus burial scenarios and study the impact of phosphorus alone and in combination with nitrogen in a global 3D ocean biogeochemistry model. Again, we analyse the effects on near-coastal and open-ocean biogeochemistry. We find that riverine export of bioavailable phosphorus alone or in conjunction with nitrogen affects marine biology on millennial timescales more than riverine nitrogen alone. Biogeochemical feedbacks in the marine nitrogen cycle are strongly influenced by additional phosphorus. Where bioavailable phosphorus increases with river input, nitrogen concentration increases as well, except for in regions with diminishing oxygen concentrations. High phosphorus burial rates decrease biological production significantly. Globally, the addition of riverine phosphorus in the modelled ocean leads to elevated primary production rates in the coastal and open oceans.

Reduced Arctic Ocean CO2 uptake due to coastal permafrost erosion

Arctic coastal permafrost erosion is projected to increase by a factor of 2–3 by 2100. However, organic matter fluxes from the coastal permafrost into the ocean have not been considered in Earth system models so far. Here we represent coastal permafrost erosion in an Earth system model and perform simulations with varying permafrost organic matter properties, such as sinking fraction and nutrient content. We find that coastal erosion reduces the Arctic Ocean CO2 uptake from the atmosphere in all simulations: by 4.6–13.2 TgC yr−1 by 2100, which is ~7–14% of the Inner Arctic Ocean uptake. We show that coastal permafrost erosion exerts a positive biogeochemical feedback on climate, increasing atmospheric CO2 by 1–2 TgC yr−1 per °C of increase in global surface air temperature. Our work will allow coastal permafrost erosion to be considered in future climate change assessments.

Large Igneous Province Control on Ocean Anoxia and Eutrophication in the North Sea at the Paleocene–Eocene Thermal Maximum

AbstractThe Paleocene–Eocene Thermal Maximum (PETM) was a global hyperthermal event ∼56 Ma characterized by massive input of carbon into the ocean–atmosphere system and global warming. A leading hypothesis for its trigger is the emplacement of the North Atlantic Igneous Province (NAIP), with extensive extrusion/intrusion of igneous material into nearby sedimentary basins, forcing local uplift and warming‐inducing carbon emissions. It remains unclear if oceanographic changes in the North Sea–Norwegian Sea–Arctic basins, such as anoxia and productivity, were causally linked to local NAIP uplift/activity, and at what time scales these perturbations occurred. To test mechanisms and time scales, we present geochemical proxies (XRF analysis, clay mineralogy, molybdenum isotopes, and pyrite framboid size distribution) in undisrupted marine sediment core E−8X located in the central North Sea. We find evidence for a rapid onset of anoxia/euxinia at the negative carbon isotope excursion from redox proxies, followed by a gradual drawdown of molybdenum/total organic carbon (Mo/TOC) during the PETM main phase indicative of tectonically‐restricted basin likely from NAIP uplift. A short‐lived increase in Mo, pyrite and TOC occurred during a precursor event associated with a sedimentary mercury pulse indicative of volcanic activity. We suggest thermal uplift and flood basalt volcanism tectonically restricted the North Sea and tipped it into an euxinic state via volcanic emission–oceanographic feedbacks inducing eutrophication. This fine temporal separation of tectonic versus climatic geochemical proxies, combined with pulsed NAIP volcanism, demonstrates that Large Igneous Province emplacements can, at least locally, result in ocean biogeochemical feedbacks operating on relatively short timescales.

Seasonal drivers and patterns of sediment temperatures in a New England salt marsh

AbstractSalt marshes are globally important ecosystems and thus their resilience to climate change holds societal importance. To date, studies addressing salt marsh responses to climate change have focused on sea‐level rise and biogeochemical feedbacks with increasing inundation. Less is known about how salt marsh sediment temperatures, which impact physical, biological, and chemical ecosystem processes, will respond to climate change. In this study, we present multi‐depth sediment temperature and porewater level data from low‐, mid‐, and high‐marsh sites at a New England salt marsh for a 1‐year period and investigate how salt marsh sediment temperatures respond to atmospheric and oceanic forcing. We use spectral analyses to identify the frequencies at which sediment temperatures vary and link the temperature variations to specific forcing mechanisms. We find that all sites across the marsh responded to air temperature with roughly equal amplitude whereas the responses to radiative heating and ocean tides varied spatially. The high‐marsh site is more sensitive to radiative heating than the mid‐ and low‐marsh sites. The low‐marsh is affected by tidal processes and inundation whereas the high‐ and mid‐marsh sites are not. In addition, we find that the bulk thermal diffusivity of the saturated sediments decreases with distance from the tidal channel. These factors contribute to considerable temporal and spatial variability in sediment temperatures with elevation, distance from the tidal channel, and time of year (season) being most important.

Toward generalized Milankovitch theory (GMT)

Abstract. In recent decades, numerous paleoclimate records and results of model simulations have provided strong support for the astronomical theory of Quaternary glacial cycles formulated in its modern form by Milutin Milankovitch. At the same time, new findings have revealed that the classical Milankovitch theory is unable to explain a number of important facts, such as the change in the dominant periodicity of glacial cycles from 41 to 100 kyr about 1 million years ago. This transition was also accompanied by an increase in the amplitude and asymmetry of the glacial cycles. Here, based on the results of a hierarchy of models and data analysis, a framework of the extended (generalized) version of the Milankovitch theory is presented. To illustrate the main elements of this theory, a simple conceptual model of glacial cycles was developed using the results of an Earth system model, CLIMBER-2. This conceptual model explicitly assumes the multistability of the climate–cryosphere system and the instability of the “supercritical” ice sheets. Using this model, it is shown that Quaternary glacial cycles can be successfully reproduced as the strongly nonlinear response of the Earth system to the orbital forcing, where 100 kyr cyclicity originates from the phase locking of the precession and obliquity-forced glacial cycles to the corresponding eccentricity cycle. The eccentricity influences glacial cycles solely through its amplitude modulation of the precession component of orbital forcing, while the long timescale of the late Quaternary glacial cycles is determined by the time required for ice sheets to reach their critical size. The postulates used to construct this conceptual model were justified using analysis of relevant physical and biogeochemical processes and feedbacks. In particular, the role of climate–ice sheet–carbon cycle feedback in shaping and globalization of glacial cycles is discussed. The reasons for the instability of the large northern ice sheets and the mechanisms of the Earth system escape from the “glacial trap” via a set of strongly nonlinear processes are presented. It is also shown that the transition from the 41 to the 100 kyr world about 1 million years ago can be explained by a gradual increase in the critical size of ice sheets, which in turn is related to the gradual removal of terrestrial sediments from the northern continents. The implications of this nonlinear paradigm for understanding Quaternary climate dynamics and the remaining knowledge gaps are finally discussed.

An Overview of Anthropogenic Changes in Land Use and Land Cover, with Specific Attention to Climate Change and Unsustainable Agriculture

The worldwide need for food and the alterations in bioenergy demands associated with shifts in land use and land cover changes (LULCC) have given rise to environmental apprehensions, including global warming and climate change. There is substantial evidence indicating that we are currently witnessing human-induced transformations that could potentially result in a sixth mass extinction event. This study seeks to explore the connections between land use and land cover change (LULCC), forests, biodiversity, the carbon cycle, hydrology, climate, and various biogeochemical feedback mechanisms. These interrelationships play a crucial role in shaping future environmental outcomes and influencing climate change. Additionally, the review critically assesses existing research on the correlation between LULCC and climate, along with associated processes. The paper delves into ongoing investigations surrounding LULCC, its ecological repercussions, and its influence on the regulatory mechanisms of terrestrial ecosystems, ecosystem services, atmospheric chemistry, and the climate. It underscores the importance of developing land-use strategies rooted in a scientific evaluation of LULCC to promote sustainability and wider economic development, with a particular emphasis on addressing climate change. Global agreements and protocols also call for land use strategies that consider the implications of environmental degradation.

Coastal inundation regime moderates the short-term effects of sediment and soil additions on seawater oxygen and greenhouse gas dynamics: a microcosm experiment

The frequency and persistence of tidal inundation varies along the coastal terrestrial-aquatic interface, from frequently inundated wetlands to rarely inundated upland forests. This inundation gradient controls soil and sediment biogeochemistry and influence the exchange of soils and sediments from terrestrial to aquatic domains. Although a rich literature exist on studies of the influence of tidal waters on the biogeochemistry of coastal ecosystem soils, few studies have experimentally addressed the reverse question: How do soils (or sediments) from different coastal ecosystems influence the biogeochemistry of the tidal waters that inundate them? To better understand initial responses of coastal waters that flood coastal wetlands and uplands, we conducted short-term laboratory experiments where seawater was amended with sediments and soils collected across regional gradients of inundation exposure (i.e., frequently to rarely inundated) for 14 sites across the Mid-Atlantic, USA. Measured changes in dissolved oxygen and greenhouse gas concentrations were used to calculate gas consumption or production rates occurring during seawater exposure to terrestrial materials. We also measured soil and water physical and chemical properties to explore potential drivers. We observed higher oxygen consumption rates for seawater incubated with soils/sediments from frequently inundated locations and higher carbon dioxide production for seawater incubated with soils from rarely inundated transect locations. Incubations with soil from rarely inundated sites produced the highest global warming potential, primarily driven by carbon dioxide and secondarily by nitrous oxide. We also found environmental drivers of gas rates varied notably between transect locations. Our findings indicate that seawater responses to soil and sediment inputs across coastal terrestrial-aquatic interfaces exhibit some consistent patterns and high intra- and inter-site variability, suggesting potential biogeochemical feedback loops as inundation regimes shift inland.

Streamflow decreases in response to acid deposition in a subtropical forest watershed in China

Streamflow reductions have been attributed to the impacts of soil nutrient availability on plant transpiration, connecting soil biogeochemical and hydrological processes. Here we conducted a plot-scale acid addition experiment and monitored long-term hydrology in a subtropical watershed to provide direct evidence for the underlying mechanisms of these connections. These results showed that acid deposition enhanced plant growth and thus increased plant transpiration in the early treatment period. It indicates that plants can increase their transport of water and nutrients to satisfy physiological demands under continuous acid deposition. Acid deposition mainly contributed to increased evapotranspiration and decreased streamflow at the watershed scale. These results provide complementary evidence of plants adjusting to acid deposition-induced changes in soil nutrient availability and these acclimations result in streamflow reductions at the watershed scale. Our results call for integrating forest biogeochemical feedback into watershed hydrology.

Cascade damming impacts on microbial mediated nitrogen cycling in rivers

Rivers display essential role in nitrogen (N) cycling in terrestrial and aquatic ecosystems, but now they are suffering from damming worldwide, especially from cascade damming. Despite of the importance of microorganisms in biogeochemical nutrient cycling, little attention has been paid to microbial functional biogeography under damming disturbances. Here, the Geochip microarray was applied to investigate the microbial mediated N cycling across the single-dammed Yarlung Tsangpo-Brahmaputra River and the cascade-dammed Lancang-Mekong River in southwest China. Our results showed that the N cycling processes (nitrogen fixation, ammonification, denitrification, nitrification and anammox) were stimulated in reservoirs in both rivers and the enhancement was inversely coupled with hydraulic retention time, but the recovery of N-cycling gene abundance in downstream of dam was intervened by cascade damming. Moreover, N-cycling gene composition was significantly altered in the single-dammed river, while no remarkable change was found in the cascade-dammed reaches. However, different from the unvaried gene composition, cascade damming intervened the recovery of N-cycling gene flow connectivity and resulted in the continuous decrease of connectivity in cascade damming reaches. In addition, in the single-dammed river, nutrients were the important drivers for variation in gene abundance, while they did not influence gene composition. Meanwhile, the abundance and composition of N-cycling genes in the cascade-dammed river were both significantly correlated to geographical parameters and water physical characteristics. Therefore, our study has vital implications for anticipating microbial functional response and biogeochemical feedback to ongoing cascade damming, contributing to the protection of river ecosystems under river regulation.

Widespread and complex drought effects on vegetation physiology inferred from space

The response of vegetation physiology to drought at large spatial scales is poorly understood due to a lack of direct observations. Here, we study vegetation drought responses related to photosynthesis, evaporation, and vegetation water content using remotely sensed data, and we isolate physiological responses using a machine learning technique. We find that vegetation functional decreases are largely driven by the downregulation of vegetation physiology such as stomatal conductance and light use efficiency, with the strongest downregulation in water-limited regions. Vegetation physiological decreases in wet regions also result in a discrepancy between functional and structural changes under severe drought. We find similar patterns of physiological drought response using simulations from a soil–plant–atmosphere continuum model coupled with a radiative transfer model. Observation-derived vegetation physiological responses to drought across space are mainly controlled by aridity and additionally modulated by abnormal hydro-meteorological conditions and vegetation types. Hence, isolating and quantifying vegetation physiological responses to drought enables a better understanding of ecosystem biogeochemical and biophysical feedback in modulating climate change.

Biodiversity and Stoichiometric Plasticity Increase Pico‐Phytoplankton Contributions to Marine Net Primary Productivity and the Biological Pump

AbstractEarth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low‐latitude oligotrophic systems harboring pico‐sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon‐climate system. We have added three pico‐phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico‐eukaryote populations that contribute significantly to productivity and export within the subtropical to mid‐latitude transition zone, with the western subtropical regions of each basin supporting the most P‐poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C‐rich/P‐poor stoichiometry, preferentially inhabited by pico‐eukaryotes and diatoms. Collectively, pico‐phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico‐phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.

Effects of Elevated Ozone Exposure on Regional Meteorology and Air Quality in China Through Ozone‐Vegetation Coupling

AbstractOzone is a phytotoxic pollutant that could damage the vegetation growth and lead to complicated impacts on air quality through meteorological and biogeochemical feedbacks. This study implements a semi‐empirical parameterization regarding the impacts of ozone exposure on photosynthesis rate and stomatal resistance into the Noah‐Multi‐parameterization (Noah‐MP) dynamic vegetation module of Weather Research and Forecasting with Chemistry (WRF/Chem) model. The gaseous dry deposition and biogenic emission algorithms are also coupled with Noah‐MP to enable the ozone‐vegetation coupling. This model reproduces the near‐surface meteorology, air pollutants and vegetation physiology in China, with a spatial correlation coefficient more than 0.9 and normalized mean bias from −0.19 to 0.42. The optimized model also improves the simulations of vegetation physiology (e.g., a reduction of model error by 18%–32%) and ozone dry deposition velocity. The elevated ozone damages plant photosynthesis, and decreases the national gross primary productivity (−28.85%) and leaf area index (−17.41%). The plant transpiration and surface heat flux, as well as air temperature (e.g., up to +0.16°C in summer) and other associated meteorological variables are also altered, finally contributing to 0.49 μg m−3 increase of surface ozone. Otherwise, the suppressed vegetation LAI and biogenic emissions, as well as the lower dry deposition velocity in response to the ozone‐vegetation coupling contribute to the remaining ozone changes by −1.07 μg m−3 and 1.18 μg m−3, jointly constituting the complicated ozone‐vegetation feedbacks on air quality. Our results highlight the necessity of including the ozone‐vegetation coupling in models for reliable prediction of regional climate and air quality.

Deoxygenation turns the coastal Red Sea lagoons into sources of nitrous oxide

Direct measurements of dissolved N2O concentrations, fluxes and saturation percentages undertaken for the first time in two coastal lagoons — Al-Shabab and Al-Arbaeen, along the east coast of the Red Sea, revealed the region as a significant source of N2O to the atmosphere. The exacerbated dissolved inorganic nitrogen (DIN) from various anthropogenic sources led to substantial oxygen depletion in both the lagoons, which turned to bottom anoxia at Al-Arbaeen lagoon during the spring season. We assume that the accumulation of N2O is caused by nitrifier-denitrification in the hypoxic/anoxic boundaries. In fact, the results indicated that oxygen-depleted bottom waters favoured denitrification when the oxygenated surface waters recorded nitrification signals. Overall, the N2O concentration ranged from 109.4 to 788.6 nM (40.6–325.6 nM) in spring and 58.7 to 209.8 nM (35.8–89.9 nM) in winter in the Al-Arbaeen (Al-Shabab) lagoon. The N2O flux ranged from 647.1 to 1763.2 μmol m−2 day−1 (85.9 to 160.2 μmol m−2 day−1) and 112.5 to 150.8 μmol m−2 day−1 (76.1 to 88.7 μmol m−2 day−1) in the spring and winter respectively, in the Al-Arbaeen (Al-Shabab) lagoons. The ongoing developmental activities may worsen the current situation of hypoxia and associated biogeochemical feedbacks; therefore, the present results underline the need for continuous monitoring of both lagoons to restrict more severe oxygen depletion in future.

Comparable biophysical and biogeochemical feedbacks on warming from tropical moist forest degradation

Comparable biophysical and biogeochemical feedbacks on warming from tropical moist forest degradation

Protecting nature's most effective carbon sink: an important role for fungi in peatlands.

This article is a Commentary on Shao et al. (2023), 238: 80–95.

Source data for 'Comparable biophysical and biogeochemical feedbacks on warming from tropical moist forest degradation'

Source data for 'Comparable biophysical and biogeochemical feedbacks on warming from tropical moist forest degradation'

Drip irrigation improves spring wheat water productivity by reducing leaf area while increasing yield

To mitigate the climate change-induced water shortage and realize the sustainable development of agriculture, drip irrigation, a more efficient water-saving irrigation method, has been intensively implemented in most arid agricultural regions in the world. However, compared to traditional border irrigation, how drip irrigation affects the biophysical conditions in the cropland and how crops physiologically respond to changes in biophysical conditions in terms of water, heat and carbon exchange remain largely unknown. In view of the above situation, to reveal the mechanism of drip irrigation in improving spring wheat water productivity, paired field experiments based on drip irrigation and border irrigation were conducted to extensively monitor water and heat fluxes at a typical spring wheat field (Triticum aestivum L.) in Northwest China during 2017–2020. The results showed that drip irrigation improved yield by 10.3 % and crop water productivity (i.e., yield-to-evapotranspiration-ratio) by 15.6 %, but reduced LAI by 16.9 % in contrast with border irrigation. Under drip irrigation, the lateral development of spring wheat roots was promoted by higher soil temperature combined with frequent dry-wet alternation in the shallow soil layer (0–20 cm), which was the basis for efficient absorption of water and fertilizer, as well as efficient formation of photosynthate. Meanwhile, drip irrigation increased net radiation and decreased latent heat flux by inhibiting leaf growth, thereby increased sensible heat, causing a higher soil temperature (+1.10 ℃) and canopy temperature (+1.11 ℃). Further analysis proved that soil temperature was the key factor affecting yield formation. Based on the above conditions, the decrease in leaf distribution coefficient (−0.030) led to the decrease in evapotranspiration (−5.7 %) and the increase in ear distribution coefficient (+0.029). Therefore, drip irrigation emphasized the role of soil moisture in the soil-plant-atmosphere continuum, enhanced crop activity by increasing field temperature, especially soil temperature, and finally improved yield and water productivity via carbon reallocation. The study revealed the mechanism of drip irrigation for improving spring wheat yield, and would contribute to improving Earth system models in representing agricultural cropland ecosystems with drip irrigation and predicting the subsequent biophysical and biogeochemical feedbacks to climate change.

Modeling demographic-driven vegetation dynamics and ecosystem biogeochemical cycling in NASA GISS's Earth system model (ModelE-BiomeE v.1.0)

Abstract. We developed a demographic vegetation model, BiomeE, to improve the modeling of vegetation dynamics and ecosystem biogeochemical cycles in the NASA Goddard Institute of Space Studies' ModelE Earth system model. This model includes the processes of plant growth, mortality, reproduction, vegetation structural dynamics, and soil carbon and nitrogen storage and transformations. The model combines the plant physiological processes of ModelE's original vegetation model, Ent, with the plant demographic and ecosystem nitrogen processes that have been represented in the Geophysical Fluid Dynamics Laboratory's LM3-PPA. We used nine plant functional types to represent global natural vegetation functional diversity, including trees, shrubs, and grasses, and a new phenology model to simulate vegetation seasonal changes with temperature and precipitation fluctuations. Competition for light and soil resources is individual based, which makes the modeling of transient compositional dynamics and vegetation succession possible. Overall, the BiomeE model simulates, with fidelity comparable to other models, the dynamics of vegetation and soil biogeochemistry, including leaf area index, vegetation structure (e.g., height, tree density, size distribution, and crown organization), and ecosystem carbon and nitrogen storage and fluxes. This model allows ModelE to simulate transient and long-term biogeophysical and biogeochemical feedbacks between the climate system and land ecosystems. Furthermore, BiomeE also allows for the eco-evolutionary modeling of community assemblage in response to past and future climate changes with its individual-based competition and demographic processes.