Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reanalysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Low ozone dry deposition rates to sea ice during the MOSAiC field campaign: Implications for the Arctic boundary layer ozone budget

Dry deposition to the surface is one of the main removal pathways of tropospheric ozone (O 3 ). We quantified for the first time the impact of O 3 deposition to the Arctic sea ice on the planetary boundary layer (PBL) O 3 concentration and budget using year-round flux and concentration observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign and simulations with a single-column atmospheric chemistry and meteorological model (SCM). Based on eddy-covariance O 3 surface flux observations, we find a median surface resistance on the order of 20,000 s m –1 , resulting in a dry deposition velocity of approximately 0.005 cm s –1 . This surface resistance is up to an order of magnitude larger than traditionally used values in many atmospheric chemistry and transport models. The SCM is able to accurately represent the yearly cycle, with maxima above 40 ppb in the winter and minima around 15 ppb at the end of summer. However, the observed springtime ozone depletion events are not captured by the SCM. In winter, the modelled PBL O 3 budget is governed by dry deposition at the surface mostly compensated by downward turbulent transport of O 3 towards the surface. Advection, which is accounted for implicitly by nudging to reanalysis data, poses a substantial, mostly negative, contribution to the simulated PBL O 3 budget in summer. During episodes with low wind speed (<5 m s –1 ) and shallow PBL (<50 m), the 7-day mean dry deposition removal rate can reach up to 1.0 ppb h –1 . Our study highlights the importance of an accurate description of dry deposition to Arctic sea ice in models to quantify the current and future O 3 sink in the Arctic, impacting the tropospheric O 3 budget, which has been modified in the last century largely due to anthropogenic activities.

54 ENVIRONMENTAL SCIENCES↗

Wave and Offshore Wind Resource in the U.S. Pacific Ocean Minor Outlying Islands

Coastal environments such as islands have unique opportunities for renewable energy resources. This work explores the wave and offshore wind energy potential for the U.S. Pacific Ocean Minor Outlying Islands, including Baker Island, Howland Island, Jarvis Island, Johnston Atoll, Kingman Reef, Palmyra Atoll, and Wake Island. A numerical wave model based on WAVEWATCH III and SWAN (both are NWS / NOAA models) was developed, validated, and executed to generate a 32-year hindcast dataset suitable for resource characterization for each island. A complementary offshore wind resource characterization is provided for the same time period based on the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5). Each island has five netcdf files: 1. wave_[island]_annual_J_map.nc: the data used to generate the map of annual averaged omnidirectional wave power (J [kW/m]) 2. wave_[island]_100m_IEC: six IEC parameters at 100m isobath from 1979/01/01 to 2020/12/31, output hourly 3. wave_[island]_5km_IEC: six IEC parameters at 5km from the shoreline 1979/01/01 to 2020/12/31, output hourly 4. wind_island_100mASL_windspeed_timeseries.nc: wind speed timeseries from ERA5 at 100 m above sea level within a 60 km radius of the island (m/s) 1979/01/01 to 2010/12/31, output hourly 5. wind_island_100mASL_winddirection_timeseries.nc: wind direction timeseries from ERA5 at 100 m above sea level within a 60 km radius of the island (degrees) 1979/01/01 to 2010/12/31, output hourly

16 TIDAL AND WAVE POWER↗

Dynamically Downscaled (WRF) 1km, Hourly Meteorological Conditions 1987-2020. East/Taylor Watersheds

This dataset contains meteorological output from the Weather Research and Forecasting (WRF) version 3.8.1. This dataset has been created to 1) investigate hydrometeorological processes impacting the East River and water-delivery to the Critical Zone, and 2) provide meteorological forcing data for distributed Earth-science modeling applications in the East River watershed. Variables have a 1 kilometer spatial resolution and hourly temporal resolution and encompass a rectangular region encompassing the East and Taylor River watersheds, Colorado, near the town of Crested Butte. WRF was forced using Climate Forecast System Reanalysis (CFSR) lateral boundary conditions. Each .zip file contains one "water year" of data (October 1 -- September 30; i.e. water year 2017 starts October 1, 2016 and ends September 30, 2017). Each zip folder contains 12 netcdf (.nc) files containing one month of hourly data each and are approximately 250mb. Model timestamps are in UTC time.The files contain the following data variables:EAST_MASK: binary mask of the watershed regionTAYLOR_MASK: binary mask of the watershed regionGLW downwelling longwave radiation (w/m2)HR_PRCP: Hourly Precipitation Rate (mm/hr). Includes all hydrometeors (solid+liquid). HFX NoahMP LSM total grid-cell modelled sensible heat flux (w/m2) [positive towards atmosphere]LH NoahMP LSM total grid-cell modelled latent heat flux (w/m2) [positive towards atmosphere; can be converted to ET]PSFC Surface Barometric Pressure (hPa)Q2 Two-meter specific humidity (kg/kg)SWDOWN Downwelling shortwave solar radiation (w/m2)SWNORM Terrain-normal downwelling shortwave radiation (w/m2)T2 Two-meter air temperature (deg K)U10 10-m U-component of wind velocity (m/s)V10 10m V-component of wind velocity (m/s)XLAT Latitude of grid-center point XLONG Longitude of grid-center point XTIME Model timestamp, **in UTC**

54 ENVIRONMENTAL SCIENCES↗

Recent streamflow trends across permafrost basins of North America: Datasets

Climate change impacts, including changing temperatures, precipitation, and vegetation, are widely anticipated to cause major shifts to the permafrost with resulting impacts to hydro-ecosystems across the high latitudes of the globe. However, it is challenging to examine streamflow shifts in these regions owing to a paucity of data, discontinuity of records, and other issues related to data consistency and accuracy. We looked at recent changes in streamflow over 1976-2021 in watersheds affected by varying degrees of permafrost coverage to characterize trends and drivers for a range of watersheds across North America. Data sets are described in detail in the paper associated with this data set, Bennett et al. 2023, Front. Water - Water and Critical Zone, DOI: 10.3389/frwa.2023.1099660.These data contain CSV files of the streamflow, climate, and land surface characteristics for several sites located across the high latitude regions of North America. Both observed and reanalysis data products are provided. These files can be opened using Excel or a text editor, or they can be read, and analyzed in software tools such as Python or R. A brief description of the files is below, and more details can be found in the Methods section.rabpro_stats_north_select_74_55m.csv - This file describes the observed gages used in the analysis.GF31_23_metadat.csv - This file describes the 23 permafrost systems. Columns are as described in rabpro_stats_north_select_74_55m.csv above, with rabpro_id, the id used for the timeseries file mapping in GF31_23_time_series.csv.GF31_23_time_series.csv - This file contains the time series data for the stations described in GF31_23_metadat.csv.GF31_random_reaches_1583.csv - This file describes the 1583 randomly selected permafrost-dominant sites for machine learning analysis.era5_GF31_monthly_vars_random_reaches_1583.csv - This file contains the monthly ERA5 land data for the 1583 randomly selected permafrost-dominant sites.observed.zip: USGS and Hydat station data for the 74 gages analyzed in this study. Monthly, seasonal, and annual streamflow observations for minimum streamflow, mean streamflow, and maximum streamflow. Units are m3/sec. 1975-2022.streamflow_daily_GF31_infilled_1979_2022.csv - infilled daily streamflow data (infilled using GloFAS v 3.1) for 55 gages. Units are m3/sec. 1979-2022.glofas_23.zip - Glofas v3.1 file for the 23 permafrost-dominant gages in the study. Monthly, seasonal, and annual streamflow observations for minimum streamflow, mean streamflow, and maximum streamflow. Units are m3/sec. 1979-2021.glofas_1583.zip - Glofas v3.1 file for the 1583 randomly selected permafrost-dominant gages in the study. Monthly, seasonal, and annual streamflow observations for minimum streamflow, mean streamflow, and maximum streamflow. Units are m3/sec. 1979-2021.

54 ENVIRONMENTAL SCIENCES↗

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗

Integrated Hourly Meteorological Database of 20 Meteorological Stations (1981-2022) for Watershed Function SFA Hydrological Modeling

This dataset contains (a) a script “R_met_integrated_for_modeling.R”, and (b) associated input CSV files: 3 CSV files per location to create a 5-variable integrated meteorological dataset file (air temperature, precipitation, wind speed, relative humidity, and solar radiation) for 19 meteorological stations and 1 location within Trail Creek from the modeling team within the East River Community Observatory as part of the Watershed Function Scientific Focus Area (SFA). As meteorological forcings varied across the watershed, a high-frequency database is needed to ensure consistency in the data analysis and modeling. We evaluated several data sources, including gridded meteorological products and field data from meteorological stations. We determined that our modeling efforts required multiple data sources to meet all their needs. As output, this dataset contains (c) a single CSV data file (*_1981-2022.csv) for each location (20 CSV output files total) containing hourly time series data for 1981 to 2022 and (d) five PNG files of time series and density plots for each variable per location (100 PNG files). Detailed location metadata is contained within the Integrated_Met_Database_Locations.csv file for each point location included within this dataset, obtained from Varadharajan et al., 2023 doi:10.15485/1660962. This dataset also includes (e) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata and (f) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. Review the (g) ReadMe_Integrated_Met_Database.pdf file for additional details on the script, methods, and structure of the dataset.The script integrates Northwest Alliance for Computational Science and Engineering’s PRISM gridded data product, National Oceanic and Atmospheric Administration’s NCEP-NCAR Reanalysis 1 gridded data product (through the `RCNEP` R package, Kemp et al., doi:10.32614/CRAN.package.RNCEP), and analytical-based calculations. Further, this script downscales the input data into hourly frequency, which is necessary for the modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

A model for $F_L$ structure function at low values of $Q^2$ and $x$ - revisited

A reanalysis of the model for the longitudinal structure function F_L (x,Q^2) F L ( x , Q 2 ) at low x x and low Q^2 Q 2 was undertaken, in view of the advent of the EIC. The model includes the kinematic constraint F_L\sim Q^4 F L ∼ Q 4 as Q^2\rightarrow Q 2 → 0. It is based on the photon-gluon fusion mechanism suitably extrapolated to the region of low Q^2 Q 2 . Revised model was critically updated, extended to the EIC kinematic region, and e.g. contains new parameterisations of the parton distribution functions.

Badelek, Barbara↗

Comparison of the Three NIF Ablators

Indirect drive implosion experiments on NIF have now been performed using three different ablator materials: glow discharge polymer (GDP) or CH, high density carbon (HDC, which we also refer to as diamond), and sputtered beryllium (Be). It has been appreciated for some time that each of these materials has specific advantages and disadvantages as an ICF ablator.[1-4] In light of experiments conducted on NIF in the last few years, how do these ablators compare? Given current understanding, is any ablator more or less likely to reach ignition on NIF? Has the understanding of their respective strengths and weaknesses changed since NIF experiments began? How are those strengths and weaknesses highlighted by implosion designs currently being tested or planned for testing soon? This document aims to address these questions by combining modern simulation results with a survey of the current experimental data base. More particularly, this document is meant to fulfill an L2 Milestone for FY17 to “Document our understanding of the relative advantages and disadvantages of CH, HDC, and Be designs.” Note that this document does not aim to recommend a down-selection of the current three ablator choices. It is intended only to gather and document the current understanding of the differences between these ablators and thereby inform the choices made in planning future implosion experiments. This document has two themes: (i) We report on a reanalysis project in which post-shot simulations were done on a common basis for layered shots using each ablator. This included data from keyholes, 2D ConA, and so forth, from each campaign, leading up to the layered shots. (“Keyholes” are shots dedicated to measuring the shock timing in a NIF target, as described in Ref. 5. “2DConAs” are backlit implosions in which the symmetry of the implosion is measured between about half and full convergence, as described in Ref. 6.) This set of common-basis postshot simulations is compared to the respective shots. Each was then scaled to a “full NIF” experiment that could be done using the respective ablators at full NIF power and/or energy, and these scaled-up designs were simulated in detail. (ii) The report also contains a general survey of experimental and simulated results as pertinent to comparing and evaluating the three ablators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Methodology for Yield Determination of Surface Shots Over Water

Shockwave analysis is one of the few methods used to determine the yield of a nuclear detonation. It has been noted that the shockwave of these surface shots over water are asymmetric making accurate yield determination difficult. To calculate the yield of these shots more accurately, reanalysis of these films is required. Armed with the digitized replicas of the original nuclear test films and the latest image processing tools, a methodology for segmenting the shockwave from the digitized film and measuring these shockwaves is devised to quickly and accurately quantify the shape of these asymmetric shockwaves. With a two-dimensional result from these asymmetric shockwaves, several methods are discussed to transform the final shockwave contour to a scalar value which can be used in existing yield determination equations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DE-NE0008277_PROTEUS final technical report 2018

This project details re-evaluations of experiments of gas-cooled fast reactor (GCFR) core designs performed in the 1970s at the PROTEUS reactor and create a series of International Reactor Physics Experiment Evaluation Project (IRPhEP) benchmarks. Currently there are no gas-cooled fast reactor (GCFR) experiments available in the International Handbook of Evaluated Reactor Physics Benchmark Experiments (IRPhEP Handbook). These experiments are excellent candidates for reanalysis and development of multiple benchmarks because these experiments provide high-quality integral nuclear data relevant to the validation and refinement of thorium, neptunium, uranium, plutonium, iron, and graphite cross sections. It would be cost prohibitive to reproduce such a comprehensive suite of experimental data to support any future GCFR endeavors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Building and Running TDAAPS Models: WRF Postdictions

This document serves to guide a researcher through the process of running the Weather Research and Forecasting (WRF) model and incorporating observations into coarse resolution reanalysis products to model atmospheric conditions at high (50 m) resolution. This documentation is specific to WRF and the WRF Preprocessing System (WPS) version 3.8.1 and the Objective Analysis (OBSGRID) code released on April 8, 2016. Output from WRF serves as an input into the Time-Domain Atmospheric Acoustic Propagation Suite (TDAAPS) which performs staggered-grid finite difference modeling of the acoustic velocity pressure system to produce Green's functions through these atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Off-Gas Condensate and Cast Stone Analysis Results

A sample of off-gas condensate produced during a demonstration test of the direct-feed low-activity waste (DFLAW) Radioactive Waste Test Platform using waste retrieved from Hanford storage tank AP-105 waste was analyzed for major cations, anions, total cyanide, Resource Conservation and Recovery Act (RCRA) metals, selected radionuclides, and selected organic compounds. The remaining off-gas condensate was concentrated through evaporation and then used to produce Cast Stone waste forms. After the solid waste forms cured for 28 days, U.S. Environmental Protection Agency (EPA) Method 1311 (EPA 1992) Toxicity Characteristic Leach Procedure (TCLP) tests were conducted on two samples. The TCLP leachates were analyzed for RCRA metals, fluoride, and total cyanide. These chemical analyses were conducted to collect data on the waste form performance and to determine if the waste form will meet the disposal requirements of the Waste Control Specialists LLC (WCS) Federal Waste Disposal Facility (FWF) in Texas (WCS 2015) and the Integrated Disposal Facility at Hanford. The radionuclide results in the off-gas condensate were used to estimate concentrations in the Cast Stone waste forms; these estimates were then compared to Class C limits (30 TAC 336.362 Appendix E). With the possible exception of U-232, all measured radionuclides were determined to be below their respective Class C limits. Inductively coupled plasma-mass spectrometry (ICP-MS) was used to measure most of the radionuclides. It is not possible to differentiate the isotope measured at mass 232 between U 232 and Th-232 by mass spectroscopy. Class C limit comparisons were made by assuming two maximum possible scenarios; one in which all the isotope at mass 232 occurred as Th-232, and another assigning the entire mass contribution to U-232. If the isotope at mass 232 is all Th-232 then the waste form is well below its Class C limit; however, if the isotope at mass 232 is all U-232, then the Class C limit would be exceeded by more than a factor of 8. A chemical procedure to separate uranium from thorium followed by reanalysis by ICP-MS would be required to resolve this issue. A similar issue occurs at mass 238, which could be U-238 or Pu-238; however, in this case Pu-238 can be ruled out as a contributor to the isotope measured at mass 238 for the following reasons. Pu-238 has a half-life of 87.7 years and its daughter product (U-234) has a half-life of 2.5 × 105 years and a measured concentration of < 0.002 µg/L. If significant concentrations of Pu-238 were present in this waste, then the amount of measurable U 234 would have increased over the decades of storage in the Hanford tanks. Because this was not observed, it is concluded that the isotope measured at mass 238 is all U-238. A total of 68 organic compounds were targeted for analysis in the off-gas condensate. Of these compounds, only nine were measured above their respective minimum reportable concentration. These compounds were acenaphthene, acetone, n-butyl alcohol, p,p’-DDT, methyl ethyl ketone, nitrobenzene, o nitrophenol, pyridine, and toluene. Of the detectable compounds, only acenaphthene, acetone, o nitrophenol, and pyridine were above wastewater standards (40 CFR 268.48 - Universal Treatment Standards). None of the organic compounds exceeded the non-wastewater standards (40 CFR 268.48 - Universal Treatment Standards). Analysis of the TCLP leachates from the Cast Stone waste forms indicated that none of the analytes exceeded the Universal Treatment Standard and, in general, all were well below the standards.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Quasi-Biennial Oscillation in Model-Simulated and Observed Equatorial Temperature of the Lower Stratosphere

In an otherwise thorough survey of metrics for models that simulate the Quasi-Biennial Oscillation, Schenzinger et al. say “there is no comparable observational dataset” to compare with model-simulated stratospheric temperatures (Geosci. Model Dev. 10, 2157-2168, 2017). They compare model temperatures only with reanalysis datasets, which come from numerical weather prediction models and therefore are not direct observations. But stratospheric brightness temperatures have been directly observed by microwave radiometers on board a series of satellites since 1979. Christy and Drouilhet (J. Climate 7, 106-120, 1994) found a QBO signal in 1979-1992 data from the Microwave Sounding Unit’s lower stratospheric Channel 4 (with peak weighting function at the 75hPa pressure level and half-power values at 120 and 40hPa, about 20±5 km altitude). To detect the QBO signal, Christy and Drouilhet used a latitude weighting that matches the QBO’s anti-correlation between the Equator and ±15°. Fernandez et al. also found a QBO signal in 1979-2000 MSU data, using EOF analysis (J. Climate 17, 3934-3946, 2004). Recently, Christy and Covey extended Christy and Drouilhet’s analysis through 2017 (“Stratospheric Temperature,” Section 2b5 in State of the Climate in 2017. Bull. Amer. Meteor. Soc. 99, S150–S152; Figure 2.10c). The QBO signal was again clearly visible. Although temperature trends for the lower atmosphere are sensitive to inter-satellite calibration and other processing details, giving rise to some controversy, data for the lower stratosphere are more robust. Christy and Covey’s Figure 2.9 shows that data processed separately by the National Oceanic and Atmospheric Administration (NOAA) and by Remote Sensing Systems (RSS) gives a QBO index that is virtually identical to that of Christy’s group at the University of Alabama at Huntsville (UAH).

58 GEOSCIENCES↗

Storm Surge Modeling in Puget Sound

In this study, a storm surge model of the Salish Sea was developed to simulate storm surge in Puget Sound using the unstructured-grid, finite-volume community ocean model (FVCOM). The model was based on Pacific Northwest National Laboratory’s existing Puget Sound hydrodynamic model that has been validated extensively in previous studies. To identify representative historical storm events, a systematic approach was developed based on the non-tidal residual (NTR) method. A total of 34 storm surge events between 1980 and 2016 were identified and simulated. To accurately simulate storm surge inside Puget Sound, the model was driven by reanalyzed Climate Forecast System Reanalysis meteorological forcing and observed water levels at the entrance of Strait of Juan de Fuca and southern entrance of Johnstone Strait. Model results for total water level, storm surge, and currents were compared with field measurements. A series of error statistics parameters were also calculated to quantify the model’s skills in simulating storm surge in the Salish Sea. Error statistics of model performance metrics demonstrated that the Salish Sea storm surge model was able to simulate storm surge inside Puget Sound with high accuracy. Storm surge propagation into Puget Sound is a nonlinear process that cannot be simply determined by a static approach for the entire domain. The maximum NTR distribution inside Puget Sound suggests that storm surge tends to be amplified in several sub-basins and inlets of the Puget Sound, such as Bellingham Bay, Hood Canal, south Puget Sound, and the multi-inlet basin behind Agate and Rich Passages. Therefore, future study is necessary to further quantify the storm surge level in these high-risk areas through field measurements and refined model simulations.

16 TIDAL AND WAVE POWER↗

Ocean and Sea Ice and their Interactions around Greenland and the West Antarctic Peninsula in Forced Fine-Resolution Global Simulations

Anomalously warm ocean waters have been implicated as the cause of accelerated ablation along the margin of the Greenland ice sheet, as well as the cause of increasing basal melt rate, mass loss, and grounding line retreat of many Antarctic ice shelves. The overarching science objective of this project is to investigate how the interplay of regional processes and decadal changes in local and remote forcing impact the delivery and end member composition of waters over the continental shelves of Greenland and Antarctica. To address this goal, an existing global eddy-active ocean general circulation model (OGCM), coupled to a thermodynamic/dynamic ice model and forced with atmospheric reanalysis fluxes from 1948-2009, was analyzed to understand these processes around the continental margins of Greenland and Antarctica. This simulation did not contain a representation of ice sheet melt or resolve mesoscale eddies over these continental margins, so two new coupled OGCM/sea-ice models were constructed and analyzed to address the roles of small scale features and processes (< 5km) and freshwater flux surpluses.

54 ENVIRONMENTAL SCIENCES↗

High-Resolution Regional Wave Hindcast for Hawaii

This report summarizes modeling efforts for hindcasting of wave climate within the exclusive economic zone (EEZ) around the state of Hawaii. Specifically, the report provides detailed description of model mesh development, sensitivity analysis of the configuration of model parameters, data used for model forcing and validation, model skill assessment, and model outputs for resource characterization. The unstructured, nested-grid modeling approach incorporates a regional high-resolution SWAN (Simulating WAves Nearshore) model within the global-regional WAVEWATCH III® (WW3) wave model. The models were configured to comply with the International Electrotechnical Commission (IEC) Technical Standards for wave energy resource characterization. The SWAN model domain covers the entire EEZ with a spatial resolution of approximately 300 m in the nearshore region. Two levels of WW3 nested grids, from 0.5 arc-degree at global to 7.5 arc-minute at regional scale, were used to provide boundary conditions that drive the high-resolution nearshore SWAN model. The WW3 and SWAN models were forced by 0.5 arc-degree Climate Forecast System Reanalysis wind and 5 km regional Weather Research and Forecasting wind, respectively. Wave resource and bulk wave parameters were simulated for a 32-year period from 1979 to 2010, and were subsequently validated using wave buoy data within the model domain. Model performance was evaluated using a set of error statistics. Overall, the model shows good skills in reproducing the wave climate in the Hawaii EEZ in both spatial and temporal domains. Comparisons of time series and scatter plots for the six IEC resource parameters and the associated error statistics tables, and horizontal two-dimensional plots of annually averaged and monthly averaged IEC parameters are provided.

16 TIDAL AND WAVE POWER↗

Extraction and Parametrization of Isobaric Trinucleon Elastic Cross Sections and Form Factors

By mining data from Jefferson Lab Hall A experiment E08-014 a new 3 He elastic cross section was extracted from a large quasielastic background. This measurement was taken with an initial beam energy of 3.356 GeV and an angle of 20.51°. The cross section was found to be 1.345 × 10 -6 µb/sr ± 0.086 × 10 -6 µb/sr at Q 2 = 34.19 fm -2 . This new data point falls approximately halfway between the first and second diffractive minima of the 3 He form factors. When combined with recent high Q 2 3 He elastic cross section measurements from Jefferson Lab and MIT-Bates this new data point improves our knowledge of the cross section and form factors at large momentum transfers. The new high Q 2 data motivate a reanalysis of the 3 He elastic cross section world data and provide an improved understanding of the magnetic form factor in particular. For this analysis the elastic cross section world data for 3 He, and its mirror nuclei 3 H, were collected. The world data spans a time frame from 1965 to 2016. The dataset contains electron energy ranges from tens of MeV to above 12 GeV for measurements performed at many different facilities. The world data were then fit using a sum of Gaussians parametrization which allowed for the extraction of both targets’ magnetic and electric form factors which were then used to calculate charge densities and radii. The new charge and magnetic form factors for 3 H and the charge form factor for 3 He are in good agreement with the 1994 fits from Amroun et al . However, the addition of the new high Q 2 data has caused the 3 He magnetic form factor’s first diffractive minimum to shift up in Q 2 by 1-3 fm -2 while also decreasing the magnitude of the magnetic form factor above Q2 ≈ 20 fm -2 . The first diffractive minima for 3H are located at Q2 ≈ 13 fm -2 and Q 2 ≈ 23-24 fm -2 for the charge and magnetic form factors respectively. The first diffractive minima for 3 He are located at Q 2 ≈ 11 fm -2 and Q 2 ≈ 17-20 fm -2 for the charge and magnetic form factors respectively. The charge radius for 3 He was found to be 1.90 fm ± 0.00144 fm in reasonable agreement with past measurements, and the charge radius for 3 H was found to be 2.02 fm ± 0.0133 fm which is much larger than past measurements. However, each of these charge radii has an additional uncertainty that must be applied to them due to allowing all parameters to float freely during the sum of Gaussians fitting procedure. This additional uncertainty should be small for 3 He, but it is likely quite significant for 3 H and would help bring this charge radius closer to agreement with past measurements that made different fitting choices. Unfortunately, this analysis was unable to quantify this additional uncertainty. The new form factor fits were compared to modern theoretical predictions from the 2016 paper of Marcucci et al . The ‘conventional’ theoretical approach applied in this paper modelled two and three-body nucleon interactions with relativistic corrections and was reasonably successful at predicting the charge form factors of 3 H and 3 He. χEFT predictions were also often successful. However, while the ‘conventional’ approach still performed best, theory failed to accurately predict the magnetic form factors of either 3 H or 3 He. The first diffractive minimum of the new 3 He magnetic form factor fits actually moved further away from theory. This disagreement between theory and experiment provides motivation for new asymmetry measurements using polarized 3 He and a polarized electron beam. When the beam is scanned in Q 2 on the target the sign of the asymmetry will flip at the form factor minima pinning down their true location.

Barcus, Scott↗