Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Light-quark connected intermediate-window contributions to the muon 𝑔 − 2 hadronic vacuum polarization from lattice QCD

We present a lattice-QCD calculation of the light-quark connected contribution to window observables associated with the leading-order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, $a^{HVP,LO}_μ$. We employ the MILC Collaboration’s isospin-symmetric QCD gauge-field ensembles, which contain four flavors of dynamical highly improved staggered quarks with four lattice spacings between $a ≈ 0.06-0.15$ fm and close-to-physical quark masses. We consider several effective-field-theory-based schemes for finite volume and other lattice corrections and combine the results via Bayesian model averaging to obtain robust estimates of the associated systematic uncertainties. After unblinding, our final results for the intermediate and "W2" windows are $a^{U,W}_μ$(conn) = 206.6(1.0) x 10 -10 and $a^{U,W2}_μ$ (conn) = 100.7(3.2) x 10 -10 , respectively

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Effect of competition of spherical and deformed shells in quasifission of superheavy nuclei

Quasifission and fusion-fission represent the two most likely reaction outcomes to occur postcapture in collisions leading to superheavy nuclei. As such, understanding these mechanisms and how they relate to one another is key to understanding the intricate dynamics that drive the formation (or dissociation) of the nascent compound nuclei formed in fusion reactions. This understanding directly translates to a more informed picture of suitable reaction partners and can provide vital information for experimental efforts to study the physics and chemistry of superheavy elements. In this work we report results from time-dependent simulations of 48 Ca + 238 U and 50 Ti + 236 Th reactions at incident energies just above the Coulomb barrier with a focus on the quasifission process that prevent the formation of a fully equilibrated 286 Cn compound nucleus. We study these reactions systematically and consider a wide range of initial configurations to extract a robust estimate of primary fragment yields for the quasifission process. Multiple preferred exit channels are observed, with both spherical and deformed shell effects in the heavy and light fragments driving contributions to the production yields depending on the initial configuration of the system. Orientation effects of the deformed actinide targets are found to be a primary driver of which exit channels are populated. Here, the impact of moving away from a doubly magic projectile is explored with implications towards the reactions considered for current and future superheavy searches.

Gumbel, Richard [Michigan State University, East L↗
DOE OSTI2020

Stochastic Gradients for Large-Scale Tensor Decomposition

Tensor decomposition is a well-known tool for multiway data analysis. This work proposes using stochastic gradients for efficient generalized canonical polyadic (GCP) tensor decomposition of large-scale tensors. GCP tensor decomposition is a recently proposed version of tensor decomposition that allows for a variety of loss functions such as Bernoulli loss for binary data or Huber loss for robust estimation. Here, the stochastic gradient is formed from randomly sampled elements of the tensor and is efficient because it can be computed using the sparse matricized-tensor times Khatri--Rao product tensor kernel. For dense tensors, we simply use uniform sampling. For sparse tensors, we propose two types of stratified sampling that give precedence to sampling nonzeros. Numerical results demonstrate the advantages of the proposed approach and its scalability to large-scale problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantile-Dependent Expressivity of Serum Uric Acid Concentrations

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., serum uric acid) is high or low relative to its distribution. Analyses were performed to test whether serum uric acid heritability is quantile-specific and whether this could explain some reported gene-environment interactions. Methods. Serum uric acid concentrations were analyzed from 2151 sibships and 12,068 offspring-parent pairs from the Framingham Heart Study. Quantile-specific heritability from offspring-parent regression slopes (β OP , h 2 = 2β OP /(1 + r spouse )) and full-sib regression slopes (β FS ,h 2 = {(1 + 8r spouse β FS ) 0.5 - 1}/(2r spouse ) was robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. Results. Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s sex- and age-adjusted uric acid distribution when estimated from β OP (P trend = 0.001): 0.34 ± 0.03 at the 10th, 0.36 ± 0.03 at the 25th, 0.41 ± 0.03 at the 50th, 0.46 ± 0.04 at the 75th, and 0.49 ± 0.05 at the 90th percentile and when estimated from β FS (P trend = 0.006). This is consistent with the larger genetic effect size of (1) the SLC2A9 rs11722228 polymorphism in gout patients vs. controls, (2) the ABCG2 rs2231142 polymorphism in men vs. women, (3) the SLC2A9 rs13113918 polymorphism in obese patients prior to bariatric surgery vs. two-year postsurgery following 29 kg weight loss, (4) the ABCG2 rs6855911 polymorphism in obese vs. nonobese women, and (5) the LRP2 rs2544390 polymorphism in heavier drinkers vs. abstainers. Quantile-dependent expressivity may also explain the larger genetic effect size of an SLC2A9/PKD2/ABCG2 haplotype for high vs. low intakes of alcohol, chicken, or processed meats. Conclusions. Heritability of serum uric acid concentrations is quantile-specific.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Changes in Annual Extremes of Daily Temperature and Precipitation in CMIP6 Models

This study presents an analysis of daily temperature and precipitation extremes with return periods ranging from 2 to 50 years in phase 6 of the Coupled Model Intercomparison Project (CMIP6) multimodel ensemble of simulations. Judged by similarity with reanalyses, the new-generation models simulate the present-day temperature and precipitation extremes reasonably well. In line with previous CMIP simulations, the new simulations continue to project a large-scale picture of more frequent and more intense hot temperature extremes and precipitation extremes and vanishing cold extremes under continued global warming. Changes in temperature extremes outpace changes in global annual mean surface air temperature (GSAT) over most landmasses, while changes in precipitation extremes follow changes in GSAT globally at roughly the Clausius–Clapeyron rate of ~7% °C –1 . Changes in temperature and precipitation extremes normalized with respect to GSAT do not depend strongly on the choice of forcing scenario or model climate sensitivity, and do not vary strongly over time, but with notable regional variations. Over the majority of land regions, the projected intensity increases and relative frequency increases tend to be larger for more extreme hot temperature and precipitation events than for weaker events. To obtain robust estimates of these changes at local scales, large initial-condition ensemble simulations are needed. Appropriate spatial pooling of data from neighboring grid cells within individual simulations can, to some extent, reduce the needed ensemble size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantile-specific heritability of plasma fibrinogen concentrations

Fibrinogen is a moderately heritable blood protein showing different genetic effects by sex, race, smoking status, pollution exposure, and disease status. These interactions may be explained in part by “quantile-dependent expressivity”, where the effect size of a genetic variant depends upon whether the phenotype (e.g. plasma fibrinogen concentration) is high or low relative to its distribution. Determine whether fibrinogen heritability (h 2 ) is quantile-specific, and whether quantile-specific h 2 could account for fibrinogen gene-environment interactions. Plasma fibrinogen concentrations from 5689 offspring-parent pairs and 1932 sibships from the Framingham Heart Study were analyzed. Quantile-specific heritability from offspring-parent (β OP , h 2 = 2β OP /(1+r spouse )) and full-sib regression slopes (β FS , h 2 = {(1+8r spouse β FS ) 0.05 –1}/(2r spouse )) were robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s age- and sex-adjusted fibrinogen distribution when estimated from β OP (P trend = 5.5x10 -6 ): 0.30±0.05 at the 10 th , 0.37±0.04 at the 25 th , 0.48±0.05 at the 50 th , 0.61±0.06 at the 75 th , and 0.65±0.08 at the 90 th percentile, and when estimated from β FS (P trend = 0.008): 0.28±0.04 at the 10 th , 0.31±0.04 at the 25 th , 0.36±0.03 at the 50 th , 0.41±0.05 at the 75 th , and 0.50±0.06 at the 90 th percentile. The larger genetic effect at higher average fibrinogen concentrations may contribute to fibrinogen’s greater heritability in women than men and in Blacks than Whites, and greater increase from smoking and air pollution for the FGB -455G>A A-allele. It may also explain greater fibrinogen differences between: 1) FGB -455G>A genotypes during acute phase reactions than usual conditions, 2) GTSM1 and IL-6 -572C>G genotypes in smokers than nonsmokers, 3) FGB -148C>T genotypes in untreated than treated diabetics, and LPL PvuII genotypes in macroalbuminuric than normoalbuminuric patients. Fibrinogen heritability is quantile specific, which may explain or contribute to its gene-environment interactions. The analyses do not disprove the traditional gene-environment interpretations of these examples, rather quantile-dependent expressivity provides an alternative explanation that warrants consideration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗
DOE OSTI2021

Using machine learning to improve land use/cover characterization and projection for scenario-based global modeling

The characterization of the land surface in models is critical for robust estimation of the water cycle and associated extremes. Land use determines water demand, and land cover affects water availability. The interaction between water demand and availability is the key determinant of whether systems are resilient to variability in the overall water cycle. Furthermore, activities such as irrigation and deforestation or afforestation may influence regional to global precipitation patterns, depending on the extent of the activity. Thus this proposed work aims to address the following question: How does integrated land use/cover data and improved land use/cover projection, as informed by machine learning approaches, better resolve human-earth system resilience to water cycle variability and extremes?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving parameterization of ice microphysical processes in Arctic clouds using a synergistic modeling and observational approach (Final Report)

The primary goals of the project were to: 1) Improve the fundamental quantitative models of ice growth processes using ARM radar observations. 2) Better understand the growth of ice particle populations through vapor deposition, riming, and aggregation in Arctic clouds using a new particle-property microphysical model. 3) Robustly estimate and constrain key unknown model parameters using ARM radar observations combined with the particle-property model to develop an improved parameterization of ice growth processes for traditional bulk parameterizations used in cloud and climate models.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing the Behavior of Moderately Volatile Elements on the Moon in Order to Constrain Processes of Magmatic Evolution and Planetary Accretion

Heterogeneous isotopic compositions of Zn and Ga in lunar basalts and highlands rocks indicate that secondary processing affected the stable isotopic composition of major lunar reservoirs. It follows that these fractionating processes must be identified and characterized in order to make robust estimates for the isotopic composition of the mare basalt source regions or the bulk silicate Moon (BSM). Contrary to previous assertions (e.g., Kato & Moynier, 2017), MVEs such as Zn and Ga do not behave similarly on the Moon, in line with differences in their chemical behavior and volatility. There is clear evidence that impact reworking fractionates zinc isotopes but, in addition, zinc isotopes are also highly sensitive to the formation of Ti-bearing oxides such as ilmenite. In contrast, gallium is far less susceptible to volatile reworking, but its isotopic composition is fractionated during incorporation into plagioclase (Wimpenny et al., 2022; Render et al., 2023). Based on these findings we conclude that current isotopic compositions of Zn and Ga in lunar rocks, such as those measured for the mare basalt suite, are unlikely to be representative of the bulk Moon when it accreted. As such, estimates for the MVE isotopic composition of the BSM derived from mare basalts may not be accurate, which may have implications for the incorporation of these data into dynamical models of the Moon’s formation. In the following, we discuss the main findings from this project in more detail, focusing on the behaviour of zinc and gallium isotopic systematics in lunar rocks and the processes that control their isotopic fractionation on the Moon.

58 GEOSCIENCES↗
DOE OSTI2024

Peer-to-peer communication control for resilient operations of networked cyberphysical systems

This report includes two main accomplishments of the peer-to-peer communication control for resilient operation of networked microgrids project in FY24, which include a scheme for cyberattack-aware coordination of networked microgrids for supporting voltages of bulk power systems and a scheme for price signal-based operations of EV-rich networked microgrids with mixed ownership. First, the cyberattack-aware scheme enables networked microgrids to distributedly determine the amount of reactive power injection to support the voltage of bulk power system (BPS) in a fair manner. In this scheme, a risk-informed algorithm is presented to generate the peer-to- peer (P2P) communication graph with minimal risk of attack on communication links. To deal with cyberattacks on MG controllers, the resilient consensus algorithm (CA) is utilized for MG controllers to robustly estimate the total reactive power headroom, from which the MGs can accurately provide the needed amount of reactive power injection for supporting the voltage of BPS. The CA implementation and performance within the P2P communication framework are demonstrated on the IEEE 39-bus system with 6 microgrids contained in the distribution feeder under different cyberattack scenarios. Second, the price-based scheme enables the usage of the real-time price signal for the operations of electric vehicle (EV)-rich networked-microgrids with mixed ownership, in which not all the microgrids can communicate with the distribution system operator (DSO). In this scheme, a max consensus is introduced to enable the real-time price signal to be propagated from the DSO to all the microgrids, from which each microgrid controller will manage the DERs to balance the load demand and the power injection from the EV charging stations within its microgrid. Numerical results over one day with 288 slots of 5-minute intervals on the modified 123-node test feeder including 3 microgrids with high penetration of EV are presented to evaluate how the price signal affects the operations of networked microgrids under different charging strategies of the EV charging stations. The result indicates that our proposed EVCS (dis)charging strategy, which leverages the flexibility of EVs to support the grid through discharging during peak demand, proves to be a cost-effective solution that reduces operational costs while improving the social welfare of EV charging.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Untangling the Galaxy. II. Structure within 3 kpc

We present the results of the hierarchical clustering analysis of the Gaia DR2 data to search for clusters, comoving groups, and other stellar structures. The current paper builds on the sample from the previous work, extending it in distance from 1 to 3 kpc and increasing the number of identified structures up to 8292. To aid in the analysis of the population properties, we developed a neural network called Auriga to robustly estimate the age, extinction, and distance of a stellar group based on the input photometry and parallaxes of the individual members. We apply Auriga to derive the properties of not only the structures found in this paper, but also previously identified open clusters. Through this work, we examine the temporal structure of the spiral arms. Specifically, we find that the Sagittarius Arm has moved by >500 pc in the last 100 Myr and the Perseus Arm has been experiencing a relative lull in star formation activity over the last 25 Myr. We confirm the findings of the previous paper on the transient nature of the spiral arms, with the timescale of transition of a few 100 Myr. Finally, we find a peculiar ~1 Gyr old stream of stars that appears to be heliocentric. Its origin is unclear.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The z ∼ 1.03 Merging Cluster SPT-CL J0356–5337: New Strong Lensing Analysis with HST and MUSE

We present a strong lensing analysis and reconstruct the mass distribution of SPT-CL J0356−5337, a galaxy cluster at redshift z = 1.034. Our model supersedes previous models by making use of new multiband Hubble Space Telescope data and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. We identify two additional lensed galaxies to inform a more well-constrained model using 12 sets of multiple images in five separate lensed sources. The three previously known sources were spectroscopically confirmed by G. Mahler et al. at redshifts of z = 2.363, z = 2.364, and z = 3.048. We measured the spectroscopic redshifts of two of the newly discovered arcs using MUSE data, at z = 3.0205 and z = 5.3288. We increase the number of cluster member galaxies by a factor of 3 compared to previous work. We also report the detection of extended Lyα emission from several background galaxies. We measure the total projected mass density of the two major subcluster components, one dominated by the brightest cluster galaxy and the other by a compact group of luminous red galaxies. We find ${M}_{{\rm{B}}{\rm{C}}{\rm{G}}}(\lt 80\,\rm{kpc})=3.9{3}_{-0.14}^{+0.21}\times 1{0}^{13}$ M⊙ and ${M}_{{\rm{LRG}}}(\lt 80\,\rm{kpc})=2.9{2}_{-0.23}^{+0.16}\times 1{0}^{13}$ M⊙, yielding a mass ratio of $1.3{5}_{-0.08}^{+0.16}$ . The strong lensing constraints offer a robust estimate of the projected mass density regardless of modeling assumptions; allowing more substructure in this line of sight does not change the results or conclusions. Our results corroborate the conclusion that SPT-CL J0356−5337 is dominated by two mass components and is likely undergoing a major merger on the plane of the sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Retrieving microphysical properties of concurrent pristine ice and snow using polarimetric radar observations

Abstract. Ice and mixed-phase clouds play a key role in our climate system because of their strong controls on global precipitation and radiation budget. Their microphysical properties have been characterized commonly by polarimetric radar measurements. However, there remains a lack of robust estimates of microphysical properties of concurrent pristine ice and aggregates because larger snow aggregates often dominate the radar signal and mask contributions of smaller pristine ice crystals. This paper presents a new method that separates the scattering signals of pristine ice embedded in snow aggregates in scanning polarimetric radar observations and retrieves their respective abundances and sizes for the first time. This method, dubbed ENCORE-ice, is built on an iterative stochastic ensemble retrieval framework. It provides the number concentration, ice water content, and effective mean diameter of pristine ice and snow aggregates with uncertainty estimates. Evaluations against synthetic observations show that the overall retrieval biases in the combined total microphysical properties are within 5 % and that the errors with respect to the truth are well within the retrieval uncertainty. The partitioning between pristine ice and snow aggregates also agrees well with the truth. Additional evaluations against in situ cloud probe measurements from a recent campaign for a stratiform cloud system are promising. Our median retrievals have a bias of 98 % in the total ice number concentration and 44 % in the total ice water content. This performance is generally better than the retrieval from empirical relationships. The ability to separate signals of different ice species and to provide their quantitative microphysical properties will open up many research opportunities, such as secondary ice production studies and model evaluations for ice microphysical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sizing ice hydrometeor populations using the dual-wavelength radar ratio

Dual-wavelength (3.2 and 0.32 cm, i.e., X- and W-radar bands) radar ratio (DWR) measurements in ice clouds and precipitation using Canada's National Research Council Institute for Aerospace Research airborne radar are compared to closely collocated particle microphysical in situ sampling data in order to develop relations between DWR and characteristic hydrometeor size. This study uses the radar and in situ data sets collected during the In-Cloud ICing and Large-drop Experiment (ICICLE) campaign in midlatitude frontal clouds. Since atmospheric particle scattering at X band is predominantly in the Rayleigh regime and the W-band frequency is the highest frequency usually used for hydrometeor remote sensing, the X–W-band combination provides a relatively strong dual-wavelength reflectivity difference. This study considers radar and in situ measurements conducted in relatively homogeneous cloud and precipitation conditions. Measurements show that under these conditions, the difference between the X-band radar reflectivities observed with vertical and horizontal pointing of the radar beam are generally small and often negligible. However, W-band reflectivities at vertical beam pointing are, on average, larger than those for horizontal beam pointing by about 4 dB, which is a non-Rayleigh scattering effect from preferentially oriented non-spherical particles. A horizontal radar beam DWR–mean volume particle size relation, D v , provides robust estimates of this characteristic size for populations of particles with different habits. Uncertainties of D v retrievals using DWR are around 0.6 mm when D v is greater than approximately 1 mm. Size estimates using vertical radar beam DWRs have larger uncertainties due to smaller dual-wavelength signals and stronger influences of hydrometeor habits and orientations at this geometry of beam pointing. Mean relations among different characteristic sizes, which describe the entire particle size distribution (PSD), such as D v , and other sizes used in various applications (e.g., the mean, effective, and median sizes) are derived, so the results of this study can be used for estimating different PSD characteristic sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Planetary Boundary Layer Height (PBLH) over SGP from 1998 to 2023

PBLH is a critical parameter influencing weather phenomena, air quality, and various meteorological processes. However, accurately retrieving PBLH has been a challenging task due to limitations such as coarse temporal resolution and measurement drift in traditional radiosonde observations. To address these limitations, we have devised a lidar-based methodology that capitalizes on a newly developed algorithm for PBLH retrieval. This algorithm demonstrates enhanced capabilities in capturing diurnal fluctuations in PBLH compared to existing lidar-based methods (Su et al. 2020). In addition, we have refined this algorithm specifically for PBLH retrieval under cloudy conditions through a novel scheme (Su et al. 2022). To ensure data reliability, a quality-control process has been implemented to filter out questionable data points. Accompanying the data set is a quality-control flag for ease of reference. It should be noted that we have assimilated all available radiosonde observations to provide a more robust estimate of PBLH, making the data set valuable for a variety of related studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Robust Distributed State Estimator for Interconnected Transmission and Distribution Networks (Final Report RPPR-1)

This project’s objective is to develop a combined transmission and distribution state estimator which accounts for very large system size and model complexity (by way of distributing the computations) and large number of solar PV units connected to the distribution system on multiple feeders. The project not only provides a robust formulation and solution to this problem but also tests the solution by implementing it on a well-established large utility system. It introduces several improvements with respect to the state of the art in existing state estimation software: (a) The developed state estimator (SE) allows robust and accurate monitoring of bidirectional flows in distribution systems which result due to the distributed energy sources which are not observable and thus not incorporated in generation dispatch; (b) Large utility systems with tens of thousands of transmission buses and hundreds of thousands of distribution nodes are difficult to model as a single integrated system. This shortcoming is addressed by developing a “scalable distributed computational framework” which allows splitting the ultra large system models into several small subsystems and coordinating their solution by a robust and practical state estimation formulation; (c) Measurement errors irrespective of their locations are detected and removed by the developed state estimator. Historically, transmission and distribution systems were analyzed and operated as two independent systems. Given the non-transposed short feeder sections, unevenly loaded phases, strictly radial configuration and unidirectional power flows in the absence of remote generation, distribution system analysis was customized to account for these characteristics. However, some of these assumptions are no longer valid (non-radial configuration, bidirectional power flows) and thus distribution system analysis should be revisited. Furthermore, in the past, the interaction between the transmission and distribution systems was quite passive, where distribution substations were modeled as lumped loads in the transmission system model. With substantial generation injected by renewable generation located in the distribution systems, such modeling will no longer be accurate. The developed state estimator facilitates proper monitoring of the interactions between the transmission and distribution systems and enables smart dispatch of these units which are made observable by the state estimator.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗