Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regional climate modeling”

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 451 records · Page 25

Relative Influence of Initial Surface and Atmospheric Conditions on Seasonal Water and Energy Balances

We constructed and analyzed wet and dry soil moisture composites for the mid-latitude GCIP region of the central US using long climate model simulations made with the NCAR CCM3 and reanalysis products from NCEP. Using the diagnostic composites as a guide, we have completed a series of predictability experiments in which we imposed soil water initial conditions in CCM3 for the GCIP region for June 1 from anomalously wet and dry years, with atmospheric initial conditions taken from June 1 of a year with 'near-normal' soil water, and initial soil water from the near-normal year and atmospheric initial conditions from the wet and dry years. Preliminary results indicate that the initial state of the atmosphere is more important than the initial state of soil water determining the subsequent late spring and summer evolution of sod water over the GCIP region. Surprisingly, neither the composites or the predictability experiments yielded a strong influence of soil moisture on the atmosphere. To explore this further, we have made runs with extreme dry soil moisture initial anomalies imposed over the GCIP region (the soil close to being completely dry). These runs did yield a very strong effect on the atmosphere that persisted for at least three months. We conclude that the magnitude of the initial soil moisture anomaly is crucial, at least in CCM3, and are currently investigating whether a threshold exists, below which little impact is seen. In a complementary study, we compared the impact of the initial condition of snow cover versus the initial atmospheric state over the western US (corresponding to the westward extension of the GAPP program follow-on to GCIP). In this case, the initial prescription of snow cover is far more important than the initial atmospheric state in determining the subsequent evolution of snow cover. We are currently working to understand the very different soil water and snow cover results.

Oglesby, Robert J.↗

Glacial meltwater cooling of the Gulf of Mexico - GCM implications for Holocene and present-day climates

The NCAR Community Climate Model GCM is presently used to investigate the possible effects on regional and hemispheric climates of reduced SSTs in the Gulf of Mexico, in view of delta-O-18 records and terrestrial evidence for at least two major glacial meltwater discharges after the last glacial maximum. Three numerical experiments have been conducted with imposed gulfwide SST coolings of 3, 6, and 12 C; in all cases, significant reductions arise in the North Atlantic storm-track intensity, together with a strong decrease in transient eddy water vapor transport out of the Gulf of Mexico. Other statistically significant changes occur across the Northern Hemisphere.

Oglesby, Robert J.↗

Using CESM-RESFire to understand climate–fire–ecosystem interactions and the implications for decadal climate variability

Large wildfires exert strong disturbance on regional and global climate systems and ecosystems by perturbing radiative forcing as well as the carbon and water balance between the atmosphere and land surface, while short- and long-term variations in fire weather, terrestrial ecosystems, and human activity modulate fire intensity and reshape fire regimes. The complex climate–fire–ecosystem interactions were not fully integrated in previous climate model studies, and the resulting effects on the projections of future climate change are not well understood. Here we use the fully interactive REgion-Specific ecosystem feedback Fire model (RESFire) that was developed in the Community Earth System Model (CESM) to investigate these interactions and their impacts on climate systems and fire activity. We designed two sets of decadal simulations using CESM-RESFire for present-day (2001–2010) and future (2051–2060) scenarios, respectively, and conducted a series of sensitivity experiments to assess the effects of individual feedback pathways among climate, fire, and ecosystems. Our implementation of RESFire, which includes online land–atmosphere coupling of fire emissions and fire-induced land cover change (LCC), reproduces the observed aerosol optical depth (AOD) from space-based Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products and ground-based AErosol RObotic NETwork (AERONET) data; it agrees well with carbon budget benchmarks from previous studies. We estimate the global averaged net radiative effect of both fire aerosols and fire-induced LCC at –0.59±0.52 W m –2 , which is dominated by fire aerosol–cloud interactions (–0.82±0.19 W m –2 ), in the present-day scenario under climatological conditions of the 2000s. The fire-related net cooling effect increases by ~170 % to –1.60±0.27 W m –2 in the 2050s under the conditions of the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Such considerably enhanced radiative effect is attributed to the largely increased global burned area (+19 %) and fire carbon emissions (+100 %) from the 2000s to the 2050s driven by climate change. The net ecosystem exchange (NEE) of carbon between the land and atmosphere components in the simulations increases by 33 % accordingly, implying that biomass burning is an increasing carbon source at short-term timescales in the future. High-latitude regions with prevalent peatlands would be more vulnerable to increased fire threats due to climate change, and the increase in fire aerosols could counter the projected decrease in anthropogenic aerosols due to air pollution control policies in many regions. We also evaluate two distinct feedback mechanisms that are associated with fire aerosols and fire-induced LCC, respectively. On a global scale, the first mechanism imposes positive feedbacks to fire activity through enhanced droughts with suppressed precipitation by fire aerosol–cloud interactions, while the second one manifests as negative feedbacks due to reduced fuel loads by fire consumption and post-fire tree mortality and recovery processes. These two feedback pathways with opposite effects compete at regional to global scales and increase the complexity of climate–fire–ecosystem interactions and their climatic impacts.

54 ENVIRONMENTAL SCIENCES↗

Influences of North Pacific Ocean Domain Extent on the Western U.S. Winter Hydroclimatology in Variable-Resolution CESM

Variable-resolution global climate models (VRGCMs) are a dynamical downscaling method that can reach spatiotemporal scales needed for regional climate assessments. Over the years, several users of VRGCMs have assumed where the location and extent of the refinement domain should be based on knowledge of the prevailing storm tracks and resolution dependence of important regional climate processes (e.g., atmospheric rivers [ARs] and orographic uplift), but the effect of high-resolution domain size and extent on the simulation of downstream hydroclimatic phenomena has not been systematically evaluated. Here, we use variable resolution in the Community Earth System Model (VR-CESM) to perform such a test. To do this, three VR-CESM grids were generated that span the entire, two thirds, and one third of the North Pacific and evaluated for a 30-year climatology using Atmospheric Model Intercomparison Project protocols. Simulations are compared with reanalysis products offshore (fifth-generation of the European Centre for Medium-Range Weather Forecasts atmospheric reanalysis [ERA5]) and onshore (Livneh, 2015, https://doi.org/10.1038/sdata.2015.42, and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) of the western United States. The westward expansion of refinement domain influenced integrated vapor transport (IVT), which was generally high biased but minimally impacted AR characteristics. Due to slight differences in landfalling AR counts in the western United States, California winter precipitation generally improved with westward expansion of the refinement domain. Western U.S. mountain snowpack and surface temperatures were insensitive to refinement domain size and were more influenced by changes in topographic resolution and/or land surface model version. Here, given minimal dependence of simulated western U.S. hydroclimate on refinement domain size over the North Pacific, we advise future VR-CESM studies to focus grid resolution on better resolving land surface heterogeneity.

54 ENVIRONMENTAL SCIENCES↗

Warming induces unexpectedly high soil respiration in a wet tropical forest

Abstract Tropical forests are a dominant regulator of the global carbon cycle, exchanging more carbon dioxide with the atmosphere than any other terrestrial biome. Climate models predict unprecedented climatic warming in tropical regions in the coming decades; however, in situ field warming studies are severely lacking in tropical forests. Here we present results from an in situ warming experiment in Puerto Rico, where soil respiration responses to + 4 o C warming were assessed half-hourly for a year. Soil respiration rates were 42-204% higher in warmed relative to ambient plots, representing some of the highest soil respiration rates reported for any terrestrial ecosystem. While respiration rates were significantly higher in the warmed plots, the temperature sensitivity (Q 10 ) was 71.7% lower, pointing to a mechanistic shift. Even with reduced Q 10 , if observed soil respiration rates persist in a warmer world, the feedback to future climate could be considerably greater than previously predicted or observed.

Science & Technology - Other Topics↗

AeDES: a next-generation monitoring and forecasting system for environmental suitability of Aedes-borne disease transmission

Aedes-borne diseases, such as dengue and chikungunya, are responsible for more than 50 million infections worldwide every year, with an overall increase of 30-fold in the last 50 years, mainly due to city population growth, more frequent travels and ecological changes. In the United States of America, the vast majority of Aedes-borne infections are imported from endemic regions by travelers, who can become new sources of mosquito infection upon their return home if the exposed population is susceptible to the disease, and if suitable environmental conditions for the mosquitoes and the virus are present. Since the susceptibility of the human population can be determined via periodic monitoring campaigns, the environmental suitability for the presence of mosquitoes and viruses becomes one of the most important pieces of information for decision makers in the health sector. We present a next-generation monitoring and forecasting system for Aedes-borne diseases’ environmental suitability (AeDES) of transmission in the conterminous United States and transboundary regions, using calibrated ento-epidemiological models, climate models and temperature observations. After analyzing the seasonal predictive skill of AeDES, we briefly consider the recent Zika epidemic, and the compound effects of the current Central American dengue outbreak happening during the SARS-CoV-2 pandemic, to illustrate how a combination of tailored deterministic and probabilistic forecasts can inform key prevention and control strategies .

54 ENVIRONMENTAL SCIENCES↗

Spatial pattern of lake evaporation increases under global warming linked to regional hydroclimate change

Abstract Lakes are critical natural resources that are vulnerable to climate change. In a warmer climate, lake evaporation is projected to increase globally, but with substantial variation between regions. Here, based on ensemble projections of climate and lake models and an attribution method, we show that future lake evaporation increase is strongly modulated by regional hydroclimate change. Specifically, a drying hydroclimate will amplify evaporation increase by enlarging surface vapor pressure deficit and reducing cloud shortwave reflection. Future lake evaporation increase is amplified in tropical America, the Mediterranean and Southeast China with drier future hydroclimates, and dampened in high latitudes and the Tibetan Plateau with wetter future hydroclimates. Such spatially coupled changes in lake evaporation and hydroclimate have important implications on regional lake water balance and volume change, which can aggravate water scarcity and flood risks.

54 ENVIRONMENTAL SCIENCES↗

Development of an Unbiased Future Solar Dataset for Solar Resource Adequacy Research Over CONUS

A high-resolution, long-term solar dataset is essential for capturing the variability of solar energy resources and informing strategies to ensure grid reliability and resilience in systems with high levels of solar energy integration. This study focuses on generating unbiased, high-resolution projections of solar irradiance through a statistical downscaling framework, using Earth system model (ESM) simulations obtained from the North American Coordinated Regional Climate Downscaling Experiment (NA-CORDEX). The National Solar Radiation Database (NSRDB) is used to calibrate statistical downscaling models. The newly developed dataset provides solar irradiance, surface air temperature, and surface wind speed at 4-km and hourly resolutions across the contiguous United States (CONUS), based on two future scenarios (RCP4.5 and RCP8.5). This study outlines key steps in developing the high-resolution future solar dataset, including (1) regridding ESM data to a common 20-km resolution grid, (2) correcting ESM biases using the NSRDB, and (3) applying temporal and spatial downscaling methods to generate high-resolution (4-km, hourly) solar projections. Preliminary results indicate that downscaled projections (4-km) captured reasonable spatial patterns when compared to observations across CONUS for four variables. On average across all pixels, 4-km daily-total GHI and DNI projections showed normalized bias (nBias) less than 1% and 6% for GHI and DNI against NSRDB, respectively (nBias less than 1% and 5% for daily-average surface air temperature and surface wind speed). In terms of long-term trend for GHI and DNI, there was no strong increasing or decreasing trend (when compared to surface air temperature), but it showed a very weak decreasing trend.

14 SOLAR ENERGY↗

Thermal remote sensing of surface soil water content with partial vegetation cover for incorporation into climate models

This study outlines a method for the estimation of regional patterns of surface moisture availability (M(sub 0)) and fractional vegetation (Fr) in the presence of spatially variable vegetation cover. The method requires relating variations in satellite-derived (NOAA, Advanced Very High Resolution Radiometer (AVHRR)) surface radiant temperature to a vegetation index (computed from satellite visible and near-infrared data) while coupling this association to an inverse modeling scheme. More than merely furnishing surface soil moisture values, the method constitues a new conceptual and practical approach for combining thermal infrared and vegetation index measurements for incorporating the derived values of M(sub 0) into hydrologic and atmospheric prediction models. Application of the technique is demonstrated for a region in and around the city of Newcastle upon Tyne situated in the northeast of England. A regional estimate of M(sub 0) is derived and is probabbly good for fractional vegetation cover up to 80% before errors in the estimated soil water content become unacceptably large. Moreover, a normalization scheme is suggested from which a nomogram, `universal triangle,' is constructed and is seen to fit the observed data well. The universal triangle also simplifies the inclusion of remotely derived M(sub 0) in hydrology and meteorological models and is perhaps a practicable step toward integrating derived data from satellite measurements in weather forecasting.

Gillies, Robert R.↗

BOREAS TGB-5 CO2, CH4 and CO Chamber Flux Data Over the NSA

The BOReal Ecosystem-Atmosphere Study Trace Gas Biogeochemistry (BOREAS TGB-5) team collected a variety of trace gas concentration and flux measurements at several NSA sites. This data set contains carbon dioxide (CO2), methane (CH4), and carbon monoxide (CO) chamber flux measurements conducted in 1994 at upland forest sites that experienced stand-replacement fires. These measurements were acquired to understand the impact of fires on soil biogeochemistry and related changes in trace gas exchange in boreal forest soils. Relevant ancillary data, including data concerning the soil temperature, solar irradiance, and information from nearby un-burned control sites, are included to provide a basis for modeling the regional impacts of fire and climate changes on trace gas biogeochemistry. The data are provided in tabular ASCII files.

Burke, Roger↗

Daytime Cloud Property Retrievals Over the Arctic from Multispectral MODIS Data

Improving climate model predictions over Earth's polar regions requires a complete understanding of polar clouds properties. Passive satellite remote sensing techniques can be used to retrieve macro and microphysical properties of polar cloud systems. However, over the Arctic, there is minimal contrast between clouds and the background snow surface observed in satellite data, especially for visible wavelengths. This makes it difficult to identify clouds and retrieve their properties from space. Variable snow and ice cover, temperature inversions, and the predominance of mixed-phase clouds further complicate cloud property identification. For this study, the operational Clouds and the Earth s Radiant Energy System (CERES) cloud mask is first used to discriminate clouds from the background surface in Terra Moderate Resolution Imaging Spectroradiometer (MODIS) data. A solar-infrared infrared nearinfrared technique (SINT) first used by Platnick et al. (2001) is used here to retrieve cloud properties over snow and ice covered regions.

Spangenberg, Douglas A.↗

The Salience of Nonlinearities in the Boreal Winter Response to ENSO: North Pacific and North America

The prominence of nonlinearities in the response to El Niño as compared to La Niña, to moderate El Niño events as compared to extreme El Ninño events, and to different flavors of El Niño events, are analyzed using the NASA Goddard Earth Observing System Chemistry-Climate Model. In the Central North Pacific region where the sea level pressure response to El Niño-Southern Oscillation (ENSO) peaks, nonlinearities are relatively muted. In contrast, changes to the east of this region (i.e. the far-Northeastern Pacific) and to the north of this region (over Alaska) in response to different ENSO phases are more clearly nonlinear, and become statistically robust after more than 15 events are considered. The relative prominence of these nonlinearities is related to the zonal wavenumber of the tropical precipitation response. Associated with these nonlinearities over the far-Northeastern Pacific are nonlinearities in precipitation over Western United States and surface temperature over Northwest North America and Midwestern United States. In all regions at least 15 events of each type are necessary before nonlinearities can be identified as statistically significant at the 95% confidence level due to the presence of internal atmospheric variability. As there have only been a similar number of ENSO events to the total needed for significance since 1920, it is not surprising that it has been difficult to establish statistically significant nonlinearities using observational data.

Nonlinearities↗

Seasonal Glacier and Snow Loading in Svalbard Recovered from Geodetic Observations

We processed time series from seven Global Navigation Satellite System (GNSS) stations and one Very Long Baseline Interferometry (VLBI) station in Svalbard. The goal was to capture the seasonal vertical displacements caused by elastic response of variable mass load due to ice and snow accumulation. We found that estimates of the annual signal in different GNSS solutions disagree by more than 3 mm which makes geophysical interpretation of raw GNSS time series problematic. To overcome this problem, we have used an enhanced Common Mode (CM) filtering technique. The time series are differentiated by the time series from remote station BJOS with known mass loading signals removed a priori. Using this technique, we have achieved a substantial reduction of the differences between the GNSS solutions. We have computed mass loading time series from a regional Climatic Mass Balance (CMB) and snow model that provides the amount of water equivalent at a 1 km resolution with a time step of 7 days. We found that the entire vertical loading signal is present in data of two totally independent techniques at a statistically significant level of 95%. This allowed us to conclude that the remaining errors in vertical signal derived from the CMB model are less than 0.2 mm at that significance level. Refining the land water storage loading model with a CMB model resulted in a reduction of the annual amplitude from 2:1 mm to 1:1 mm in the CM filtered time series, while it had only a marginal impact on raw time series. This provides a strong evidence that CM filtering is essential for revealing local periodic signals when a millimetre level of accuracy is required.

H P Kierulf↗

The Impact of Decarbonization on Particulate Soiling of Solar Panels

Climate researchers have examined many impacts of climate change on energy supply and demand under various scenarios. However, the effect of changing particulate deposition onto solar panel surfaces on solar power production efficiency (i.e., soiling) has not been studied. We therefore characterize probabilistic outcomes across multiple climate models and scenarios. We find large current regional losses (up to 40% without manual cleaning, up to 20% with monthly cleaning and rain removal) in generation that grow slightly under a high-emission scenario, largely due to regional increases in windblown dust. In contrast, under a low-emissions scenario, potential production increases significantly (2–8% interquartile range with only rain removal) due to reduced soiling, especially in regions of Asia and Africa where anthropogenic aerosols are major contributors to soiling. Projected changes vary widely across models in many dusty areas outside of the Sahara and Arabia. Differences can also be large in regions dominated by anthropogenic aerosols, such as Nigeria, eastern China, and northern India, where the full range across modeled potential power production changes extends from −1 to +11% for the end of the century (without manual cleaning), underscoring the need to consider multiple climate models. With large increases in projected solar power deployment, the relatively small potential production increases reported here could nevertheless represent a large dividend in additional energy production. Hence, reductions in air pollution attributable to decarbonization could provide positive feedback under which a greater deployment of solar power (or other renewables) increases the production of solar power, facilitating the transition to renewable energy.

renewable energy↗

Surface Radiation Budget and Climate Classification

The surface radiation budget of a region is strongly tied to its climate. An 8-yr climatology of surface radiation budget components for 2.5 deg over the earth is examined in order to learn how the regional climate and surface radiation are related. The yearly cycles of a few individual regions were studied by plotting monthly mean net longwave flux as a function of net shortwave flux at the surface. These plots show trajectories that are characteristic of the climate class. The behavior of the trajectories of surface radiation and their relation to the regional climate can be understood with simple conceptual models for many cases. From an examination of these trajectories, a set of parameters is developed, such as mean net longwave flux and range of net shortwave flux, which distinguish various climate classes on the basis of the surface radiation. These criteria are applied to produce a map of regional climate classes based on surface radiation, similar to those of Koeppen or Trewartha and Horn, which were based on vegetation, temperature, and precipitation. The current maps can be used to explore the relationships between surface radiation and regional climate.

Smith, G. Louis↗

IN11B-1621: Quantifying How Climate Affects Vegetation in the Amazon Rainforest

Amazon droughts in 2005 and 2010 have raised serious concern about the future of the rainforest. Amazon forests are crucial because of their role as the largest carbon sink in the world which would effect the global warming phenomena with decreased photosynthesis activity. Especially, after a decline in plant growth in 1.68 million km2 forest area during the once-in-a-century severe drought in 2010, it is of primary importance to understand the relationship between different climatic variables and vegetation. In an earlier study, we have shown that non-linear models are better at capturing the relation dynamics of vegetation and climate variables such as temperature and precipitation, compared to linear models. In this research, we learn precise models between vegetation and climatic variables (temperature, precipitation) for normal conditions in the Amazon region using genetic programming based symbolic regression. This is done by removing high elevation and drought affected areas and also considering the slope of the region as one of the important factors while building the model. The model learned reveals new and interesting ways historical and current climate variables affect the vegetation at any location. MAIAC data has been used as a vegetation surrogate in our study. For temperature and precipitation, we have used TRMM and MODIS Land Surface Temperature data sets while learning the non-linear regression model. However, to generalize the model to make it independent of the data source, we perform transfer learning where we regress a regularized least squares to learn the parameters of the non-linear model using other data sources such as the precipitation and temperature from the Climatic Research Center (CRU). This new model is very similar in structure and performance compared to the original learned model and verifies the same claims about the nature of dependency between these climate variables and the vegetation in the Amazon region. As a result of this study, we are able to learn, for the very first time how exactly different climate factors influence vegetation at any location in the Amazon rainforests, independent of the specific sources from which the data has been obtained.

global warming↗

Precipitation recycling in the Amazon basin

Precipitation recycling is the contribution of evaporation within a region to precipitation in that same region. The recycling rate is a diagnostic measure of the potential for interactions between land surface hydrology and regional climate. In this paper we present a model for describing the seasonal and spatial variability of the recycling process. The precipitation recycling ratio, rho, is the basic variable in describing the recycling process. Rho is the fraction of precipitation at a certain location and time which is contributed by evaporation within the region under study. The recycling model is applied in studyiing the hydrologic cycle in the Amazon basin. It is estimated that about 25% of all the rain that falls in the Amazon basin is contributed by evaporation within the basin. This estimate is based on analysis of a data set supplied by the European Centre for Medium-range Weather Forecasts (ECMWF). The same analysis is repeated using a different data set from the Geophysical Fluid Dynamics Laboratory (GFDL). Based on this data set, the recycling ratio is estimated to be 35%. The seasonal variability of the recycling ratio is small compared with the yearly average. The new estimates of the recycling ratio are compared with results of previous studies, and the differences are explained.

Eltahir, E. A. B.↗

Modelled land use and land cover change emissions – a spatio-temporal comparison of different approaches

Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC/yr for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC/yr). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

land use and land cover change↗