Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 469 records · Page 26

Determining the kinetic rate constants of Fe 3 O 4 -to-Fe and FeO-to-Fe reduction by H 2

Steel production using coal accounts for ~ 8% of global carbon emissions. “Green Steel” is a new grand concept proposed recently to make steel from iron ores using renewable derived “Green Hydrogen” to achieve zero carbon emission. The kinetics and rate-limiting steps of iron ore reduction into iron with H 2 as a reducing agent is critically important to the success of this new technology. While reduction of Fe 2 O 3 into Fe by H 2 follows multiple steps, the past research on this topic mainly deals with the overall averaged kinetics, giving little information on the elemental and rate-limiting steps. In this work, we report a kinetic study specifically design to attain kinetic rate constants of one-step reduction of Fe 3 O 4 -to-Fe and FeO-to-Fe. Guided by thermodynamics, we show first how to create in situ the desirable starting oxide phases, i.e., Fe 3 O 4 and FeO, with precisely controlled the ratio of partial pressures of H 2 O and H 2 . We then show time-dependent raw H 2 O content data collected by a mass spectrometer and the processed reduction data to extract kinetic rate constants. We found that the kinetics of the two one-step reduction reactions follows nicely the Johnson-Mehl-Avrami (JMA) phase transformation model. The one-step reduction mechanisms and activation energy are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simplifying activations with linear approximations in neural networks

A key step in Neural Networks is activation. Among the different types of activation functions, sigmoid, tanh, and others involve the usage of exponents for calculation. From a hardware perspective, exponential implementation implies the usage of Taylor series or repeated methods involving many addition, multiplication, and division steps, and as a result are power-hungry and consume many clock cycles. We implement a piecewise linear approximation of the sigmoid function as a replacement for standard sigmoid activation libraries. This approach provides a practical alternative by leveraging piecewise segmentation, which simplifies hardware implementation and improves computational efficiency. In this paper, we detail piecewise functions that can be implemented using linear approximations and their implications for overall model accuracy and performance gain. Our results show that for the DenseNet, ResNet, and GoogLeNet architectures, the piecewise linear approximation of the sigmoid function provides faster execution times compared to the standard TensorFlow sigmoid implementation while maintaining comparable accuracy. Specifically, for MNIST with DenseNet, accuracy reaches 99.91% (Piecewise) vs. 99.97% (Base) with up to 1.31x speedup in execution time. For CIFAR-10 with DenseNet, accuracy improves to 98.97% (Piecewise) vs. 99.40% (Base) while achieving 1.24x faster execution. Similarly, for CIFAR-100 with DenseNet, the accuracy is 97.93% (Piecewise) vs. 98.39% (Base), with a 1.18x execution time reduction. These results confirm the proposed method’s capability to efficiently process large-scale datasets and computationally demanding tasks, offering a practical means to accelerate deep learning models, including LSTMs, without compromising accuracy.

Activation function↗

STRIDES: automated uniform models for 30 quadruply imaged quasars

ABSTRACT Gravitational time delays provide a powerful one-step measurement of H0, independent of all other probes. One key ingredient in time-delay cosmography are high-accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (105–106 CPU hours and ∼0.5–1 yr, respectively). Major improvements in modelling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, we develop an automated modelling pipeline and apply it to a sample of 31 lens systems, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in the Fermat potential (proportional to time delay) with respect to modelling choices. We find that for 10/30 lenses, our models are cosmography or nearly cosmography grade (<3 per cent and 3–5 per cent variations). For 6/30 lenses, the models are close to cosmography grade (5–10 per cent). These results utilize informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modelling sequence and options. In conclusion, we show that uniform cosmography grade modelling of large strong lens samples is within reach.

79 ASTRONOMY AND ASTROPHYSICS↗

Surface elevation, snow depth, vegetation height, and color imagery from multiple UAV surveys from 2018 to 2023 across a watershed near Teller road mile marker 27, Seward Peninsula, Alaska

Multiple photogrammetric surveys using an unoccupied aerial system (UAS) were performed across a watershed along Teller road on the Seward Peninsula (Alaska) to investigate terrain, vegetation and snowpack properties. This work was led by the environmental geophysics team from Lawrence Berkeley National Laboratory as part of the Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic). The multiple photogrammetric surveys used mostly a DJI Matrice 210 UAS with a DJI X5S camera, Ground Control Points (GCPs) surveyed with a Real Time Kinematic (RTK)-GPS, and a structure from motion (SfM) technique for photogrammetric reconstruction. Digital Surface elevation Models (DSMs) were inferred at various times of year, including near peak in plant and leaf density (July 19 2017), close to peak in snow depth (April 1 2019 and April 11 2022), and at low plant and leaf density of tall after the first bare-ground date (June 9 2019). Snow and canopy height were obtained by subtracting DSMs from a Digital Terrain Model (DTM) proxy inferred from the June 9 2019 survey. The derived products include ortho-mosaiced RGB map, DSMs, a DTM proxy, snow depth (April 2019 and 2022) and canopy height (July 2017). Unprocessed and processed data products are included in this package (processing levels 0-2). Data and metadata are provided as text (*.csv), image (*.jpg), and GeoTiff (*.tif) formats. This metadata document contains a description of the survey and processing steps and the inferred products.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Common Analysis Tool Being Developed for Aeropropulsion: The National Cycle Program Within the Numerical Propulsion System Simulation Environment

The NASA Lewis Research Center is developing an environment for analyzing and designing aircraft engines-the Numerical Propulsion System Simulation (NPSS). NPSS will integrate multiple disciplines, such as aerodynamics, structure, and heat transfer, and will make use of numerical "zooming" on component codes. Zooming is the coupling of analyses at various levels of detail. NPSS uses the latest computing and communication technologies to capture complex physical processes in a timely, cost-effective manner. The vision of NPSS is to create a "numerical test cell" enabling full engine simulations overnight on cost-effective computing platforms. Through the NASA/Industry Cooperative Effort agreement, NASA Lewis and industry partners are developing a new engine simulation called the National Cycle Program (NCP). NCP, which is the first step toward NPSS and is its initial framework, supports the aerothermodynamic system simulation process for the full life cycle of an engine. U.S. aircraft and airframe companies recognize NCP as the future industry standard common analysis tool for aeropropulsion system modeling. The estimated potential payoff for NCP is a $50 million/yr savings to industry through improved engineering productivity.

Follen, Gregory J.↗

Extending substructure based iterative solvers to multiple load and repeated analyses

Direct solvers currently dominate commercial finite element structural software, but do not scale well in the fine granularity regime targeted by emerging parallel processors. Substructure based iterative solvers--often called also domain decomposition algorithms--lend themselves better to parallel processing, but must overcome several obstacles before earning their place in general purpose structural analysis programs. One such obstacle is the solution of systems with many or repeated right hand sides. Such systems arise, for example, in multiple load static analyses and in implicit linear dynamics computations. Direct solvers are well-suited for these problems because after the system matrix has been factored, the multiple or repeated solutions can be obtained through relatively inexpensive forward and backward substitutions. On the other hand, iterative solvers in general are ill-suited for these problems because they often must restart from scratch for every different right hand side. In this paper, we present a methodology for extending the range of applications of domain decomposition methods to problems with multiple or repeated right hand sides. Basically, we formulate the overall problem as a series of minimization problems over K-orthogonal and supplementary subspaces, and tailor the preconditioned conjugate gradient algorithm to solve them efficiently. The resulting solution method is scalable, whereas direct factorization schemes and forward and backward substitution algorithms are not. We illustrate the proposed methodology with the solution of static and dynamic structural problems, and highlight its potential to outperform forward and backward substitutions on parallel computers. As an example, we show that for a linear structural dynamics problem with 11640 degrees of freedom, every time-step beyond time-step 15 is solved in a single iteration and consumes 1.0 second on a 32 processor iPSC-860 system; for the same problem and the same parallel processor, a pair of forward/backward substitutions at each step consumes 15.0 seconds.

Farhat, Charbel↗

Ghent Generating Station Nuclear Study: Siting Evaluation

The purpose of this initial siting evaluation is to assess the suitability of the PPL Corporation (PPL) owned land at/near the Kentucky Utilities Ghent Generating Station (GGS) for the deployment of advanced nuclear reactors and to determine if there are any exclusionary or avoidance factors associated with siting nuclear power at GGS. The formal siting process for a nuclear reactor requires a great level of time (i.e., multiple years), effort, and detail. The Electric Power Research Institute (EPRI) estimates the formal siting process could take a dedicated team between one and three years to complete. This initial siting evaluation considers several characteristics (e.g., environmental conditions, seismic concerns, footprint, water use, etc.) and highlights favorable/preferred characteristics as well as potential risks to feasibility. Potential risks were explored in greater detail to provide PPL with additional insights to support decision planning. The results indicate the strengths and weaknesses associated with the GGS site as a location for a nuclear generating station and are intended to inform PPL’s next steps. The results of this initial siting evaluation will also inform the selection of candidate nuclear technologies for the GGS site that align with PPL’s mission and business objectives.

99 GENERAL AND MISCELLANEOUS↗

Rapid Corner Detection Using FPGAs

In order to perform precision landings for space missions, a control system must be accurate to within ten meters. Feature detection applied against images taken during descent and correlated against the provided base image is computationally expensive and requires tens of seconds of processing time to do just one image while the goal is to process multiple images per second. To solve this problem, this algorithm takes that processing load from the central processing unit (CPU) and gives it to a reconfigurable field programmable gate array (FPGA), which is able to compute data in parallel at very high clock speeds. The workload of the processor then becomes simpler; to read an image from a camera, it is transferred into the FPGA, and the results are read back from the FPGA. The Harris Corner Detector uses the determinant and trace to find a corner score, with each step of the computation occurring on independent clock cycles. Essentially, the image is converted into an x and y derivative map. Once three lines of pixel information have been queued up, valid pixel derivatives are clocked into the product and averaging phase of the pipeline. Each x and y derivative is squared against itself, as well as the product of the ix and iy derivative, and each value is stored in a WxN size buffer, where W represents the size of the integration window and N is the width of the image. In this particular case, a window size of 5 was chosen, and the image is 640 480. Over a WxN size window, an equidistance Gaussian is applied (to bring out the stronger corners), and then each value in the entire window is summed and stored. The required components of the equation are in place, and it is just a matter of taking the determinant and trace. It should be noted that the trace is being weighted by a constant k, a value that is found empirically to be within 0.04 to 0.15 (and in this implementation is 0.05). The constant k determines the number of corners available to be compared against a threshold sigma to mark a valid corner. After a fixed delay from when the first pixel is clocked in (to fill the pipeline), a score is achieved after each successive clock. This score corresponds with an (x,y) location within the image. If the score is higher than the predetermined threshold sigma, then a flag is set high and the location is recorded.

Morfopoulos, Arin C.↗

A terminology for scientific workflow systems

The term “scientific workflow” has evolved over the last two decades to encompass a broad range of compositions of interdependent compute tasks and data movements. It has also become an umbrella term for processing in modern scientific applications. Today, many scientific applications can be considered as workflows made of multiple dependent steps, and hundreds of workflow systems have been developed to manage and run these scientific workflows. However, no turnkey solution has emerged from the field to address the diversity of scientific processes and the infrastructure on which they are supposed to be implemented. Instead, new research problems requiring the execution of scientific workflows with some novel feature often lead to the development of an entirely new workflow system. A direct consequence of this situation is that many existing workflow management systems (WMSs) share some salient features, offer similar functionalities, and can manage the same categories of workflows but at the same time also have some distinct capabilities that can be important for specific applications. This situation makes researchers who develop workflows face the complex question of selecting a WMS. This selection can be driven by technical considerations, to find the system that is the most appropriate for their application and for the computing and storage resources available to them, or other factors such as reputation, adoption, strong community support, or long-term sustainability. To address this problem, a group of WMS developers and practitioners joined their efforts to produce a community-based terminology of WMSs. This paper summarizes their findings and introduces this new terminology to characterize WMSs. Furthermore, this terminology is composed of fives axes: workflow structure and characteristics, composition, orchestration, data management, and metadata capture. Each axis comprises several concepts that capture the prominent features of WMSs. Based on this terminology, this paper also presents a classification of 23 existing WMSs according to the proposed axes and terms.

Community-based terminology↗

Constructing high-fidelity halo merger trees in abacussummit

ABSTRACT Tracking the formation and evolution of dark matter haloes is a critical aspect of any analysis of cosmological N-body simulations. In particular, the mass assembly of a halo and its progenitors, encapsulated in the form of its merger tree, serves as a fundamental input for constructing semi-analytic models of galaxy formation and, more generally, for building mock catalogues that emulate galaxy surveys. We present an algorithm for constructing halo merger trees from abacussummit, the largest suite of cosmological N-body simulations performed to date consisting of nearly 60 trillion particles, and which has been designed to meet the Cosmological Simulation Requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. Our method tracks the cores of haloes to determine associations between objects across multiple time slices, yielding lists of halo progenitors and descendants for the several tens of billions of haloes identified across the entire suite. We present an application of these merger trees as a means to enhance the fidelity of abacussummit halo catalogues by flagging and ‘merging’ haloes deemed to exhibit non-monotonic past merger histories. We show that this cleaning technique identifies portions of the halo population that have been deblended due to choices made by the halo finder, but which could have feasibly been part of larger aggregate systems. We demonstrate that by cleaning halo catalogues in this post-processing step, we remove potentially unphysical features in the default halo catalogues, leaving behind a more robust halo population that can be used to create highly accurate mock galaxy realizations from abacussummit.

79 ASTRONOMY AND ASTROPHYSICS↗

Mars dust storms - Interannual variability and chaos

The hypothesis is that the global climate system, consisting of atmospheric dust interacting with the circulation, produces its own interannual variability when forced at the annual frequency. The model has two time-dependent variables representing the amount of atmospheric dust in the northern and southern hemispheres, respectively. Absorption of sunlight by the dust drives a cross-equatorial Hadley cell that brings more dust into the heated hemisphere. The circulation decays when the dust storm covers the globe. Interannual variability manifests itself either as a periodic solution in which the period is a multiple of the Martian year, or as an aperiodic (chaotic) solution that never repeats. Both kinds of solution are found in the model, lending support to the idea that interannual variability is an intrinsic property of the global climate system. The next step is to develop a hierarchy of dust-circulation models capable of being integrated for many years.

Ingersoll, Andrew P.↗

Comparative Magma Oceanography

The question of whether the Earth ever passed through a magma ocean stage is of considerable interest. Geochemical evidence strongly suggests that the Moon had a magma ocean and the evidence is mounting that the same was true for Mars. Analyses of martian (SNC) meteorites have yielded insights into the differentiation history of Mars, and consequently, it is interesting to compare that planet to the Earth. Three primary features of Mars contrast strongly to those of the Earth: (i) the extremely ancient ages of the martian core, mantle, and crust (about 4.55 b.y.); (ii) the highly depleted nature of the martian mantle; and (iii) the extreme ranges of Nd isotopic compositions that arise within the crust and depleted mantle. The easiest way to explain the ages and diverse isotopic compositions of martian basalts is to postulate that Mars had an early magma ocean. Cumulates of this magma ocean were later remelted to form the SNC meteorite suite and some of these melts assimilated crustal materials enriched in incompatible elements. The REE pattern of the crust assimilated by these SNC magmas was LREE enriched. If this pattern is typical of the crust as a whole, the martian crust is probably similar in composition to melts generated by small degrees of partial melting (about 5%) of a primitive source. Higher degrees of partial melting would cause the crustal LREE pattern to be essentially flat. In the context of a magma ocean model, where large degrees of partial melting presumably prevailed, the crust would have to be dominated by late-stage, LREE-enriched residual liquids. Regardless of the exact physical setting, Nd and W isotopic evidence indicates that martian geochemical reservoirs must have formed early and that they have not been efficiently remixed since. The important point is that in both the Moon and Mars we see evidence of a magma ocean phase and that we recognize it as such. Several lines of theoretical inference point to an early Earth that was also hot and, perhaps, mostly molten. The Giant Impact hypothesis for the origin of the Moon offers a tremendous input of thermal energy and the same could be true for core formation. And current solar system models favor the formation of a limited number of large (about 1000 km) planetesimals that, upon accreting to Earth, would cause great heating, being lesser versions of the Giant Impact. Several lines of geochemical evidence do not favor this hot early Earth scenario. (i) Terrestrial man-tle xenoliths are sometimes nearly chondritic in their major element compositions, suggesting that these rocks have never been much molten. Large degrees of partial melting probably promote differentiation rather than homogenization. (ii) Unlike the case of Mars, the continental crust probably did not form as a highly fractionated residual liquid from a magma ocean (about 99% crystallization), but, rather, formed in multiple steps. [The simplest model for the formation of continental crust is complicated: (a) about 10% melting of a primitive mantle, making basalt; (b) hydrothermal alteration of that basalt, converting it to greenstone; and (c) 10% partial melting of that greenstone, producing tonalite.] This model is reinforced by the recent observation from old (about 4.1 b.y.) zircons that the early crust formed from an undepleted mantle having a chondritic Lu/Hf ratio. (iii) If the mantle were once differentiated by a magma ocean, the mantle xenolith suite requires that it subsequently be homogenized. The Os isotopic compositions of fertile spinel lherzolites place constraints on the timing of that homogenization. The Os isotopic composition of spinel lherzolites approaches that of chondrites and correlates with elements such as Lu and Al. As Lu and Al concentrations approach those of the primitive mantle, Os isotopic compositions approach chondritic. The Re and Os in these xenoliths were probably added as a late veneer. Thus, the mantle that received the late veneer must have been nearly chondritic in terms of its major elements (excluding Fe). If the mantle that the veneer was mixed into was not al-ready homogenized, then Os isotopes should not correlate with incompatible elements such as Al. Consequently, either early differentiation of the mantle did not occur or the homogenization of this differentiation must have occurred before the late veneer was added. The timing of the late veneer is itself uncertain but presumably postdated core formation at about 4.45 b.y. and did not postdate the 3.8-3.9 b.y. late bombardment of the Moon. This timing based on siderophile elements is consistent with the Hf isotopic evidence cited above. If the Earth, Moon and Mars had magma oceans, the Earth subsequently rehomogenized whereas the Moon and Mars did not. The simplest solution to this observation is that homogenization of igneous differentiates was never necessary on Earth, either because the hypothetical magma ocean never occurred or because this event did not produce mantle differentiation.

Jones, J. H.↗

Directing Nanoparticle Organization in Response to Diverse Chemical Inputs

Signaling cascades are crucial for transducing stimuli in biological systems, enabling multiple stimuli to regulate a downstream target with precisely controlled timing and amplifying signals through a series of intermediary reactions. Developing a robust signaling system with such capabilities would be pivotal for programming complex behaviors in synthetic DNA-based molecular devices. However, although “software” such as nucleic acid circuits could potentially be harnessed to relay signals to DNA-based nanostructure hardware, such explorations have been limited. Here, in this study, we develop a platform for transducing a variety of stimuli via messenger-mediated reactions to regulate the release and reloading of gold nanoparticles (AuNPs) in a 3D DNA framework. In the first step, an in vitro transcription circuit is engineered to sense and amplify chemical stimuli, including arbitrary DNA sequences and proteins, producing RNA. In the second step, the RNA releases the DNA-coated AuNPs from the DNA framework via a strand displacement reaction. AuNP reloading is controlled by a separate step driven by degradation of the RNA. Our platform holds promise for applications requiring dynamic multiagent control over DNA-based devices, offering a versatile tool for advanced molecular device engineering.

36 MATERIALS SCIENCE↗

Signal processing applications of massively parallel charge domain computing devices

The present invention is embodied in a charge coupled device (CCD)/charge injection device (CID) architecture capable of performing a Fourier transform by simultaneous matrix vector multiplication (MVM) operations in respective plural CCD/CID arrays in parallel in O(1) steps. For example, in one embodiment, a first CCD/CID array stores charge packets representing a first matrix operator based upon permutations of a Hartley transform and computes the Fourier transform of an incoming vector. A second CCD/CID array stores charge packets representing a second matrix operator based upon different permutations of a Hartley transform and computes the Fourier transform of an incoming vector. The incoming vector is applied to the inputs of the two CCD/CID arrays simultaneously, and the real and imaginary parts of the Fourier transform are produced simultaneously in the time required to perform a single MVM operation in a CCD/CID array.

Fijany, Amir↗

A design study for the use of a multiple aperture deployable antenna for soil moisture remote sensing satellite applications

The instrumentation problems associated with the measurement of soil moisture with a meaningful spatial and temperature resolution at a global scale are addressed. For this goal only medium term available affordable technology will be considered. The study while limited in scope, will utilize a large scale antenna structure, which is being developed presently as an experimental model. The interface constraints presented by a singel Space Transportation System (STS) flight will be assumed. Methodology consists of the following steps: review of science requirements; analyze effects of these requirements; present basic system engineering considerations and trade-offs related to orbit parameters, number of spacecraft and their lifetime, observation angles, beamwidth, crossover and swath, coverage percentage, beam quality and resolution, instrument quantities, and integration time; bracket the key system characteristics and develop an electromagnetic design of the antenna-passive radiometer system. Several aperture division combinations and feed array concepts are investigated to achieve maximum feasible performacne within the stated STS constraints.

Foldes, P.↗

Correlated Trajectory Uncertainty for Adaptive Sequential Decision Making

One of the great challenges with decision making tasks on real world systems is the fact that data is sparse and acquiring additional data is expensive. In these cases, it is often crucial to make a model of the environment to assist in making decisions. At the same time, limited data means that learned models are erroneous, making it just as important to equip the model with good predictive uncertainties. In the context of learning sequential decision making policies, these uncertainties can prove useful for informing which data to collect for the greatest improvement in policy performance \citep{mehta2021experimental, mehta2022exploration} or informing the policy about unsure regions of state and action space to avoid during test time \citep{yu2020mopo}. Additionally, assuming that realistic samples of the environment can be drawn, an adaptable policy can be trained that attempts to make optimal decisions for any given possible instance of the environment \citep{ghosh2022offline, chen2021offline}. In this work, we examine the so-called ``probabilistic neural network'' (PNN) model that is ubiquitous in model-based reinforcement learning (MBRL) works. We argue that while PNN models may have good marginal uncertainties, they form a distribution of non-smooth transition functions. Not only are these samples unrealistic and may hamper adaptability, but we also assert that this leads to poor uncertainty estimates when predicting multiple step trajectory estimates. To address this issue, we propose a simple sampling method that can be implemented on top of pre-existing models.We evaluate our sampling technique on a number of environments, including a realistic nuclear fusion task, and find that, not only do smooth transition function samples produce more calibrated uncertainties, but they also lead to better downstream performance for an adaptive policy.

Offline Reinforcement Learning↗

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗

Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds

Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.

Farrugia, C. J.↗