Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “event data models”

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 253 records · Page 14

Modeling the Contribution of Precipitation Loss to a Radiation Belt Electron Dropout Observed by Van Allen Probes

A drift-diffusion model is used to simulate the low-altitude electron distribution, accounting for azimuthal drift, pitch angle diffusion, and atmospheric backscattering effects during a rapid electron dropout event on 21st August 2013, at L = 4.5. Additional external loss effects are introduced during times when the low-altitude electron distribution cannot be reproduced by diffusion alone. Here, the model utilizes low-altitude electron count rate data from five POES/MetOp satellites to quantify pitch angle diffusion rates. Low-altitude data provides critical constraint on the model because it includes the drift loss cone region where the electron distribution in longitude is highly dependent on the balance between azimuthal drift and pitch angle diffusion. Furthermore, a newly derived angular response function for the detectors onboard POES/MetOp is employed to accurately incorporate the bounce loss cone measurements, which have been previously contaminated by electrons from outside the nominal field-of-view. While constrained by low-altitude data, the model also shows reasonable agreement with high-altitude data. Pitch angle diffusion rates during the event are quantified and are faster at lower energies. Precipitation is determined to account for all of the total loss observed for 450 keV electrons, 88% for 600 keV and 38% for 900 keV. Predictions made in the MeV range are deemed unreliable as the integral energy channels E3 and P6 fail to provide the necessary constraint at relativistic energies.

79 ASTRONOMY AND ASTROPHYSICS↗

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

A Convolution Neural Network for Voltage Event Classification at a Photovoltaic Inverter

This paper presents a convolutional neural network (CNN) developed to identify voltage events in photovoltaic (PV) inverters. The CNN is trained on synthetic data generated using the IEEE 13-bus distribution feeder model and evaluated on field measured data collected from Energy Northwest’s Horn Rapids Solar, Storage, and Training (HRSST) facility. The study focuses on two common voltage events: faults and voltage sags. The CNN is configured to analyze voltage and current waveforms from three-phase PV systems, demonstrating excellent accuracy during training. Field data from the HRSST facility is employed to assess its real-world performance, where the CNN achieves perfect identification of faults and voltage sags in a sample of nine events. This work highlights the potential of the proposed method to enhance PV protection schemes, providing a robust foundation for improved voltage event detection and grid reliability.

Cornachione, Matthew A.↗

Multi-scale Multi-physics Scientific Machine Learning for Water Cycle Extreme Events Identification, Labelling, Representation, and Characterization

Impacts of climate are usually felt through extreme events such as droughts, floods, thunderstorms, windstorms, wildfires, and so on, that are intimately tied to the water cycle. Predicting the frequency and severity of extreme events under climate change remains a significant challenge; meanwhile, the mechanisms and impacts of these extremes are far from well understood. There are several major science challenges: (1) Lack of labelled extreme events data and missing standards in defining extremes; (2) Computational demand of high-resolution ensemble climate modeling; (3) Modeling the multiscale multi-physics hierarchical structure of compound extremes; (4) Lack of understanding of mechanisms of extreme events; (5) Large uncertainty in extreme events impacts on infrastructure; (6) Subjective assessment of weather-related risk from seasonal to multi-decadal time scales and lack of metrics for risk assessment and mitigation control.

54 ENVIRONMENTAL SCIENCES↗

DEMOS (Demographic Microsimulator Tool for Longitudinal Synthetic Population) [SWR-25-135] related to NLR SWR-26-076

The Demographic Microsimulator (DEMOS) is an agent-based simulation framework used to model the evolution of population demographic characteristics and lifecycle events, such as education attainment, marital status, and other key transitions. DEMOS modules are designed to capture the interdependencies between short-term and long-term lifecycle events, which are often influential in downstream transportation and land-use modeling. A key feature of DEMOS is its ability to track changes in an agent’s demographic status from year t to year t + 1. This structure allows the model to evolve populations over any user-defined time horizon. As a result, DEMOS is well suited for analyzing medium- and long-term transportation-related decisions, including household vehicle transactions (e.g., purchasing, selling, or replacing vehicles) and work location choices. Core features of DEMOS include the modeling of more than ten lifecycle events, behaviorally realistic patterns informed by long-running panel data, explicit representation of interdependencies among lifecycle processes, and a flexible, modular simulation architecture. A technical memorandum describing DEMOS is available here. The memorandum provides an overview of the framework’s functionality, model structure, input and output data, and its applications in transportation planning and broader policy analysis contexts. Interested readers are also encouraged to consult the paper listed below for additional details on the DEMOS methodology. Sun, Bingrong, Shivam Sharda, Venu M. Garikapati, Mohamed Amine Bouzaghrane, Juan Caicedo, Srinath Ravulaparthy, Isabel Viegas de Lima, Ling Jin, C. Anna Spurlock, and Paul Waddell. "Demographic Microsimulator for Integrated Urban Systems: Adapting Panel Survey of Income Dynamics to Capture the Continuum of Life." Transportation Research Record (2025): 03611981251333339.

Sun, Bingrong [National Laboratory of the Rockies ↗

The 2024 “Hacking Limnology” Workshop Series and Virtual Summit: Increasing Inclusion, Participation, and Representation in the Aquatic Sciences

The 4th Aquatic Ecosystem MOdeling Network—Junior (AEMON-J) Hacking Limnology Workshop and 5th Virtual Summit: Incorporating Data Science and Open Science in the Aquatic Sciences (DSOS) convened 15–19 July 2024. During the week, these joint communities engaged in activities at the intersection of big data, open science, modeling, remote sensing, and the aquatic sciences. The weeklong event, with over 100 aquatic science practitioners and enthusiasts, followed a similar structure to previous years, comprising three days of workshops followed by two days of the virtual summit.

54 ENVIRONMENTAL SCIENCES↗

Deciphering supramolecular and polymer-like behavior in metallogels: real-time insights into temperature-modulated gelation and rapid self-assembly dynamics

Bis(pyridyl) urea-based gelators, namely L2 and its isomeric mixture ( L1 + L2 ), are known to self-assemble into 1D architectures capable of inducing supramolecular gelation. Coordination with metal ions such as Ag( I ), Cu( II ), and Fe( III ) introduces structural reinforcement, enabling the formation of distinct 3D networks governed by metal-specific coordination geometries. Here, we present a comprehensive investigation into the temperature-responsive behavior (20–60 °C) of L2 and L1 + L2 , both in the absence and presence of Ag( I ), Dy( III ), Fe( III ), Cu( II ), and Ho( III ), using real-time small-angle neutron scattering (SANS). To probe long-term structural evolution/kinetics of self-assembly, real-time small-angle X-ray scattering (SAXS) was employed on L2 + Ag gels, complemented by differential scanning calorimetry (DSC) to evaluate thermal transitions. Our results reveal strikingly divergent gelation behaviors: L2 forms a highly rigid, covalent polymer-like network, while L1 + L2 exhibits remarkable thermal adaptability. Upon metal coordination, the assemblies exhibit pronounced crystallinity and exceptional thermal stability, as evidenced by persistent Bragg reflections and invariant d-spacings. Intriguingly, L2 : Fe (2 : 1) and L1 : L2 : Fe (0.5 : 0.5 : 1) in acetonitrile-d 3 (ACN-d 3 ) deviate from this trend, forming thermally labile amorphous gels. These systems show a complete loss of crystalline order, reduced Porod exponents—indicative of collapsed or branched fiber morphologies—and prominent melting and glass transition events in DSC. Fitting SANS and SAXS data to the correlation length model unveiled insightful nanostructural features. While most systems displayed minimal temperature-induced variation in mesh size or surface morphology, L2 : Ag in dimethyl sulfoxide-d 6 (DMSO-d 6 )/D 2 O and L2 : Fe (1 : 1) in ACN-d 3 exhibited a rare combination of thermally stable correlation lengths and increasing high- q exponents—strongly suggesting progressive fiber densification or surface smoothing within a robust gel framework. These findings highlight the tunability and structural resilience of supramolecular gels through precise control of ligand architecture, metal coordination, and temperature, offering valuable design principles for functional soft materials.

Pajoubpong, Jinnipha [Univ. of Cincinnati, OH (Uni↗

Expected ${\boldsymbol H \xrightarrow{} {\boldsymbol{\mu^+\mu^-}}} $ measurement precision with ${\boldsymbol { e^{+}e^{-}}{\bf \xrightarrow{}}{ \boldsymbol{Z(q\bar{q})H}}}$ production at the CEPC *

Here, a search for the dimuon decay of the Standard Model Higgs boson is performed using Monte Carlo simulated events to mimic data corresponding to an integrated luminosity of 5.6 ab -1 collected with the Circular Electron-Positron Collider detector in e + e - collisions at $\sqrt{s}$ =240 GeV. This study investigates the e + e - →ZH, Z → $q\bar{q}$, H → μ + μ - process, and the expected significance considering only the statistical uncertainty in the data for a background-only hypothesis for a Higgs boson with a mass of 125 GeV is found to be 6.1 σ, corresponding to a precision of 19%. The systematic impacts from the background Monte Carlo statistical fluctuations are estimated to be negligible. Moreover, the dependence of the measurement accuracy on the muon momentum resolution of the CEPC detector is investigated. It is found that the muon momentum resolution must be better than 204 MeV to discover the H →μμ process at the nominal integrated luminosity. If the resolution is 100% worse than the designed parameter, the integrated luminosity must be greater than 7.2 ab -1 to reach 5σ significance.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data-Based Resilience Enhancement Strategies for Electric-Gas Systems Against Sequential Extreme Weather Events

Some extreme weather events, such as the hurricane, pass through an area sequentially and thus are called sequential extreme weather events (SEWEs). This paper proposes a data-based robust optimization (RO) model to enhance the resilience of the integrated electricity and gas system (IEGS) against SEWEs. Specifically, the SEWE strikes the IEGS sequentially. After each attack, the system state is adjusted immediately to minimize the maximized expected system cost caused by the SEWE. The attack-defense procedures are repeated alternatively during the SEWE. Preventive measures, hardening, are made in advance to reduce the impact of sequential attacks. The entire process is formulated as a multi-period RO model. Furthermore, it is proved that the most effective resilience enhancement strategies for this model are the same as those for a two-stage RO model, which can be solved by the nested column-and-constraint generation (C&CG) algorithm. In addition, the property of SEWEs, sequentially endangering limited regions of the IEGS, is incorporated to build a data-based uncertainty set and reduce its conservativeness. Simulation results on two IEGSs validate the effectiveness of the proposed model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Factors affecting powerhouse passage of spring migrant smolts at federally operated hydroelectric dams of the Snake and Columbia rivers

From 2008 to 2018, acoustic telemetry studies were conducted to evaluate dam passage survival of spring migrant Chinook salmon and steelhead smolts at seven of the eight federally operated dams on the lower Snake and Columbia rivers. Data from over 87 000 dam passage events were evaluated using regression modeling to identify the effect of spill operations, environmental conditions, and fish characteristics on powerhouse passage probability. In general, powerhouse passage was positively correlated with discharge, negatively correlated with forebay temperature and fish size, and higher for fish that passed the dam at night and for those that approached from the powerhouse side of the river, suggesting powerhouse passage is largely a function of smolt activity level and swimming ability. As such, spilling large volumes of water to reduce powerhouse passage is likely to be most effective during times of reduced activity and swimming ability (e.g., at night, high flows, and cold temperatures). This information can be used to develop dam- and time-specific spill operations that optimize smolt passage, power generation, and other competing demands, such as adult passage.

60 APPLIED LIFE SCIENCES↗

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗

Study of CC2p0$\pi$ Event Topologies in the MicroBooNE Detector

The precise measurement of neutrino oscillation parameters is a critical component of understanding neutrino oscillations. As future accelerator neutrino experiments will utilize 40 Ar as the target nucleus it is important to understand the influence that nuclear effects have on neutrino-nucleon cross-sections. Processes that lead to 2 particle-2 hole (2p2h) states are of particular interest, such as Meson Exchange Currents (MEC) and Short Range Nucleon-Nucleon Correlations. In the MicroBooNE detector, a charged-current interaction that produced a 2p2h state would create a 3 track topology with 1 muon and 2 protons connected to a single interaction vertex (CC2p0$\pi$). We present here a study of MicroBooNE’s potential sensitivity to a variety of MEC models as well as CC2p0$\pi$ events selected from a combination of data collected from the first 3 years of running of the MicroBooNE experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

AGN-201 Digital Twin Concept of Operations

The AGN-201 is a nuclear reactor at Idaho State University (ISU) and is currently being used in the development of a digital twin (DT) with Idaho National Laboratory (INL). The goal is to create a DT (called the AGN-201 DT) which can monitor an operating nuclear reactor to determine when the reactor is being operated normally; normal operations are any operations that is declared by the ISU staff. To accomplish this, researchers at INL developed a DT ecosystem which can ingest data from the AGN-201 reactor. This data is fed into a series of machine learning (ML) and reactor physics models. The ML and reactor physics models assess the data and determine if the reactor is operating normally. Event(s) that are flagged as anomalies are investigated by the INL staff to determine if any undeclared experiments were conducted. The AGN-201 DT was verified in July/August 2023, when the ISU staff performed undeclared experiments and the INL staff were able to assess the anomalous events and determine what likely caused each event. This project was a steppingstone to promote the use of DTs for international safeguards and marks the first time a DT was able to monitor a nuclear reactor for this purpose. This document provides the concept of operations (CONOPS) for the AGN-201 digital twin (DT). The CONOPS describes how the AGN-201 and AGN-201 DT are expected to operate and details how each are established.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

BROOD: Bilevel and Robust Optimization and Outlier Detection for Efficient Tuning of High-Energy Physics Event Generators

The parameters in Monte Carlo (MC) event generators are tuned on experimental measurements by evaluating the goodness of fit between the data and the MC predictions. The relative importance of each measurement is adjusted manually in an often time-consuming, iterative process to meet different experimental needs. In this work, we introduce several optimization formulations and algorithms with new decision criteria for streamlining and automating this process. These algorithms are designed for two formulations: bilevel optimization and robust optimization. Both formulations are applied to the datasets used in the ATLAS A14 tune and to the dedicated hadronization datasets generated by the SHERPA generator, respectively. The corresponding tuned generator parameters are compared using three metrics. We compare the quality of our automatic tunes to the published ATLAS A14 tune. Moreover, we analyze the impact of a pre-processing step that excludes data that cannot be described by the physics models used in the MC event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗

Search for stealth supersymmetry in final states with two photons, jets, and low missing transverse momentum in proton-proton collisions at s = 13 TeV

The results of a search for stealth supersymmetry in final states with two photons and jets, targeting a phase space region with low missing transverse momentum ( p T miss ), are reported. The study is based on a sample of proton-proton collisions at s = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb - 1 . As LHC results continue to constrain the parameter space of the minimal supersymmetric standard model, the low p T miss regime is increasingly valuable to explore. To estimate the backgrounds due to standard model processes in such events, we apply corrections derived from simulation to an estimate based on a control selection in data. The results are interpreted in the context of simplified stealth supersymmetry models with gluino and squark pair production. The observed data are consistent with the standard model predictions, and gluino (squark) masses of up to 2150 (1850) GeV are excluded at the 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First measurement of the forward rapidity gap distribution in pPb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV

For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_{NN}}$ = 8.16 TeV , almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of 2 below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC , QGSJET II , and HIJING predictions are all at least a factor of 5 lower than the data. The latter effect might be explained by a significant contribution of ultraperipheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗