Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPACE LAW”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Study of structural, transport and magneto-crystalline anisotropy in La 1-x Sr x MnO 3 (0.30 ≤ x ≤ 0.40) perovskite manganites

In this paper we present structural, transport, magnetic and magneto-crystalline anisotropy of La 1–x Sr x MnO 3 (х = 0.30, 0.33, 0.36, and 0.40) synthesized using the solid-state reaction method. The X-ray diffraction pattern of the samples is well indexed to the rhombohedral structure with R 3 - c space group. Lattice parameters and unit cell volume are found to decrease monotonically upon increasing Sr 2+ concentration at La 3+ site. X-ray photoelectron spectroscopy confirms the presence of Mn 4+ and Mn 3+ , the ratio of which (Mn 4+ /Mn 3+ ) increases from 0.44 to 0.64 with increasing Sr 2+ concentration. From the ρ(T) and M(T) data analysis, we find that T MI and M S are larger for x = 0.33 composition, suggesting a weakening of the double-exchange interaction due to competing super-exchange interactions in x = 0.30, 0.36 and 0.40 samples. We calculated MCA constant (K 1 ) by fitting the saturation magnetisation data at a higher magnetic field using the law of approach model, M = M S [1-a/H-b/H 2 ]-χH with K 1 2 = (105/8)bMs 2 and is in the order of 10 5 erg/cm 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational framework for behind-the-meter DER techno-economic modeling and optimization: REopt Lite

The energy system is undergoing a major transformation with the global emphasis on decarbonization. Distributed generation is projected to play a significant role in the new energy system, and energy models are informing how distributed generation can be integrated reliably and economically. In this work, we present an end-to-end computational framework for distributed energy resource (DER) modeling, REopt Lite™, which captures the interface of technology, economics, and policy in the energy modeling process. We describe the problem space, the building blocks of the model, the scaling capabilities of the design, the optimization formulation, and the extensibility of the model. We present a framework for accelerating the techno-economic analysis of behind-the-meter distributed energy resources to enable rapid planning and decision-making, thereby enabling greater renewable energy deployment. This computation framework is open-sourced to facilitate transparency, flexibility, and wider collaboration opportunities within the worldwide energy modeling community.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Hybrid Simulation and Quasi-linear Theory of Bi-Kappa Proton Instabilities

The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropy T ⊥ > T ∥ (where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κ power-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Continuum-mediated self-interacting dark matter

Dark matter may self-interact through a continuum of low-mass states. This happens if dark matter couples to a strongly-coupled nearly-conformal hidden sector. This type of theory is holographically described by brane-localized dark matter interacting with bulk fields in a slice of 5D anti-de Sitter space. The long-range potential in this scenario depends on a non-integer power of the spatial separation, in contrast to the Yukawa potential generated by the exchange of a single 4D mediator. The resulting self-interaction cross section scales like a non-integer power of velocity. We identify the Born, classical and resonant regimes and investigate them using state-of-the-art numerical methods. We demonstrate the viability of our continuum-mediated framework to address the astrophysical small-scale structure anomalies. Investigating the continuum-mediated Sommerfeld enhancement, we demonstrate that a pattern of resonances can occur depending on the non-integer power. We conclude that continuum mediators introduce novel power-law scalings which open new possibilities for dark matter self-interaction phenomenology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

We measure the effective opacity (τ{sub eff}) of the intergalactic medium from the composite spectra of 281 Lyman-break galaxies in the redshift range 2 ≲ z ≲ 3. Our spectra are taken from the COSMOS Lyα Mapping And Tomographic Observations survey derived from the Low Resolution Imaging Spectrometer on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Lyα forest, we measure the effective Lyα opacity (τ{sub eff}) in the 2.02 ≤ z ≤ 2.44 range. At z = 2.22, we estimate τ{sub eff} =0.159±0.001 from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at z < 2.5 indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian processes model of our results and previous τ{sub eff} measurements that describes the steep redshift evolution in τ{sub eff} from z = 1.5–4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

FEgen (v.1): Field Emission Distribution Generator Freeware Based on Fowler–Nordheim Equation

This work presents an open-source code, FEgen, that generates initial field emission distributions based on the time-dependent Fowler–Nordheim equation. Current in-demand and open-source beam tracking software—such as ASTRA, IMPACT-T, and GPT—do not contain native distribution generator functions suitable for field emitted electron tracking and phase space analysis. Furthermore, the intuitive graphical user interface of FEgen generates input files shown to be fully compatible with both ASTRA and GPT. Here, detailed application examples demonstrating a subset of FEgen’s capabilities are given by creating a customized emitter pattern with application to transverse beam shaping and longitudinal beam dynamic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Galaxy morphologies revealed with Subaru HSC and super-resolution techniques. I. Major merger fractions of L UV ∼ 3–15 L*UV dropout galaxies at z ∼ 4–7

Abstract We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z ∼ 4–7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ∼1″ to $\lesssim \!\!{0{^{\prime \prime }_{.}}1}$, which is comparable to that of the Hubble Space Telescope, allowing us to identify z ∼ 4–7 bright major mergers at a high completeness value of $\gtrsim \!\!90\%$. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z ∼ 4, 5, 6, and 7, respectively, in a UV luminosity range of LUV ∼ 3–$15\, L_{\rm UV}^*$ corresponding to −24 ≲ MUV ≲ −22. The major merger fractions are estimated to be $f_{\rm merger}\sim 10\%$–$20\%$ at z ∼ 4 and $\sim 50\%$–$70\%$ at z ∼ 5–7, which shows no fmerger difference compared to those of a control faint galaxy sample. Based on the fmerger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z ∼ 4–7 galaxy UV LFs. While these two effects partly explain the DPL shape at LUV ∼ 3–$10\, L_{\rm UV}^*$, the DPL shape cannot be explained at the very bright end of $L_{\rm UV}\gtrsim 10\, L_{\rm UV}^*$, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs.

Shibuya, Takatoshi↗
DOE OSTI2025

Photocurrent Nanoscopy of Quantum Hall Bulk

Understanding nanoscale electronic and thermal transport of two-dimensional (2D) electron systems in the quantum Hall regime, particularly in the bulk insulating state, poses considerable challenges. One of the primary difficulties arises from the presence of chiral edge channels, whose transport behavior obscures the investigation of the insulating bulk. Using near-field optical and photocurrent nanoscopy, we probe real-space variations of the optical and thermal dynamics of graphene in the quantum Hall regime without relying on complex sample or electrode geometries. Near the charge-neutrality point, we detect strong optical and photothermal signals from resonant inter-Landau-level (LL) magnetoexciton excitations between the zeroth and ± 1 st LLs, which gradually weaken with increasing doping due to Pauli blocking. Interestingly, at higher doping levels and full-integer LL fillings, photothermal signals reappear across the entire sample over an approximately 10 − μ m scale, indicating unexpectedly long cooling lengths and nonlocal photothermal heating through the insulating bulk. This observation suggests that thermal conductivity persists for the localized states even as electronic transport is suppressed—a clear violation of the Wiedemann-Franz law. Our experiments provide novel insights into nanoscale thermal and electronic transport in incompressible 2D gases, highlighting the roles of magnetoexcitons and chiral edge states in the thermo-optoelectric dynamics of the Dirac quantum Hall state. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uncertainty quantification in scientific machine learning: Methods, metrics, and comparisons

Neural networks (NNs) are currently changing the computational paradigm on how to combine data with mathematical laws in physics and engineering in a profound way, tackling challenging inverse and ill-posed problems not solvable with traditional methods. However, quantifying errors and uncertainties in NN-based inference is more complicated than in traditional methods. This is because in addition to aleatoric uncertainty associated with noisy data, there is also uncertainty due to limited data, but also due to NN hyperparameters, overparametrization, optimization and sampling errors as well as model misspecification. Although there are some recent works on uncertainty quantification (UQ) in NNs, there is no systematic investigation of suitable methods towards quantifying the total uncertainty effectively and efficiently even for function approximation, and there is even less work on solving partial differential equations and learning operator mappings between infinite-dimensional function spaces using NNs. In this work, we present a comprehensive framework that includes uncertainty modeling, new and existing solution methods, as well as evaluation metrics and post-hoc improvement approaches. Further, to demonstrate the applicability and reliability of our framework, we present an extensive comparative study in which various methods are tested on prototype problems, including problems with mixed input-output data, and stochastic problems in high dimensions. In the Appendix, we include a comprehensive description of all the UQ methods employed. Further, to help facilitate the deployment of UQ in Scientific Machine Learning research and practice, we present and develop in [1] an open-source Python library (github.com/Crunch-UQ4MI/neuraluq), termed NeuralUQ, that is accompanied by an educational tutorial and additional computational experiments.

11 physics-informed neural networks↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Feynman diagrams for matter wave interferometry

We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.

Glick, Jonah [Northwestern U.; Fermilab] (ORCID:00↗
DOE OSTI2026

Illinois Storage Corridor - CarbonSAFE Phase III: Policy, Regulatory, Legal and Permitting Characteristics; Subtask 5.5

The Illinois Storage Corridor (ISC) project aims to advance the commercial development and implementation of carbon capture and storage (CCS) technologies within a region in Illinois of suitable geology for carbon dioxide (CO 2 ) storage in deep saline aquifers. The project partners have annual emissions greater than 6.5 million tonnes per year; storage hubs are being explored at two sites—one near the One Earth Energy facility in east-central Illinois, and one near the Prairie State Generating Company campus in southwest Illinois. To implement the technology, legal, policy, and economic considerations must be explored. The United States Environmental Protection Agency (US EPA) administers the Underground Injection Control (UIC) program and is responsible for issuing Class VI permits to construct and operate CO 2 injection wells, i.e., a UIC Class VI well permit is required to inject CO 2 into the subsurface for geologic sequestration. The largest consideration of the Class VI well requirements is to protect underground sources of drinking water (USDWs). Additionally discussed are considerations for Class VI permits relating to public engagement. Besides permitting, property rights to storage sites, subsurface pore spaces, and areas for pipeline transportation must be secured. Pore space rights is still a novel concept being explored and not yet addressed by the Illinois legislature. It is believed that surface property rights are required for the subsurface pore space below, so long as there is not a separated mineral estate in the subsurface. Illinois legislature has addressed securing rights-of-way for CO 2 transportation, allowing for easements and the exercise of eminent domain to secure such rights. Economically, the incentives for CCS are ever-expanding. Recently, the federal government broadened the availability and increased the dollar-amounts for the § 45Q tax credits for geological storage of carbon oxides. Congress has also authorized the Advanced Industrial Facilities Development Program which allocates billions of dollars in funds to installing technology at industrial facilities to reduce greenhouse gas emissions. On the state level, Illinois has had its own incentives for CCS projects since 2009 and expanded its emission-related goals again in 2021. The law enacted in 2021 specifically creates a commission to explore implementing CCS at Prairie State Generating Company, a partner on the Illinois Storage Corridor project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Minimizing exposure to legacy wells and avoiding conflict between storage projects: Exploring area of review as a screening tool

Elevated pressure from large-volume injection is a key driver of risk and project cost. If transmissive features (e.g., non-isolating wells or fracture systems) are present, increased injection-zone pressure can drive fluids from depth toward protected freshwater resources. In US Carbon Capture and Storage (CCS) law, the area at risk is known as the Area of Review (AoR). The size and number of potentially transmissive features to be evaluated and possibly remediated or managed is a function of the size and location of the AoR. The size of the AoR depends on several variables, including properties of the injection zone, properties of protected resources, and injection rate and duration. Evaluation of the intersection of these variables across a portfolio of sites highlights the injection zone depth and boundary conditions as top-level controls. Deep injection, use of multiple stacked injection zones, reduced injection rate and choice of injection well location can all be used to minimize AoR and the number of potentially transmissive features within it. Here, we introduce the concept of pressure space (defined as connected pore volume times pressure) as the key subsurface commodity for CO 2 storage and we suggest that it forms a more robust basis for leasing and regulation than pore space alone.

58 GEOSCIENCES↗
DOE OSTI2021

Structural, vibrational, thermal, and magnetic properties of mullite-type NdMnTiO 5 ceramic

Mullite-type $R$Mn 2 O 5 ($R$ = Y, rare-earth element) ceramics are of ongoing research attention because of their interesting crystal-chemical and magnetic properties. We report nuclear and magnetic structures of NdMnTiO 5 together with its spectroscopic, thermogravimetric, and magnetic properties. The polycrystalline sample is prepared by solid-state synthesis and characterized from neutron and X-ray powder diffraction data Rietveld refinements. NdMnTiO 5 crystallizes in the orthorhombic space group $Pbam$ with metric parameter $a$ = 755.20(1) pm, $b$ = 869.91(1) pm, $c$ = 582.42(1) pm, and $V$ = 382.62(1) 10 6 pm 3 . The Mn 3+ and Ti 4+ cations are observed to be located in the octahedral and pyramidal sites, respectively. The vibrational features in these polyhedral sites are characterized by Raman and Fourier transform infrared spectroscopes. The higher decomposition temperature of NdMnTiO 5 , compared to other $R$Mn 2 O 5 phases, is explained in terms of the higher bond strength of Ti-O bonds than those of Mn-O bonds. Temperature-dependent DC magnetic susceptibility suggests a paramagnetic to antiferromagnetic phase transition at 43(1) K. Inverse susceptibility in the paramagnetic region above 120 K follows the Curie-Weiss law, resulting in a magnetic moment of 6.33(1) μ B per formula unit. Neutron diffraction data collected at 7.5 K reveal that the magnetic moments of Nd 3+ and Mn 3+ in NdMnTiO 5 are incommensurately ordered with a propagation vector $k$ = (0, 0.238, 0.117).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Seamless Hybrid-integrated Interconnect NEtwork (SHINE)

Large-scale date centers are becoming increasingly limited by the capacity of interconnects due to the rapid increasing power consumption and bandwidth demands driven by the rigorous down scaling of CMOS critical dimensions in according with the Moore’s Law. As a result of the continuing bandwidth scaling, optical interconnects based on parallel optics, where data are simultaneously transmitted along more than one lane, have becoming a preferred interconnected solution over copper cables. To building an efficient optical interconnect system, optical coupling between different elements, such as optical fibers, photonics chips, is a critical step, which is challenging due to the mode mismatch of varies optical ports. In this project, we proposed and demonstrated an optical interfacing scheme, which provided a light coupling approach with low energy loss, broadband operation, high alignment tolerance and compatible with scalable fabrication and testing. Here, the mode transformation is achieved by a total internal reflection from a chip integrated microscale freefrom surface. Due to the enormous light wavefront manipulating ability associated with the gigantic design space, the freeform couplers platform is universally applicable to a wide range of optical coupling scenarios, such as the interfacing between a waveguide to a fiber, another waveguide, surface normal device, or free space.

42 ENGINEERING↗
DOE OSTI2024

New Measurements of the Lyα Forest Continuum and Effective Optical Depth with LyCAN and DESI Y1 Data

Abstract We present the Ly α Continuum Analysis Network (LyCAN), a convolutional neural network that predicts the unabsorbed quasar continuum within the rest-frame wavelength range of 1040–1600 Å based on the red side of the Ly α emission line (1216–1600 Å). We developed synthetic spectra based on a Gaussian mixture model representation of nonnegative matrix factorization (NMF) coefficients. These coefficients were derived from high-resolution, low-redshift ( z < 0.2) Hubble Space Telescope/Cosmic Origins Spectrograph (COS) quasar spectra. We supplemented this COS-based synthetic sample with an equal number of DESI Year 5 mock spectra. LyCAN performs extremely well on testing sets, achieving a median error in the forest region of 1.5% on the DESI mock sample, 2.0% on the COS-based synthetic sample, and 4.1% on the original COS spectra. LyCAN outperforms principal component analysis (PCA) and NMF-based prediction methods using the same training set by 40% or more. We predict the intrinsic continua of 83,635 DESI Year 1 spectra in the redshift range of 2.1 ≤ z ≤ 4.2 and perform an absolute measurement of the evolution of the effective optical depth. This is the largest sample employed to measure the optical depth evolution to date. We fit a power law of the form τ ( z ) = τ 0 ( 1 + z ) γ to our measurements and find τ 0 = (2.46 ± 0.14) × 10 −3 and γ = 3.62 ± 0.04. Our results show particular agreement with high-resolution, ground-based observations around z = 2, indicating that LyCAN is able to predict the quasar continuum in the forest region with only spectral information outside the forest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Thermal transport in compensated semimetals: Effect of electron-electron scattering on Lorenz ratio

It is well known that the electronic thermal conductivity of clean compensated semimetals can be greatly enhanced over the electric conductivity by the availability of an ambipolar mechanism of conduction, whereby electrons and holes flow in the same direction experiencing negligible Coulomb scattering as well as negligible impurity scattering. This enhancement—resulting in a breakdown of the Wiedemann-Franz law with an anomalously large Lorenz ratio—has been recently observed in two-dimensional monolayer and bilayer graphene near the charge neutrality point. In contrast to this, three-dimensional compensated semimetals such as WP 2 and Sb are typically found to show a reduced Lorenz ratio. Here, we investigate the reasons for this difference, focusing on the low-temperature regime where the electron-electron scattering is expected to dominate over other scattering mechanisms. We show that the different regimes of Fermi statistics (nondegenerate electron-hole liquid in graphene versus degenerate electron-hole liquid in compensated semimetals) are not sufficient to explain the reduction of the Lorenz ratio in the latter. We propose that the solution of the puzzle lies in the large separation of electron and hole pockets in momentum space, which allows compensated semimetals to sustain sizable regions of electron-hole accumulation near the contacts. These accumulations suppress the ambipolar conduction mechanism and effectively split the system into two independent electron and hole conductors. We present a quantitative theory of the crossover from ambipolar to unipolar conduction as a function of the size of the electron-hole accumulation regions, and show that it naturally leads to a sample-size-dependent thermal conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗