Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “event study”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Single Entity Electrochemistry in Nanopore Electrode Arrays: Ion Transport Meets Electron Transfer in Confined Geometries

Electrochemical measurements conducted in confined volumes provide a powerful and direct means to address scientific questions at the nexus of nanoscience, biotechnology, and chemical analysis. How are electron transfer and ion transport coupled in confined volumes and how does understanding them require moving beyond macroscopic theories? Also, how do these coupled processes impact electrochemical detection and processing? We address these questions by studying a special type of confined-volume architecture, the nanopore electrode array, or NEA, which is designed to be commensurate in size with physical scaling lengths, such as the Debye length, a concordance that offers performance characteristics not available in larger scale structures. The experiments described here depend critically on carefully constructed nanoscale architectures that can usefully control molecular transport and electrochemical reactivity. We begin by considering the experimental constraints that guide the design and fabrication of zero-dimensional nanopore arrays with multiple embedded electrodes. These zero-dimensional structures are nearly ideal for exploring how permselectivity and unscreened ion migration can be combined to amplify signals and improve selectivity by enabling highly efficient redox cycling. Our studies also highlight the benefits of arrays, in that molecules escaping from a single nanopore are efficiently captured by neighboring pores and returned to the population of active redox species being measured, benefits that arise from coupling ion accumulation and migration. These tools for manipulating redox species are well-positioned to explore single molecule and single particle electron transfer events through spectroelectrochemistry, studies which are enabled by the electrochemical zero-mode waveguide (ZMW), a special hybrid nanophotonic/nanoelectronic architecture in which the lower ring electrode of an NEA nanopore functions both as a working electrode to initiate electron transfer reactions and as the optical cladding layer of a ZMW. Furthermore, while the work described here is largely exploratory and fundamental, we believe that the development of NEAs will enable important applications that emerge directly from the unique coupled transport and electron-transfer capabilities of NEAs, including in situ molecular separation and detection with external stimuli, redox-based electrochemical rectification in individually encapsulated nanopores, and coupled sorters and analyzers for nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Landslides in West Coast Metropolitan Areas: The Role of Extreme Weather Events

Rainfall-induced landslides represent a pervasive issue in areas where extreme rainfall intersects complex terrain. A farsighted management of landslide risk requires assessing how landslide hazard will change in coming decades and thus requires, inter alia, that we understand what rainfall events are most likely to trigger landslides and how global warming will affect the frequency of such weather events. We take advantage of 9 years of landslide occurrence data compiled by collating Google news reports and of a high-resolution satellite-based daily rainfall data to investigate what weather triggers landslide along the West Coast US. We show that, while this landslide compilation cannot provide consistent and widespread monitoring everywhere, it captures enough of the events in the major urban areas that it can be used to identify the relevant relationships between landslides and rainfall events in Puget Sound, the Bay Area, and greater Los Angeles. In all these regions, days that recorded landslides have rainfall distributions that are skewed away from dry and low-rainfall accumulations and towards heavy intensities. However, large daily accumulation is the main driver of enhanced hazard of landslides only in Puget Sound. There, landslide are often clustered in space and time and major events are primarily driven by synoptic scale variability, namely "atmospheric rivers" of high humidity air hitting anywhere along the West Coast, and the interaction of frontal system with the coastal orography. The relationship between landslide occurrences and daily rainfall is less robust in California, where antecedent precipitation (in the case of the Bay area) and the peak intensity of localized downpours at sub-daily time scales (in the case of Los Angeles) are key factors not captured by the same-day accumulations. Accordingly, we suggest that the assessment of future changes in landslide hazard for the entire the West Coast requires consideration of future changes in the occurrence and intensity of atmospheric rivers, in their duration and clustering, and in the occurrence of short-duration (sub-daily) extreme rainfall as well. Major regional landslide events, in which multiple occurrences are recorded in the catalog for the same day, are too rare to allow a statistical characterization of their triggering events, but a case study analysis indicates that a variety of synoptic-scale events can be involved, including not only atmospheric rivers but also broader cold- and warm-front precipitation. That a news-based catalog of landslides is accurate enough to allow the identification of different landslide/ rainfall relationships in the major urban areas along the US West Coast suggests that this technology can potentially be used for other English-language cities and could become an even more powerful tool if expanded to other languages and non-traditional news sources, such as social media.

landslides↗
DOE OSTI2022

Turbulent Flux Measurements of the Near‐Surface and Residual‐Layer Small Particle Events

Abstract According to recent field studies, almost half of the New Particle Formation (NPF) events occur aloft, in a residual layer, near the top of the boundary layer. Therefore, measurements of the meteorological parameters, precursor gas concentrations, and aerosol loadings conducted at the ground level are often not representative of the conditions where the NPFs take place. This paper presents new measurements obtained during the Turbulent Flux Measurements of the Residual Layer Nucleation Particles, conducted at the Southern Great Plains research site. Vertical turbulent fluxes of 3–10 nm‐sized particles were measured using a sonic anemometer and two condensation particle counters with nominal cutoff diameters of 3 nm and 10 nm mounted at the top of the 10‐m telescoping tower. Aerosol number size distribution (5–300 nm) was determined through the ground‐based Scanning Mobility Particle Sizers. The size selected (15–50 nm) particle hygroscopicity was derived with the Humidified Tandem Differential Mobility Analyzer. The ground level observations were supplemented by vertically‐resolved measurements of horizontal and vertical wind speeds and aerosol backscatter. The data analysis suggests that (a) turbulent flux measurements of 3–10 nm particles can distinguish between near‐surface and residual‐layer small particle events; (b) sub‐50 nm particles had a hygroscopicity value of 0.2, suggesting that organic compounds dominate atmospheric nanoparticle chemical composition at the site; and (c) current methodologies are inadequate for estimating the dry deposition velocity of sub‐10 nm particles because it is not feasible to measure particle concentration very near the surface, in the diffusion sublayer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Choosing Transport Events for Initiating Splitting and Rouletting

A study was performed to determine which transport events should be used to initiate a weight window lookup to achieve the best variance reduction performance. A weight window lookup potentially triggers particle splitting (in important regions of phase space) or rouletting (in unimportant regions), thereby optimizing computational effort. Potential initiating transport events include collisions (both pre- and post-collision), geometry surface crossings, traversing a mean-free path, and streaming across a weight window boundary. Permutations of these initiating events were tested on an urban model with background radiation sources and a spent fuel cask with a neutron dose mesh tally. Generally, all methods perform better with finer weight window meshes. Tracking on weight windows performs well for coarse weight window meshes, while a combination of splitting each mean-free path, geometric surface crossing, and before collisions performs well for fine weight window meshes.

42 ENGINEERING↗
DOE OSTI2025

Studies for the selection of fully contained muon neutrino CC events in the ICARUS T600 detector at Fermilab

The ICARUS T600 LAr-TPC detector has started its new physics runs in June 2022 at Fermilab within the SBN program. This detector is recording neutrino interactions from the Booster BNB neutrino beam in order to definitively clarify the open questions related to the possible existence of sterile neutrinos, as suggested by numerous observed experimental anomalies. In addition the neutrinos events recorded from the NuMI off-axis beam are under study in order to perform neutrino-Argon cross section measurements. At Fermilab, ICARUS is facing a challenging experimental condition: the detector is presently installed essentially at Earth surface where cosmic ray particles can become a serious source of background for the neutrino event search. This condition makes it necessary to deploy suitable automatic tools for the identification, selection, and measurement of the neutrino events among the millions of events triggered by cosmics. In this contribution, two automatic selection procedures devoted to the identification of fully contained muon neutrino interactions in the BNB neutrino beam will be presented, together with the first results of their application on simulated and on recently recorded events in the T600 detector.

Farnese, Christian [INFN, Padua] (ORCID:0000000163↗
DOE OSTI2025

Studies for the Selection of Fully Contained Muon Neutrino CC Events in the ICARUS T600 Detector at Fermilab

The ICARUS T600 LAr-TPC detector has started its new physics runs in June 2022 at Fermilab within the SBN program. This detector is recording neutrino interactions from the Booster BNB neutrino beam in order to definitively clarify the open questions related to the possible existence of sterile neutrinos, as suggested by numerous observed experimental anomalies. In addition the neutrinos events recorded from the NuMI off-axis beam are under study in order to perform neutrino-Argon cross section measurements. At Fermilab, ICARUS is facing a challenging experimental condition: the detector is presently installed essentially at Earth surface where cosmic ray particles can become a serious source of background for the neutrino event search. This condition makes it necessary to deploy suitable automatic tools for the identification, selection, and measurement of the neutrino events among the millions of events triggered by cosmics. In this contribution, two automatic selection procedures devoted to the identification of fully contained muon neutrino interactions in the BNB neutrino beam will be presented, together with the first results of their application on simulated and on recently recorded events in the T600 detector.

Farnese, Christian [INFN, Padua]↗
DOE OSTI2020

Analysis of the Frequency of Extreme Precipitation Events at Los Alamos National Laboratory

The Environmental Protection and Compliance (EPC) Storm Water Permitting and Compliance Program is currently updating its stormwater compliance standards for the Los Alamos National Laboratory (LANL) Engineering Standards Manual (2015) and is addressing the required current meteorological data needed for storm water applications. Information on the frequency of extreme precipitation events at LANL is needed. In response to this request, LANL Meteorology & Air Quality (MAQ) was tasked to estimate the frequency of extreme precipitation events. An earlier study by Foley (2012) which estimated extreme precipitation rates at LANL for a set of duration periods and average recurrence intervals (return periods) was evaluated. However, it was determined that this study would not be helpful since its technical support documentation was unavailable and the meteorological data and the methodology employed were not well understood. This led to the need to develop a new peer-reviewed methodology. This report documents the development of the new methodology and the application of a temporally representative precipitation dataset from February 1990 through December 2018, monitored at the official LANL meteorological station, Technical Area (TA)-6, to develop the frequency of extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Transform fault earthquakes in the North Atlantic: Source mechanisms and depth of faulting

The centroid depths and source mechanisms of 12 large earthquakes on transform faults of the northern Mid-Atlantic Ridge were determined from an inversion of long-period body waveforms. The earthquakes occurred on the Gibbs, Oceanographer, Hayes, Kane, 15 deg 20 min, and Vema transforms. The depth extent of faulting during each earthquake was estimated from the centroid depth and the fault width. The source mechanisms for all events in this study display the strike slip motion expected for transform fault earthquakes; slip vector azimuths agree to 2 to 3 deg of the local strike of the zone of active faulting. The only anomalies in mechanism were for two earthquakes near the western end of the Vema transform which occurred on significantly nonvertical fault planes. Secondary faulting, occurring either precursory to or near the end of the main episode of strike-slip rupture, was observed for 5 of the 12 earthquakes. For three events the secondary faulting was characterized by reverse motion on fault planes striking oblique to the trend of the transform. In all three cases, the site of secondary reverse faulting is near a compression jog in the current trace of the active transform fault zone. No evidence was found to support the conclusions of Engeln, Wiens, and Stein that oceanic transform faults in general are either hotter than expected from current thermal models or weaker than normal oceanic lithosphere.

Bergman, Eric A.↗
NASA NTRS1988

Transform fault earthquakes in the North Atlantic - Source mechanisms and depth of faulting

The centroid depths and source mechanisms of 12 large earthquakes on transform faults of the northern Mid-Atlantic Ridge were determined from an inversion of long-period body waveforms. The earthquakes occurred on the Gibbs, Oceanographer, Hayes, Kane, 15 deg 20 min, and Vema transforms. The depth extent of faulting during each earthquake was estimated from the centroid depth and the fault width. The source mechanisms for all events in this study display the strike slip motion expected for transform fault earthquakes; slip vector azimuths agree to 2 to 3 deg of the local strike of the zone of active faulting. The only anomalies in mechanism were for two earthquakes near the western end of the Vema transform which occurred on significantly nonvertical fault planes. Secondary faulting, occurring either precursory to or near the end of the main episode of strike-slip rupture, was observed for 5 of the 12 earthquakes. For three events the secondary faulting was characterized by reverse motion on fault planes striking oblique to the trend of the transform. In all three cases, the site of secondary reverse faulting is near a compression jog in the current trace of the active transform fault zone. No evidence was found to support the conclusions of Engeln, Wiens, and Stein that oceanic transform faults in general are either hotter than expected from current thermal models or weaker than normal oceanic lithosphere.

Bergman, Eric A.↗
NASA NTRS1994

Cosmic ray decreases and shock structure: A multispacecraft study

We examine greater than 60-MeV/amu ion data from three spacecraft (IMP 8 and Helios 1 and 2) at the time of a number of short term (less than 20-day duration) cosmic ray decreases (greater than 1 GeV) detected by ground-based neutron monitors in the years 1976 to 1979. The multispacecraft data allow us to investigate the structure of the modulation region and in particular the relative importance, as a function of location, of the shock and shock driver (ejecta) in causing the reduction in particle densities. Although the shocks contributing to cosmic ray decrease often have particle enhancements associated with them in the greater than 60-meV/amu data, this is not the case for three of the events discussed in this paper where a shock-associated decrease is also evident. Whereas the shock can cause an increase or decrease at low (i.e., less than neutron monitor) energies, the reduction of particle densities in the driver, if it is intercepted, is usually evident at all energies. Thus the overall shape of a decrease at greater than 60 MeV/amu depends primarily on whether the ejecta is intercepted. We find that the particle density inside ejecta increases with increasing radical distance from the Sun. In many of the events in this study, entry and exit of ejecta are accompanied by abrupt changes in the decrease and recovery rates which indicate that the effect of the ejecta is local. In contrast, the effect of the shock lasts many days after the shock has passed by and is evident at large angular distances from the longitude of the solar source, i.e., the effect of the shock is nonlocal. Within 1 AU there seems to be no radial dependence of the shock effect. One cosmic ray decrease seen at Earth, which had an unusual profile, can be understood if the median plane of the ejecta was inclined to the ecliptic.

Cane, H. V.↗
DOE OSTI2024

Sharp decline of dust events induces regional wetting over arid and semi-arid Northwest China in the NCAR Community atmosphere model

Abstract Multiple lines of observational evidence have indicated a significant wetting over the arid and semi-arid Northwest China (NWC) during recent decades, coinciding with a simultaneous sharp decline of dust events. Although recent studies have attributed NWC wetting to different anthropogenic and natural forcings, the mechanisms are not definitive and the regional wetting has been greatly underestimated in the Coupled Model Intercomparison Project historical simulations. Based on sensitivity experiments with different dust emission amounts using the NCAR Community Atmospheric Model version 5 (CAM5), here we find that decreasing dusts exert significant impacts on mixed-phase clouds through reducing the concentration of ice nucleating particles, increase the NWC precipitation and thus induce regional wetting through enhancing convection precipitation. A possible convection invigoration mechanism whereby the atmospheric vertical temperature gradient and convective instability are strengthened by reduced dusts, leading to convection invigoration and increased precipitation. These results are reinforced by simulations over the dust region in North Africa where mixed-phase and ice clouds are rare and reduced dusts do not increase precipitation. This study highlights the possible mechanism of dust-ice cloud interactions in recent NWC wetting and future regional climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Leveraging Data Intensive Computing to Support Automated Event Services

A large portion of Earth Science investigations is phenomenon- or event-based, such as the studies of Rossby waves, mesoscale convective systems, and tropical cyclones. However, except for a few high-impact phenomena, e.g. tropical cyclones, comprehensive records are absent for the occurrences or events of these phenomena. Phenomenon-based studies therefore often focus on a few prominent cases while the lesser ones are overlooked. Without an automated means to gather the events, comprehensive investigation of a phenomenon is at least time-consuming if not impossible. An Earth Science event (ES event) is defined here as an episode of an Earth Science phenomenon. A cumulus cloud, a thunderstorm shower, a rogue wave, a tornado, an earthquake, a tsunami, a hurricane, or an EI Nino, is each an episode of a named ES phenomenon," and, from the small and insignificant to the large and potent, all are examples of ES events. An ES event has a finite duration and an associated geolocation as a function of time; its therefore an entity in four-dimensional . (4D) spatiotemporal space. The interests of Earth scientists typically rivet on Earth Science phenomena with potential to cause massive economic disruption or loss of life, but broader scientific curiosity also drives the study of phenomena that pose no immediate danger. We generally gain understanding of a given phenomenon by observing and studying individual events - usually beginning by identifying the occurrences of these events. Once representative events are identified or found, we must locate associated observed or simulated data prior to commencing analysis and concerted studies of the phenomenon. Knowledge concerning the phenomenon can accumulate only after analysis has started. However, except for a few high-impact phenomena. such as tropical cyclones and tornadoes, finding events and locating associated data currently may take a prohibitive amount of time and effort on the part of an individual investigator. And even for these high-impact phenomena, the availability of comprehensive records is still only a recent development. A major reason for the lack of comprehensive ,records for the majority of the ES phenomena is the perception that they do not pose immediate and/or severe threat to life and property and are thus not consistently tracked. monitored, and catalogued. Many phenomena even lack commonly accepted criteria for definitions. However. the lack of comprehensive records is also due to the increasingly prohibitive volume of observations and model data that must be examined. NASA Earth Observing System Data Information System (EOSDIS) alone archives several petabytes (PB) of satellite remote sensing data and steadily increases. All of these factors contribute to the difficulty of methodically identifying events corresponding to a given phenomenon and significantly impede systematic investigations. In the following we present a couple motivating scenarios, demonstrating the issues faced by Earth scientists studying ES phenomena.

Clune, Thomas L.↗
NASA NTRS1999

A mechanical energy analysis of gait initiation

The analysis of gait initiation (the transient state between standing and walking) is an important diagnostic tool to study pathologic gait and to evaluate prosthetic devices. While past studies have quantified mechanical energy of the body during steady-state gait, to date no one has computed the mechanical energy of the body during gait initiation. In this study, gait initiation in seven normal male subjects was studied using a mechanical energy analysis to compute total body energy. The data showed three separate states: quiet standing, gait initiation, and steady-state gait. During gait initiation, the trends in the energy data for the individual segments were similar to those seen during steady-state gait (and in Winter DA, Quanbury AO, Reimer GD. Analysis of instantaneous energy of normal gait. J Biochem 1976;9:253-257), but diminished in amplitude. However, these amplitudes increased to those seen in steady-state during the gait initiation event (GIE), with the greatest increase occurring in the second step due to the push-off of the foundation leg. The baseline level of mechanical energy was due to the potential energy of the individual segments, while the cyclic nature of the data was indicative of the kinetic energy of the particular leg in swing phase during that step. The data presented showed differences in energy trends during gait initiation from those of steady state, thereby demonstrating the importance of this event in the study of locomotion.

Leg/physiology↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

Data Albums: An Event Driven Search, Aggregation and Curation Tool for Earth Science

One of the largest continuing challenges in any Earth science investigation is the discovery and access of useful science content from the increasingly large volumes of Earth science data and related information available. Approaches used in Earth science research such as case study analysis and climatology studies involve gathering discovering and gathering diverse data sets and information to support the research goals. Research based on case studies involves a detailed description of specific weather events using data from different sources, to characterize physical processes in play for a specific event. Climatology-based research tends to focus on the representativeness of a given event, by studying the characteristics and distribution of a large number of events. This allows researchers to generalize characteristics such as spatio-temporal distribution, intensity, annual cycle, duration, etc. To gather relevant data and information for case studies and climatology analysis is both tedious and time consuming. Current Earth science data systems are designed with the assumption that researchers access data primarily by instrument or geophysical parameter. Those who know exactly the datasets of interest can obtain the specific files they need using these systems. However, in cases where researchers are interested in studying a significant event, they have to manually assemble a variety of datasets relevant to it by searching the different distributed data systems. In these cases, a search process needs to be organized around the event rather than observing instruments. In addition, the existing data systems assume users have sufficient knowledge regarding the domain vocabulary to be able to effectively utilize their catalogs. These systems do not support new or interdisciplinary researchers who may be unfamiliar with the domain terminology. This paper presents a specialized search, aggregation and curation tool for Earth science to address these existing challenges. The search tool automatically creates curated "Data Albums", aggregated collections of information related to a specific science topic or event, containing links to relevant data files (granules) from different instruments; tools and services for visualization and analysis; and information about the event contained in news reports, images or videos to supplement research analysis. Curation in the tool is driven via an ontology based relevancy ranking algorithm to filter out non-relevant information and data.

Ramachandran, Rahul↗
DOE OSTI2023

Adaptive cold-load pickup considerations in 2-stage microgrid unit commitment for enhancing microgrid resilience

In an extended main grid outage spanning multiple days, load shedding serves as a critical mechanism for islanded microgrids to maintain essential power and energy reserves that are indispensable for fulfilling reliability and resiliency mandates. However, using load shedding for such purposes leads to increasing occurrence of cold load pickup (CLPU) events. Here, this study presents an innovative adaptive CLPU model that introduces a method for determining and incorporating parameters related to CLPU power and energy requirements into a two-stage microgrid unit commitment (MGUC) algorithm. In contrast to the traditional fixed-CLPU-curve approach, this model calculates CLPU duration, power, and energy demands by considering outage durations and ambient temperature variations within the MGUC process. By integrating the adaptive CLPU model into the MGUC problem formulation, it allows for the optimal allocation of energy resources throughout the entire scheduling horizon to fulfill the CLPU requirements when scheduling multiple CLPU events. The performance of the enhanced MGUC algorithm considering CLPU needs is assessed using actual load and photovoltaic (PV) data. Simulation results demonstrate significant improvements in dispatch optimality evaluated by the amount of load served, customer comfort, energy storage operation, and adherence to energy schedules. These enhancements collectively contribute to reliable and resilient microgrid operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1985

Comparison of big event with calculations of the air shower development

The incidence of high energy hadrons and electron-photons in air showers at various stages of development is calculated. Numerical calculation is used to solve the diffusion equation for a nuclear cascade and analytical calculation for cascade shower induced gamma rays. From these calculations, one can get the longitudinal development of the high energy hadron and electron-photon components, and the energy spectra of these components at various depths of air shower development. The total number of hadrons (N sub H) and electron-photon components (N sub gamma) are related according to stages of the air shower development and primary energy. The relation of the total energy of hadron and electron-photon component above the threshold energy is given. The energy balance between both components is also a useful parameter to study high energy events accompanying air showers. The relation of N sub H and fractional hadronic energy E (sum E sub H sup gamma/sum E sub H sup gamma + Sum E sub gamma) is calculated. This relation is helpful to understand the stage of air shower development(t) and primary energy (E sub p).

Niwa, M.↗