Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth”

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 1,459 records · Page 81

Modeled and observed properties related to the direct aerosol radiative effect of biomass burning aerosol over the southeastern Atlantic

Biomass burning smoke is advected over the southeastern Atlantic Ocean between July and October of each year. This smoke plume overlies and mixes into a region of persistent low marine clouds. Model calculations of climate forcing by this plume vary significantly in both magnitude and sign. NASA EVS-2 (Earth Venture Suborbital-2) ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) had deployments for field campaigns off the west coast of Africa in 3 consecutive years (September 2016, August 2017, and October 2018) with the goal of better characterizing this plume as a function of the monthly evolution by measuring the parameters necessary to calculate the direct aerosol radiative effect. Here, this dataset and satellite retrievals of cloud properties are used to test the representation of the smoke plume and the underlying cloud layer in two regional models (WRF-CAM5 and CNRM-ALADIN) and two global models (GEOS and UM-UKCA). The focus is on the comparisons of those aerosol and cloud properties that are the primary determinants of the direct aerosol radiative effect and on the vertical distribution of the plume and its properties. The representativeness of the observations to monthly averages are tested for each field campaign, with the sampled mean aerosol light extinction generally found to be within 20 % of the monthly mean at plume altitudes. When compared to the observations, in all models, the simulated plume is too vertically diffuse and has smaller vertical gradients, and in two of the models (GEOS and UM-UKCA), the plume core is displaced lower than in the observations. Plume carbon monoxide, black carbon, and organic aerosol masses indicate underestimates in modeled plume concentrations, leading, in general, to underestimates in mid-visible aerosol extinction and optical depth. Biases in mid-visible single scatter albedo are both positive and negative across the models. Observed vertical gradients in single scatter albedo are not captured by the models, but the models do capture the coarse temporal evolution, correctly simulating higher values in October (2018) than in August (2017) and September (2016). Uncertainties in the measured absorption Ångstrom exponent were large but propagate into a negligible (<4 %) uncertainty in integrated solar absorption by the aerosol and, therefore, in the aerosol direct radiative effect. Model biases in cloud fraction, and, therefore, the scene albedo below the plume, vary significantly across the four models. The optical thickness of clouds is, on average, well simulated in the WRF-CAM5 and ALADIN models in the stratocumulus region and is underestimated in the GEOS model; UM-UKCA simulates cloud optical thickness that is significantly too high. Overall, the study demonstrates the utility of repeated, semi-random sampling across multiple years that can give insights into model biases and how these biases affect modeled climate forcing. The combined impact of these aerosol and cloud biases on the direct aerosol radiative effect (DARE) is estimated using a first-order approximation for a subset of five comparison grid boxes. A significant finding is that the observed grid box average aerosol and cloud properties yield a positive (warming) aerosol direct radiative effect for all five grid boxes, whereas DARE using the grid-box-averaged modeled properties ranges from much larger positive values to small, negative values. It is shown quantitatively how model biases can offset each other, so that model improvements that reduce biases in only one property (e.g., single scatter albedo but not cloud fraction) would lead to even greater biases in DARE. Across the models, biases in aerosol extinction and in cloud fraction and optical depth contribute the largest biases in DARE, with aerosol single scatter albedo also making a significant contribution.

54 ENVIRONMENTAL SCIENCES↗

An approach to sulfate geoengineering with surface emissions of carbonyl sulfide

Sulfate geoengineering (SG) methods based on lower stratospheric tropical injection of sulfur dioxide (SO 2 ) have been widely discussed in recent years, focusing on the direct and indirect effects they would have on the climate system. Here a potential alternative method is discussed, where sulfur emissions are located at the surface or in the troposphere in the form of carbonyl sulfide (COS) gas. There are two time-dependent chemistry–climate model experiments designed from the years 2021 to 2055, assuming a 40 Tg-Syr -1 artificial global flux of COS, which is geographically distributed following the present-day anthropogenic COS surface emissions (SG-COS-SRF) or a 6 Tg-Syr -1 injection of COS in the tropical upper troposphere (SG-COS-TTL). The budget of COS and sulfur species is discussed, as are the effects of both SG-COS strategies on the stratospheric sulfate aerosol optical depth (~Δτ=0.080 in the years 2046–2055), aerosol effective radius (0.46 µm), surface SO x deposition (+8.9 % for SG-COS-SRF; +3.3 % for SG-COS-TTL), and tropopause radiative forcing (RF; ~-1.5 W m –2 in all-sky conditions in both SG-COS experiments). Indirect effects on ozone, methane and stratospheric water vapour are also considered, along with the COS direct contribution. According to our model results, the resulting net RF is –1.3 W m –2 , for SG-COS-SRF, and –1.5 W m –2 , for SG-COS-TTL, and it is comparable to the corresponding RF of –1.7 W m –2 obtained with a sustained injection of 4 Tg-Syr -1 in the tropical lower stratosphere in the form of SO 2 (SG-SO2, which is able to produce a comparable increase of the sulfate aerosol optical depth). Significant changes in the stratospheric ozone response are found in both SG-COS experiments with respect to SG-SO2 (~5 DU versus +1.4 DU globally). According to the model results, the resulting ultraviolet B (UVB) perturbation at the surface accounts for –4.3 % as a global and annual average (versus –2.4 % in the SG-SO2 case), with a springtime Antarctic decrease of –2.7 % (versus a +5.8 % increase in the SG-SO2 experiment). Overall, we find that an increase in COS emissions may be feasible and produce a more latitudinally uniform forcing without the need for the deployment of stratospheric aircraft. However, our assumption that the rate of COS uptake by soils and plants does not vary with increasing COS concentrations will need to be investigated in future work, and more studies are needed on the prolonged exposure effects to higher COS values in humans and ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Shallow- and deep-convection characteristics in the greater Houston, Texas, area using cell tracking methodology

The convective lifecycle, from initiation to maturity and dissipation, is driven by a combination of kinematic, thermodynamic, microphysical, and radiative processes that are strongly coupled and variable in time and space. Weather radars have been traditionally used to provide various convective-cloud characteristics. Here, we analyzed climatological convective-cell radar characteristics to obtain and assess the diurnal cycles of three convective-cell types – shallow, modest deep, and vigorous deep convective cells – that formed in the greater Houston area, using the National Weather Service radar from Houston, Texas, and a multi-cell identification and tracking algorithm. The examined dataset spans 4 years (2018–2021) and covers the warm-season months (June to September) in those years. The analysis showed clear diurnal cycles in cell initiation (CI) consistent with the sea breeze circulation and showed diurnal and normalized lifetime relationships in cell evolution parameters (e.g., maximum reflectivity, echo-top height, Geostationary Operational Environmental Satellite-16 (GOES-16) channel 13 brightness temperature, and the height of maximum reflectivity). The cell evolution is well represented by relationships between (1) the height and value of the maximum radar reflectivity, (2) the minimum GOES-16 channel 13 brightness temperature and the maximum vertically integrated liquid, (3) the maximum reflectivity and columnar-average reflectivity, and (4) the echo-top ascent rate and cell lifetime. The relationships presented herein help to identify the cell lifecycle stages such as early shallow convection, vigorous vertical development, anvil development, and convective core dissipation. GOES-16 Aerosol Optical Depth values are also used as a proxy for cell initiation aerosol concentrations to investigate any potential relationships between initiation location and aerosol concentration. Overall, no significant relationships between initiation location and aerosol concentration were found for the three cell types investigated, but there are some minor differences in the pre-CI aerosol optical depth for vigorous deep convective cells.

54 ENVIRONMENTAL SCIENCES↗

Impact of host climate model on contrail cirrus effective radiative forcing estimates

Estimates of aviation effective radiative forcing (ERF) indicate that contrail cirrus is currently its largest contributor, although with a substantial associated uncertainty of ∼ 70 %. Here, we implement the contrail parameterisation developed for the Community Atmosphere Model (CAM) in the UK Met Office Unified Model (UM), allowing us to compare, for the first time, the impact of key features of the host climate model on contrail cirrus ERF. We find that differences in background humidity between the models result in the UM-simulated contrail fractions being 2 to 3 times larger than in CAM. Additionally, the models show contrasting responses in overall global cloud fraction, with contrails increasing the total cloud fraction in the UM and decreasing it in CAM. Differences in the complexity of the cloud microphysics schemes lead to significant differences in simulated changes to cloud ice water content due to aviation. After compensating for the unrealistically low contrail optical depth in the UM, we estimate the 2018 contrail cirrus ERF to be 40.8 mW m−2 in the UM, compared to 60.1 mW m−2 in CAM. These values highlight the substantial uncertainty in contrail cirrus ERF due to differences in microphysics and radiation schemes between the two models. We also find a factor-of-8 uncertainty in contrail cirrus ERF due to existing uncertainty in contrail cirrus optical depth. Future research should focus on better representing microphysical and radiative contrail characteristics in climate models and on improved observational constraints.

Zhang, Weiyu [University of Leeds (United Kingdom)↗

Impact on cloud properties of reduced-sulphur shipping fuel in the Eastern North Atlantic

The global reduction in shipping fuel sulphur that culminated in 2020 with an ∼ 80 % reduction has created a large-scale natural experiment on the role of aerosol-cloud interaction (ACI) in the climate system. We compare observations from the Atmospheric Radiation Measurement program's Eastern North Atlantic site (ARM-ENA; 39.1° N, 28.0° W) during two June to September periods: 2016–2018 (pre-2020) and 2021–2023 (post-2020). We find a significant (∼ 15 %) decrease in cloud condensation nuclei concentrations post-2020, which resulted in a decrease in cloud droplet number (N d ) and an increase in effective radius (r e ) of marine boundary layer clouds. However, cloud liquid water path (LWP) increased post-2020. The increase in LWP offset the increase in r e , resulting in insignificant changes to optical depth. MODIS and CERES data in the vicinity of ENA during these periods produce similar results also with negligible change in albedo and optical depth. Regional cloud occurrence declined in line with changes in the large-scale meteorology. Our results highlight the complex interplay of factors that modulate cloud feedbacks in the Eastern North Atlantic.

Mace, Gerald G. [Univ. of Utah, Salt Lake City, UT↗

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗

Shortwave Array Spectroradiometer-Hemispheric (SAS-He): design and evaluation

A novel ground-based radiometer, referred to as the Shortwave Array Spectroradiometer-Hemispheric (SAS-He), is introduced. This radiometer uses the shadow-band technique to report total irradiance and its direct and diffuse components frequently (every 30 s) with continuous spectral coverage (350–1700 nm) and moderate spectral (~2.5 nm ultraviolet–visible and ~6 nm shortwave-infrared) resolution. The SAS-He's performance is evaluated using integrated datasets collected over coastal regions during three field campaigns supported by the US Department of Energy's Atmospheric Radiation Measurement (ARM) program, namely the (1) Two-Column Aerosol Project (TCAP; Cape Cod, Massachusetts), (2) Tracking Aerosol Convection Interactions Experiment (TRACER; in and around Houston, Texas), and (3) Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE; La Jolla, California). We compare (i) aerosol optical depth (AOD) and total optical depth (TOD) derived from the direct irradiance, as well as (ii) the diffuse irradiance and direct-to-diffuse ratio (DDR) calculated from two components of the total irradiance. As part of the evaluation, both AOD and TOD derived from the SAS-He direct irradiance are compared to those provided by a collocated Cimel sunphotometer (CSPHOT) at five (380, 440, 500, 675, 870 nm) and two (1020, 1640 nm) wavelengths, respectively. Additionally, the SAS-He diffuse irradiance and DDR are contrasted with their counterparts offered by a collocated multifilter rotating shadowband radiometer (MFRSR) at six (415, 500, 615, 675, 870, 1625 nm) wavelengths. Overall, reasonable agreement is demonstrated between the compared products despite the challenging observational conditions associated with varying aerosol loadings and diverse types of aerosols and clouds. For example, the AOD- and TOD-related values of root mean square error remain within 0.021 at 380, 440, 500, 675, 870, 1020, and 1640 nm wavelengths during the three field campaigns.

47 OTHER INSTRUMENTATION↗

A large-eddy simulation exploration of the assumptions used in retrieving entrainment from a mixing diagram approach with ground-based remote sensors

Entrainment is a crucial component of the atmospheric boundary layer (BL) moisture and heat budget. While usually thought of as only entrainment flux, entrainment within the mixed layer budget equation is really composed of two terms: the flux of a property across the boundary separating the BL from the free troposphere and the change in the concentration of a property as the depth of the BL changes. In a recent study, Wakefield et al. (2023) used ground-based remote-sensing observations to estimate entrainment flux as the residual of a mixing diagram framework that was applied to the daytime convective boundary layer. This present work uses large-eddy simulation (LES) to examine how well this residual assumption for entrainment fluxes alone compares to the actual sum of those two entrainment terms derived from spatial averages of the LES output. We highlight the importance of the second entrainment term in closing the mixed layer budget and show that the residual assumption does not represent entrainment flux only but rather a total entrainment term when the boundary layer depth is changing.

Rosenberger, Tessa E. [Cleveland State Univ., Clev↗

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗

Organic carbon burial by river meandering partially offsets bank erosion carbon fluxes in a discontinuous permafrost floodplain

Abstract. Arctic river systems erode permafrost in their banks and mobilize particulate organic carbon (OC). Meandering rivers can entrain particulate OC from permafrost many meters below the depth of annual thaw, potentially enabling the production of greenhouse gases. However, the amount and fate of permafrost OC that is mobilized by river erosion is uncertain. To constrain OC fluxes due to riverbank erosion and deposition, we collected riverbank and floodplain sediment samples along the Koyukuk River, which meanders through discontinuous permafrost in the Yukon River watershed, Alaska, USA, with an average migration rate of 0.52 m yr−1. We measured sediment total OC (TOC) content, radiocarbon activity, water content, bulk density, grain size, and floodplain stratigraphy. Radiocarbon activity and TOC content were higher in samples dominated by silt as compared to sand, which we used to map OC content onto floodplain stratigraphy and estimate carbon fluxes due to river meandering. Results showed that the Koyukuk River erodes and re-deposits a substantial flux of OC each year due to its depth and high migration rate, generating a combined OC flux of a similar magnitude to the floodplain net ecological productivity. However, sediment being eroded from cutbanks and deposited as point bars had similar OC stocks (mean ± 1 SD of 125.3±13.1 kg OC m−2 in cutbanks versus 114.0±15.7 kg OC m−2 in point bars) whether or not the banks contained permafrost. We also observed radiocarbon-depleted biospheric OC in both cutbanks and permafrost-free point bars. These results indicate that a substantial fraction of aged biospheric OC that is liberated from floodplains by bank erosion is subsequently re-deposited in point bars rather than being oxidized. The process of aging, erosion, and re-deposition of floodplain organic material may be intrinsic to river–floodplain dynamics, regardless of permafrost content.

54 ENVIRONMENTAL SCIENCES↗

Cosmogenic nuclide and solute flux data from central Cuban rivers emphasize the importance of both physical and chemical mass loss from tropical landscapes

We use 25 new measurements of in situ produced cosmogenic 26 Al and 10 Be in river sand, paired with estimates of dissolved load flux in river water, to characterize the processes and pace of landscape change in central Cuba. Long-term erosion rates inferred from 10 Be concentrations in quartz extracted from central Cuban river sand range from 3.4–189 Mg km –2 yr –1 (mean 59, median 45). Dissolved loads (10–176 Mg km –2 yr –1 ; mean 92, median 97), calculated from stream solute concentrations and modeled runoff, exceed measured cosmogenic- 10 Be-derived erosion rates in 18 of 23 basins. This disparity mandates that in this environment landscape-scale mass loss is not fully represented by the cosmogenic nuclide measurements. The 26 Al / 10 Be ratios are lower than expected for steadystate exposure or erosion in 16 of 24 samples. Depressed 26 Al/ 10 Be ratios occur in many of the basins that have the greatest disparity between dissolved loads (high) and erosion rates inferred from cosmogenic nuclide concentrations (low). Depressed 26 Al/ 10 Be ratios are consistent with the presence of a deep, mixed, regolith layer providing extended storage times on slopes and/or burial and extended storage during fluvial transport. River water chemical analyses indicate that many basins with lower 26 Al/ 10 Be ratios and high 10 Be concentrations are underlain at least in part by evaporitic rocks that rapidly dissolve. Our data show that when assessing mass loss in humid tropical landscapes, accounting for the contribution of rock dissolution at depth is particularly important. In such warm, wet climates, mineral dissolution can occur many meters below the surface, beyond the penetration depth of most cosmic rays and thus the production of most cosmogenic nuclides. Our data suggest the importance of estimating solute fluxes and measuring paired cosmogenic nuclides to better understand the processes and rates of mass transfer at a basin scale.

54 ENVIRONMENTAL SCIENCES↗

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗

The DOE E3SM version 2.1: overview and assessment of the impacts of parameterized ocean submesoscales

The U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2.1 builds on E3SMv2 with several changes, with the most notable being the addition of the Fox-Kemper et al. (2011) mixed-layer eddy parameterization. This parameterization captures the effect of finite-amplitude, mixed-layer eddies as an overturning streamfunction and has the primary function of restratification. Herein, we outline the changes to the mean climate state of E3SM that were introduced by the addition of this parameterization. Overall, the presence of the submesoscale parameterization improves the fidelity of the v2.1 simulation by reducing the ocean surface biases in the North Atlantic present in v2, as illustrated by changes in the climatological sea surface temperature and salinity and the Arctic sea-ice extent. Other impacts include a slight shoaling of the mixed-layer depths in the North Atlantic and a small improvement in the Atlantic Meridional Overturning Circulation (AMOC). We note that the expected shoaling due to the parameterization is regionally dependent in our coupled configuration. In addition, we investigate why the parameterization and its impacts on mixed-layer depth have little impact on the simulated AMOC: despite increased dense-water formation in the Norwegian Sea, only a small fraction of the water formed makes its way south into the North Atlantic basin. Version 2.1 also exhibits small improvements in the atmospheric climatology, with smaller biases in many notable quantities and modes of variability.

54 ENVIRONMENTAL SCIENCES↗