Search NASA⌕ Search

SEARCH · Search NASA

Results for “local significance”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Thermal Modeling of Mars Ground for Surface Missions

Thermal analysis packages are capable of including ground temperature effects for orbital and surface thermal analyses. In particular, Thermal Desktop® offers the option of specifying ground temperatures as a function of time for planetary surface modeling. While suitable for many cases, this approach is not sufficient if an object has local interactions with the ground that could significantly affect the ground temperature. Ground modeling is necessary for the Mars 2020 rover thermal design and analysis since shadowing and heat dissipation from the rover’s Multi-Mission Radioisotope Thermoelectric Generator (MMRTG) can result in significant temperature deviations of the local ground. The Mars 2020 thermal team is explicitly modeling the Martian ground so that these local temperature effects can be captured. The upper portion of the ground is modeled, and material and optical properties of the ground are varied in order to match data collected from orbiting science instruments. Atmospheric surface temperature, sky temperatures, and solar loads from a Mars General Circulation Model (GCM) are used as boundary conditions, resulting in a ground surface temperature profile consistent with the GCM predictions. The rover model is then placed on this modeled ground so that the effects of shadowing and MMRTG dissipation on the ground temperature can be captured.

Farias, Edgardo↗
DOE OSTI2020

Rapid Memory Footprint Access Diagnostics

Footprint and reuse distance measure temporal locality and therefore do not capture the significance of access patterns (spacial locality). A strided access pattern has the largest possible footprint but usually has the best performance. To highlight exposed memory latency, we separate footprint into strided (prefetchable) and irregular (non-prefetchable) access components and calculate the growth rate of each. To rapidly compute these footprint access diagnostics, we present two methods, whole-program and precise. Current footprint analyses can cause 200× or more slowdown with realistic inputs and are therefore impractical. Our whole-program method reduces the overhead to 10% by computing upper bounds, but still yields inter-procedural insight through a call path profile. Our precise method uses additional static analysis and profiling to refine the upper bounds for intra-procedural loop nests. We evaluate our approaches using benchmarks that vary access patterns (strided \vs unpredictable), sparsity (all words in a cache line \vs some), and reuse (varying and repeated accesses per element). Notably, for loop nests with unpredictable accesses, the precise method's accuracy is within 10% of ideal. The whole-program method has sufficient accuracy to diagnose bottlenecks.

data locality, memory footprint, footprint access ↗
DOE OSTI2020

Uncertainty analysis of correlated parameters in automated reaction mechanism generation

Abstract Uncertainty analysis is a useful tool for inspecting and improving detailed kinetic mechanisms because it can identify the greatest sources of model output error. Owing to the very nonlinear relationship between kinetic and thermodynamic parameters and computed concentrations, model predictions can be extremely sensitive to uncertainties in some parameters while uncertainties in other parameters can be irrelevant. Error propagation becomes even more convoluted in automatically generated kinetic models, where input uncertainties are correlated through kinetic rate rules and thermodynamic group values. Local and global uncertainty analyses were implemented and used to analyze error propagation in Reaction Mechanism Generator (RMG), an open‐source software for generating kinetic models. A framework for automatically assigning parameter uncertainties to estimated thermodynamics and kinetics was created, enabling tracking of correlated uncertainties. Local first‐order uncertainty propagation was implemented using sensitivities computed natively within RMG. Global uncertainty analysis was implemented using adaptive Smolyak pseudospectral approximations as implemented in the MIT Uncertainty Quantification Library to efficiently compute and construct polynomial chaos expansions to approximate the dependence of outputs on a subset of uncertain inputs. Cantera was used as a backend for simulating the reactor system in the global analysis. Analyses were performed for a phenyldodecane pyrolysis model. Local and global methods demonstrated similar trends; however, many uncertainties were significantly overestimated by the local analysis. Both local and global analyses show that correlated uncertainties based on kinetic rate rules and thermochemical groups drastically reduce a model's degrees of freedom and have a large impact on the determination of the most influential input parameters. These results highlight the necessity of incorporating uncertainty analysis in the mechanism generation workflow.

Gao, Connie W.↗
DOE OSTI2021

Local prediction of Laser Powder Bed Fusion porosity by short-wave infrared imaging thermal feature porosity probability maps

We report that local thermal history can significantly vary in parts during metal Additive Manufacturing (AM), leading to local defects. However, the sequential layer-by-layer nature of AM facilitates in-situ part voxelmetric observations that can be used to detect and correct these defects for part qualification and quality control. The challenge is to relate this local radiometric data with local defect information to estimate process error likelihood in future builds. This paper uses a Short-Wave Infrared (SWIR) camera to record the temperature history for parts manufactured with Laser Powder Bed Fusion (LPBF) processes. The porosity from a cylindrical specimen is measured by ex-situ micro-computed tomography (μCT). Specimen data from the SWIR camera, combined with the μCT data, are used to generate thermal feature-based porosity probability maps. The porosity predictions made by various SWIR thermal feature-porosity probability maps of a specimen with a complex geometry are scored against the true porosity obtained via μCT. The receiver operating characteristic curves constructed from the predictions for the complex sample demonstrate the porosity probability mapping methodology’s potential for in-situ based porosity detection.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Effect of simulated hyperemia on the flow field in a mildly atherosclerotic coronary artery casting of man

Changes in an arterial flow field due to mild atherosclerosis were determined using a main coronary artery casting with a maximum obstruction of about 50 percent by area. Local pressure changes were measured using six pressure tap holes along the wall of the casting. The test-fluid was a 33 percent sugar-water solution of approximately the same viscosity as human blood. Flow visualization results were obtained by injecting blue-dye through the pressure tap holes. Measurement of local pressure demonstrated a significant Reynolds number effect. At Reynolds numbers of 80-710, a local pressure rise was observed downstream of the mild atherosclerotic constriction due to momentum changes. The Reynolds number necessary for flow separation in the divergent region of the coronary casting was about 330. The experimental results can be used to obtain a quantitative relation between coronary morphology and the fluid dynamic consequences of mild diffuse disease under conditions of maximum cardiac demand i.e., higher coronary flow rates and Reynolds numbers associated with space and atmospheric flight.

Cho, Y. I.↗
NASA NTRS2020

Constraining Fossil Fuel CO2 Emissions from Urban Area Using OCO-2 Observations of Total Column CO2

Satellite observations of the total column dry-air CO2(XCO2) are expected to support the quantification and monitoring of fossil fuel CO2 (ffCO2) emissions from urban areas. We evaluate the utility of the Orbiting Carbon Observatory 2 (OCO-2) XCO2 retrievals to optimize whole-city emissions, using a Bayesian inversion system and high-resolution transport modeling. The uncertainties of constrained emissions related to transport model, satellite measurements, and local biospheric fluxes are quantified. For the first two uncertainty sources, we examine cities of different landscapes: “plume city” located in relatively flat terrain, represented by Riyadh and Cairo; “basin city” located in basin terrain, represented by Los Angeles (LA). The retrieved scaling factors of emissions and their uncertainties show prominent variabilities from track to track, due to the varying meteorological conditions and relative locations of the tracks transecting plumes. To explore the performance of multiple tracks in retrieving emissions, pseudo data experiments are carried out. The estimated least number of tracks required to constrain the total emissions for Riyadh (<10% uncertainty), Cairo (<10%), and LA (<5%) are 5, 8, and 7, respectively. Additionally, to evaluate the impact of biospheric fluxes on derivation of the ffXCO2 enhancements, we conduct simulations for Pearl River Delta metropolitan area. Significant fractions of local XCO2 enhancements associated with local biospheric XCO2 variations are shown, which potentially lead to biased estimates of ffCO2 emissions. We demonstrate that satellite measurements can be used to improve urban ffCO2 emissions with a sufficient amount of measurements and appropriate representations of the uncertainty components.

Xinxin Ye↗
DOE OSTI2022

Projecting the Thermal Response in a HTGR-Type System during Conduction Cooldown Using Graph-Laplacian Based Machine Learning

Accurate prediction of an off-normal event in a nuclear reactor is dependent upon the availability of sensory data, reactor core physical condition, and understanding of the underlying phenomenon. This work presents a method to project the data from some discrete sensory locations to the overall reactor domain during conduction cooldown scenarios similar to High Temperature Gas-cooled Reactors (HTGRs). The existing models for conductive cooldown in a heterogeneous multi-body system, such as an assembly of prismatic blocks or pebble beds relies on knowledge of the thermal contact conductance, requiring significant knowledge of local thermal contacts and heat transport possibilities across those contacts. With a priori knowledge of bulk geometry features and some discrete sensors, a machine learning approach was devised. The presented work uses an experimental facility to mimic conduction cooldown with an assembly of 68 cylindrical rods initially heated to 1200 K. High-fidelity temperature data were collected using an infrared (IR) camera to provide training data to the model and validate the predicted temperature data. The machine learning approach used here first converts the macroscopic bulk geometry information into Graph-Laplacian, and then uses the eigenvectors of the Graph-Laplacian to develop Kernel functions. Support vector regression (SVR) was implemented on the obtained Kernels and used to predict the thermal response in a packed rod assembly during a conduction cooldown experiment. The usage of SVR modeling differs from most models today because of its representation of thermal coupling between rods in the core. When trained with thermographic data, the average normalized error is less than 2% over 400 s, during which temperatures of the assembly have dropped by more than 500 K. The rod temperature prediction performance was significantly better for rods in the interior of the assembly compared to those near the exterior, likely due to the model simplification of the surroundings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Multiphysics simulation of recent experiments on alkali‐silica reaction expansion in reinforced concrete members

Alkali‐silica reaction (ASR) is an important degradation process that causes volumetric expansion and damage in concrete, and is affected significantly by the local temperature, moisture and stress conditions that often vary across the regions of a structure. Numerical simulation is essential to predict the progression and effects of ASR on the performance of structures. Because of the interactions between thermal and moisture transport and mechanical deformation, it is important for numerical models to represent all these physical phenomena and the coupling between them. Simulations of ASR in reinforced concrete (RC) structures are further complicated by the need to capture interactions between concrete and embedded reinforcing bars. Here, this paper describes the implementation of a scalable, coupled‐physics ASR model for simulating RC structures and assesses the ability of that model to predict ASR‐induced expansion in recent laboratory tests on RC block and beam specimens. These laboratory tests and the simulation approach were selected because of their applicability to RC structural‐scale simulations. This validation study helps builds confidence the ability of this approach to model ASR expansion in large, complex RC structures, which is a current high‐priority need.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dual carbonate clumped isotopes (Δ 47 -Δ 48 ) constrains kinetic effects and timescales in peridotite-associated springs at the Cedars, Northern California

The Cedars is an area in Northern California with a chain of highly alkaline springs resulting from CO 2 -charged meteorological water interacting with a peridotite body. Serpentinization resulting from this interaction at depth leads to the sequestration of various carbonate minerals into veins accompanied by a release of Ca 2+ and OH – enriched water to the surface, creating an environment which promotes rapid precipitation of CaCO 3 at surface springs. This environment enables us to apply the recently developed Δ 47 -Δ 48 dual clumped isotope analysis to probe kinetic isotope effects (KIEs) and timescales of CO 2 transformation in a region with the potential for geological CO 2 sequestration. Here, we analyzed CaCO 3 recovered from various localities and identified significant kinetic fractionations associated with CO 2 absorption in a majority of samples, characterized by enrichment in Δ 47 values and depletion in Δ 48 values relative to equilibrium. Surface floes exhibited the largest KIEs (ΔΔ 47 : 0.163‰, ΔΔ 48 : -0.761‰). Surface floe samples begin to precipitate out of solution within the first hour of CO 2 absorption, and the dissolved inorganic carbon (DIC) pool requires a residence time of >100 h to achieve isotopic equilibria. The Δ 48 /Δ 47 slope of samples from the Cedars (-3.223 ± 0.519) is within the range of published theoretical values designed to constrain CO 2 hydrolysis-related kinetic fractionation (-1.724 to -8.330). The Δ 47 /δ 18 O slope (-0.009 ± 0.001) and Δ 47 /δ 13 C slope (-0.009 ± 0.001) are roughly consistent with literature values reported from a peridotite in Oman of -0.006 ± 0.002 and -0.005 ± 0.002, respectively. The consistency of slopes in the multi-isotope space suggests the Δ 47 -Δ 48 dual carbonate clumped isotope framework can be applied to study CO 2 -absorption processes in applied systems, including sites of interest for geological sequestration.

58 GEOSCIENCES↗
DOE OSTI2023

Role of grain boundary crystallography on void growth in FCC metals

The growth of an isolated spherical void at and near a grain boundary or a triple junction is studied by means of full-field crystal plasticity simulations with explicit representation of the void. To examine the fundamental aspects of void/grain boundary interactions in these cases, a three-dimensional large-strain elasto-viscoplastic fast-Fourier transform (LS-EVP-FFT) model with axisymmetric tensile periodic boundary conditions was developed. We employ this model to reveal the role of crystallography (grain orientations, grain boundary misorientation and inclination angle with respect to the largest component of the applied stress) and type of loading on void growth. It was found that grain orientation was the primary crystallographic feature governing void growth in strain-rate controlled simulations. For strains small enough to not significantly alter the local geometry, the overall growth of the void in bicrystal boundaries was not strongly affected by grain boundary inclination angle. The relative amount of void growth experienced within each grain in a bicrystal or tri-crystal was strongly dependent on the grain boundary inclination angle. We found a subset of orientations that exhibit faster void growth in bicrystal boundaries than inside the constituent single crystals. Last, we show that when placed primarily in one grain but near the grain boundary, void growth was similar to void growth when placed at the grain boundary. Furthermore, the growing void remained in the original grain and did not grow into the other grain, regardless of crystallographic character.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Study of growth twins produced through heat treatment of fine-grained magnesium

This study investigated the microstructural features of magnesium after the fine-grained magnesium sample was heat treated at 400 °C for 2 h. Both fine-grained and heat-treated magnesium samples were observed using the electron backscatter diffraction technique in scanning electron microscopy. The results show that the grains were substantially coarsened, the local strain and the local misorientation was significantly reduced, and a large number of twins were observed in the heat-treated magnesium sample. The misorientation angle analysis determined that the twins were oriented about 86.3° away from the orientation of their corresponding matrices, and the twins were the variants of $\{1 0 \bar{1} 2\}$ tensile twins. The relative orientation analysis also identified that the fast growth directions of the twins were closely along the $<2 \bar{1} \bar{1} 0>$ or $<1 0 \bar{1} 0>$ directions in the basal plane. The twins can merge or impinge onto each other. Most boundaries between the twins have the misorientation angle of about 7.4°.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of the temperature-dependent failure processes in PVD Cr-coated ZIRLO nuclear fuel cladding using in situ X-ray micro-tomography imaging

In this work, an accident tolerant fuel cladding system with Cr coating deposited using physical vapour deposition (PVD) method onto commercial Optimized ZIRLO™ was studied. The cladding tubes were machined into C-rings and loaded to failure under compression with real-time synchrotron micro-X-ray computed tomography (XCT) imaging at room temperature (RT), 345 °C, 650 °C and 950 °C in argon (Ar) atmosphere. The mechanical behaviour and failure processes were found strongly temperature-dependent where the Cr coating showed brittle fracture at RT and 345 °C, ductile fracture at 650 °C and a reversion to brittle fracture at 950 °C. Nanoindentation measurements and scanning electron microscopy (SEM) imaging were conducted on the materials after high temperature testing. It was found that recrystallisation of the Cr coating occurred at elevated temperatures, which significantly affected its local properties hence the failure behaviour at different temperatures. This work represents the first in situ 3D XCT observation of progressive failure processes in PVD Cr-coated ZIRLO claddings up to 950 °C providing critical insights into its brittle-to-ductile transition (BDTT) behaviour and subsequent ductile-to-brittle reversion with the increase in temperature. Moreover, results are compared with other PVD Cr-coated Zircaloy materials from open literature; the influences of the coating microstructure and local properties on the failure stress/strain and fracture processes are discussed.

Yuan, Guanjie↗
DOE OSTI2024

Analysis of the Structure of Heavy Ion Irradiated LaFeO 3 Using Grazing Angle X-ray Absorption Spectroscopy

Crystalline ceramics are candidate materials for the immobilization of radionuclides, particularly transuranics (such as U, Pu, and Am), arising from the nuclear fuel cycle. Due to the α-decay of transuranics and the associated recoil of the parent nucleus, crystalline materials may eventually be rendered amorphous through changes to the crystal lattice caused by these recoil events. Previous work has shown irradiation of titanate-based ceramics to change the local cation environment significantly, particularly in the case of Ti which was shown to change from 6- to 5-fold coordination. Here, this work expands the Ti-based study to investigate the behavior in Fe-based materials, using LaFeO 3 as an example material. Irradiation was simulated by heavy ion implantation of the bulk LaFeO 3 ceramic, with the resulting amorphous layer characterized with grazing angle X-ray absorption spectroscopy (GA-XAS). Insights into the Fe speciation changes exhibited by the amorphized surface layer were provided through quantitative analysis, including pre-edge analysis, and modeling of the extended X-ray absorption fine structure (EXAFS), of the GA-XAS data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Layer-Dependent Magnetism and Spin Fluctuations in Atomically Thin van der Waals Magnet CrPS 4

van der Waals (vdW) magnets, an emerging family of two-dimensional (2D) materials, have received tremendous attention due to their rich fundamental physics and significant potential for cutting-edge technological applications. In contrast to the conventional bulk counterparts, vdW magnets exhibit significant tunability of local material properties, such as stacking engineered interlayer coupling and layer-number dependent magnetic and electronic interactions, which promise to deliver previously unavailable merits to develop multifunctional microelectronic devices. As a further ingredient of this emerging topic, here we report nanoscale quantum sensing and imaging of atomically thin vdW magnet chromium thiophosphate CrPS 4 , revealing its characteristic layer-dependent 2D static magnetism and dynamic spin fluctuations. We also show a large tunneling magnetoresistance in CrPS 4 -based spin filter vdW heterostructures. Here, the excellent material stability, robust strategy against environmental degradation, in combination with tailored magnetic properties highlight the potential of CrPS 4 in developing state-of-the-art 2D spintronic devices for next-generation information technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High-Mass Loading of Flower-like Ni-MoS 2 microspheres Toward Efficient Intercalation pseudocapacitive Electrode

This work reports the exploration of intercalation pseudocapacitance in a thicker electrode of flowerlike Ni-doped MoS 2 microspheres that features a mass loading of ~10 mg/cm 2 without sacrificing the gravimetric capacitance (~425 F/g at 5 mV/s). Integration of Ni atoms into MoS 2 microspheres not only stabilized the structural integrity but also ameliorated the rapid intercalation and deintercalation of electrolyte ions even at a commercial-level mass loading. The energy instability by Ni doping significantly changed the local bonding behavior and the overall electronic structure of MoS 2 , facilitating the breaking of the MoS 2 layer and generation of more active edge sites, which are responsible for faster reaction kinetics. The experiments attribute the overall capacitance enhancement in (Mo-Ni)S 2 to the increased rate of electrolyte ion insertion and extraction, which is confirmed by b-values close to 0.5, at different potentials, indicating that the current response predominantly depends on the diffusive mechanism for both MoS 2 and Ni-MoS 2 thicker electrodes. The symmetric device constructed with Ni-MoS 2 microspheres exhibited a capacitance value of 101 F/g in 1 mV/s, for which the energy density is 9 Wh/kg, as well as attained an outstanding cycling stability of 10000 cycles with 60% retention at 2 A/g. In addition to providing insights into the development of 2D TMDs, this work explores the design of robust and highly efficient intercalation electrode material for electrochemical energy storage devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Importance of Van der Waals Interactions in Hydrogen Adsorption on a Silicon-carbide Nanotube Revisited with vdW-DFT and Quantum Monte Carlo

Density functional theory (DFT) is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the local/semi-local theory, it remains unclear as to how the consideration of van der Waals (vdW) interactions affects such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H 2 -SiCNT). Within the DFT framework, we benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H 2 -SiCNT geometry within the DFT framework and investigated the correlation between the structure and charge density. The vdW contribution to the adsorption was found to be non-negligible at similar to 1 kcal/mol per hydrogen molecule, which amounts to 9-29% of the ideal adsorption energy required for hydrogen storage applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Strong Ferromagnetic Exchange Coupling and Single-Molecule Magnetism in MoS 4 3– -Bridged Dilanthanide Complexes

We report the synthesis and characterization of the trinuclear 4d-4f compounds [Co(C 5 Me 5 ) 2 ][(C 5 Me 5 ) 2 Ln(μ-S) 2 Mo(μ-S) 2 Ln(C 5 Me 5 ) 2 ], 1-Ln (Ln = Y, Gd, Tb, Dy), containing the highly polarizable MoS 4 3 - bridging unit. UV-Vis-NIR diffuse reflectance spectra and DFT calculations of 1-Ln reveal a low-energy metal-to-metal charge transfer transition assigned to charge transfer from the singly occupied 4d z 2 orbital of Mo V to the empty 5d orbitals of the lanthanides (4d in the case of 1-Y ), mediated by sulfur-based 3p orbitals. Electron paramagnetic resonance spectra collected for 1-Y in a tetrahydrofuran solution show large 89 Y hyperfine coupling constants of A ⊥ = 23 MHz and A || = 26 MHz, indicating the presence of significant yttrium-localized unpaired electron density. Magnetic susceptibility data support similar electron delocalization and ferromagnetic Ln-Mo exchange for 1-Gd , 1-Tb , and 1-Dy . This ferromagnetic exchange gives rise to an S = 15/2 ground state for 1-Gd and one of the largest magnetic exchange constants involving Gd III observed to date, with J Gd-Mo = +16.1(2) cm -1 . Additional characterization of 1-Tb and 1-Dy by ac magnetic susceptibility measurements reveals that both compounds exhibit slow magnetic relaxation. Although a Raman magnetic relaxation process is dominant for both 1-Tb and 1-Dy , an extracted thermal relaxation barrier of U eff = 68 cm -1 for 1-Dy is the largest yet reported for a complex containing a paramagnetic 4d metal center. Furthermore, these results provide a potentially generalizable route to enhanced n d-4f magnetic exchange, revealing opportunities for the design of new n d-4f single-molecule magnets and bulk magnetic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Seasonal Seismicity in the Lake Biwa Region of Central Japan Moderately Modulated by Lake Water Storage Changes

Seasonal seismicity in Japan has been observed and attributed to variations induced by hydrological loading, as well as atmospheric and pore-fluid pressure. We examine the seasonal variations of crustal (<15 km depth) and low magnitude (M < 4.5) earthquakes in central Japan near the Biwako-Seigan Fault Zone (BSFZ), and analyze their correlation with the annual hydrological cycle of nearby Lake Biwa, the largest freshwater lake in Japan. The 2002–2018 lake water storage is estimated using water-level gauge and satellite altimetry data. Regional surface mass loading from soil moisture and snow predicted by the Global Land Data Assimilation System is used to augment the lake mass variation. We compute the stress and stressing-rate changes induced by the surface mass loading from lake storage, snow, and soil moisture at seismogenic depth (10 km) along the BSFZ using a finite element model incorporating a 3D Earth structure. The results show that seasonal surface mass loading generates 1–5 kPa peak-to-peak Coulomb stress change and 2.5–15 kPa/year stressing-rate variations on the fault segments along the BSFZ. Two of the four fault segments present significant correlation between local seismicity rate and loading-derived stress and/or stressing rate changes. Among the hydrological load components, lake water, snow, and soil moisture exhibit different annual phases, and the lake water storage dominates the total annual stress changes. Our analysis indicates that seasonal surface hydrological loads, especially lake water variations, modulate the seismicity rate along the BSFZ in central Japan.

58 GEOSCIENCES↗