SEARCH · Engineering Papers
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.
T2Well/ECO2M simulates coupled well-reservoir systems (where the reservoir is a porous medium) consisting of water, salt, and CO2. T2Well/ECO2M is an enhancement of the T2Well/ECO2N because it can simulate transitions between gaseous and liquid CO2, and model three-phase flow in the well.
Powersheds is an open scientific software project for simulating river–reservoir cascades. It combines the performance of Rust with a friendly Python interface to model storage, pool elevation, head, releases, spills, routing lags, and power generation at hourly resolution. Designed for coupling with power-system models, simulations are driven by plant-level target power schedules and report realized generation after accounting for hydrologic and operational constraints.
The microphysical quasi-equilibrium in an ascending adiabatic parcel is elucidated by an analytical theory in 0D with drastic simplifications. The theory predicts how in-cloud activation is most likely to be triggered by the onset of precipitation during sufficient ascent, with the ascent only needing to approach almost twice the cloud-base updraft speed aloft. The precipitation and cloud condensate mass fields are coupled in a closed system in a simplified version of the model. The initial state of no precipitation is unstable with respect to a perturbation. In the 2D phase space of both mass fields, there is a neutral line. Unstable growth of precipitation mass drives the system to cross the line into a regime of stability. An attractor is then approached consisting of precipitation mass balanced by accretion of cloud mass and fallout. The cloud-particle number concentration also approaches a stable equilibrium governed by the balance between in-cloud activation aloft from the inexorably increasing supersaturation during vertical acceleration and accretion of cloud droplets by precipitation. Dimensionless numbers characterizing the microphysical equilibria and their stability are derived mathematically, including a condensation–precipitation efficiency and an in-cloud activation efficiency. The theory explains common observations of the orders of magnitude of liquid water content in convective and stratiform clouds. Finally, sensitivity tests of the numerically integrated theoretical equations are documented with respect to variations in cloud condensation nucleus (CCN) aerosols and updraft speed. In conclusion, it is shown that this theory of in-cloud activation, with an increasing supersaturation during ascent, applies to both ice-only and liquid-only cloud.
Plasmon-phonon coupling between metamaterials and molecular vibrations provides a new path for studying mid-infrared light-matter interactions and molecular detection. So far, the coupling between the plasmonic resonances of metamaterials and the phonon vibrational modes of molecules has been realized under linearly polarized light. Here, mid-infrared chiral plasmonic metasurfaces with high circular dichroism (CD) in absorption over 0.65 in the frequency range of 50 to 60 THz are demonstrated to strongly interact with the phonon vibrational resonance of polymethyl methacrylate (PMMA) molecules at 52 THz, under both left-handed and right-handed circularly polarized (LCP and RCP) light. The mode splitting features in the absorption spectra of the coupled metasurface-PMMA systems under both circular polarizations are studied in PMMA layers with different thicknesses. The relation between the mode splitting gap and the PMMA thickness is also revealed. The demonstrated results can be applied in areas of chiral molecular sensing, thermal emission, and thermal energy harvesting.
We demonstrate an optical waveguide device, capable of supporting the high, in-vacuum, optical power necessary for trapping a single atom or a cold atom ensemble with evanescent fields. Our photonic integrated platform, with suspended membrane waveguides, successfully manages optical powers of 6 mW (500 μ m span) to nearly 30 mW (125 μ m span) over an un-tethered waveguide span. This platform is compatible with laser cooling and magneto-optical traps (MOTs) in the vicinity of the suspended waveguide, called the membrane MOT and the needle MOT, a key ingredient for efficient trap loading. We evaluate two novel designs that explore critical thermal management features that enable this large power handling. This work represents a significant step toward an integrated platform for coupling neutral atom quantum systems to photonic and electronic integrated circuits on silicon.
A quantum transducer converts an input signal to an output probe at a distant frequency band while maintaining the quantum information with high fidelity, which is crucial for quantum networking and distributed quantum sensing and computing. In terms of microwave–optical quantum transduction, the state-of-the-art quantum transducers suffer low transduction efficiency from weak nonlinear coupling, wherein increasing pump power to enhance efficiency inevitably leads to thermal noise from heating. Moreover, we reveal that the efficiency-bandwidth product of a cavity electro-optical or electro-optomechanical transducer is fundamentally limited by pump power and nonlinear coupling coefficient, irrespective of cavity engineering efforts. To overcome this fundamental limit, we propose to noiselessly boost the transduction efficiency by consuming intraband entanglement (e.g., microwave–microwave or optical–optical entanglement in the case of microwave–optical transduction). Via a squeezer–coupler–antisqueezer sandwich structure, the protocol enhances the transduction efficiency to unity in the ideal lossless case, given an arbitrarily weak pump and nonlinear coupling. In practical cavity systems, our entanglement-assisted protocol surpasses the non-assisted fundamental limit of the efficiency-bandwidth product and reduces the threshold cooperativity for positive quantum capacity by a factor proportional to two-mode squeezing gain. Given a fixed cooperativity, our approach increases the broadband quantum capacity by orders of magnitude. The entanglement-assisted advantage is robust to ancilla loss and cavity detuning.
We propose a novel quantum transduction hybrid system based on the coupling of long-coherence time superconducting cavities with electro-optic resonators to achieve high-efficiency and high-fidelity in quantum communication protocols and quantum sensing.
Year-round observations of the physical snow and ice properties and processes that govern the ice pack evolution and its interaction with the atmosphere and the ocean were conducted during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to September 2020. This work was embedded into the interdisciplinary design of the 5 MOSAiC teams, studying the atmosphere, the sea ice, the ocean, the ecosystem, and biogeochemical processes. The overall aim of the snow and sea ice observations during MOSAiC was to characterize the physical properties of the snow and ice cover comprehensively in the central Arctic over an entire annual cycle. This objective was achieved by detailed observations of physical properties and of energy and mass balance of snow and ice. By studying snow and sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the variability across scales can be considered. On-ice observations of in situ and remote sensing properties of the different surface types over all seasons will help to improve numerical process and climate models and to establish and validate novel satellite remote sensing methods; the linkages to accompanying airborne measurements, satellite observations, and results of numerical models are discussed. We found large spatial variabilities of snow metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly variable snow cover needs to be considered in more detail (in observations, remote sensing, and models) to better understand snow-related feedback processes. The ice pack revealed rapid transformations and motions along the drift in all seasons. The number of coupled ice–ocean interface processes observed in detail are expected to guide upcoming research with respect to the changing Arctic sea ice.
This data package contains mean values for cation concentration for water samples taken from the East River Watershed in Colorado. Inductively coupled plasma mass spectrometry (ICP-MS) has been used to measure the concentrations of elements of interest simultaneously for the East River Watershed, Colorado groundwater and surface water samples to inform insights on the biogeochemistry processes within the watershed. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States. For samples collected prior to 06-16-2021, the instrumentation, Elan DRC II, PerkinElmer SCIEX, automatically switches among the three models necessary to analyze all 37 elements. These 37 elements include: (1) Lithium (Li), Beryllium (Be), Boron (B), Sodium (Na), Magnesium (Mg), Aluminium (Al), Silicon (Si), Phosphorus (P), Titanium (Ti), Cobalt (Co), Nickel (Ni), Copper (Cu), Zinc (Zn), Germanium (Ge), Arsenic (As), Rubidium (Rb), Strontium (Sr), Zirconium (Zr), Molybdenum (Mo), Silver (Ag), Cadmium (Cd), Tin (Sn), Antimony (Sb), Caesium (Cs), Barium (Ba), Europium (Eu), Lead (Pb), Thorium (Th), Uranium (U) using standard model, argon Ar as reaction gas, (2) Potassium (K), Calcium (Ca), Vanadium (V), Chromium (Cr), Manganese (Mn), Iron (Fe) using dynamic reaction cell (DRC) model, ammonia NH3 as reaction gas, and (3) Phosphorus (P) and Selenium (Se) using DRC model, oxygen O2 as reaction gas. Note for the samples with higher concentrations of chloride (Cl-), asenic (As) concentrations were analysed with DRC model (oxygen O2 as reaction gas) to avoid the interference of chloride. For samples collected on and after 06-16-2021, an advanced Agilent 8900 triple quadrupole inductively coupled plasma mass spectrometry system (Agilent 8900 QQQ ICP-MS, Agilent Technologies) has been used to measure the concentrations of interested 36 elements simultaneously for environmental samples, including (1) Lithium (Li), Beryllium (Be) and Boron (B) using standard no gas mode, (2) Sodium (Na), Magnesium (Mg), Aluminium (Al) Phosphorus (P), Potassium (K), Chromium (Cr), Manganese (Mn), Iron (Fe), Cobalt (Co), Nickel (Ni), Copper (Cu), Zinc (Zn), Germanium (Ge), Arsenic (As), Rubidium (Rb), Strontium (Sr), Zirconium (Zr), Molybdenum (Mo), Silver (Ag), Cadmium (Cd), Tin (Sn), Antimony (Sb), Cesium (Cs), Barium (Ba), Europium (Eu), Lead (Pb), Thorium (Th) and Uranium (U) using standard helium (He) collision mode, (3) Titanium (Ti) and Vanadium (V) using high Energy (HEHe) helium (He) collision mode, and (4) Silicon (Si), Calcium (Ca) and Selenium (Se) using standard H2 reaction mode. All samples were prepared/diluted with 2% (v/v) ultrapure nitric acid in Milli-Q water (18.2 mega ohm-cm), and analyzed under a rigorous quality assurance and quality control (QA/QC) process. This data package contains (1) a zip file (cation_data_2014_2025.zip) containing a total of 5,849 files: 5.848 data files of cation data from across the Lawrence Berkeley National Laboratory (LBNL) Watershed Function Scientific Focus Area (SFA) which is reported in .csv files per location and a locations.csv (1 file) with latitude and longitude for each location; (2) a file-level metadata (v6_20260901_flmd.csv) file that lists each file contained in the dataset with associated metadata; (3) a data dictionary (v6_20260901_dd.csv) file that contains terms/column_headers used throughout the files along with a definition, units, and data type; (4) PDF and docx files for the detemination of Method Detection Limits (MDLs) for ICP-MS PerkinElmer DRC II instrumentation (Detemination_of_Method_Detection_Limits__MDLs__for_ICP_MS__PerkinElmer_Elan_DRC_II__LBL_Bldg74_Lab214D) for samples before November 2021; (5) PDF and docx files for the determination of MDLs for ICP-MS Agilent 8900 QQQ instrumentation (ICP_MS_Analysis_detection_limits_and_QA_QC_WenmingDong_updated_2026-08-06) for samples November 2021 and onward. Missing values within the anion data files are noted as either "-9999" or "0.0" for not detectable (N.D.) data. There are a total of 113 locations containing cation data. Update on 2021-04-11: Added Detemination of Method Detection Limits (MDLs) for ICP-MS document, which can be accessed as a PDF or with Microsoft Word. Update on 2022-06-10: versioned updates to this dataset was made along with these changes: (1) updated cation data for all locations up to 2021-12-31, (2) removal of units from column headers in datafiles, (3) added row underneath headers to contain units of variables, (4) removed suffix and prefix on two variables (“aqberylliumion_asberyllium” and “aqlithiumion_aslithium”), (5) added -9999 for empty numerical cells, and (6) the addition of the file-level metadata (flmd.csv) and data dictionary (dd.csv) were added to comply with the File-Level Metadata Reporting Format. Update on 2022-09-09: Updates were made to reporting format specific files (file-level metadata and data dictionary) to correct swapped file names, add additional details on metadata descriptions on both files, add a header_row column to enable parsing, and add version number and date to file names (v2_20220909_flmd.csv and v2_20220909_dd.csv). Update on 2023-08-08: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-01-05. The file level metadata and data dictionary files were updated to reflect the additional data added. Update on 2024-03-11: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-10-16. Further, revisions to the data files were made to remove incorrect data points (from 1970 and 2001). The reporting format specific files were updated to reflect the additional data added. Updated versions of the PDF and docx files for determination of MDLs for ICP-MS data were added to this dataset for samples starting in November 2021. Update on 2025-05-15: Updates were made to both the data files and reporting format specific files. New available cation data was added, up until the end of WY2024 (September 30, 2024). International Generic Sample Numbers (IGSNs), when registered, were added to the data files. The reporting format specific files were updated to reflect the additional data added. Update on 2026-09-01: Updates were made to both the data files and reporting format specific files. New available cation data was added, up until the end of WY2025 (September 30, 2025). Updated versions, as of 2026-08-06, of the PDF and docx files for determination of MDLs for ICP-MS data were added to this dataset for samples starting in November 2021.
This data package contains raw leaf gas exchange data used to support the publication "Leaf isoprene and monoterpene emission distribution across hyperdominant tree genera in the Amazon basin," Phytochemistry, in-press. Data was collected from 162 individuals, distributed among 25 botanical families and 113 species, from four sites in the Amazon. The data was collected using a modified portable photosynthesis system for environmental control coupled to a thermal desorption tube for volatile isoprenoid collections. Collected TD tube samples were analysed in Manaus, Brazil, using thermal desorption-gas chromatography-mass spectrometry (TD-GC-MS). The dataset includes both photosynthesis and isoprene/monoterpene emission rates from leaves under standard environmental conditions as well as controlled light response curves. For detailed methods and raw data, refer to the publication and related material cited in Dataset References.
plutus is an R package that post-processes outputs from the Global Change Analysis Model (GCAM) to calculate the power sector capital investment costs and stranded asset costs associated with GCAM projections of electricity generation by technology. GCAM is a market equilibrium model used to examine the dynamics of the coupled human-Earth system and the response of this system to global change, including socioeconomics, technology, climate, and policy. GCAM tracks electricity generation by technology and vintage over 32 geopolitical regions throughout the lifetime of each generating technology. plutus extends GCAM functionality by (1) estimating the foregone economic value of prematurely retired power plants as a result of economic-induced retirements compared to scheduled lifetimes; and (2) estimating the new installations and capital investments driven by future changes in economic, energy, agriculture, and land-use systems in GCAM.
The Colorado River Basin (CRB) is experiencing an unprecedented water shortage crisis brought upon by a combination of factors arising from interactions across the region’s coupled human and natural systems. Allocation of water to the seven states that rely on the Colorado River was settled in the Colorado River Compact of 1922 during a period now known to be characterized by atypically high flows (Christensen et al., 2004). Since then, aridification due to anthropogenic-driven warming has steadily reduced the overall water supply available in the basin, with a 10% decrease in the river’s flow occurring over just the past two decades (Bass et al., 2023). The river is further strained by increasing demands associated with a growing population and diverse multi-sectoral demands. Navigating these challenges also requires accounting for the complex prior appropriation water rights system governing water allocation across the region’s diverse users.
In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.
Information and data from radiolysis testing to measure hydrogen (H₂) generated from hydrated oxides on aluminum exposed to ⁶⁰Co radiation were reviewed to evaluate hydrogen generation with radiation dose. Hydrogen generation rate is a primary input to the coupled thermal-chemical system model of the performance of aluminum-clad spent nuclear fuel in a dry storage canister (ASNF-in-canister) over its storage life. Hydrogen generation data and hydrogen generation rate (hydrogen generated per absorbed radiation dose, $\Delta$(H₂)/$\Delta$(dose)) were reported in two reports for radiolysis testing [1, 2]. The testing involved first-time data¹ for radiolytic yield of hydrogen from hydrated-oxides-on-aluminum substrates, and from pristine (non-corroded) aluminum substrates. The laboratory-grown hydrated-oxides-on-aluminum substrates contained a ~ 5 μm film consisting of the trihydroxide bayerite (Al(OH)₃) with amounts of the oxyhydroxide boehmite (AlOOH). These specimens were assumed to also contain physisorbed water (unquantified). The pristine substrates were flat coupons of aluminum that were ground to a 600-grit finish but were not immersed in water to grow a hydrated oxide. The pristine substrates were tested as companion specimens to provide information on hydrogen generation from material without chemisorbed water and were assumed to contain physisorbed water (unquantified).
The industrial manufacturing industry is looking to implement new methods of energy production to ensure a consistent energy supply while reducing economic costs and environmental impacts. Much of the manufacturing industry relies on fossil fuels—primarily coal and natural gas—of which there are finite resources subject to price volatility due to an inelastic demand. These resources also come with significant negative environmental impacts related to emissions. While complete independence from fossil fuels is not immediately realistic, options are available to significantly reduce dependency on fossil fuel supplies, including integrated energy systems (IESs) which tightly couple nuclear energy source generators with energy consumers (i.e., industrial factories) to fulfill power and energy requirements. Once realized, optimized IESs may yield significant economic and environmental benefits over traditionally isolated generator and consumer facilities. This report presents a feasibility study for siting an IES to meet the steam and electricity needs of Eastman Chemical Company’s facility in Kingsport, Tennessee. This study explored reactor technology options, evaluation, and optimization, with a focus on meeting the facility’s operational and reliability requirements. This work is part of an ongoing effort by Eastman to be good environmental stewards and meet the demands of their customers by shifting to more environmentally sustainable solutions for their energy needs. Outcomes of this report are also of general interest and value to any design or development team seeking to deploy an IES, as many topics described herein are generally applicable to any nuclear power facility.
A grand challenge in hydrologic science is to understand why signals of climate change and variability, which are often visible in precipitation extremes at aggregate scales, are not consistently observed in the case of extreme flooding. However, a solution to this challenge may prove elusive unless the water cycle is viewed in an integrative manner. Thus, for riverine flooding, while Hortonian (infiltration excess) runoff may have stronger correlation with precipitation extremes and hence perhaps to warming trends or climate oscillators, Dunne (saturation excess) runoff may have a more complex relationships with time series of precipitation and with evaporation and transpiration, but rain-on-snow and snowmelt events may depend on land-surface and atmospheric temperatures. Atmospheric rivers and tropical cyclones lead to precipitation or flooding and are impacted by climate. Flooding assessments need to consider long-term baselines, evolving risk factors, coupled natural-human systems, and novel adaptation such as nature-inspired design.
Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.
The slowdown of the Atlantic Meridional Overturning Circulation (AMOC) and associated consequences on ocean carbon uptake could have large implications for the Earth's climate system and its global carbon cycle. This study analyzes ten Earth System Models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and reveals that the regional carbon uptake in the subpolar North Atlantic under a high CO 2 emission scenario moderately correlates with the decline in AMOC at 40°N. AMOC transports warm and salty subtropical waters to the subpolar regions. Models with stronger AMOC slowdown generally exhibit weaker surface warming and larger decline of surface salinity and alkalinity. We consider two plausible mechanisms linking the AMOC slowdown to the decline of regional CO 2 uptake: the reduction in surface alkalinity and diminished subduction. The decline of surface salinity and alkalinity reduces the ocean's capacity to buffer acids leading to a reduced CO 2 uptake. This important contribution is unique to the North Atlantic. Diminished convective mixing and subduction of surface water can further decrease the downward transport of anthropogenic carbon, as also shown in previous research. The centennial trends of p CO 2 are decomposed into four components driven by temperature, salinity, alkalinity and dissolved inorganic carbon, revealing that alkalinity and dissolved inorganic carbon are both significant contributors. The alkalinity-driven p CO 2 essentially follows surface salinity, establishing the linkage between AMOC slowdown and alkalinity decline. Our results indicate that alkalinity changes are important for the interplay between AMOC and the regional carbon sequestration ability across the late 20 th and the entirety of the 21 st century in the subpolar North Atlantic.