Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “climate dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

Data-driven global ocean modeling for seasonal to decadal prediction

Accurate modeling of ocean dynamics is crucial for enhancing our understanding of complex ocean circulation processes, predicting climate variability, and tackling challenges posed by climate change. Although great efforts have been made to improve traditional numerical models, predicting global ocean variability over multiyear scales remains challenging. Here, we propose ORCA-DL (Oceanic Reliable foreCAst via Deep Learning), a data-driven three-dimensional ocean model for seasonal to decadal prediction of global ocean dynamics. ORCA-DL accurately simulates the three-dimensional structure of global ocean dynamics with high physical consistency and outperforms state-of-the-art numerical models in capturing extreme events, including El Niño–Southern Oscillation and upper ocean heat waves. Moreover, ORCA-DL stably emulates ocean dynamics at decadal timescales, demonstrating its potential even for skillful decadal predictions and climate projections. Our results demonstrate the high potential of data-driven models for providing efficient and accurate global ocean modeling and prediction.

Science & Technology - Other Topics↗

Response of the Current Climate to Land‐Ocean Contrasts in Parameterized Cumulus Entrainment

Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.

54 ENVIRONMENTAL SCIENCES↗

w Mel replacement of dengue-competent mosquitoes is robust to near-term climate change

Rising temperatures are impacting the range and prevalence of mosquito-borne diseases. A promising biocontrol technology replaces wild mosquitoes with those carrying the virus-blocking Wolbachia bacterium. Because the most widely used strain, wMel, is adversely affected by heat stress, we examined how global warming may influence wMel-based replacement. We simulated interventions in two locations with successful field trials using Coupled Model Intercomparison Project Phase 5 climate projections and historical temperature records, integrating empirical data on wMel’s thermal sensitivity into a model of Aedes aegypti population dynamics to evaluate introgression and persistence over one year. We show that in Cairns, Australia, climatic futures necessitate operational adaptations for heatwaves exceeding two weeks. In Nha Trang, Vietnam, projected heatwaves of three weeks and longer eliminate wMel under the most stringent assumptions of that symbiont’s thermal limits. We conclude that this technology is generally robust to near-term (2030s) climate change. Accelerated warming may challenge this in the 2050s and beyond.

54 ENVIRONMENTAL SCIENCES↗

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗

Global maximum carboxylation rate trait scaling in the Sheffield Dynamic Global Vegetation Model 1901-2012

This data package contains model output data and climate input data associated with the New Phytologist publication: Walker, A.P., Quaife, T., van Bodegom, P.M., De Kauwe, M.G., Keenan, T.F., Joiner, J., Lomas, M.R., MacBean, N., Xu, C., Yang, X., Woodward, F.I., 2017. The impact of alternative trait-scaling hypotheses for the maximum photosynthetic carboxylation rate (Vcmax) on global gross primary production. New Phytologist 215, 1370–1386. https://doi.org/10.1111/nph.14623. This dataset contains maximum carboxylation rate at 25 ºC (Vcmax,25), Gross Primary Production (GPP), Net Biome Production (NBP), vegetation carbon, and soil carbon data from Sheffield Dynamic Global Vegetation Model (SDGVM) for multiple trait-scaling methods for scaling Vcmax,25 across Earth, and alternative instantaneous temperature scaling functions. Vcmax,25 values are mean top-leaf values taken during the growing seasons, which were defined as periods during which leaf area index (LAI) was greater than one. Vcmax,25 values were reported prior to scaling of Vcmax by water-stress or leaf-age. The attached metadata file (SDGVM_vcmax_dataset.docx) describes the methods used and explains the organization of data within the included .zip file. This metadata file is also included in PDF format.

54 ENVIRONMENTAL SCIENCES↗

Late Holocene history of Tokotan Lake (Kuril Archipelago, Russian Far East): The use of lacustrine records for paleoclimatic reconstructions from geologically dynamic settings

A high-resolution, multiproxy analysis of a core from Tokotan Lake, Urup Island, illustrates the complex paleoenvironmental signals that can be preserved in sedimentary deposits from geologically active areas, such as the Kuril Archipelago. Diatom and geochemical analyses of the Tokotan core suggest a shallowing and lessening of the lake area over the last ~3450 cal BP. Without other proxies, the record could be interpreted as showing a general decrease in precipitation over the Late Holocene. However, palynological data indicate no change in vegetation, thereby implying no significant climatic variation. Neither of these paleoclimatic interpretations agrees with Late Holocene reconstructions from the Kuril Archipelago or the Russian mainland. Paleomagnetic data imply variable sediment sources during the lake's history, indicating alterations within the Tokotan catchment (e.g., shifts in gradients or sediment loads in the inflowing streams). Such changes would be particularly important in the Shabalina River, which today has formed a partial dam across the lake's outlet. Here, the results of the Tokotan study indicate that other factors (e.g., volcanism, seismic activity, sea-level changes) rather than climate were more probable determinants of the Late Holocene lake record. Given the geologic dynamism of the region, it is perhaps not surprising that paleo-records from the Kuril Islands display an inconsistency in the numbers and/or duration of proposed Holocene climatic oscillations, underscoring the importance of multiproxy investigations in areas that are or have been geologically dynamic.

58 GEOSCIENCES↗

Recent Trends in Transport of Surface Carbonaceous Aerosols to the Upper-Troposphere-Lower-Stratosphere Linked to Expansion of the Asian Summer Monsoon Anticyclone

Using Modern-Era Retrospective analysis for Research and Applications Version-2 (MERRA-2) re-analyses, we have examined recent trends (2000–2019) in transport of surface carbonaceous aerosols (CAs), that is, organic carbon and black carbon, to the upper troposphere and lower stratosphere (UTLS) during the Asian summer monsoon (ASM). We find a significant increased concentration of UTLS CA, linked to a linear trend in an expansion of the Asian Summer Monsoon Anticyclone (ASMA). Over core monsoon latitudes (25°–35°N), the trends in UTLS CA are enabled by increased upward transport from surface sources via an intensified monsoon meridional circulation, with increased latent heating over the southern Tibetan Plateau and foothill regions, enhanced by feedback processes associated with radiative heating by UTLS CA. In the extra-tropics, increased UTLS CA stems primarily from an extended source of increased wildfire emissions over eastern Siberia and northern Asia, coincident with a large-scale anomalous anticyclone with enhanced surface warming and drying near (80°–120°E, 55°–70°N). Here, CAs are transported upward from the surface to the ULS, likely by increased pyro-convection associated with enhanced wildfires, and enter the tropical UTLS via increased equatorward transport on the eastern flanks of anomalous upper-level anticyclones, coupled to the expanded ASMA. Overall, the increasing UTLS CA trends associated with expansion of the ASMA are consistent with a hydrostatic expansion of a warming troposphere, reflected in a rise in the tropical tropopause and consequential dynamical adjustments of the ASM subtropical jetstream, modulating the climate of the greater ASM region.

54 ENVIRONMENTAL SCIENCES↗

Field experiments show no consistent reductions in soil microbial carbon in response to warming

Soil microbes play an essential role in maintaining soil functions and services, but the dynamics of soil microbial biomass carbon (MBC) under global climate change remain unclear. Recently, Patoine et al. combined a global MBC data set with Random Forest modeling and reported that global MBC decreased over 1992–2013, mainly driven by increasing temperatures. Contrarily, using MBC field observations from soil warming manipulation experiments and in-situ long-term measurements across the globe, we found that MBC showed no significant changes under soil warming. Our findings indicate that soil MBC is unlikely to have decreased significantly due to the global warming of 0.28 °C during 1992–2013, and that further mechanistic studies are needed to understand potential changes in MBC under climate change.

54 ENVIRONMENTAL SCIENCES↗

Boreal–Arctic wetland methane emissions modulated by warming and vegetation activity

Abstract Wetland methane (CH 4 ) emissions over the Boreal–Arctic region are vulnerable to climate change and linked to climate feedbacks, yet understanding of their long-term dynamics remains uncertain. Here, we upscaled and analysed two decades (2002–2021) of Boreal–Arctic wetland CH 4 emissions, representing an unprecedented compilation of eddy covariance and chamber observations. We found a robust increasing trend of CH 4 emissions (+8.9%) with strong inter-annual variability. The majority of emission increases occurred in early summer (June and July) and were mainly driven by warming (52.3%) and ecosystem productivity (40.7%). Moreover, a 2 °C temperature anomaly in 2016 led to the highest recorded annual CH 4 emissions (22.3 Tg CH 4 yr −1 ) over this region, driven primarily by high emissions over Western Siberian lowlands. However, current-generation models from the Global Carbon Project failed to capture the emission magnitude and trend, and may bias the estimates in future wetland CH 4 emission driven by amplified Boreal–Arctic warming and greening.

54 ENVIRONMENTAL SCIENCES↗

Using Sentinel-1 and GRACE satellite data to monitor the hydrological variations within the Tulare Basin, California

Abstract Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth’s surface. California’s Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018–2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

54 ENVIRONMENTAL SCIENCES↗

Rapid organic carbon spiraling in a headwater stream linked with streamflow, biogeochemistry, and canopy phenology

Headwater streams are abundant worldwide and important to global biogeochemical cycles, serving as critical processors and transporters of C. C spiraling is a useful way to understand the retention and mineralization of organic C (OC) in streams. However, analyses of seasonal and interannual variability in OC spiraling are currently limited. In this study, we aimed to understand the temporal patterns and driving mechanisms of OC spiraling, which will inform our understanding of future OC changes under climate change. We used 7 y of daily data in a small headwater stream (Walker Branch, Tennessee, USA) to assess seasonal and interannual variability in OC spiraling length (S OC ) and mineralization velocity (v fOC ), as well as their potential related variables. On average, S OC in Walker Branch was ~10× shorter than in previously studied small streams, indicating strong connections between the water column and the benthic environment where OC mineralization mostly takes place. OC spiraling was faster during the more biologically active periods of spring and autumn compared with more elongated OC spiraling in summer and winter, when OC retention was lower and downstream transport was higher. Gross primary production (GPP) was most strongly related to S OC and v fOC . Photosynthetically active radiation (PAR) and NO 3 − were also positively and negatively related to v fOC , respectively. Trends toward earlier and longer canopy cover and reduced GPP and PAR may result in longer S OC and slower v fOC , reducing localized instream processing of OC and potentially shunting more OC downstream. However, long-term observations indicate reduced NO 3 − at Walker Branch, suggesting opposing effects to those of GPP and PAR, leading to faster v fOC and greater OC retention. Time-series analyses of OC spiraling in streams can enhance our understanding of current and future responses of OC processing and downstream transport to climate change, as well as implications for downstream OC dynamics.

biological activity↗

Thermal acclimation of stem respiration implies a weaker carbon-climate feedback

The efflux of carbon dioxide (CO2) from woody stems, a proxy for stem respiration, is a critical carbon flux from ecosystems to the atmosphere, which increases with temperature on short timescales. However, plants acclimate their respiratory response to temperature on longer timescales, potentially weakening the carbon-climate feedback. The magnitude of this acclimation is uncertain despite its importance for predicting future climate change. We develop an optimality-based theory dynamically linking stem respiration with leaf water supply to predict its thermal acclimation. We show that the theory accurately reproduces observations of spatial and seasonal change. We estimate the global value for current annual stem CO2 efflux as 27.4 ± 5.9 PgC. By 2100, incorporating thermal acclimation reduces projected stem respiration without considering acclimation by 24 to 46%, thus reducing land ecosystem carbon emissions.

Zhang, Han↗

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Environmental sciences↗

CLM simulated data at three FLUXNET sites and three SNOTEL sites in the Colorado River Basin

This data package includes 6 csv files that include Community Land Model (CLM) simulations for three Fluxnet sites (i.e., US-NR1, US-GLE and US-VCM) and three SNOTEL sites (i.e., ER-BT, ER-SP and ER-PK) in the Upper Colorado River Basin area. Data inputs include meteorological forcing data (obtained from Fluxnet and SNOTEL databases) and CLM-simulated evapotranspiration (ET) data from 2005 to 2016. These datasets support the development of subseasonal regime approach based on statistics driven Hidden Markov Model to assess intra- and inter-annual variability of ET, which improves our understanding of climate change's impacts on water resources and ecosystem dynamics in the Colorado River Basin.

54 ENVIRONMENTAL SCIENCES↗

Exploring Trait Trade-Offs for Fungal Decomposers in a Southern California Grassland

Fungi are important decomposers in terrestrial ecosystems, so their responses to climate change might influence carbon (C) and nitrogen (N) dynamics. We investigated whether growth and activity of fungi under drought conditions were structured by trade-offs among traits in 15 fungal isolates from a Mediterranean Southern California grassland. We inoculated fungi onto sterilized litter that was incubated at three moisture levels (4, 27, and 50% water holding capacity, WHC). For each isolate, we characterized traits that described three potential lifestyles within the newly proposed “YAS” framework: growth yield, resource acquisition, and stress tolerance. Specifically, we measured fungal hyphal length per unit litter decomposition for growth yield; the potential activities of the extracellular enzymes cellobiohydrolase (CBH), β -glucosidase (BG), β -xylosidase (BX), and N-acetyl- β - D -glucosaminidase (NAG) for resource acquisition; and ability to grow in drought vs. higher moisture levels for drought stress tolerance. Although, we had hypothesized that evolutionary and physiological trade-offs would elicit negative relationships among traits, we found no supporting evidence for this hypothesis. Across isolates, growth yield, drought stress tolerance, and extracellular enzyme activities were not significantly related to each other. Thus, it is possible that drought-induced shifts in fungal community composition may not necessarily lead to changes in fungal biomass or decomposer ability in this arid grassland.

59 BASIC BIOLOGICAL SCIENCES↗

Snowfall and snow accumulation during the MOSAiC winter and spring seasons

Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 mm SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio 2 pluviometer and the OTT Parsivel 2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based K a -band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the K a -band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 mm and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98–114 mm.

54 ENVIRONMENTAL SCIENCES↗

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗

Belowground cross-trophic networks impact CH4 and CO2 emissions in degraded alpine peatlands

Belowground organisms forming complex cross-trophic ecological networks are essential for maintaining peatland carbon stability and energy flow. However, how peatland degradation affects the biodiversity and cross-trophic ecological networks of soil communities remains poorly understood. Here, we examined the degradation effects on soil prokaryotes (i.e., bacteria, archaea), fungi and nematodes in alpine peatlands on the eastern Tibetan Plateau, characterized by varying water table depths (indicating degradation levels). We found that peatland degradation, accompanied by significant shifts in soil moisture and pH (P < 0.05), reduced the taxonomic richness and phylogenetic diversity of prokaryotes, fungi, and nematodes, particularly in deeper soil layers (20–50 cm). Crucially, peatland degradation weakened potential cross-trophic interactions within bipartite networks of prokaryotes-nematodes and fungi-nematodes, resulting in less than 6.5 %–28.8 % of unchanged modules. Degradation-induced changes in soil moisture and pH were identified as primary drivers of biodiversity loss and network restructuring. Furthermore, such changes of belowground cross-trophic networks (particularly prokaryote-nematode) were significantly correlated with greenhouse gas emissions, such as decreased CO2 emissions, maintained CH4 emissions (leading to a higher CH4/CO2 ratio in deep layers), and reduced temperature sensitivity (Q10) of soil respiration. These findings underscore the critical need to protect soil biodiversity and cross-trophic networks in peatlands, particularly under the threat of climate change, to preserve peatland carbon stocks and maintain ecosystem stability. Our findings highlight that belowground cross-trophic networks are pivotal to decipher soil carbon dynamics of degraded peatlands and project the fate of peatland carbon stocks under future climate change scenarios.

Cross-trophic networks↗