Understanding and reducing the uncertainties of land surface energy flux partitioning within CMIP6 land models
Not Available
54 ENVIRONMENTAL SCIENCES↗
SEARCH · Engineering Papers
Results for “land surface model”
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.
Changes in the Regional Water Cycle and Their Impact on Societies
ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.
54 ENVIRONMENTAL SCIENCES↗
Implementing belowground controls on nutrient uptake in ELMv2-SPRUCE improves representation of a boreal peatland ecosystem
Boreal peatlands store 13 %–32 % of the global soil carbon (C) stock, a service dependent on plant-mycorrhizal fungi associations. In these nutrient poor systems, ectomycorrhizal and ericoid mycorrhizal fungi supply up to >80 % of the nutrient requirements of their plant hosts, partly with mined nitrogen (N) and phosphorus (P) from soil organic matter that are otherwise inaccessible to plants. Despite the ecological significance, mycorrhizal associations are only represented in a few land surface or ecosystem models. We modify the peatland branch of version 2 of the Energy Exascale Earth System Land Model (ELMv2-SPRUCE) to replace the default photosynthesis-driven inorganic N and P (NP) uptake process with a more realistic representation of the process via three pathways: (1) direct inorganic NP uptake by uncolonized fine roots, (2) indirect inorganic NP acquisition and (3) indirect NP acquisition from organic sources by mycorrhizal roots. We systematically evaluated the performance of the default and modified models with field observations from a whole ecosystem warming and carbon dioxide fertilization experimental site: Spruce and Peatland Responses Under Changing Environment (SPRUCE), in northern Minnesota, USA. The modified model reduces the underestimation of the growth response of shrubs in the default model to warming from 40 %–80 % to 17 %–35 % and reduces the overall relative absolute error on C fluxes from 1.61 to 1.54 in calibration. Improvements on modeled shrub growths and shrub-moss community net ecosystem exchanges are also seen in validation. The improved growth response of shrubs to warming is accompanied by several-fold increase in direct inorganic NP uptake and decrease in fungal colonization rate. The modified model simulates a smaller magnitude of transition of the ecosystem from C sink to C source under warming due to alleviation of plant nutrient limitation. Equifinality analysis shows the newly added parameters in the modified model can be constrained by the observed C fluxes. Sensitivity analysis shows the newly added parameters have stronger statistical interactions than the preexisting parameters in the default model. Overall, the modified model is an improvement over the default ELMv2-SPRUCE and will be a useful tool for understanding boreal peatland change.
Wang, Yaoping [Oak Ridge National Laboratory (ORNL↗
The Terrestrial Biosphere Model Farm
Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901–2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.
54 ENVIRONMENTAL SCIENCES↗
Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw: Modeling Archive
This Modeling Archive is in support of a NGEE Arctic publication in review "Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw". The study used a cryohydrology model to assess the potential risk of abrupt permafrost thaw triggered by melting ground ice, a key open question associated with permafrost response to a warming Arctic. The spatially resolved simulations are for a small catchment 465 ice-wedge polygons in polygonal tundra near Utqiagvik, Alaska in the high-emissions RCP8.5 climate scenario. The simulations are compared to runoff, evapotranspiration and subsidence in the current climate and agree well. The study used the ATS code configured as an intermediate-scale cryohydrology model (Advanced Terrestrial Simulator, Version v2). The archive includes input files for spinup (1985 to 2005) and projections (2006 to 2100 - the manuscript reports 2006 to 2098). Three of the projections include the effects of subsidence and microtopography change with different depth profiles of ice content corresponding to the median, 20th percentile and 80th percentile. The fourth projection has subsidence turnoff and uses the reference case (median) ice content. Spatially averaged or aggregated output variables are included for the projections. Selected checkpoint files for the projections and postprocessing scripts are also included. Files included are *.xml; *.h5; *.exo; *.py; *.sh, *. nb, mesh files, and one file in the original Excel format plus *.csv and *.pdf to conserve formatting contained in the original file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
54 ENVIRONMENTAL SCIENCES↗
Root Foraging Alters Global Patterns of Ecosystem Legacy From Climate Perturbations
The response of terrestrial ecosystems to climate perturbations typically persist longer than the timescale of the forcing, a phenomenon broadly referred to as legacy. Understanding the strength of legacy is critical for predicting ecosystem sensitivity to climate extremes and the extent that persistent changes in surface-atmosphere exchange might feedback onto climate. The cause of ecosystem legacy has been associated with myriad factors such as changes in aboveground biomass, however, few studies have tested how changes in the depth distribution of fine roots in response to perturbation might alter an ecosystem's recovery time. We explore this question using an Earth System Model that includes a dynamic root module where vegetation can forage for water and nutrients by altering their root profiles. The global simulations presented here show that in response to water stress events most ecosystems deepen their root profiles. In semi-arid ecosystems, this response expedites recovery (i.e., less legacy) relative to simulations without dynamics roots because access to deeper water pools after the initial event remains favorable. In wetter ecosystems, the development of deeper root profiles slows down the recovery timescale (i.e., more legacy) because the deeper root profile reduces access to nutrients and is thus unfavorable. Here, the simulations show that while the inclusion of dynamic roots might only minimally affect global patterns of Gross Primary Productivity and Evapotranspiration, the shift in root profile alters the timescale of recovery. Studies interested in the sustained response of land surfaces fluxes to climate disturbances should consider models that include dynamic root capability.
54 ENVIRONMENTAL SCIENCES↗
Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model
Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.
58 GEOSCIENCES↗
Global Methane Budget 2000–2020
Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. CH4 is the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2), and both emissions and atmospheric concentrations of CH4 have continued to increase since 2007 after a temporary pause. The relative importance of CH4 emissions compared to those of CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in quantifying the factors responsible for the observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise, and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in situ and Greenhouse Gases Observing SATellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full data sets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this 2025 edition benefits from important progress in estimating inland freshwater emissions, with better counting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double counting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double counting that may exist (average of 23 Tg CH4 yr−1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr−1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr−1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr−1 or ∼ 65 % is attributed to direct anthropogenic sources in the fossil, agriculture, and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr−1 or 63 %–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr−1 (range 9–40). The 2020 emission rate is the highest of the period and reaches 608 Tg CH4 yr−1 (range 581–627), which is 12 % higher than the average emissions in the 2000s. Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr−1) larger global emissions (669 Tg CH4 yr−1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr−1 in Saunois et al. (2016, 2020) respectively), and for the first time uncertainties in bottom-up and top-down budgets overlap. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. The tropospheric loss of methane, as the main contributor to methane lifetime, has been estimated at 563 [510–663] Tg CH4 yr−1 based on chemistry–climate models. These values are slightly larger than for 2000–2009 due to the impact of the rise in atmospheric methane and remaining large uncertainty (∼ 25 %). The total sink of CH4 is estimated at 633 [507–796] Tg CH4 yr−1 by the bottom-up approaches and at 554 [550–567] Tg CH4 yr−1 by top-down approaches. However, most of the top-down models use the same OH distribution, which introduces less uncertainty to the global budget than is likely justified. For 2010–2019, agriculture and waste contributed an estimated 228 [213–242] Tg CH4 yr−1 in the top-down budget and 211 [195–231] Tg CH4 yr−1 in the bottom-up budget. Fossil fuel emissions contributed 115 [100–124] Tg CH4 yr−1 in the top-down budget and 120 [117–125] Tg CH4 yr−1 in the bottom-up budget. Biomass and biofuel burning contributed 27 [26–27] Tg CH4 yr−1 in the top-down budget and 28 [21–39] Tg CH4 yr−1 in the bottom-up budget. We identify five major priorities for improving the CH4 budget: (i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; (ii) further development of process-based models for inland-water emissions; (iii) intensification of CH4 observations at local (e.g. FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture, and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).
54 ENVIRONMENTAL SCIENCES↗
Influence of plateau, slope, and valley on soil hydrology during the dry season in a Central Amazon old‐growth forest
Soil moisture regulates plant water supply and drought sensitivity in tropical forests, yet its vertical and topographic variation remains poorly characterized. We combined high-frequency time-domain reflectometry measurements from 5 to 100 cm across plateau, slope, and valley landforms at the Zona Florestal 2 research site north of Manaus, Central Amazonia, to quantify how soil moisture memory, timing of responses to rainfall, dry-down rates (τ), and soil–water depletion vary across these contrasting landforms. Landform-specific soil moisture calibration curves ensured accurate volumetric water content estimates in these highly weathered soils. During the 2023 dry-to-wet transition (August–November), soil moisture memory showed strong topographic contrasts, with valley profiles increasing from ∼47 h at 5 cm to ∼154 h at 100 cm, while plateaus exhibited higher near-surface persistence (∼124 h at 5 cm) but weaker memory at depth. Dry-down behavior reinforced these differences as valley soils exhibited τ values exceeding ∼200 h, more than double the characteristic τ of plateau soils (∼90 h). Rainfall–soil moisture correlations indicated immediate responses at shallow depths in valleys and progressively longer lags with depth on plateaus and slopes. These hydrologic patterns were mirrored in depletion profiles, which declined sharply below 30 cm on plateaus but remained high and sustained throughout the upper meter in slopes and valleys. Together, these findings provide the first depth-resolved field measurements of soil moisture memory, rainfall coupling, dry-down constants, and depletion dynamics across major upland landforms in Central Amazonia and offer clear observational benchmarks for improving land-surface and ecosystem model representations of soil–water processes.
Hillslope↗
The Global Methane Budget 2000–2017
Understanding and quantifying the global methane (CH 4 ) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH 4 are continuing to increase, making CH 4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO 2 ). Assessing the relative importance of CH 4 in comparison to CO 2 is complicated by its shorter atmospheric lifetime, stronger warming potential, and atmospheric growth rate variations over the past decade, the causes of which are still debated. Two major difficulties in reducing uncertainties arise from the variety of geographically overlapping CH 4 sources and from the destruction of CH 4 by short-lived hydroxyl radicals (OH). To address these difficulties, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations).
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Transformation rate maps of dissolved organic carbon in the contiguous US
Riverine dissolved organic carbon (DOC) plays a vital role in regional and global carbon cycles. However, the processes of DOC conversion from soil organic carbon (SOC) and leaching into rivers are insufficiently understood, inconsistently represented, and poorly parameterized, particularly in land surface and Earth system models. As a first attempt to fill this gap, we propose a generic formula that directly connects SOC concentration with DOC concentration in headwater streams, where a single parameter, the transformation rate from SOC in the soil to DOC leaching flux (P r ), accounts for the overall processes governing SOC conversion to DOC and leaching from soils (along with runoff) into headwater streams. We then derive high-resolution P r maps over the contiguous US (CONUS) using SOC data from two different sources: the Harmonized World Soil Database v1.2 (HWSD) and SoilGrids 2.0. Both maps are developed following the same five major steps: (1) selecting independent catchments where observed riverine DOC data are available with reasonable quality; (2) estimating catchment-average SOC for the independent catchments; (3) estimating the P r values for these catchments based on the generic formula and catchment-average SOC; (4) developing a predictive model of P r with machine learning (ML) techniques and catchment-scale climate, hydrology, geology, and other attributes; and (5) deriving a national map of P r based on the ML model. For evaluation, we compare the DOC concentration derived using the P r map and the observed DOC concentration values at evaluation catchments. The resulting mean absolute scaled error and coefficient of determination are 0.73 and 0.47 for the HWSD-based model and 0.58 and 0.72 for the SoilGrids-based model, respectively, suggesting the effectiveness of the overall methodology. Efforts to constrain uncertainty and evaluate sensitivity of P r to different factors are discussed. To illustrate the use of such maps, we derive a riverine DOC concentration reanalysis dataset over CONUS. The two P r maps, robustly derived and empirically validated, lay a critical cornerstone for better simulating the terrestrial carbon cycle in land surface and Earth system models. Our findings not only set a foundation for improving our predictive understanding of the terrestrial carbon cycle at the regional and global scales, but also hold promises for informing policy decisions related to decarbonization and climate change mitigation. The data presented in this study are publicly available at https://doi.org/10.5281/zenodo.14563816 (Li et al., 2024).
54 ENVIRONMENTAL SCIENCES↗
Deforestation reshapes land-surface energy-flux partitioning
Land-use and land-cover change significantly modify local land-surface characteristics and water/energy exchanges, which can lead to atmospheric circulation and regional climate changes. In particular, deforestation accounts for a large portion of global land-use changes, which transforms forests into other land cover types, such as croplands and grazing lands. Many previous efforts have focused on observing and modeling land–atmosphere–water/energy fluxes to investigate land–atmosphere coupling induced by deforestation. However, interpreting land–atmosphere–water/energy-flux responses to deforestation is often complicated by the concurrent impacts from shifts in land-surface properties versus background atmospheric forcings. In this study, we used 29 paired FLUXNET sites, to improve understanding of how deforested land surfaces drive changes in surface-energy-flux partitioning. Each paired sites included an intact forested and non-forested site that had similar background climate. We employed transfer entropy, a method based on information theory, to diagnose directional controls between coupling variables, and identify nonlinear cause–effect relationships. Transfer entropy is a powerful tool to detective causal relationships in nonlinear and asynchronous systems. The paired eddy covariance flux measurements showed consistent and strong information flows from vegetation activity (gross primary productivity (GPP)) and physical climate (e.g. shortwave radiation, air temperature) to evaporative fraction (EF) over both non-forested and forested land surfaces. More importantly, the information transfers from radiation, precipitation, and GPP to EF were significantly reduced at non-forested sites, compared to forested sites. We then applied these observationally constrained metrics as benchmarks to evaluate the Energy Exascale Earth System Model (E3SM) land model (ELM). ELM predicted vegetation controls on EF relatively well, but underpredicted climate factors on EF, indicating model deficiencies in describing the relationships between atmospheric state and surface fluxes. Moreover, changes in controls on surface energy flux partitioning due to deforestation were not detected in the model. We highlight the need for benchmarking model simulated surface-energy fluxes and the corresponding causal relationships against those of observations, to improve our understanding of model predictability on how deforestation reshapes land surface energy fluxes.
54 ENVIRONMENTAL SCIENCES↗
Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality
Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.
54 ENVIRONMENTAL SCIENCES↗
Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon
Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.
Wang, Chen↗
STAIR 2.0: A Generic and Automatic Algorithm to Fuse Modis, Landsat, and Sentinel-2 to Generate 10 m, Daily, and Cloud-/Gap-Free Surface Reflectance Product
Remote sensing datasets with both high spatial and high temporal resolution are critical for monitoring and modeling the dynamics of land surfaces. However, no current satellite sensor could simultaneously achieve both high spatial resolution and high revisiting frequency. Therefore, the integration of different sources of satellite data to produce a fusion product has become a popular solution to address this challenge. Many methods have been proposed to generate synthetic images with rich spatial details and high temporal frequency by combining two types of satellite datasets—usually frequent coarse-resolution images (e.g., MODIS) and sparse fine-resolution images (e.g., Landsat). In this paper, we introduce STAIR 2.0, a new fusion method that extends the previous STAIR fusion framework, to fuse three types of satellite datasets, including MODIS, Landsat, and Sentinel-2. In STAIR 2.0, input images are first processed to impute missing-value pixels that are due to clouds or sensor mechanical issues using a gap-filling algorithm. The multiple refined time series are then integrated stepwisely, from coarse- to fine- and high-resolution, ultimately providing a synthetic daily, high-resolution surface reflectance observations. We applied STAIR 2.0 to generate a 10-m, daily, cloud-/gap-free time series that covers the 2017 growing season of Saunders County, Nebraska. Moreover, the framework is generic and can be extended to integrate more types of satellite data sources, further improving the quality of the fusion product. View Full-Text
47 OTHER INSTRUMENTATION↗
The Modeled Seasonal Cycles of Surface N 2 O Fluxes and Atmospheric N 2 O
Nitrous oxide (N 2 O) is a greenhouse gas and stratospheric ozone-depleting substance with large and growing anthropogenic emissions. Previous studies identified the influx of N 2 O-depleted air from the stratosphere to partly cause the seasonality in tropospheric N 2 O (aN 2 O), but other contributions remain unclear. Here, we combine surface fluxes from eight land and four ocean models from phase 2 of the Nitrogen/N 2 O Model Intercomparison Project with tropospheric transport modeling to simulate aN 2 O at eight remote air sampling sites for modern and pre-industrial periods. Models show general agreement on the seasonal phasing of zonal-average N 2 O fluxes for most sites, but seasonal peak-to-peak amplitudes differ several-fold across models. The modeled seasonal amplitude of surface aN 2 O ranges from 0.25 to 0.80 ppb (interquartile ranges 21%–52% of median) for land, 0.14–0.25 ppb (17%–68%) for ocean, and 0.28–0.77 ppb (23%–52%) for combined flux contributions. The observed seasonal amplitude ranges from 0.34 to 1.08 ppb for these sites. The stratospheric contributions to aN 2 O, inferred by the difference between the surface-troposphere model and observations, show 16%–126% larger amplitudes and minima delayed by ~1 month compared to Northern Hemisphere site observations. Land fluxes and their seasonal amplitude have increased since the pre-industrial era and are projected to grow further under anthropogenic activities. Our results demonstrate the increasing importance of land fluxes for aN 2 O seasonality. Considering the large model spread, in situ aN 2 O observations and atmospheric transport-chemistry models will provide opportunities for constraining terrestrial and oceanic biosphere models, critical for projecting carbon-nitrogen cycles under ongoing global warming.
54 ENVIRONMENTAL SCIENCES↗
Improved subseasonal-to-seasonal precipitation prediction of climate models with nudging approach for better initialization of Tibetan Plateau-Rocky Mountain Circumglobal wave train and land surface conditions
Abstract Reliable subseasonal-to-seasonal (S2S) precipitation prediction is highly desired due to the great socioeconomical implications, yet it remains one of the most challenging topics in the weather/climate prediction research area. As part of the Impact of Initialized Land Temperature and Snowpack on Sub-seasonal to Seasonal Prediction (LS4P) project of the Global Energy and Water Exchanges (GEWEX) program, twenty-one climate models follow the LS4P protocol to quantify the impact of the Tibetan Plateau (TP) land surface temperature/subsurface temperature (LST/SUBT) springtime anomalies on the global summertime precipitation. We find that nudging towards reanalysis winds is crucial for climate models to generate atmosphere and land surface initial conditions close to observations, which is necessary for meaningful S2S applications. Simulations with nudged initial conditions can better capture the summer precipitation responses to the imposed TP LST/SUBT spring anomalies at hotspot regions all over the world. Further analyses show that the enhanced S2S prediction skill is largely attributable to the substantially improved initialization of the Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train pattern in the atmosphere. This study highlights the important role that initial condition plays in the S2S prediction and suggests that data assimilation technique (e.g., nudging) should be adopted to initialize climate models to improve their S2S prediction.
54 ENVIRONMENTAL SCIENCES↗
Upscaling Methane Flux From Plot Level to Eddy Covariance Tower Domains in Five Alaskan Tundra Ecosystems
Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2 s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2 s −1 and 0.758 using HF and 0.007 μmol m −2 s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems.
54 ENVIRONMENTAL SCIENCES↗