Life Cycle Analysis Tools to Evaluate Environmental Footprint of Carbon Conversion Projects
Description of NETL CO2U LCA Guidance Toolkit that includes goals, key points, overview of addenda and toolkit utilities.
SEARCH · Search NASA
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.
Description of NETL CO2U LCA Guidance Toolkit that includes goals, key points, overview of addenda and toolkit utilities.
A feasibility investigation into using optical fiber for interrogating arrays of piezoelectric sensors for Electromagnetic Surveying Applications.
This dataset provides projections of urban building morphologies for the Los Angeles urban area at 30-meter spatial resolution. It contains 192 raster files that detail two primary building attributes: building footprint fractions (ranging from 0 to 1) and average building heights (ranging from 0 to 75 meters). The projections account for a wide range of future pathways, covering two Shared Socioeconomic Pathway (SSP) scenarios (SSP3 and SSP5), two population scenarios, two developed land intensification scenarios, and four distinct levels of intensification. The dataset was created using dual Generative Adversarial Networks (GANs) trained on 2015 land cover and building properties from the National Land Cover Database (NLCD) and Model America datasets. Supporting information on the dataset has been described in the LAUrbanAreaMorphologyProjections2100_README.txt file.
Hydrogen Fueling Infrastructure Research and Station Technology (H2FIRST) is a project initiated by the DOE in 2015 and executed by Sandia National Laboratories and the National Renewable Energy Laboratory to address R&D barriers to the deployment of hydrogen fueling infrastructure. One key barrier to the deployment of fueling stations is the land area they require (i.e. "footprint"). Space is particularly a constraint in dense urban areas where hydrogen demand is high but space for fueling stations is limited. This work presents current fire code requirements that inform station footprint, then identifies and quantifies opportunities to reduce footprint without altering the safety profile of fueling stations. Opportunities analyzed include potential new methods of hydrogen delivery, as well as alternative placements of station technologies (i.e. rooftop/underground fuel storage). As interest in heavy-duty fueling stations and other markets for hydrogen grows, this study can inform techniques to reduce the footprint of heavy-duty stations as well. This work characterizes generic designs for stations with a capacity of 600 kg/day hydrogen dispensed and 4 dispenser hoses. Three base case designs (delivered gas, delivered liquid, and on-site electrolysis production) have been modified in 5 different ways to study the impacts of recently released fire code changes, colocation with gasoline refueling, alternate delivery assumptions, underground storage of hydrogen, and rooftop storage of hydrogen, resulting in a total of 32 different station designs. The footprints of the base case stations range from 13,000 to 21,000 ft 2 . A significant focus of this study is the NFPA 2 requirements, especially the prescribed setback distances for bulk gaseous or liquid hydrogen storage. While the prescribed distances are large in some cases, these setback distances are found to have a nuanced impact on station lot size; considerations of the delivery truck path, traffic flow, parking, and convenience store location are also important. Station designs that utilize underground and rooftop storage can reduce footprint but may not be practical or economical. For example, burying hydrogen storage tanks underground can reduce footprint, but the cost savings they enable depend on the cost of burial and the cost land. Siting and economic analysis of station lot sizes illustrate the benefit of smaller station footprints in the flexibility and cost savings they can provide. This study can be used as a reference that provides examples of the key design differences that fueling stations can incorporate, the approximate sizes of generic station lots, and considerations that might be unique to particular designs.
Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2 s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2 s −1 and 0.758 using HF and 0.007 μmol m −2 s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems.
As photometric surveys reach unprecedented statistical precision, systematic uncertainties increasingly dominate large-scale structure probes relying on galaxy number density. Defining the final survey footprint is critical, as it excludes regions affected by artefacts or suboptimal observing conditions. For galaxy clustering, spatially varying observational systematics, such as seeing, are a leading source of bias. Template maps of contaminants are used to derive spatially dependent corrections, but extreme values may fall outside the applicability range of mitigation methods, compromising correction reliability. The complexity and accuracy of systematics modelling depend on footprint conservativeness, with aggressive masking enabling simpler, robust mitigation. We present a unified approach to define the DES Year 6 joint footprint, integrating observational systematics templates and artefact indicators that degrade mitigation performance. This removes extreme values from an initial seed footprint, leading to the final joint footprint. By evaluating the DES Year 6 lens sample MagLim++ plus plus on this footprint, we enhance the Iterative Systematics Decontamination (ISD) method, detecting non-linear systematic contamination and improving correction accuracy. While the mask's impact on clustering is less significant than systematics decontamination, it remains non-negligible, comparable to statistical uncertainties in certain w(theta) scales and redshift bins. Supporting coherent analyses of galaxy clustering and cosmic shear, the final footprint spans 4031.04 deg2, setting the basis for DES Year 6 1x2pt, 2x2pt, and 3x2pt analyses. This work highlights how targeted masking strategies optimise the balance between statistical power and systematic control in Stage-III and -IV surveys.
Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.
Hydrogen (H 2 ), as a high-energy-density molecule, offers a clean solution to carry energy. However, the high diffusivity and low volumetric density of H 2 pose a challenge for long-term storage and transportation. Liquid organic hydrogen carriers (LOHCs) have been suggested as a strategic way to store and transport hydrogen in stable molecules. More so, electrochemical LOHC cycling renders an opportunity to utilize renewable energy for hydrogen storage and transportation toward the goal of eliminating carbon emissions. In this Perspective, examples of electrochemical reactions of organic molecules and their suitability for LOHC couples are examined. A comparative carbon footprint assessment of electrochemical LOHC cycling processes against thermochemical and hybrid LOHC cycling processes was performed. The electrochemical LOHC cycling process had the lowest relative carbon footprint only when highly concentrated LOHCs were used as the feed or when purification of the LOHC product was not required. The carbon footprint in electrochemical cycling of diluted LOHC was primarily contributed to by the LOHC distillation separation process. A sensitivity analysis showed the carbon footprint LOHC concentration dependence during the electrochemical cycling process. Moreover, the electrolyte composition significantly affects the carbon footprint during electrochemical LOHC cycling. Energy utilization, water usage, and toxicity for electrochemical LOHC cycling are discussed to provide an overview for better economic and environmental practices. There are significant opportunities in the electrochemical cycling of LOHCs if appropriate conditions such as high concentrations of reactant, reversible redox cycling ability, high Faradaic efficiencies, and catalyst stabilities are achieved.
Background Accurate estimation of net ecosystem exchange of CO 2 fluxes (Fc) is essential for understanding carbon cycle processes and assessing ecosystem carbon budgets. However, conventional modeling approaches often emphasize temporal dynamics while overlooking the pronounced spatial heterogeneity within the footprint of eddy covariance (EC) towers, potentially limiting predictive accuracy and interpretability of Fc estimates. To address this challenge, we developed a spatiotemporal model that integrates high-resolution footprint-weighted spatial information with sequential environmental drivers. Results The integrated model combines a deeper graph convolutional network to characterize fine-scale spatial variability within EC footprints and a gated recurrent unit network to capture temporal dependencies in biophysical conditions. Using multi-year flux tower observations, remote sensing vegetation indices and footprint modeling, we evaluate the proposed method across three land cover types. This spatiotemporal model consistently outperforms temporal-only and spatial-only baselines, achieving the highest overall accuracy (R 2 = 0.9569) and the lowest RMSE (1.8128 μmol m −2 s −1 ) and MAE (1.1939 μmol m −2 s −1 ). Performance gains are particularly evident in ecosystems with strong vegetation heterogeneity, where spatial structure substantially modulates Fc variability. Conclusions This study demonstrates the importance of joint modeling spatial heterogeneity and temporal dynamics for improving Fc estimation and provides a robust method for advancing footprint-based Fc estimates across diverse ecosystems, supporting refined assessments of terrestrial carbon fluxes, and enhancing scientific foundations for carbon studies.
On average, one quarter of the global carbon footprint is embodied in traded goods. These emissions are a growing issue for global efforts to decarbonize the world economy. Embodied emissions in trade are not accounted for by existing greenhouse gas accounting systems. For example, countries only report their domestic carbon emissions (also known as production-based or territorial accounting) to the Intergovernmental Panel on Climate Change (IPCC). If the embodied carbon in trade were accounted for and reported, the promising climate trends depicted by many countries would be negated or reversed. For example, many achievements of reducing emissions by developed countries under the Kyoto Protocol would actually appear as emissions outsourced to developing countries. Like many other developed countries, the United States is a net importer of embodied greenhouse gas (GHG) emissions in trade. The goal of this study was to calculate the carbon footprint of various sectors of the U.S. economy using the latest available data and based on that, estimate the embodied carbon in certain manufacturing sector products imported and exported by the U.S. First, we conducted an input-output analysis to calculate the carbon footprint of 401 aggregate level sectors of the U.S. economy. After that, for a selected number of manufacturing subsectors /products(hereafter referred to as products),we collected their trade(import/export) data and combined that with their computed carbon footprint to estimate the embodied carbon in import and export of these products for the U.S. While we conducted the carbon footprint analysis for sectors covering the entire U.S. economy (401 sectors), we only conducted the embodied carbon in trade analysis for a few products within specific aggregate level manufacturing sectors,such as computer and electronic product manufacturing (NAICS 334), transport equipment manufacturing (NAICS 336), machinery manufacturing (NAICS 333), and chemical manufacturing (NAICS 325). It is possible to do the detailed embodied carbon analysis for all other products within the various aggregate level manufacturing sectors of the U.S. economy, but given the scope of this study and space constraint of this report, we only present the results for a few products. The results of our analysis for the selected products show the embodied carbon in trade for each product by country and also show the rankings of countries in terms of top importers and exporters of embodied carbon in trade with the U.S. Our results show that in majority of cases, the U.S. is a net importer of embodied carbon in trade. Unless consumption-based accounting is used, the U.S. and other net embodied carbon importing countries may continue to outsource their emissions to meet their climate change mitigation targets under the Paris Agreement, as observed previously with the Kyoto Protocol. Some countries reported reductions that exceeded their Kyoto targets, however, the changes in emissions embodied in imports were comparable to or larger than changes in domestic emissions. Traded emissions undermined emissions reductions in the Kyoto Protocol, and threaten to continue to do so for the Paris Agreement. New climate policies such as California’s Buy Clean act, EU’s proposed border carbon tax adjustment, as well as efforts in the private sector can help address the embodied carbon in traded products.
A large sample of ground-based evapotranspiration (ET) measurements made in the United States, primarily from eddy covariance systems, were post-processed to produce a benchmark ET dataset. The dataset was produced primarily to support the intercomparison and evaluation of the OpenET satellite-based remote sensing ET (RSET) models and could also be used to evaluate ET data from other models and approaches. OpenET is a web-based service that makes field-delineated and pixel-level ET estimates from well-established RSET models readily available to water managers, agricultural producers, and the public. The benchmark dataset is composed of flux and meteorological data from a variety of providers covering native vegetation and agricultural settings. Flux footprint predictions were developed for each station and included static flux footprints developed based on average wind direction and speed, as well as dynamic hourly footprints that were generated with a physically based model of upwind source area. The two footprint prediction methods were rigorously compared to evaluate their relative spatial coverage. Data from all sources were post-processed in a consistent and reproducible manner including data handling, gap-filling, temporal aggregation, and energy balance closure correction. The resulting dataset included 243,048 daily and 5,284 monthly ET values from 194 stations, with all data falling between 1995 and 2021. We assessed average daily energy imbalance using 172 flux sites with a total of 193,021 days of data, finding that overall turbulent fluxes were understated by about 12% on average relative to available energy. Multiple linear regression analyses indicated that daily average latent energy flux may be typically understated slightly more than sensible heat flux. This dataset was developed to provide a consistent reference to support evaluation of RSET data being developed for a wide range of applications related to water accounting and water resources management at field to watershed scales.
Population data are normally collected at various census administrative levels, and areal interpolation of population is often required to transform population to the desired spatial resolution. Building footprint datasets, such as Microsoft building footprints, have proven to be useful in estimating population distribution and can therefore be used for areal interpolation of population. In addition to Microsoft building footprints, the recently released USA Structures dataset provides additional information such as building type and building height for some regions, which may provide valuable information for a better depiction of population distribution and improved population areal interpolation accuracy. In this study, we have conducted areal interpolation of population projections consistent with three different Shared Socioeconomic Pathways (SSP2, SSP3, and SSP5) from 1-km grid cells to block level in Washington state for every ten years from 2020 to 2040 based on USA Structures. We assessed USA Structures-based population downscaling accuracy using U.S. decennial survey data in 2020 under three different downscaling schemes, including population downscaling from census tracts to block groups, from census tracts to blocks, and from block groups to blocks. The resulting accuracies were compared with those based on Microsoft building footprints. The comparison showed that USA Structures achieved higher accuracies across different population density regions and areas with different urbanization extent within our study area.
Surface-atmosphere fluxes are known to vary at multiple time scales, but uncertainty is high as to how fluxes change spatially within regions. With an increase in the number of eddy covariance towers, we are now able to examine the geospatial coherence of ecosystem fluxes, using time-series correlation. Eighteen sites from Michigan and Wisconsin were used in this study, ranging from 100 m to 600 km apart. Surface-atmosphere fluxes from a six-month period were used to quantify spatial coherence on a pair-wise basis. Using geospatial statistics, carbon and sensible heat (H) fluxes were found to be 95% correlated directly outside of their flux footprint and 56% correlated up to a distance of ~35 km. Latent (LE) and momentum (τ) fluxes were less correlated, 83% directly outside of their flux footprint and 40% at a distance of ~130 km albeit, at a much larger spatial distance than for the carbon and sensible heat fluxes. All fluxes showed strong spectral resonance at diel and seasonal timescales, with 1-, 2- and 3-month periods being common modes of variability among H, LE, and τ fluxes. Results based on Empirical Orthogonal Function show distinct transitions of net ecosystem exchange from fall to winter before photosynthesis or respiration while H and τ do not exhibit coherent trends. Here, this work demonstrates the potential of quantifying geospatial coherence of surface-atmosphere fluxes in the Midwestern United States, with the ability to predict fluxes beyond the spatial limit of a single flux tower footprint. Ultimately, expanding the flux measurements to larger scales would allow better spatial scaling of terrestrial surface-atmosphere fluxes between tower footprint and modeling or remote sensing scales.
A seamless chip-to-world photonic interface enables broad advancements in optical ranging, display, communication, computation and quantum information science. The ideal solution enables two-dimensional scanning of a diffraction-limited beam from anywhere on a photonic integrated circuit to a large number of resolvable spots. Current beam-scanning technologies are limited by a fundamental trade-off: photonic-integrated-circuits with diffractive optics offer scalability but have poor mode quality, whereas inertially limited micromechanical scanners provide high-quality beams but lack scalable integration. Here we report a photonic ski-jump—a nanoscale waveguide monolithically integrated on a piezoelectric cantilever—to overcome these limitations. It passively curls ~90° out-of-plane within a less-than-0.1 mm 2 footprint, emits a submicrometre, broadband diffraction-limited beam, and exhibits kilohertz-rate mechanical resonances with quality factors of over 10,000. Fabricated in a volume complementary metal–oxide–semiconductor (CMOS) foundry, our device enables scalable two-dimensional beam scanning. Driven on-resonance at CMOS-level voltages, it achieves a footprint-adjusted spot rate of 68.6 mega spots s –1 mm–², exceeding state-of-the-art micro-electro-mechanical systems mirrors by more than 50-fold, which is sufficient for one million pixels at 100 Hz from an approximately 1.5 mm diameter footprint. We demonstrate full-colour image and video projection, and single-photon initialization and readout from silicon vacancy centres in diamond. Finally, by demonstrating uniformity across a 64 ski-jump array, we establish a pathway to achieving greater than one gigaspot resolution at kilohertz rates within a sub-5-cm-diameter footprint, creating a seamless optical pipeline between integrated photonic processors and the free-space world.
Large datasets of carbon dioxide, energy, and water fluxes were measured with the eddy-covariance (EC) technique, such as FLUXNET2015. These datasets are widely used to validate remote-sensing products and benchmark models. One of the major challenges in utilizing EC-flux data is determining the spatial extent to which measurements taken at individual EC towers reflect model-grid or remote sensing pixels. To minimize the potential biases caused by the footprint-to-target area mismatch, it is important to use flux datasets with awareness of the footprint. This study analyze the spatial representativeness of global EC measurements based on the open-source FLUXNET2015 data, using the published flux footprint model (SAFE-f). The calculated annual cumulative footprint climatology (ACFC) was overlaid on land cover and vegetation index maps to create a spatial representativeness dataset of global flux towers. The dataset includes the following components: (1) the ACFC contour (ACFCC) data and areas representing 50%, 60%, 70%, and 80% ACFCC of each site, (2) the proportion of each land cover type weighted by the 80% ACFC (ACFCW), (3) the semivariogram calculated using Normalized Difference Vegetation Index (NDVI) considering the 80% ACFCW, and (4) the sensor location bias (SLB) between the 80% ACFCW and designated areas (e.g. 80% ACFCC and window sizes) proxied by NDVI. Finally, we conducted a comprehensive evaluation of the representativeness of each site from three aspects: (1) the underlying surface cover, (2) the semivariogram, and (3) the SLB between 80% ACFCW and 80% ACFCC, and categorized them into 3 levels. The goal of creating this dataset is to provide data quality guidance for international researchers to effectively utilize the FLUXNET2015 dataset in the future.
Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.
Asymmetric double-crystal monochromators (aDCMs) and inclined DCMs (iDCMs) can significantly expand the X-ray beam footprint and consequently reduce the heat load density and gradient. Based on rigorous dynamical theory calculations, the major principles and properties of aDCMs and iDCMs are presented to guide their design and development, particularly for fourth-generation synchrotrons. In addition to the large beam footprint, aDCMs have very large bandwidths (up to ∼10 eV) and angular acceptance, but the narrow angular acceptance of the second crystal requires precise control of the relative orientations and strains. Based on Fourier coupled-wave diffraction theory calculations, it is rigorously proved that the iDCM has almost the same properties as the conventional symmetric DCM, including the efficiency, angular acceptance, bandwidth, tuning energy range and sensitivity to misalignment. The exception is that, for the extremely inclined geometry that can achieve very large footprint expansion, the iDCM has (beneficially) a larger bandwidth and wider angular acceptance. Inclined diffraction has the `rho-kick effect' that can be cancelled by the second reflection of the iDCM (even with misalignment), except that inhomogeneous strains may cause non-uniform rho-kick angles. At present, fabrication/mounting-induced strains pose low risk since they can be controlled to <0.5 µrad over large areas. The only uncertain challenge is the thermally induced strains, yet it is estimated that these strains are naturally lowered by the large footprint and may be further mitigated by optimized cryogenic cooling to the 1–2 µrad level. Overall, aDCMs and iDCMs have more stringent requirements than normal DCMs, but they are feasible schemes in practice.
Optical delay lines control the flow of light in time, introducing phase and group delays for engineering interferences and ultrashort pulses. Photonic integration of such optical delay lines is essential for chip-scale lightwave signal processing and pulse control. However, typical photonic delay lines based on long spiral waveguides require extensively large chip footprints, ranging from mm 2 to cm 2 scales. Here we present a scalable, high-density integrated delay line using a skin-depth engineered subwavelength grating waveguide, i.e., an extreme skin-depth (eskid) waveguide. The eskid waveguide suppresses the crosstalk between closely spaced waveguides, significantly saving the chip footprint area. Our eskid-based photonic delay line is easily scalable by increasing the number of turns and should improve the photonic chip integration density.