SEARCH · Engineering Papers
Results for “PP”
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.
Evolution_and_engineering_of_pathways_for_aromatic_O-demethylation_in_Pseudomonas_putida_KT2440
Biological conversion of lignin from biomass offers a promising strategy for sustainable production of fuels and chemicals. However,aromatic compoundsderived from lignin commonly contain methoxy groups, andO-demethylation of these substrates is often a rate-limiting reaction that influences catabolic efficiency. Severalenzymefamilies catalyze aromaticO-demethylation, but they are rarely comparedin vivoto determine an optimal biocatalytic strategy. Here, two pathways for aromaticO-demethylation were compared inPseudomonas putidaKT2440. The native Rieske non-heme ironmonooxygenase(VanAB) and, separately, a heterologous tetrahydrofolate-dependentdemethylase(LigM) were constitutively expressed inP. putida, and the strains were optimized via adaptive laboratory evolution (ALE) withvanillateas a model substrate. All evolved strains displayed improved growth phenotypes, with the evolved strains harboring the native VanAB pathway exhibiting growth rates ∼1.8x faster than those harboring the heterologous LigM pathway.Enzyme kineticsandtranscriptomicsstudies investigated the contribution of selected mutations toward enhanced utilization of vanillate. The VanAB-overexpressing strains contained the most impactful mutations, including those in VanB, thereductasefor vanillateO-demethylase, PP_3494, a global regulator of vanillate catabolism, andfghA, involved in formaldehydedetoxification. These three mutations were combined into a single strain, which exhibited approximately 5x faster vanillate consumption than the wild-type strain in the first 8 h of cultivation. Overall, this study illuminates the details of vanillate catabolism in the context of two distinct enzymatic mechanisms,yielding a platform strain for efficientO-demethylation of lignin-related aromatic compounds to value-added products. This DOI contains supplementary material associated with the published manuscript.
Biofilm mitigation in hybrid chemical-biological upcycling of waste polymers
Accumulation of plastic waste in the environment is a serious global issue. To deal with this, there is a need for improved and more efficient methods for plastic waste recycling. One approach is to depolymerize plastic using pyrolysis or chemical deconstruction followed by microbial-upcycling of the monomers into more valuable products. Microbial consortia may be able to increase stability in response to process perturbations and adapt to diverse carbon sources, but may be more likely to form biofilms that foul process equipment, increasing the challenge of harvesting the cell biomass. To better understand the relationship between bioprocess conditions, biofilm formation, and ecology within the bioreactor, in this study a previously-enriched microbial consortium (LS1_Calumet) was grown on (1) ammonium hydroxide-depolymerized polyethylene terephthalate (PET) monomers and (2) the pyrolysis products of polyethylene (PE) and polypropylene (PP). Bioreactor temperature, pH, agitation speed, and aeration were varied to determine the conditions that led to the highest production of planktonic biomass and minimal formation of biofilm. The community makeup and diversity in the planktonic and biofilm states were evaluated using 16S rRNA gene amplicon sequencing. Results showed that there was very little microbial growth on the liquid product from pyrolysis under all fermentation conditions. When grown on the chemically-deconstructed PET the highest cell density (0.69 g/L) with minimal biofilm formation was produced at 30°C, pH 7, 100 rpm agitation, and 10 sL/hr airflow. Results from 16S rRNAsequencing showed that the planktonic phase had higher observed diversity than the biofilm, and that Rhodococcus, Paracoccus, and Chelatococcus were the most abundant genera for all process conditions. Biofilm formation by Rhodococcus sp. And Paracoccus sp. Isolates was typically lower than the full microbial community and varied based on the carbon source. Ultimately, the results indicate that biofilm formation within the bioreactor can be significantly reduced by optimizing process conditions and using pure cultures or a less diverse community, while maintaining high biomass productivity. The results of this study provide insight into methods for upcycling plastic waste and how process conditions can be used to control the formation of biofilm in bioreactors.
Seismic Characterization of the Blue Mountain Geothermal Field
Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.
Testing Various Cement Formulations under Temperature Cycles and Drying Shrinkage for Low-Temperature Geothermal Wells
Low-enthalpy geothermal wells are considered a sustainable energy source, particularly for district heating in the Netherlands. The cement sheath in these wells experiences thermal cycles. The stability of cement recipes under such conditions is not well understood. In this work, thermal cycling experiments for intermediate- and low-temperature geothermal well cements have been conducted. The samples were cured either under ambient conditions or under realistic pressure and temperature for 7 days. The samples did not show any signs of failure after performing 10 cycles of thermal treatment between 100 °C and 18 °C. We also tested cement formulations under drying conditions. Drying shrinkage is caused by a reduction in the water content of cement, which leads to capillary forces that can damage cement. Such circumstances lead to tensile stresses causing radial cracks. Most samples exhibited cracks under low humidity conditions (drying). Fiber reinforcement, especially using short PP fibers, improved the cement’s resilience to temperature and humidity changes. Such additives can improve the longevity of cement sheaths in geothermal wells.
Comparing the Use of High- to Low-Cost Black Carbon and Carbon Dioxide Sensors for Characterizing On-Road Diesel Truck Emissions
The exhaust plume capture method is a commonly used approach to measure pollutants emitted by in-use heavy-duty diesel trucks. Lower cost sensors, if used in place of traditional research-grade analyzers, could enable wider application of this method, including use as a monitoring tool to identify high-emitting trucks that may warrant inspection and maintenance. However, low-cost sensors have for the most part only been evaluated under ambient conditions as opposed to source-influenced environments with rapidly changing pollutant concentrations. This study compared black carbon (BC) emission factors determined using different BC and carbon dioxide (CO 2 ) sensors that range in cost from $200 to $20,000. Controlled laboratory experiments show that traditional zero and span steady-state calibration checks are not robust indicators of sensor performance when sampling short duration concentration peaks. Fleet BC emission factor distributions measured at two locations at the Port of Oakland in California with 16 BC/CO 2 sensor pairs were similar, but unique sensor pairs identified different high-emitting trucks. At one location, the low-cost PP Systems SBA-5 agreed on the classification of 90% of the high emitters identified by the LI-COR LI-7000 when both were paired with the Magee Scientific AE33. Conversely, lower cost BC sensors when paired with the LI-7000 misclassified more than 50% of high emitters when compared to the AE33/LI-7000. Confidence in emission factor quantification and high-emitter identification improves with larger integrated peak areas of CO 2 and especially BC. This work highlights that sensor evaluation should be conducted under application-specific conditions, whether that be for ambient air monitoring or source characterization.
Prospects for Heavy Neutral SUSY HIGGS Scalars in the hMSSM and Natural SUSY at LHC Upgrades
We examine production and decay of heavy neutral SUSY Higgs bosons pp→H,A→ττ¯ within the hMSSM and compare against a perhaps more plausible natural supersymmetry scenario dubbed mh125(nat) which allows for a natural explanation for mweak≃mW,Z,h∼100 GeV while maintaining mh≃125 GeV. We evaluate signal against various Standard Model backgrounds from γ,Z→ττ¯, tt¯ and vector boson pair production VV. We combine the transverse mass method for back-to-back (BtB) taus along with the ditau mass peak mττ method for acollinear taus as our signal channels. This technique ultimately gives a boost to the signal significance over the standard technique of using just the BtB signal channel. We evaluate both the 95% CL exclusion and 5σ discovery reach in the mA vs. tanβ plane for present LHC with 139 fb−1, Run 3 with 300 fb−1 and high luminosity LHC (HL-LHC) with 3000 fb−1 of integrated luminosity. For tanβ=10, the exclusion limits range up to mA∼1, 1.1 and 1.4 TeV, respectively. These may be compared to the range of mA values gleaned from a statistical analysis of the string landscape wherein mA can range up to ∼8 TeV.
Charmonium Transport in Heavy-Ion Collisions at the LHC
We provide an update on our semi-classical transport approach for quarkonium production in high-energy heavy-ion collisions, focusing on 𝐽/𝜓 and 𝜓(2𝑆) mesons in 5.02 TeV Pb-Pb collisions at the Large Hadron Collider (LHC) at both forward and mid-rapidity. In particular, we employ the most recent charm-production cross sections reported in pp collisions, which are pivotal for the magnitude of the regeneration contribution, and their modifications due to cold-nuclear-matter (CNM) effects. Multi-differential observables are calculated in terms of nuclear modification factors as a function of centrality, transverse momentum, and rapidity, including the contributions from feeddown from bottom hadron decays. For our predictions for 𝜓(2𝑆) production, the mechanism of sequential regeneration relative to the more strongly bound 𝐽/𝜓 meson plays an important role in interpreting recent ALICE data.
Machine Learning Using Rapidity-Mass Matrices for Event Classification Problems in HEP
In this work, supervised artificial neural networks (ANN) with rapidity–mass matrix (RMM) inputs are studied using several Monte Carlo event samples for various pp collision processes. The study shows the usability of this approach for general event classification problems. The proposed standardization of the ANN feature space can simplify searches for signatures of new physics at the Large Hadron Collider (LHC) when using machine learning techniques. In particular, we illustrate how to improve signal-over-background ratios in the search for new physics, how to filter out Standard Model events for model-agnostic searches, and how to separate gluon and quark jets for Standard Model measurements.
Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy
Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.
Impact of Biodiesel, Renewable Diesel, 1-Octanol, Dibutoxymethane, n-Undecane, Hexyl hexanoate and 2-Nonanone with Infrastructure Plastics as Blends with Diesel
In this work the volume and hardness were measured for thermoplastics and thermosetting resins with diesel containing up to 30% of the following blend stocks: biodiesel, renewable diesel, n-undecane, dibutoxymethane, 1-octanol, hexyl hexanoate, and 2-nonanone. Thermoplastics included polyphenylene sulfide (PPS), polyethylene terephthalate (PET), polytetrafluoroethylene (PTFE), polyvinylidene fluoride (PVDF), polyoxymethylene (POM), polybutylene terephthalate (PBT), polypropylene (PP), high density polyethylene (HDPE), nylons, acetals, polyetherimide (PEI), polyetheretherketone (PEEK), a PET co-polymer, polyphthalamides (PPAs), polyarylamide (PARA) and ethylene tetrafluoroethylene (ETFE). Three thermosetting resins were also evaluated. The material specimens were exposed to the test fuels under ambient conditions for 16 weeks. In general, the volume and hardness of the specimens were relatively unaffected following exposure to the test fuels; however, n-undecane produced significant swelling in polypropylene and may not be suitable for use with this material. N-undecane, along with hexyl hexanoate and 2-nonanone also caused low levels of swelling in nylon, which was not significant enough to preclude the use of nylons in sealing applications with these fuels.
Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators
The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.
Visible to near-infrared spectral reflectance along the NGEE-Arctic tram transect, Utqiagvik (Barrow), Alaska, 2015-2018
Canopy spectral reflectance collected measured in 2015-2018 along the NGEE-Arctic tram located on the Barrow Environmental Observatory (BEO), Barrow, Alaska. Measurements were made using a PP-Systems UniSpec-DC spectrometer mounted on an automatically operated instrument cart. Spectra were collected at 137 stop points at 0.5 m intervals along the tram transect. Transect runs were collected daily, or multiple times per day, from April or May through to August or September each year. Downwelling radiance was measured using a 2 meter fiber optic cable connected to a cosine diffuser. Upwelling (i.e. reflected) radiance was measured using a 2 meter cable connected to a 12 degree field-of-view (FOV) lens. Canopy reflectance was calculated using the ratio of upwelling to downwelling radiance measured over a 99.99% reflectance Spectralon standard measured at the start of each measurement set. Unprocessed spectral data are presented in .spu files. Metadata are provided in csv files, and the instrument manual (pdf). The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) 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).
Precision Small Scattering Angle Measurements of Proton–Proton and Proton–Nucleus Analyzing Powers at the RHIC Hydrogen Jet Polarimeter
At RHIC, the hydrogen jet target polarimeter (HJET) is used to measure proton beam polarization with accuracy $σ^{syst}_{P}$ / $P$ ≲ 0.5% by counting low energy (1–10 MeV) recoil protons in left-right symmetric detectors. The HJET performance also allowed us to precisely measure $pp$ and $pA$ (where $A$ is any ion stored at RHIC) analyzing powers in the CNI region. The results of the measurements are discussed.
Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins
This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties.csv contains all processed outputs, including aquifer properties like porosity, permeability, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R: R script that processes the shapefiles to produce the aquifer_properties file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Grid_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Grid_area: Area of the grid cell (in square meters) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data Repository. https://doi.org/10.57931/2307831 Contact Reach out to Hassan Niazi or open an issue in superwell repository for questions or suggestions.
Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins
This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.
Analysis of Infrastructures for Processing Plastic Waste using Pyrolysis-Based Chemical Upcycling Pathways
Modern mechanical recycling infrastructure for plastic is capable of processing only a small subset of waste plastics, reinforcing the need for parallel disposal methods such as landfilling and incineration. Emerging pyrolysis-based chemical technologies can "upcycle" plastic waste into high-value polymer and chemical products and process a broader range of waste plastics. In this work, we study the economic and environmental benefits of deploying an upcycling infrastructure in the continental United States for producing low-density polyethylene (LDPE) and polypropylene (PP) from post-consumer mixed plastic waste. Our analysis aims to determine the market size that the infrastructure can create, the degree of circularity that it can achieve, the prices for waste and derived products it can propagate, and the environmental benefits of diverting plastic waste from landfill and incineration facilities it can produce. We apply a computational framework that integrates techno-economic analysis, life cycle assessment, and value chain optimization. Our results demonstrate that the infrastructure generates an economy of nearly 20 billion USD and positive prices for plastic waste, opening opportunities for compensation to residents who provide plastic waste. Our analysis also indicates that the infrastructure can achieve a plastic-to-plastic degree of circularity of 34% and remains viable under various external factors (including technology efficiencies, capital investment budgets, and polymer market values). Finally, we present significant environmental benefits of upcycling over alternative landfill and incineration waste disposal methods, and comment on ongoing work expanding our modeling methodology to other chemical upcycling pathway case studies, including hydroformylation of specific plastics to chemicals.
Predictive performance of multi-model ensemble forecasts of COVID-19 across European nations
Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022. Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1–4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models’ predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models’ forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models’ past predictive performance. Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models’ forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models’ forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models’ forecasts of deaths (N=763 predictions from 20 models). Across a 1–4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models. Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks. Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 ( https://www.nfdi4health.de/task-force-covid-19-2 ) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z).