Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “carbonization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Soil carbon stocks from soil organic carbon at the Teller 47, Teller 27, and Kougarok 86 field sites, Seward Peninsula, Alaska, 2021 – 2023

Carbon stored in permafrost soils represents the largest terrestrial storage and hence has key significance for the global carbon cycle, yet the site-specific estimates of permafrost soil carbon are limited. To better understand permafrost carbon storage, soil organic content (SOC) was measured from soils collected at the Teller 27, Teller 47, and Kougarok 86 field sites on the Seward Peninsula, Alaska in the summers of 2019 and 2021-2023. Soil samples were collected from locations with and without near-surface permafrost to varied depths within the permafrost, where present. In addition to sample collection for soil organic carbon analysis, bulk density samples were also collected within each soil profile to allow for calculation of carbon stocks. Carbon to nitrogen ratios are also reported. Location data were collected at each observation point using a handheld GPS unit with waypoint averaging to reduce location measurement error to within 3 m. SOC was measured using elemental analysis in the laboratory. This dataset contains a *.csv file of soil carbon stocks, a *.kml file of measurement locations, a *.pdf user guide, a *.csv data dictionary, and a *.csv file level metadata.

EARTH SCIENCE > CRYOSPHERE > FROZEN GROUND↗
NASA NTRS2018

The Ocean Carbon States Database: A Proof-of-Concept Application of Cluster Analysis in the Ocean Carbon Cycle

In this paper, we present a database of the basic regimes of the carbon cycle in the ocean, the 'ocean carbon states', as obtained using a data mining/pattern recognition technique in observation-based as well as model data. The goal of this study is to establish a new data analysis methodology, test it and assess its utility in providing more insights into the regional and temporal variability of the marine carbon cycle. This is important as advanced data mining techniques are becoming widely used in climate and Earth sciences and in particular in studies of the global carbon cycle, where the interaction of physical and biogeochemical drivers confounds our ability to accurately describe, understand, and predict CO2 concentrations and their changes in the major planetary carbon reservoirs. In this proof-of-concept study, we focus on using well-understood data that are based on observations, as well as model results from the NASA Goddard Institute for Space Studies (GISS) climate model. Our analysis shows that ocean carbon states are associated with the subtropical-subpolar gyre during the colder months of the year and the tropics during the warmer season in the North Atlantic basin. Conversely, in the Southern Ocean, the ocean carbon states can be associated with the subtropical and Antarctic convergence zones in the warmer season and the coastal Antarctic divergence zone in the colder season. With respect to model evaluation, we find that the GISS model reproduces the cold and warm season regimes more skillfully in the North Atlantic than in the Southern Ocean and matches the observed seasonality better than the spatial distribution of the regimes. Finally, the ocean carbon states provide useful information in the model error attribution. Model air-sea CO2 flux biases in the North Atlantic stem from wind speed and salinity biases in the subpolar region and nutrient and wind speed biases in the subtropics and tropics. Nutrient biases are shown to be most important in the Southern Ocean flux bias.

carbon cycle↗
NASA NTRS2020

Oxygen and carbon stable isotope composition of the weathering Mg-carbonates on the surface of the LEW 85320 ordinary chondrite: Revisited

New δ13C and δ18O values were measured to constrain the formation mechanisms of two previously identified generations (Antarctica and Houston) of terrestrial weathering carbonates, nesquehonite, which formed on the surface of the Lewis Cliff 85320 ordinary chondrite. These two nesquehonite generations exhibit a characteristic carbon and oxygen isotopic trend which starts with a linear δ13C-δ18O covariation followed by a continuous increase in δ18O with little to no change in δ13C. Based upon the newly developed nesquehonite-water oxygen isotope thermometry and the measured δ18O value of melted ice from Antarctica, the formation temperature of the Antarctic nesquehonite generation was estimated to be -1.9 ± 3℃. Houston nesquehonite generation was most likely formed from an evaporative residue of the melted Antarctic ice at 30 ± 4℃. Modeling the isotopic trend of the two generations suggests that evaporation created large δ18O variations in the nesquehonite while carbon isotopes were stabilized by exchange between the parent liquid and atmospheric CO2. Thus, we suggest that carbonates formed by terrestrial weathering in Antarctica should possess a widespread in δ18O values and a narrow range in δ13C values under dry and cold conditions. Finally, the observed wide spread of δ18O in Martian carbonates (e.g., Allan Hills 84001 and Nakhla) can be explained by evaporation near 0℃ and may have formed by a similar mechanism to that of the nesquehonite in the LEW 85320. Meanwhile, the wide spread of δ13C suggests that two carbon reservoirs with distinct isotope compositions on Mars (e.g., the atmospheric carbon and magmatic carbon) may be actively participating in the formation process.

weathering↗
NASA NTRSDate not supplied

Bosch Carbon Catalyst Regeneration in Water I: Bosch Carbon Removal from Reacted Beads Utilizing Mechanical Agitation

The ability to recover oxygen from metabolic carbon dioxide is crucial for crewed long duration missions in which the atmosphere must be recycled. Two physiochemical methods are primarily employed for this purpose: the Bosch process and the Sabatier reaction. This paper focuses on the Bosch process. It utilizes hydrogen in the presence of a catalyst to form water and solid carbon out of metabolic carbon dioxide. The solid carbon product deposits onto the catalyst and eventually causes reactor over pressurization issues. Multiple methods are currently being investigated to regenerate the catalyst’s surface and break up the pressure causing carbon clogs. This paper explores whether water and mechanical abrasion would work to remove Bosch carbon by testing multiple washing conditions. The different washing conditions involved agitation intensity, heating, duration, and volume of water. It was found that it is possible to break up the clumps of carbon formed in the reactor through agitation with room temperature water in one hour and shows promise for completely cleaning the carbon off the beads’ surface.

Agitation↗
DOE OSTI2020

Enhanced carbon-transfer and -utilization efficiencies achieved using membrane carbonation with gas sources having a range of CO 2 concentrations

The economic viability of microalgal biofuels relies on increasing productivity in a cost-effective manner. As microalgal biomass contains > 50% carbon, a high rate of CO 2 delivery is required for high productivity, and inefficient CO 2 delivery amplifies operating costs. Membrane carbonation using non-porous hollow fiber membranes can ideally deliver CO 2 without bubble formation and high carbon transfer efficiency. Because CO 2 streams from industrial resources are not 100% CO 2 , the buildup of inert gasses can significantly lower the CO 2 delivery rate when the distal end of the membrane is closed. To overcome the buildup of inert gases, we managed the distal end of the membranes with three different approaches: fully open end, restricted bleed valve, and restricted bleed valve with pH-actuated venting. For all approaches, CO 2 was delivered to membranes ondemand based on a pH set point. Evaluating a wide range of CO 2 concentrations (10% to 100%), we found that all approaches eliminated the buildup of inert gases, could maintain target pH values and gave the same biomass productivities and carbon distributions. However, carbon transfer efficiency depended on the operation of the distal end. Fully open-end operation gave a poor carbon transfer efficiency because of excessive loss of CO 2 from the distal end. However, restricting the exit flow rate to ≤4 cm 3 /min mitigated the problems of excessive CO 2 loss, but without incurring a large loss of CO 2 -delivery flux. For the continuous cultivation, combining a restricted bleed valve with pH-actuated venting improved the carbon-transfer efficiency and -utilization efficiencies up to 85% and 67%, respectively, with a sufficient CO 2 delivery flux.

42 ENGINEERING↗
DOE OSTI2025

The stable carbon isotope fractionation of methanogenesis products at complete carbon consumption

The stable carbon isotope signature (δ 13 C) of methane (CH 4 ) is used to discriminate between biological, thermogenic, and abiotic sources. Methanogens, or methane producing archaea, inhabit a broad range of chemical conditions. Many of these environments are replete in dissolved inorganic carbon (DIC), causing isotopically depleted δ 13 C biogenic CH 4 . However, some extreme environments inhabited by methanogens, such as serpentinising systems, exhibit low carbon dioxide (CO 2 ) availability, replete H 2 , and isotopically enriched δ 13 C CH 4 that is outside the known biogenic range. We measured the δ 13 C of CO 2 , biomass, lipids, and CH 4 during hydrogenotrophic methanogenesis under hydrogen replete conditions with a limited carbon pool to investigate carbon isotope dynamics at complete DIC consumption. As theory predicts, we found that the final, accumulated methane δ 13 C values closely reflect the δ 13 C of the initial DIC supply, and that methane is more 13 C enriched than biomass and lipids. This provides the first experimental evidence that methanogens can achieve complete carbon consumption and thus can produce accumulated CH 4 products that isotopically reflect the initial CO 2 . These data show that the range of possible δ 13 C values from biogenic methane needs to be expanded for natural environments impacted by extreme carbon limitation.

biomass↗
NASA NTRS2008

Ocean Carbon and Biogeochemistry Scoping Workshop on Terrestrial and Coastal Carbon Fluxes in the Gulf of Mexico, St. Petersburg, FL

Despite their relatively small surface area, ocean margins may have a significant impact on global biogeochemical cycles and, potentially, the global air-sea fluxes of carbon dioxide. Margins are characterized by intense geochemical and biological processing of carbon and other elements and exchange large amounts of matter and energy with the open ocean. The area-specific rates of productivity, biogeochemical cycling, and organic/inorganic matter sequestration are high in coastal margins, with as much as half of the global integrated new production occurring over the continental shelves and slopes (Walsh, 1991; Doney and Hood, 2002; Jahnke, in press). However, the current lack of knowledge and understanding of biogeochemical processes occurring at the ocean margins has left them largely ignored in most of the previous global assessments of the oceanic carbon cycle (Doney and Hood, 2002). A major source of North American and global uncertainty is the Gulf of Mexico, a large semi-enclosed subtropical basin bordered by the United States, Mexico, and Cuba. Like many of the marginal oceans worldwide, the Gulf of Mexico remains largely unsampled and poorly characterized in terms of its air-sea exchange of carbon dioxide and other carbon fluxes. The goal of the workshop was to bring together researchers from multiple disciplines studying terrestrial, aquatic, and marine ecosystems to discuss the state of knowledge in carbon fluxes in the Gulf of Mexico, data gaps, and overarching questions in the Gulf of Mexico system. The discussions at the workshop were intended to stimulate integrated studies of marine and terrestrial biogeochemical cycles and associated ecosystems that will help to establish the role of the Gulf of Mexico in the carbon cycle and how it might evolve in the face of environmental change.

Air sea exchanges↗
DOE OSTI2026

Wyoming Carbon Blueprint A Roadmap for Statewide Carbon Capture and Storage

Wyoming stands at a pivotal point in defining its energy future, and the Wyoming Carbon Blueprint positions the state to lead the nation in carbon management through a coordinated buildout of carbon capture, transportation, utilization, and geologic storage. Building on Wyoming’s long history as a major producer of coal, oil, and natural gas, and supported by favorable policies, skilled labor, and robust subsurface resources, the Blueprint outlines a pathway to develop one of the world’s first open-access carbon hubs capable of capturing between 10 and 25 million tonnes of CO₂ per year within five years, and up to 45 million tonnes within a decade. In doing so, it advances the concept of “No Carbon Left Behind”: a statewide strategy that prevents stranded CO₂ by aligning capture, transport, and storage systems into a single coherent infrastructure network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dual 3D Ceramic Textile Electrodes: Fast Kinetics for Carbon Oxidation Reaction and Oxygen Reduction Reaction in Direct Carbon Fuel Cells at Reduced Temperatures

Abstract Direct carbon fuel cells (DCFCs) are an efficient energy‐conversion technology capable of generating electricity with carbon‐dioxide‐capture chemistry with solid carbon as fuels. The efficiency and performance of DCFCs depend on the kinetics of the carbon oxidation reactions (COR) and the oxygen reduction reactions (ORR), each occurring at anode and cathode, respectively. The limited active sites paired with reduced temperatures greatly decrease the efficiency of the electrochemical reactions. Ultraporous dual‐3D ceramic textiles (dual‐3DCT) are integrated into electrolyte‐supported DCFCs to enhance charge and mass transfer at the electrodes. Improved COR at the anode is achieved by the synergy between the 3DCT NiO–Ce 0.8 Gd 0.2 O 1.95 (GDC) structure and optimal carbon fuel choice. In a comparative study, DCFCs using graphitic carbon (GC) as fuel show the best COR performance when compared to DCFCs utilizing alternative fuels such as carbon black (CB) and activated carbon (AC). The 3DCT Sm 0.5 Sr 0.5 CoO 3‐δ –GDC (SSC–GDC) composite cathode shows electrochemical performance superior to that of the conventional screen‐printed SSC–GDC. A peak power density of 392 mW cm −2 at 600 °C is obtained in a DCFC using the 3DCT‐anode/electrolyte/3DCT‐cathode configuration, an unprecedented value for any reported DCFC as of yet. This points toward promising applications of dual‐3DCT electrodes for reduced‐temperature DCFCs.

Bian, Wenjuan↗
DOE OSTI2024

Toward Rational Design of Nickel Catalysts for Thermocatalytic Decomposition of Methane for Carbon Dioxide-Free Hydrogen and Value-Added Carbon Co-Product: A Review

Thermocatalytic decomposition of methane provides opportunities for hydrogen (H2) production with no emission of carbon dioxide. However, high-value carbon products need to be produced for economic deployment of thermocatalytic decomposition and to achieve a minimum H 2 selling price below the U.S, Department of Energy target of $ 1/kg H 2 . In this review, we re-evaluate data on catalyst development reported in the literature and propose correlations between catalyst characteristics, catalytic stability, and properties of carbon co-products. In the first part of the review, growth mechanisms for carbon nanotubes using state-of-the-art chemical vapor deposition are reviewed to catalog the effects of catalyst characteristics, the influence of carbon sources, interactions between metal particles and supports, and metal particle sizes on carbon growth. In the second part, representative developments in mono-, bi-, and tri-metallic nickel catalysts are highlighted. We present kinetic analysis of reactions catalyzed by mono-metallic nickel catalysts, which generates a correlation between metal particle size and catalyst stability. Rational design of Ni-based catalysts for TCD of methane requires attention to the size of the metal particle and effective normalization of the reaction rates. Further attention to the distribution of the metal particle sizes may help identify catalyst properties that contribute longevity and selectivity to processes that use them. While it is tempting to focus on the highest valued carbon products (e.g., CNTs and CFs), analysis of the markets for other carbon products suggests that a more flexible approach may generate comparable returns without the risk associated with specialization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Tuning the reactivity of carbon surfaces with oxygen-containing functional groups

Oxygen-containing carbons are promising supports and metal-free catalysts for many reactions. However, distinguishing the role of various oxygen functional groups and quantifying and tuning each functionality is still difficult. Here we investigate the role of Brønsted acidic oxygen-containing functional groups by synthesizing a diverse library of materials. By combining acid-catalyzed elimination probe chemistry, comprehensive surface characterizations, 15N isotopically labeled acetonitrile adsorption coupled with magic-angle spinning nuclear magnetic resonance, machine learning, and density-functional theory calculations, we demonstrate that phenolic is the main acid site in gas-phase chemistries and unexpectedly carboxylic groups are much less acidic than phenolic groups in the graphitized mesoporous carbon due to electron density delocalization induced by the aromatic rings of graphitic carbon. The methodology can identify acidic sites in oxygenated carbon materials in solid acid catalyst-driven chemistry.

Zhou, Jiahua↗
DOE OSTI2021

Mechanistic insights into carbon–carbon coupling on NiAu and PdAu single-atom alloys

Carbon–carbon coupling is an important step in many catalytic reactions, and performing sp 3 –sp 3 carbon–carbon coupling heterogeneously is particularly challenging. It has been reported that PdAu single-atom alloy (SAA) model catalytic surfaces are able to selectively couple methyl groups, producing ethane from methyl iodide. In this work, we extend this study to NiAu SAAs and find that Ni atoms in Au are active for C–I cleavage and selective sp 3 –sp 3 carbon–carbon coupling to produce ethane. Furthermore, we perform ab initio kinetic Monte Carlo simulations that include the effect of the iodine atom, which was previously considered a bystander species. We find that model NiAu surfaces exhibit a similar chemistry to PdAu, but the reason for the similarity is due to the role the iodine atoms play in terms of blocking the Ni atom active sites. Specifically, on NiAu SAAs, the iodine atoms outcompete the methyl groups for occupancy of the Ni sites leaving the Me groups on Au, while on PdAu SAAs, the binding strengths of methyl groups and iodine atoms at the Pd atom active site are more similar. These simulations shed light on the mechanism of this important sp 3 –sp 3 carbon–carbon coupling chemistry on SAAs. Furthermore, we discuss the effect of the iodine atoms on the reaction energetics and make an analogy between the effect of iodine as an active site blocker on this model heterogeneous catalyst and homogeneous catalysts in which ligands must detach in order for the active site to be accessed by the reactants.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Total ecosystem carbon stocks at the marine-terrestrial interface: Blue carbon of the Pacific Northwest Coast, United States

The coastal ecosystems of temperate North America provide a variety of ecosystem services including high rates of carbon sequestration. Yet, little data exist for the carbon stocks of major tidal wetland types in the Pacific Northwest, United States. We quantified the total ecosystem carbon stocks (TECS) in seagrass, emergent marshes, and forested tidal wetlands, occurring along increasing elevation and decreasing salinity gradients. The TECS included the total aboveground carbon stocks and the entire soil profile (to as deep as 3 m). TECS significantly increased along the elevation and salinity gradients: 217 ± 60 Mg C/ha for seagrass (low elevation/high salinity), 417 ± 70 Mg C/ha for low marsh, 551 ± 47 Mg C/ha for high marsh, and 1,064 ± 38 Mg C/ha for tidal forest (high elevation/low salinity). Soil carbon stocks accounted for >98% of TECS in the seagrass and marsh communities and 78% in the tidal forest. Soils in the 0–100 cm portion of the profile accounted for only 48%–53% of the TECS in seagrasses and marshes and 34% of the TECS in tidal forests. Thus, the commonly applied limit defining TECS to a 100 cm depth would greatly underestimate both carbon stocks and potential greenhouse gas emissions from land-use conversion. The large carbon stocks coupled with other ecosystem services suggest value in the conservation and restoration of temperate zone tidal wetlands through climate change mitigation strategies. However, the findings suggest that long-term sea level rise effects such as tidal inundation and increased porewater salinity will likely decrease ecosystem carbon stocks in the absence of upslope wetland migration buffer zones.

Kauffman, J Boone↗
DOE OSTI2021

Global patterns of forest autotrophic carbon fluxes

Carbon (C) fixation, allocation, and metabolism by trees set the basis for energy and material flows in forest ecosystems and define their interactions with Earth’s changing climate. However, we lack a cohesive synthesis on how forest carbon fluxes vary globally with respect to climate and one another. Here, we draw upon 1,319 records from the Global Forest Carbon Database (ForC), representing all major forest types and the nine most significant autotrophic carbon fluxes, to comprehensively explore how C cycling in mature, undisturbed forests varies with latitude and climate on a global scale. We show that, across all flux variables analyzed, C cycling decreases continuously with absolute latitude – a finding that confirms multiple previous studies but contradicts the idea that net primary productivity of temperate forests rivals that of tropical forests. C flux variables generally displayed similar trends across latitude and multiple climate variables, with no differences in allocation detected at this global scale. Temperature variables in general, and mean annual temperature and temperature seasonality in particular, were the best univariate predictors of C flux, explaining 19 - 71% of variation in the C fluxes analyzed. The effects of temperature were modified by moisture availability, with C flux reduced under hot and dry conditions and sometimes under very high precipitation. C fluxes increased with growing season length, but this was never the best univariate predictor. Within the growing season, the influence of climate on C cycling was small but significant for a number of flux variables. These findings clarify how forest C flux varies with latitude and climate on a global scale. In a period of accelerating climatic change, this improved understanding of the fundamental climatic controls on forest C cycling sets a foundation for understanding patterns of change.

Banbury Morgan, Rebecca↗
DOE OSTI2026

Data for Greenhouse Gas Accounting Procedures in Low Carbon Fuel Policies Overlook the Spatial Variability of Miscanthus-Derived Sustainable Aviation Fuel

Low carbon fuel policies such as the U.S. Renewable Fuel Standard (RFS), Canada Clean Fuel Regulations (CFR), and California Low Carbon Fuel Standard (LCFS) as well as the 45Z tax credit are intended to reduce greenhouse gas (GHG) emissions from transportation. Cellulosic feedstocks, optimized biorefineries, and favorable farming locations can significantly reduce biofuel carbon intensity (CI). Despite advances in field-to-fuel GHG monitoring and flexibility in resource allocation within biorefineries (e.g., governing net electricity production), rigid CI accounting procedures in current policies may limit CI responsiveness across candidate sites and processing facilities. This work examines a hypothetical biomass-to-sustainable aviation fuel (SAF) pathway using miscanthus and alcohol-to-jet (i) to demonstrate how GHG accounting requirements drive estimates of biofuel CIs and (ii) to explore potential CI and financial implications of scenario-specific life cycle assessment (LCA). Results demonstrate that GHG accounting using the CFR/LCFS can reasonably account for distinct levels of net electricity production by a biorefinery, but only the CFR yields similar CI sensitivity to spatially explicit factors (feedstock CI, grid electricity CI) as scenario-specific LCA: most GHG accounting frameworks do not capture CI variation across candidate sites in the United States. Ultimately, this work demonstrates the importance of LCA methodological specifications in low carbon fuel policies and tax credits.

Miscanthus↗
DOE OSTI2022

Carbon Dioxide Utilization Life Cycle Analysis Guidance for the U.S. DOE Office of Fossil Energy and Carbon Management (Version 2.0)

Capturing carbon dioxide (CO 2 ) and placing it in permanent storage in geologic formations is an option for reducing CO 2 emissions, but it may not be a viable one for all CO 2 emitters. For some, the added cost of capture may be too high to implement, or the geology near the source may not be suitable for storage. In these circumstances, other options will be needed. Carbon use and reuse, or CO 2 utilization (CO2U), is an alternative approach that seeks beneficial uses for captured CO 2 , such as using it as a feedstock in the production of fuels, chemicals, and building materials. These uses would give CO 2 value that could be used by suppliers (emitters) to offset capture costs. One of the principal features and challenges associated with CO2U is that the products derived from CO 2 must have lower carbon footprints than their conventional counterparts. Previous assessments of CO2U alternatives have focused on the carbon content of utilization products as an indicator of CO 2 equivalent (CO 2 e) emissions reduction potential. However, embodied emissions are—at best—only weakly correlated with the amount of carbon contained in any physical product. Therefore, the most attractive CO2U options will both displace the carbon in an existing product and improve the overall carbon efficiency of the manufacturing process. Research to overcome barriers will include identifying existing co-feeds and available low-carbon energy sources to enable the conversion of CO 2 to value-added products under favorable processing conditions. New discoveries in the fields of nano- and bio-technology will be applied to efficiently utilize CO 2 in new applications. Development of advanced materials and processes, integrating CO 2 capture with utilization processes (e.g., algae), exploring a diverse slate of products from CO 2 to effectively offset capture costs and developing processes based on waste energy are means to overcome these barriers. The research will lead to the development of advanced catalysts, materials, and equipment that can be used to convert CO 2 into useful products. The result will be multiple flexible and adaptable technology platforms that can be used to produce suites of products spanning multiple utilization pathways.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Autonomous measurements of carbon dioxide uptake in a blue carbon habitat

A better understanding of carbon cycling dynamics in coastal ecosystems is critical to quantify “blue carbon” storage. However, obtaining the necessary measurements is complicated due to the dynamic conditions that characterize these ecosystems, and requires both spatially and temporally resolved measurements of carbon dioxide concentrations and complementary water quality parameters. In this study, we evaluate the use of a utility-class autonomous surface vehicle as a platform for measuring carbon dioxide dynamics in a shallow, intertidal, blue carbon ecosystem. We integrated a carbon dioxide sensor, a water quality sonde, and hydroacoustic sensors with an autonomous surface vehicle and performed measurements in and around a meadow of Zostera marina eelgrass. Results demonstrate the utility of this system to provide insight into complex carbon cycling dynamics in and around a shallow, intertidal blue carbon habitat.

54 ENVIRONMENTAL SCIENCES↗