Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainable”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Combined Effects of Stream Hydrology and Land Use on Basin‐Scale Hyporheic Zone Denitrification in the Columbia River Basin

Abstract Denitrification in the hyporheic zone (HZ) of river corridors is crucial to removing excess nitrogen in rivers from anthropogenic activities. However, previous modeling studies of the effectiveness of river corridors in removing excess nitrogen via denitrification were often limited to the reach‐scale and low‐order stream watersheds. We developed a basin‐scale river corridor model for the Columbia River Basin with random forest models to identify the dominant factors associated with the spatial variation of HZ denitrification. Our modeling results suggest that the combined effects of hydrologic variability in reaches and substrate availability influenced by land use are associated with the spatial variability of modeled HZ denitrification at the basin scale. Hyporheic exchange flux can explain most of spatial variation of denitrification amounts in reaches of different sizes, while among the reaches affected by different land uses, the combination of hyporheic exchange flux and stream dissolved organic carbon (DOC) concentration can explain the denitrification differences. Also, we can generalize that the most influential watershed and channel variables controlling denitrification variation are channel morphology parameters (median grain size (D50), stream slope), climate (annual precipitation and evapotranspiration), and stream DOC‐related parameters (percent of shrub area). The modeling framework in our study can serve as a valuable tool to identify the limiting factors in removing excess nitrogen pollution in large river basins where direct measurement is often infeasible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Maximum Supersaturation in the Marine Boundary Layer Clouds Over the North Atlantic

Abstract The maximum supersaturation ( S x ) in clouds is a key parameter affecting the cloud's microphysical and radiative properties. We investigate the S x of the marine boundary layer clouds by combining airborne and surface observations in the Eastern North Atlantic. The cloud droplet number concentration ( N c ) in the least diluted cloud cores agrees well with the number concentration of particles larger than the Hoppel Minimum (HM) ( N >HM ) below clouds, indicating that the HM represents the average size threshold above which particles are activated to form cloud droplets. The S x values derived from surface observations vary from 0.10% to 0.50% from June 2017 to June 2018, with a clear seasonal variation exhibiting higher values during winter. Most of the S x variance (∼60%) can be explained by the cloud condensation nuclei (CCN) concentration and updraft velocity ( w ), with the CCN concentration playing a more important role than w in explaining the variation of S x . The influence of CCN concentration on S x leads to a buffered response of N c to aerosol perturbations. The response of N c to low aerosol concentration during winter is further buffered by the high w . The global Community Earth System Model (CESM) simulated S x values in the Azores have a positive bias compared to measured S x , likely due to overestimated w and underestimated CCN concentration. The CESM simulated S x exhibits higher values further north over the North Atlantic Ocean, which is attributed to stronger w . The suppression of S x by aerosol is also evident in regions with high CCN concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Synoptic Drivers of Atmospheric River Induced Precipitation Near Dronning Maud Land, Antarctica

Atmospheric rivers (ARs) that reach the Antarctic Ice Sheet (AIS) transport anomalous moisture from lower latitudes and can impact the AIS via extreme precipitation and increased downward longwave radiation. ARs contribute significantly to the interannual variability of precipitation over the AIS and thus are likely to play a key role in understanding future changes in the surface mass balance of the AIS. Dronning Maud Land (DML) is one of four maxima in AR frequency over coastal East Antarctica, with AR precipitation explaining 77% of the interannual variability in precipitation for this region. We employ a 16-node self-organizing map (SOM) trained with MERRA-2 sea-level pressure anomalies to identify synoptic-scale environments associated with landfalling ARs in and around DML. Node composites of atmospheric variables reveal common drivers of precipitation associated with ARs reaching DML including anomalous high-low surface pressure couplets, anomalously high integrated water vapor, and coastal barrier jets. Using a quasi-geostrophic framework, we find that upward vertical motion associated with the occlusion process of attendant cyclones dominates atmospheric lift in AR environments. We further identify mechanisms that explain the variability in AR precipitation intensity across nodes, such as the lift associated with the occlusion process of attendant cyclones and the spatial coincidence of ascent induced by the occlusion process and frontogenesis. Here, the latter suggests that ARs making landfall during the mature phase of cyclogenesis result in higher precipitation intensity compared to landfalling ARs that occur during the occluded phase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multiscale Temporal Variability of the Global Air‐Sea CO 2 Flux Anomaly

Abstract The global air‐sea CO 2 flux (F) impacts and is impacted by a plethora of climate‐related processes operating at multiple time scales. In bulk mass transfer formulations, F is driven by physico‐ and bio‐chemical factors such as the air‐sea partial pressure difference (∆pCO 2 ), gas transfer velocity, sea surface temperature, and salinity–all varying at multiple time scales. To de‐convolve the impact of these factors on variability in F at different time scales, time‐resolved estimates of F were computed using a global data set assembled between 1988 and 2015. The F anomalies were defined as temporal deviations from the 28‐year time‐averaged value. Spectral analysis revealed four dominant timescales of variability in F–subseasonal, seasonal, interannual, and decadal with relative amplitude differences varying across regions. A second‐order Taylor series expansion was then conducted along these four timescales to separate drivers across differing regions. The analysis showed that on subseasonal timescales, wind speed variability explains some 66% of the global F anomaly and is the dominant driver. On seasonal, interannual, and decadal timescales, the ∆pCO 2 effect controlled by the ∆pCO 2 anomaly, explained much of the F anomaly. On decadal timescales, the F anomaly was almost entirely governed by the ∆pCO 2 effect with large contributions from high latitudes. The main drivers across timescales also dominate the regional F anomaly, particularly in the mid‐high latitude regions. Finally, the driver of the ∆pCO 2 effect was closely connected with the relative strength of atmospheric pCO 2 and the nonthermal component of oceanic pCO 2 anomaly associated with dissolved inorganic carbon and alkalinity.

Environmental Sciences & Ecology↗
DOE OSTI2023

The Zonal Patterns in Late Quaternary Tropical South American Precipitation

Speleothem oxygen isotope records (δ 18 O) of tropical South American rainfall in the late Quaternary show a zonal “South American Precipitation Dipole” (SAPD). The dipole is characterized by opposing east-west precipitation anomalies compared to the present—wetter in the east and drier in the west at the mid-Holocene (~7 ka), and drier in the east and wetter in the west at the Last Glacial Maximum (~21 ka). However, the SAPD remains enigmatic because it is expressed differently in western versus eastern δ18O records and isotope-enabled climate model simulations usually misrepresent the magnitude and/or spatial pattern of δ 18 O change. Here, we address the SAPD enigma in two parts. First, we re-interpret the δ 18 O data to account for upwind rainout effects that are known to be pervasive in tropical South America, but are not always considered in Quaternary paleoclimate studies. Our revised interpretation reconciles the δ 18 O data with cave infiltration and other proxy records, and indicates that the centroid of tropical South American rainfall has migrated zonally over time. Second, using an energy balance model of tropical atmospheric circulation, we hypothesize that zonal migration of the precipitation centroid can be explained by regional energy budget shifts, such as changing Saharan albedo associated with the African Humid Period, that have not been modeled in previous SAPD studies. This hypothesis of a migrating precipitation centroid presents a new framework for interpreting δ18O records from tropical South America and may help explain the zonal rainfall anomalies that predate the late Quaternary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaporation and Transpiration From Multiple Proximal Forests and Wetlands

Abstract Climate change is intensifying the hydrologic cycle and altering ecosystem function, including water flux to the atmosphere through evapotranspiration (ET). ET is made up of evaporation (E) via non‐stomatal surfaces, and transpiration (T) through plant stomata which are impacted by global changes in different ways. E and T are difficult to measure independently at the ecosystem scale, especially across multiple sites that represent different land use and land management strategies. To address this gap in understanding, we applied flux variance similarity (FVS) to quantify how E and T differ across 13 different ecosystems measured using eddy covariance in a 10 × 10 km area from the CHEESEHEAD19 experiment in northern Wisconsin, USA. The study sites included eight forests with a large deciduous broadleaf component, three evergreen needleleaf forests, and two wetlands. Average T/ET for the study period averaged nearly 52% in forested sites and 45% in wetlands, with larger values after excluding periods following rain events when evaporation from canopy interception may be expected. A dominance analysis revealed that environmental variables explained on average 69% of the variance of half‐hourly T, which decreased from summer to autumn. Deciduous and evergreen forests showed similar E trajectories over time despite differences in vegetation phenology, and vapor pressure deficit explained some 13% of the variance E in wetlands but only 5% or less in forests. Retrieval of E and T within a dense network of flux towers lends confidence that FVS is a promising approach for comparing ecosystem hydrology across multiple sites to improve our process‐based understanding of ecosystem water fluxes.

Environmental Sciences & Ecology↗
DOE OSTI2023

Prolonged Drought in a Northern California Coastal Region Suppresses Wildfire Impacts on Hydrology

Wildfires naturally occur in many landscapes, however they are undergoing rapid regime shifts. Despite the emphasis in the literature on the most severe hydrological responses to wildfire, there remains a knowledge gap on the thresholds of wildfire (i.e., burned area/drainage area ratio, BAR) required to initiate hydrological responses. We investigated hydrological changes in the Russian River Watershed (RRW) in California, a coastal, Mediterranean, drought-prone, wildfire-adapted ecosystem, following ten wildfires that burned 30% of the watershed. Our findings suggest that sub-watersheds of the RRW have not burned beyond an intrinsic, unknown, threshold required to initiate change. Using paired watersheds, we examined spatiotemporal patterns of pre-and-post wildfire hydrology with a rainfall-runoff hydrological model. Even though these successive wildfires burned 1%–50% of each sub-watershed (1%–30% at moderate/high severity), we found little evidence of wildfire-related shifts in hydrology. As a function of BAR, wildfire imposed limited effects on runoff ratios (runoff/precipitation) and runoff residuals (observations—model simulations). Our findings that post-wildfire runoff enhancements asymptote beyond 30% burn indicate that when a watershed is burned beyond a certain threshold, the magnitude of the hydrologic response no longer increases. Drought and storm conditions explained much of the variability observed in streamflow, whereas wildfire explained only moderate variability in streamflow even when wildfire accounted for >45% BAR. While the BAR in the RRW was sufficiently beyond previously reported minimum disturbance thresholds (>20% burned forest), the lack of hydrological response is attributed to buffering effects of wildfire adaptation and drought factors that are unique to Mediterranean ecoregions.

Newcomer, Michelle E.↗
DOE OSTI2024

A New Coupled Biogeochemical Modeling Approach Provides Accurate Predictions of Methane and Carbon Dioxide Fluxes Across Diverse Tidal Wetlands

Abstract Tidal wetlands provide valuable ecosystem services, including storing large amounts of carbon. However, the net exchanges of carbon dioxide (CO 2 ) and methane (CH 4 ) in tidal wetlands are highly uncertain. While several biogeochemical models can operate in tidal wetlands, they have yet to be parameterized and validated against high‐frequency, ecosystem‐scale CO 2 and CH 4 flux measurements across diverse sites. We paired the Cohort Marsh Equilibrium Model (CMEM) with a version of the PEPRMT model called PEPRMT‐Tidal, which considers the effects of water table height, sulfate, and nitrate availability on CO 2 and CH 4 emissions. Using a model‐data fusion approach, we parameterized the model with three sites and validated it with two independent sites, with representation from the three marine coasts of North America. Gross primary productivity (GPP) and ecosystem respiration (R eco ) modules explained, on average, 73% of the variation in CO 2 exchange with low model error (normalized root mean square error (nRMSE) <1). The CH 4 module also explained the majority of variance in CH 4 emissions in validation sites ( R 2 = 0.54; nRMSE = 1.15). The PEPRMT‐Tidal‐CMEM model coupling is a key advance toward constraining estimates of greenhouse gas emissions across diverse North American tidal wetlands. Further analyses of model error and case studies during changing salinity conditions guide future modeling efforts regarding four main processes: (a) the influence of salinity and nitrate on GPP, (b) the influence of laterally transported dissolved inorganic C on R eco , (c) heterogeneous sulfate availability and methylotrophic methanogenesis impacts on surface CH 4 emissions, and (d) CH 4 responses to non‐periodic changes in salinity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Test of Functional Balance Theory for Wetland Biomass Allocation in a Global Change Experiment

Abstract Forecasts of root growth and carbon sequestration under global change are compromised by uncertainty in how plants will allocate biomass between above and belowground pools. Here, we develop a simple model to assess whether functional balance theory can explain a complex biomass allocation response observed in a brackish marsh under experimental warming and elevated CO 2 . Our model shows how treatment‐driven changes in nitrogen supply and demand can explain divergent observations of root growth (i.e., maximum responses under intermediate warming and elevated CO 2 ). The model also reveals a surprising interaction between warming and eutrophication, where enhanced N loading to coastal marshes may reduce adverse impacts of warming on root growth. Our findings provide a mechanistic basis for incorporating biomass allocation into forecast models of marsh evolution. They also provide a general example of using ecological theory to decompose complex net responses observed in multi‐factor global change experiments into constituent processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Velocity Dependence of Rate‐And‐State Friction in Granular Fault Gouge and Implications for Slow‐Slip Events

Abstract The standard rate‐and‐state friction (RSF) has extensively captured frictional behaviors, but it fails to explain the velocity dependence of frictional stability transition and widespread slow‐slip events (SSEs) in experiments and nature adequately. An alternative microphysical Chen‐Niemeijer‐Spiers (CNS) model can well describe the velocity dependence of frictional behaviors of granular gouges. Using the original CNS model, standard RSF parameters can be quantified microphysically. However, some micro‐parameters are not easy to estimate quantitatively, making it difficult to extrapolate to natural and experimental conditions. Here, we simplify the microphysically‐derived RSF parameters including direct effect a , evolution effect b , and critical slip distance D c , as well as equivalent values ( a eq , b eq , and D eq ). The simplified friction parameters directly illustrate their velocity dependence, namely the essentially constant a , a eq , and D c , negatively velocity‐dependent b and b eq , as well as varying D eq for different laws. They are roughly consistent with experimental results in various fault gouges. A modified CNS model is further derived from the original CNS model, establishing a direct link between the standard RSF and CNS models. The modified CNS model exhibits virtually identical frictional behaviors to the original CNS, but differs from the standard RSF at large velocity perturbations. Moreover, the linearized stability analysis indicates that the critical stiffness for the modified CNS model is velocity‐dependent. Compared with the standard RSF, the modified CNS model not only explains the velocity dependence of frictional stability transition, but also exhibits a more gradual transition for SSEs with a broader range of stiffness ratios.

Geochemistry & Geophysics↗
DOE OSTI2025

Downward Convective Moisture Transport Dominated by a Few Overshooting Clouds in Marine and Continental Shallow Convection

Abstract In a previous study (Xiao et al., 2023, https://doi.org/10.1029/2022ms003526 ), we found that ignoring the moist convective downdrafts associated with overshooting clouds in parameterizations can lead to significant biases in the simulated depth and liquid water content of a shallow cloud layer. In this study, we seek to better quantify the properties of the clouds responsible for these moist downdrafts to help improve shallow convection parameterizations. We apply a 3‐D cloud‐tracking algorithm to large‐eddy simulations (LESs) of marine and continental shallow convection. We find that top 1% and 2% of the tracked cloud population ranked by lifetime‐mean cloud‐base mass flux can explain 90%–95% of the total downward moisture transport in the upper cloud layer whereas top 10%–20% is required to explain 90%–95% of the total upward moisture transport near mean cloud base. The vertical structure of the clouds in the top 1% and 2% (the overshooting “deep mode”) is also distinctively different from that of the rest of the cloud population (the “shallow mode”). Shallow convection parameterizations need to capture accurately the properties and convective transports of the clouds in both the deep and shallow modes. To do that, our results suggest that mass‐flux parameterizations need to (a) accurately predict the size and number of the deep‐mode clouds and (b) explicitly represent overshooting cloud updrafts and associated moist downdrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Understanding the Recent Increase in Landfalling Tropical Cyclones Over Florida's Gulf Coast

Unlike Florida's Atlantic Coast, the Gulf Coast of Florida has seen heightened tropical cyclone (TC) activity in recent decades with several destructive landfalls. Here, we attempt to understand this regional contrast using a suite of observations for the period 1979–2024. First, we demonstrate that while the El Niño Southern Oscillation (ENSO), the Atlantic Multidecadal Oscillation (AMO) and the North Atlantic Oscillation (NAO) can explain ~23% of the interannual variability in landfalls over the Gulf Coast, the variance explained by them for the Atlantic Coast is statistically insignificant. Next, we show that this striking difference may be attributed to the regional patterns of wind shear, steering flow and air-sea thermodynamic state excited by those modes of variability. The differential control exerted by ENSO, AMO, and NAO on landfalling Florida TCs, in combination with decadal trends in those modes, is likely responsible for the observed increases in landfalls over Florida's Gulf Coast.

Florida↗
DOE OSTI2025

Outer Radiation Belt Dynamics During the October 2012 Storm Revisited: Rapid Inward Radial Transport From a Dynamic Outer Boundary

Earth's outer radiation belt electron flux is highly variable and can be enhanced by over an order of magnitude over timescales less than one day, as observed during the October 2012 storm. Previous studies of this storm (e.g., Reeves et al., 2013, https://doi.org/10.1126/science.1237743) have invoked local acceleration to explain this. However, here, we argue that the observations can instead be explained by fast inward radial transport. One method often invoked to distinguish between these two acceleration processes is the existence of local peaks in electron phase space density (PSD) as a function of L* at fixed first, M, and second, K, adiabatic invariants. However, this method relies on the assumption that the evolution of the PSD as a function of L* occurs over timescales slower than the satellite orbital period. Here, high spatiotemporal resolution data from the Global Positioning System (GPS) spacecraft constellation is used to show that enhancements in the PSD occur during the October 2012 storm over short timescales not resolvable by the Van Allen Probes. In addition, Geostationary Operational Environmental Satellite spacecraft data also indicate that these enhancements are consistent with relativistic electron injections. A radial diffusion model is shown to reproduce the PSD dynamics observed by the Van Allen Probes, once rapid variations at the simulation outer boundary are included, consistent with GPS data. This verifies that apparently “locally growing” peaks in PSD along high apogee satellite orbits can be produced by fast inward radial transport without requiring the action of any local acceleration processes.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2020

Mountain biodiversity and ecosystem functions: interplay between geology and contemporary environments

Although biodiversity and ecosystem functions are strongly shaped by contemporary environments, such as climate and local biotic and abiotic attributes, relatively little is known about how they depend on long-term geological processes. Here, along a 3000-m elevational gradient with tectonic faults on the Tibetan Plateau (that is, Galongla Mountain in Medog County, China), we study the joint effects of geological and contemporary environments on biological communities, such as the diversity and community composition of plants and soil bacteria, and ecosystem functions. We find that these biological communities and ecosystem functions generally show consistent elevational breakpoints at 2000–2800 m, which coincide with Indus-Yalu suture zone fault and are similar to the elevational breakpoints of soil bacteria on another mountain range 1000 km away. Mean annual temperature, soil pH and moisture are the primary contemporary determinants of biodiversity and ecosystem functions, which support previous findings. However, compared with the models excluding geological processes, inclusion of geological effects, such as parent rock and weathering, increases 67.9 and 35.9% of the explained variations in plant and bacterial communities, respectively. Such inclusion increases 27.6% of the explained variations in ecosystem functions. The geological processes thus provide additional links to ecosystem properties, which are prominent but show divergent effects on biodiversity and ecosystem functions: parent rock and weathering exert considerable direct effects on biodiversity, whereas indirectly influence ecosystem functions via interactions with biodiversity and contemporary environments. Thus, the integration of geological processes with environmental gradients could enhance our understanding of biodiversity and, ultimately, ecosystem functioning across different climatic zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The temperature sensitivity of soil: microbial biodiversity, growth, and carbon mineralization

Abstract Microorganisms drive soil carbon mineralization and changes in their activity with increased temperature could feedback to climate change. Variation in microbial biodiversity and the temperature sensitivities (Q10) of individual taxa may explain differences in the Q10 of soil respiration, a possibility not previously examined due to methodological limitations. Here, we show phylogenetic and taxonomic variation in the Q10 of growth (5–35 °C) among soil bacteria from four sites, one from each of Arctic, boreal, temperate, and tropical biomes. Differences in the temperature sensitivities of taxa and the taxonomic composition of communities determined community-assembled bacterial growth Q10, which was strongly predictive of soil respiration Q10 within and across biomes. Our results suggest community-assembled traits of microbial taxa may enable enhanced prediction of carbon cycling feedbacks to climate change in ecosystems across the globe.

Wang, Chao (ORCID:0000000257567505)↗
DOE OSTI2021

Narrow-band high-lying excitons with negative-mass electrons in monolayer WSe2

Monolayer transition-metal dichalcogenides (TMDCs) show a wealth of exciton physics. Here, we report the existence of a new excitonic species, the high-lying exciton (HX), in single-layer WSe 2 with an energy of ~3.4 eV, almost twice the band-edge A-exciton energy, with a linewidth as narrow as 5.8 meV. The HX is populated through momentum-selective optical excitation in the K-valleys and is identified in upconverted photoluminescence (UPL) in the UV spectral region. Strong electron-phonon coupling results in a cascaded phonon progression with equidistant peaks in the luminescence spectrum, resolvable to ninth order. Ab initio GW-BSE calculations with full electron-hole correlations explain HX formation and unmask the admixture of upper conduction-band states to this complex many-body excitation. These calculations suggest that the HX is comprised of electrons of negative mass. The coincidence of such high-lying excitonic species at around twice the energy of band-edge excitons rationalizes the excitonic quantum-interference phenomenon recently discovered in optical second-harmonic generation (SHG) and explains the efficient Auger-like annihilation of band-edge excitons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Deep Gravity model for mobility flows generation

The movements of individuals within and among cities influence critical aspects of our society, such as well-being, the spreading of epidemics, and the quality of the environment. When information about mobility flows is not available for a particular region of interest, we must rely on mathematical models to generate them. In this work, we propose Deep Gravity, an effective model to generate flow probabilities that exploits many features (e.g., land use, road network, transport, food, health facilities) extracted from voluntary geographic data, and uses deep neural networks to discover non-linear relationships between those features and mobility flows. Our experiments, conducted on mobility flows in England, Italy, and New York State, show that Deep Gravity achieves a significant increase in performance, especially in densely populated regions of interest, with respect to the classic gravity model and models that do not use deep neural networks or geographic data. Deep Gravity has good generalization capability, generating realistic flows also for geographic areas for which there is no data availability for training. Finally, we show how flows generated by Deep Gravity may be explained in terms of the geographic features and highlight crucial differences among the three considered countries interpreting the model’s prediction with explainable AI techniques.

97 MATHEMATICS AND COMPUTING↗