Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “event location estimation”

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 181 records · Page 10

Evolution of titanium particle combustion in potassium perchlorate and air

Understanding titanium particle combustion processes is critical not only for characterizing existing pyrotechnic systems but also for creating new igniter designs. In order to characterize titanium particle combustion processes, morphologies, and temperatures, simultaneous spatially-resolved electric field holography and imaging pyrometry techniques were used to capture post-ignition data at up to 7 kHz. Due to the phase and thermal distortions present in the combustion cloud, traditional digital in-line holography techniques fail to capture accurate data. Here in this work, electric field holography techniques are used in order to cancel distortions and capture the three-dimensional spatial locations and diameters of the particles. In order to estimate the projected surface temperatures of the titanium particles, an imaging pyrometry method that ratios emission at 750 and 850 nm is utilized. Using these diagnostics, joint statistics are collected for particle size, morphology, velocity, and temperature. Results show that, early in the combustion process, the titanium particles are primarily oxidized by potassium perchlorate inside the igniter cup, resulting in projected surface temperatures near 3000 K. Later in the process, the particles interact with ambient air, resulting in lower surface temperatures around 2400 K and the formation of flame zones. These results are consistent with adiabatic flame temperature predictions as well as particle morphology observations of a titanium core with a TiO 2 surface. Late stage particle expansion, star fragmentation, and molten droplet breakup events are also observed using the time-resolved morphology and temperature diagnostics. These results illustrate the different stages of titanium particle combustion in pyrotechnic environments, which can be used to inform improvements in next-generation igniters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Tale of Two Disks: Mapping the Milky Way with the Final Data Release of APOGEE

We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), α-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-α disk, high-α disk, and total population independently. The low-α disk exhibits a negative radial metallicity gradient of -0.06 ± 0.001 dex kpc -1 , which flattens with distance from the midplane. The high-α disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ~0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Developing New Beam Constraints to Improve Sterile Neutrino Searches at NOvA

NOvA is a long-baseline neutrino oscillation experiment utilizing two functionally identical liquid scintillator detectors located 14 mrad off-axis from the NuMI beam at Fermilab. The position of the Near Detector (ND), 1 km from the neutrino production target, means it sees a range of off-axis angles (12.6-17.2 mrad), resulting in a noticeable energy gradient across the detector face. In addition, the neutrino production points vary widely along the decay pipe, which leads to a distribution of true neutrino baselines (L), and limits the sensitivity of our current sterile neutrino searches. This work presents a new "baseline reconstruction" technique designed to exploit our off-axis position to triangulate the precise origin of neutrinos within the NuMI beamline. By projecting the total momentum vector of Charged Current (CC) and Neutral Current (NC) interactions back to the central beam axis, we estimate the individual baseline for each event. We describe the technique, evaluate its performance, and demonstrate how this can improve the sensitivity of future searches for sterile neutrinos through reduced flux systematic uncertainties and improved L/E resolution.

Burns, Jessica [Cincinnati U.]↗

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

17 WIND ENERGY↗

Quantifying the impact of climate change on erosion - 20397

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha (3,338-ac) site located approximately 48 km (30 mi) south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (metric tons, or 705 short tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. One of the key processes governing erosion is extreme precipitation and it is critical that trends in the distribution of extreme precipitation events through time be represented. The evidence of climate change is increasingly well documented and projected impacts on extreme precipitation events should be incorporated in performance assessment studies when relevant. Not evaluating future climate states in a performance assessment is contradictory to good modeling practice. Specifically, excluding climate change limits development of modeling information that could aid in effective decision making. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. In this work, this information was used to assess the impacts of potential changes in climate on erosion processes. The goal was to understand how projections of future climate relate to the performance of the WNYNSC through time. In the first stage of this work multi-temporal historical aerial images were analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. The information from these analyses was then used in conjunction with simulated data from the West Valley Erosion Working Group (EWG). Analysis of EWG simulations provides estimates for the change in erosion rates through time that is driven by changes in climate. The time-varying rates of erosion change were applied to the historical aerial imagery data in order to inform time-varying rates of erosion that are driven by changes in climate. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. The analyses from the EWG were used to characterize the relative impact of climate change on the erosion rates. This relative impact was quantified by comparing the LEM based estimates of erosion that were derived using historical climate data with LEM based estimates of erosion that were derived using projections of future climate. The relative increase in the erosion rates was then applied to the historical estimates of erosion derived from the analysis of the historical aerial imagery. This approach was used since the LEM-based estimated of the historical erosion rates from the past century have a bias towards underestimating erosion. This underestimation is hypothesized to be a consequence of an inability of the LEMs to account for the impacts of land-use change (i.e. deforestation) which had been shown to significantly increase erosion in similar Northern hardwood forest ecosystems. In summary, gully head retreat rate and gully widening rate were characterized using statistical probability distributions from the historical aerial image analysis. Simulated data from the EWG was used to estimate the climatically-driven changes in erosion rates through time. The changes through time from the EWG were applied to the gully head retreat rate and gully widening rate characterized using the historical aerial image analysis. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Near-Term Electricity Requirement and Emission Implications for Sustainable Aviation Fuel Production with CO 2 -to-Fuels Technologies

Aviation contributed approximately 10% of the U.S. transportation sector's greenhouse gas (GHG) emissions and about 3% of the nation's total GHG production before the coronavirus disease (COVID-19) pandemic (EPA 2022). Unlike the ground transportation sector, which can be decarbonized by using batteries and hydrogen fuel cell powertrain technologies, the technical and economic challenges of aviation electrification open the opportunity for CO 2 utilization (CO2U) using clean power sources. The U.S. government set a goal to produce 3 billion gallons per year of sustainable aviation fuels (SAF) by 2030 and scale up the production to 35 billion gallons per year by 2050 (Bioenergy Technologies Office 2022). In this report, we examine three potential locations, in California, Iowa, and Louisiana, for SAF production using two production pathways that are expected to be available by 2030. We analyze the electricity cost to satisfy 10% of the SAF production potential for the select locations. Specifically, we consider (1) retail electricity cost from the default local utility in each location; (2) physical power purchase agreement (PPA) for renewable power and battery storage hybrid systems with preset prices; (3) financial PPA from a dedicated renewable plant; and (4) estimated real-time pricing (RTP) from the wholesale power market with utility delivery adders. Retail rates are available for the three potential locations; however, the developer should negotiate with the local utility in Iowa for new tariff riders, because the load of the proposed SAF plant (1.4 GW) is significantly higher than the current industrial tariff structure requirement (200 kW). The developers are not able to claim federal credits (i.e., the Inflation Reduction Act of 2022) when sourcing electricity from local utilities. Developers could consider prioritizing physical and financial PPAs as purchase options. While physical PPAs can be imported out of state, this structure imposes availability issues as developers must locate in the same electricity grid regions. Financial PPAs do not have the same location restrictions and provide the same cost savings as physical PPAs, but financial PPAs impose financial risks in the event of system curtailment and grid interruption. RTPs are available for customers in California with flexible load and hourly load management capability; however, utilities in our studied regions in Louisiana and Iowa only offer time-of-use and curtailment programs to incentivize lower energy usage in peak hours. These utilities do not currently offer RTP programs, and the developer may need to have further negotiations with the local utilities to access RTP programs. The results show that purchased electricity prices may range from 2.6 cents/kWh to 7.1 cents/kWh for the three studied plants.

02 PETROLEUM↗

Conceptual Design Report for the MATHUSLA Long-Lived Particle Detector near CMS

We present the Conceptual Design Report (CDR) for the MATHUSLA (MAssive Timing Hodoscope for Ultra-Stable neutraL pArticles) long-lived particle detector at the HL-LHC, covering the design, fabrication and installation at CERN Point 5. MATHUSLA is a 40 m-scale detector with an air-filled decay volume that is instrumented with scintillator tracking detectors, to be located near CMS. Its large size, close proximity to the CMS interaction point and about 100 m of rock shielding from HL-LHC backgrounds allows it to detect LLP production rates and lifetimes that are one to two orders of magnitude beyond the ultimate sensitivity of the HL-LHC main detectors for many highly motivated LLP signals. Data taking is projected to commence with the start of HL-LHC operations. We present a new 40m design for the detector: its individual scintillator bars and wavelength-shifting fibers, their organization into tracking layers, tracking modules, tower modules and the veto detector; define a high-level design for the supporting electronics, DAQ and trigger system, including supplying a hardware trigger signal to CMS to record the LLP production event; outline computing systems, civil engineering and safety considerations; and present preliminary cost estimates and timelines for the project. We also conduct detailed simulation studies of the important cosmic ray and HL-LHC muon backgrounds, implementing full track/vertex reconstruction and background rejection, to ultimately demonstrate high signal efficiency and $\ll 1$ background event in realistic LLP searches for the main physics targets at MATHUSLA. This sensitivity is robust with respect to detector design or background simulation details. Appendices provide various supplemental information.

Aitken, Branden [Victoria U.]↗

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crowd cluster data in the USA for analysis of human response to COVID-19 events and policies

We provide data on daily social contact intensity of clusters of people at different types of Points of Interest (POI) by zip code in Florida and California. This data is obtained by aggregating fine-scaled details of interactions of people at the spatial resolution of 10 m, which is then normalized as a social contact index. We also provide the distribution of cluster sizes and average time spent in a cluster by POI type. This data will help researchers perform fine-scaled, privacy-preserving analysis of human interaction patterns to understand the drivers of the COVID-19 epidemic spread and mitigation. Current mobility datasets either provide coarse-level metrics of social distancing, such as radius of gyration at the county or province level, or traffic at a finer scale, neither of which is a direct measure of contacts between people. We use anonymized, de-identified, and privacy-enhanced location-based services (LBS) data from opted-in cell phone apps, suitably reweighted to correct for geographic heterogeneities, and identify clusters of people at non-sensitive public areas to estimate fine-scaled contacts.

60 APPLIED LIFE SCIENCES↗

Tail Dependence as a Measure of Teleconnected Warm and Cold Extremes of North American Wintertime Temperatures

Current models for spatial extremes are concerned with the joint upper (or lower) tail of the distribution at two or more locations. Such models cannot account for teleconnection patterns of 2-m surface air temperature ( T 2m ) in North America, where very low temperatures in the contiguous United States may coincide with very high temperatures in Alaska in the wintertime. This dependence between warm and cold extremes motivates the need for a model with opposite-tail dependence in spatial extremes. This work develops a statistical modeling framework that has flexible behavior in all four pairings of high and low extremes at pairs of locations. In particular, we use a mixture of rotations of common Archimedean copulas to capture various combinations of four-corner tail dependence. We study teleconnected T 2m extremes using ERA5 of daily average 2-m temperature during the boreal winter. Further, the estimated mixture model quantifies the strength of opposite-tail dependence between warm temperatures in Alaska and cold temperatures in the midlatitudes of North America, as well as the reverse pattern. These dependence patterns are shown to correspond to blocked and zonal patterns of midtropospheric flow. This analysis extends the classical notion of correlation-based teleconnections to considering dependence in higher quantiles.

54 ENVIRONMENTAL SCIENCES↗

A Climatology and Extreme Value Analysis of Large Hail in China

In this paper, hail size records are analyzed at 2254 stations in China and a hail size climatology is developed based on gridded hail observations for the period 1960–2015. It is found that the annual percentiles of hail size records changed sharply and national-wide after 1980, therefore two periods, 1960–79 and 1980–2015, are studied. There are some similarities between the two periods in terms of the characteristics of hail size such as the spatial distribution patterns of mean annual maximum hail size and occurrence week of annual maximum hail size. The 1980–2015 period had higher observation density than the 1960–79 period, but showed smaller mean annual maximum hail size, especially in northern China. In the majority of grid boxes, the annual maximum hail size experienced a decreasing trend during the 1980–2015 period. A Gumbel extreme value model is fitted to each grid box to estimate the return periods of maximum hail size. The scale and location parameter of the fitted Gumbel distributions are higher in eastern China than in western China, thereby reflecting a greater likelihood of large hail in eastern China. In southern China, the maximum hail size exceeds 127 mm for a 10-yr return period, whereas in northern China maximum hail size exceeds this threshold for a 50-yr return period. Furthermore, the Gumbel model is found to potentially underestimate the maximum hail size for certain return periods, but provides a more informed picture of the spatial distribution of extreme hail size and the regional differences.

54 ENVIRONMENTAL SCIENCES↗

Oregon Roadway CSZ Liquefaction and Landslide Impact Screening Analysis

The roadway transportation routes and the surrounding environments in Oregon will be subject to ground acceleration and potential associated deformation due to the propagation of waves of energy during a Cascadia Subduction Zone (CSZ) earthquake. Linear features such as highways may be subject to partial or complete disruption of their function resulting from ground deformations at and along the roadway. A potential source of these disruptions associated with a CSZ event are landslides and liquefaction in the soils along and beneath the highways and local roadways. The construction of the Oregon roadway system and urban development have occurred largely along the banks of the rivers—those along the Pacific coast, and along the valleys throughout the Oregon Coast Range and Cascade Mountains. These locations also represent areas where landslides and liquefaction have occurred historically, or may potentially occur in the future, presenting a significant vulnerability for the roadway system and dependent post-CSZ earthquake emergency response and recovery effort. Many of the Oregon areas with landslide- and liquefaction-prone soils have the potential to impact roadway transportation routes after a major earthquake causing delays in establishing routes to transport emergency response personnel and resources to impacted areas and communities. This report presents methods used to project the impacts that CSZ earthquake-induced landslides and liquefaction may have on roadways across Oregon. First, this report presents an approach for estimating the risk of landslides to Oregon roadways, and determining approximate reopening times. Second, this report presents a similar approach for estimating the risk of liquefaction-induced ground displacements to roadway pavements, also resulting in approximate reopening times for roadway segments. The outcomes of these two ground failure analyses serve as inputs to a broader analysis of CSZ earthquake impacts to Oregon’s surface transportation system and its ability to support post-disaster response and recovery activities, as part of the Oregon Transportation Systems Regional Resiliency Assessment Program project (CISA 2021).

58 GEOSCIENCES↗

Cosmogenic surface exposure ( 10 Be) dating of raised beaches in Marguerite bay, Antarctic Peninsula: Implications for relative sea-level history

Understanding the dynamics of ice mass loss in polar regions is crucial for deciphering climate change and Glacio Isostatic Adjustment patterns. This study focuses on Marguerite Bay, located in the south-central Antarctic Peninsula. We dated raised beaches to investigate relative sea-level changes using the cosmogenic surface exposure ( 10 Be) method. Previous studies have provided valuable insights into the region's glacial history, but limitations in dating techniques and age estimates necessitate further investigation. By analysing raised shingle beaches in Gaul Cove of Horseshoe Island and the southern coast of Calmette Bay, this research aims to contribute relative sea-level change history for these areas. In Horseshoe Island's Gaul Cove, raised beaches clustered on prominent steps reveal a 15 m relative sea-level change over the last 3.31 ka. Differently, Calmette Bay exhibits a 36 m relative sea-level fall over the last 7.29 ka. Finally, these findings indicate significant and differential glacial-isostatic adjustments in both regions during the middle and late Holocene. Additionally, our data reveal accelerated sea-level fall periods corresponding to Holocene deglaciation and glacial advance events, indicating the shorelines' relative sea-level change sensitivity to climate change.

58 GEOSCIENCES↗

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

58 GEOSCIENCES↗

The Strange Case of Ground-Coupled Airwaves on Seismoacoustic Stations at Local to Near-Regional Scales

Here we investigate the application of ground-coupled airwaves observed by seismoacoustic stations at local to near-regional scales to detect signals of interest and determine back-azimuth information. Ground-coupled airwaves are created from incident pressure waves traveling through the atmosphere that couple to the earth and transmit as a seismic wave with retrograde elliptical motion. Previous studies at sub-local scales (<10 km from a source of interest) found the back-azimuth to the source could be accurately determined from seismoacoustic signals recorded by acoustic and 3-component seismic sensors spatially separated on the order of 10 to 150 m. The potential back-azimuth directions are estimated from the coherent signals between the acoustic and vertical seismic data, via a propagation-induced phase shift of the seismoacoustic signal. A unique solution is then informed by the particle motion of the 3-component seismic station, which was previously found to be less accurate than the seismoacoustic-sensor method. We investigate the applicability of this technique to greater source-receiver distances, from 50-100 km and up to 400 km, which contains pressure waves with tropospheric and stratospheric ray paths, respectively. Specifically, we analyze seismoacoustic sources with ground truth from rocket motor fuel elimination events at the Utah Test and Training Range (UTTR) as well as a 2020 rocket launch in Southern California. From these sources we observe evidence that while coherent signals can be seen from both sources on multiple seismoacoustic station pairs, the determined ground-coupled airwave back-azimuths are more complicated than results at more local scales. Our findings suggest more complex factors including incidence angle, coupling location, subsurface material, and atmospheric propagation effects need to be fully investigated before the ground-coupled airwave back-azimuth determination method can be applied or assessed at these further distances.

58 GEOSCIENCES↗

Nonlinear Photovoltaic Degradation Rates: Modeling and Comparison Against Conventional Methods

Although common practice for estimating photovoltaic (PV) degradation rate (RD) assumes a linear behavior, field data have shown that degradation rates are frequently nonlinear. This article presents a new methodology to detect and calculate nonlinear RD based on PV performance time-series from nine different systems over an eight-year period. Prior to performing the analysis and in order to adjust model parameters to reflect actual PV operation, synthetic datasets were utilized for calibration purposes. A change-point analysis is then applied to detect changes in the slopes of PV trends, which are extracted from constructed performance ratio (PR) time-series. Once the number and location of change points is found, the ordinary least squares method is applied to the different segments to compute the corresponding rates. The obtained results verified that the extracted trends from the PR time-series may not always be linear and therefore, “nonconventional” models need to be applied. All thin-film technologies demonstrated nonlinear behavior whereas nonlinearity detected in the crystalline silicon systems is thought to be due to a maintenance event. A comparative analysis between the new methodology and other conventional methods demonstrated levelized cost of energy differences of up to 6.14%, highlighting the importance of considering nonlinear degradation behavior

14 SOLAR ENERGY↗

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗