Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 865 records · Page 48

The Green Sahara: Climate Change, Hydrologic History and Human Occupation

Archaeology can provide insight into interactions of climate change and human activities in sensitive areas such as the Sahara, to the benefit of both disciplines. Such analyses can help set bounds on climate change projections, perhaps identify elements of tipping points, and provide constraints on models. The opportunity exists to more precisely constrain the relationship of natural solar and climate interactions, improving understanding of present and future anthropogenic forcing. We are beginning to explore the relationship of human occupation of the Sahara and long-term solar irradiance variations synergetic with changes in atmospheric-ocean circulation patterns. Archaeological and climate records for the last 12 K years are gaining adequate precision to make such comparisons possible. We employ a range of climate records taken over the globe (e.g. Antarctica, Greenland, Cariaco Basin, West African Ocean cores, records from caves) to identify the timing and spatial patterns affecting Saharan climate to compare with archaeological records. We see correlation in changing ocean temperature patterns approx. contemporaneous with drying of the Sahara approx. 6K years BP. The role of radar images and other remote sensing in this work includes providing a geographically comprehensive geomorphic overview of this key area. Such coverage is becoming available from the Japanese PALSAR radar system, which can guide field work to collect archaeological and climatic data to further constrain the climate change chronology and link to models. Our initial remote sensing efforts concentrate on the Gilf Kebir area of Egypt.

remote sensing↗

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗

A closer look at how symmetry constraints and the spin–orbit coupling shape the electronic structure of Bi(111)

Abstract Fully relativistic density-functional-theory calculations of Bi(111) thin films are analyzed to revisit their two metallic surface-states branches. We first contrast these metallic branches with surface states arising at gaps in the valence band opened by the spin–orbit coupling (SOC). We find that the two metallic branches along Γ M ‾ do not overlap with the bulk band at the zone boundary,M. We show that the spin texture observed in such states cannot be traced to the lifting of Kramers’ degeneracy. Instead, we track them to the m j = ± 1 / 2 – m j = ± 3 / 2 SOC splitting, the potential anisotropy for in-plane and out-of-plane states, and the coupling between the opposite surfaces of a slab occurring nearM, which is driven by a spatial redistribution of the four metallic states composing the two metallic branches. Each of these branches appears to be non-degenerate at the tested surface, yet each is degenerate with another state of opposite spin at the other surface. Nevertheless, the four metallic states bear some contribution on both surfaces of the film because of their spatial redistribution nearM. The overlapping among these states nearM, afforded by their spatial redistribution on both surfaces, causes a hybridization that perpetuates the splitting between the two branches, makes the film’s electronic structure thickness dependent nearM, extinguishes the magnetic moment of the metallic states avoiding the magnetic-moment discontinuity atM, and denies the need or expectancy of the metallic branches becoming degenerate atM. We propose that theoppositespin polarization observed for the two metallic branches occurs because the surface atoms retain their covalent bonds and thus cannot afford magnetic polarization. We show that the Rashba-splitting of the metallic states for inversion-asymmetric films does not have a fixed magnitude but can be tuned by changing the perturbation breaking inversion symmetry.

Physics↗

Comparing multi-source urban flood indicators: satellite, simulation, and citizen-reported data

Urban flooding arises from complex mechanisms, making it challenging to capture accurately with a single detection method. This study evaluates three complementary approaches to detect flooding across three Chicago neighborhoods: (i) Sentinel-1 synthetic aperture radar (SAR), offering weather-independent, high-resolution (10 m) imagery of surface inundation; (ii) the storm water management model (SWMM), simulating combined sewer overflow and drainage performance; and (iii) citizen-generated 311 service requests, capturing observed flooding impacts. By analyzing six storms ranging from severe to mild, we examine how each source uniquely contributes to identifying urban flood events. SAR imagery effectively identifies standing water but can miss brief flooding due to satellite revisit constraints. SWMM provides detailed insights into system-wide drainage behavior yet may underestimate localized street-level flooding. Meanwhile, 311 calls reflect real-world flooding impacts but are vulnerable to underreporting. Statistical overlap analysis highlights chronic flood hotspots repeatedly identified across multiple detection methods, indicating persistent infrastructure and topographic vulnerabilities. Temporal analysis further reveals that while SWMM flooding aligns closely with rainfall peaks, 311 calls typically precede or persist beyond these peaks. Our findings emphasize the value of using satellite observations, hydrological modeling, and resident-reported data in a complementary manner to better interpret patterns in flood timing, severity, and spatial distribution—providing insights that can inform targeted infrastructure improvements and contribute to urban flood resilience planning.

311↗

Delays of optical bursts in simultaneous optical and X-ray observations of MXB 1636-53

Observations of simultaneous optical and X-ray bursts from 4U/MXB 1636-53 were made using the Hakucho burst monitor system and optical telescopes at the European Southern Observatory during 1979 and 1980. The six best cases among the 10 coinciding observations are analyzed in terms of a model in which the optical emission is the result of reprocessing of X-rays (through blackbody heating). From this analysis, the temperature (spatially averaged) and size of a reprocessor, and the smearing and delay of the optical bursts are obtained. For the maximum temperatures of the optical reprocessor, the values differ from burst to burst, ranging from about 3 x 10 to the 4th to about 10 to the 5th K. The present analysis suggests that the size of the reprocessor varies by a factor of a few. For the smearing of the optical bursts an upper limit of a few seconds is derived. The most important result of this analysis is that the delay times are not the same for all bursts. The possible constraints which these results put on a low-mass binary model of this burst source are discussed.

Matsuoka, M.↗

Constraints on Decreases in Eta Carinae's Mass-loss from 3D Hydrodynamic Simulations of Its Binary Colliding Winds

Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Car's colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 × 10(exp −4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-B's wind, plus radiative inhibition of Eta-B's wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-A's spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Car's next periastron in 2014 to further test for decreases in dot-M Eta-A. If dot-M Eta-A is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Eta Carinae↗

The AWAKEN wind farm benchmark, Part 2: Modeling results

Accurately modeling wind farm performance in complex atmospheric flows remains a challenge. This paper presents the modeling results of the American WAKE experimeNt (AWAKEN) wind farm benchmark, a collaborative effort involving 16 research groups from academia and industry within the International Energy Agency Wind Technology Collaboration Programme Task 57. The study evaluates a diverse suite of simulation tools, ranging from fast-running engineering wake models to high-fidelity large-eddy simulations, against a diurnal case study observed during the AWAKEN campaign. The benchmark utilized a three-phase structure to progressively assess model performance as observational data availability increased. Initial blind predictions showed that higher-fidelity models did not uniformly outperform simpler simulation tools. A distinct spatial bias was observed where models struggled to resolve the interplay between a low-level jet, wakes, and terrain-induced flow acceleration. In subsequent phases, leveraging additional measurements for model improvement led to a reduction in mean absolute error across the model ensemble; however, this effect was most pronounced in engineering wake models, where targeted calibration reduced error by up to 40~\%. Overall, the study demonstrates that inflow characterization remains a primary prerequisite for accuracy, particularly for models relying on coarse forcing datasets. While the limited ability to resolve local terrain-flow interactions under single-day conditions represent a recognized constraint, the overall findings on wake modeling and real-world validation still provide valuable guidance for model application and for mitigating this limitation.

Bodini, Nicola↗

Performance of customized DCT quantization tables on scientific data

We show that it is desirable to use data-specific or customized quantization tables for scaling the spatial frequency coefficients obtained using the Discrete Cosine Transform (DCT). DCT is widely used for image and video compression (MP89, PM93) but applications typically use default quantization matrices. Using actual scientific data gathered from divers sources such as spacecrafts and electron-microscopes, we show that the default compression/quality tradeoffs can be significantly improved upon by using customized tables. We also show that significant improvements are possible for the standard test images Lena and Baboon. This work is part of an effort to develop a practical scheme for optimizing quantization matrices for any given image or video stream, under any given quality or compression constraints.

Ratnakar, Viresh↗

Lunar and Planetary Science XXXVI, Part 13

Contents include the following: A Fast, Non-Destructive Method for Classifying Ordinary Chondrite Falls Using Density and Magnetic Susceptibility. An Update on Results from the Magnetic Properties Experiments on the Mars Exploration Rovers, Spirit and Opportunity. Measurement Protocols for In Situ Analysis of Organic Compounds at Mars and Comets. Piping Structures on Earth and Possibly Mars: Astrobiological Implications. Uranium and Lead in the Early Planetary Core Formation: New Insights Given by High Pressure and Temperature Experiments. The Mast Cameras and Mars Descent Imager (MARDI) for the 2009 Mars Science Laboratory. MGS MOC: First Views of Mars at Sub-Meter Resolution from Orbit. Analysis of Candor Chasma Interior Layered Deposits from OMEGA/MEX Spectra. Analysis of Valley Networks on Valles Marineris Plateau Using HRSC/MEX Data. Solar Abundance of Elements from Neutron-Capture Cross Sections. Preliminary Evaluation of the Secondary Ion/Accelerator Mass Spectrometer, MegaSIMS. Equilibrium Landforms in the Dry Valleys of Antarctica: Implications for Landscape Evolution and Climate Change on Mars. Continued Study of Ba Isotopic Compositions of Presolar Silicon Carbide Grains from Supernovae. Paleoenviromental Evolution of the Holden-Uzboi Area. Stability of Magnesium Sulfate Minerals in Martian Environments. Tungsten Isotopic Constraints on the Formation and Evolution of Iron Meteorite Parent Bodies. Migration of Dust Particles and Volatiles Delivery to the Inner Planets. On the Sitting of Trapped Noble Gases in Insoluble Organic Matter of Primitive Meteorites. Trapping of Xenon Upon Evaporation-Condensation of Organic Matter Under UV Irradiation: Isotopic Fractionation and Electron Paramagnetic Resonance Analysis. Stability of Water on Mars. A Didactic Activity. Analysis of Coronae in the Parga Chasma Region, Venus. Photometric and Compositional Surface Properties of the Gusev Crater Region, Mars, as Derived from Multi-Angle, Multi-Spectral Investigation of Mars Express HRSC Data. Mapping Compositional Diversity on Mars: Spatial Distribution and Geological Implications. A New Simulation Chamber for Studying Planetary Environments. Folded Structure in Terra Sirenum. Mars. Nitrogen-Noble Gas Static Mass Spectrometry of Genesis Collector Materials. Neon Isotope Heterogeneity in the Terrestrial Mantle: Implication for the Acquisition of Volatile Elements in Terrestrial Planets. The Cosmic Clock, the Cycle of Terrestrial Mass Extinctions.

Source record↗

Ground-based Network and Supersite Measurements for Studying Aerosol Properties and Aerosol-Cloud Interactions

From radiometric principles, it is expected that the retrieved properties of extensive aerosols and clouds from reflected/emitted measurements by satellite (and/or aircraft) should be consistent with those retrieved from transmitted/emitted radiance observed at the surface. Although space-borne remote sensing observations contain large spatial domain, they are often plagued by contamination of surface signatures. Thus, ground-based in-situ and remote-sensing measurements, where signals come directly from atmospheric constituents, the sun, and the Earth-atmosphere interactions, provide additional information content for comparisons that confirm quantitatively the usefulness of the integrated surface, aircraft, and satellite datasets. The development and deployment of AERONET (AErosol RObotic NETwork) sunphotometer network and SMART-COMMIT (Surface-sensing Measurements for Atmospheric Radiative Transfer - Chemical, Optical & Microphysical Measurements of In-situ Troposphere) mobile supersite are aimed for the optimal utilization of collocated ground-based observations as constraints to yield higher fidelity satellite retrievals and to determine any sampling bias due to target conditions. To characterize the regional natural and anthropogenic aerosols, AERONET is an internationally federated network of unique sunphotometry that contains more than 250 permanent sites worldwide. Since 1993, there are more than 480 million aerosol optical depth observations and about 15 sites have continuous records longer than 10 years for annual/seasonal trend analyses. To quantify the energetics of the surface-atmosphere system and the atmospheric processes, SMART-COMMIT instrument into three categories: flux radiometer, radiance sensor and in-situ probe. Through participation in many satellite remote-sensing/retrieval and validation projects over eight years, SMART-COMMIT have gradually refine( and been proven vital for field deployment. In this paper, we will demonstrate the capability of AERONET SMART-COMMIT in current Asian Monsoon Year-2008 campaigns that are designed and being executed to study the compelling variability in temporal scale of both anthropogenic and natural aerosols (e.g., airborne dust, smoke, mega-city pollutant). Feedback mechanisms between aerosol radiative effects and monsoon dynamics have been recently proposed, however there is a lack of consensus on whether aerosol forcing would be more likely to enhance or reduce the strength of the monsoon circulation. We envision robust approaches which well-collocated ground-based measurements and space-borne observations will greatly advance our understanding of absorbing aerosols (e.g., "Global Dimming" vs. "Elevated Heat-Pump" effects) on aerosol cloud water cycle interactions.

Tsay, Si-Chee↗

Management of Operations under Visual Flight Rules in UTM for Disaster Response Missions

The disaster response domain has experienced an increased focus in recent years due to the rise in number and scale of events, lessons learned from past experience, and emerging technologies that make possible a more coordinated and effective response. As part of this focus, JAXA and NASA have been collaborating on the integration of manned and unmanned aircraft in support of disaster response operations through integrated testing of their respective mission planning and optimization system (Disaster Relief Aircraft Information Sharing Network, or D-NET) and an automated UAS traffic management system (e.g., UTM). In 2018, JAXA and NASA jointly participated in a large-scale disaster drill in Japan where the integration of systems was successfully demonstrated through real-time data exchanges, visualization, and decision making as part of the coordinated airspace management of a manned helicopter in VFR conditions and unmanned small UAS operating in common areas. This work details a flight test consisting of two flights that were conducted December 2019 near the Chofu Aerodrome in Tokyo, which focused on the evaluation of pilots operating under Visual Flight Rules (VFR) communicating through D-NET and sharing intent and position information within UTM. This work contributes to defining the necessary requirements for digital coordination between manned and unmanned operations. UTM requires the use of operation volumes, which are spatial and temporal volumes that encompass UAS flight trajectories and account for technical performance errors and deviations due to disturbances (e.g., wind). A series of flights, representing different missions, used landmark-based operation volumes and conformance of the aircraft to those operation volumes were tracked within UTM. Experienced disaster response helicopter pilots provided insight on the development of the operation plans and their usability during disaster response operations. Results from the flight test supported the suggested benefits of using landmarks for planning and positional awareness and highlighted the need for future research in advanced visualization capabilities to support operations that consider both system constraints and flight deck/airspace management interaction.

disaster response↗

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.°2 away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on 2021 December 8. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV γ -ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the γ -ray data from Fermi-LAT, VERITAS, and H.E.S.S. require a spectral cutoff near 100 GeV. Both the X-ray and γ -ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed γ -ray spectral cutoff in both the leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power.

79 ASTRONOMY AND ASTROPHYSICS↗

Electromagnetic deep-probing (100-1000 kms) of the Earth's interior from artificial satellites: Constraints on the regional emplacement of crustal resources

A spherical harmonic analysis program is being tested which takes magnetic data in universal time from a set of arbitrarily space observatories and calculates a value for the instantaneous magnetic field at any point on the globe. The calculation is done as a least mean-squares value fit to a set of spherical harmonics up to any desired order. The program accepts as a set of input the orbit position of a satellite coordinates it with ground-based magnetic data for a given time. The output is a predicted time series for the magnetic field on the Earth's surface at the (r, theta) position directly under the hypothetically orbiting satellite for the duration of the time period of the input data set. By tracking the surface magnetic field beneath the satellite, narrow-band averages crosspowers between the spatially coordinated satellite and the ground-based data sets are computed. These crosspowers are used to calculate field transfer coefficients with minimum noise distortion. The application of this technique to calculating the vector response function W is discussed.

Hermance, J. F.↗

Comparison of temporal resolution selection approaches in energy systems models

Capacity expansion models for the power sector are used to project future decisions over the coming decades by simulating investment and operation decisions for the use of electricity. Due to model performance constraints, these models typically do not explicitly simulate every hour within a year, but instead simulate representative time segments (groups of hours). This paper evaluates different approaches for selecting time segments across three methods: sequential, categorical, and clustering, across a wide range of time-segment quantities, for a total of 204 temporal profiles. To measure the performance of each profile's ability to accurately represent data, the root-mean-square-error of each profile's time segments are compared to the data's original hourly data. The temporal alignment across regions is also measured (i.e., how often windy days align across regions). Different spatial resolutions were applied for a subset of the temporal selection methods to investigate the impact spatial resolution has on performance. This paper provides a framework for measuring the value of different temporal selection methods and of adding more granular data to energy system models. Overall, multi-criteria clustering yields the lowest root-mean-square-error across all datasets evaluated and provides a holistic view of the intertwined relationships between renewable generation and electricity demand.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Optimizing the location and configuration of disaster resilience hubs under transportation and electric power network failures

Natural disasters often result in failures of transportation network components and blackouts that imperil the wellbeing of vulnerable populations. In response to these events, resilience hubs have been proposed as a pre-disaster planning strategy to improve access to critical services. This paper introduces an optimization-based approach to locate and configure electric power-generating resilience hubs considering the possibility of failures in transportation and electric power systems. The model's objective is to identify hub locations and configurations that maximize transportation accessibility to the hubs and maximize the satisfaction of basic energy needs through hub-generated electric power. Besides a budget constraint, the model accounts for limits on the levels of hub energy generation vis-à-vis community energy demands, and on the transportation network distance of communities to hubs. Three heuristics are presented for the proposed planning problem. The first heuristic is a genetic algorithm (GA) with problem-specific solution generation procedures. The other two heuristics implement greedy search techniques. Numerical experiments were conducted, using data from rural Puerto Rico, to illustrate the application of the proposed model and heuristics, and examine their performance. In the numerical experiments, the GA heuristic found better solutions than the greedy heuristics. Additionally, design solutions consisting of spatially dispersed hubs with low energy generation capacity were better than solutions with spatially concentrated high-capacity hubs. Lastly, across a wide range of hub demand scenarios, only a small number of candidate hub locations consistently ranked among the best locations for establishing a hub.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

On the Diffusivity of Moist Static Energy and Implications for the Polar Amplification Response to Climate Warming

Energy balance models (EBMs) have been widely used in a range of climate problems, but the assumption of constant diffusivity in the parameterization of the moist static energy (MSE) flux can be hardly justified. We demonstrate in this study that the diffusive MSE flux can be derived from the basic energy balance equation with a few tolerable assumptions. The estimated diffusivity is both spatially and seasonally dependent, and its midlatitude average is then tested against several scaling theories for the midlatitude eddy diffusivity. The result supports the diffusivity theory of Held and Larichev (1996) modified for the moist atmosphere, affording a dynamics-based parameterization of MSE diffusivity. The implementation of the parameterization in an EBM leads to an interactive MSE diffusivity that accounts for the midlatitude eddy response to climate forcing perturbations. Under a uniform radiative forcing, the EBM with a diffusivity so parameterized produces a weakening of the midlatitude diffusivity and a modestly polar-amplified surface temperature response as an inevitable outcome under the dual constraints of the nonlinear Clausius-Clapeyron relation and the temperature gradient-dependent diffusivity, even in the absence of any poleward amplifying radiative feedbacks. As the consequence of more isothermal temperature and reduced diffusivity, the variance of the midlatitude surface temperature also decreases with warming.

54 ENVIRONMENTAL SCIENCES↗

Big Data Meets Geothermal Exploration (CRADA Final Report)

As part of the Cyclotron Road program, Zanskar Geothermal & Minerals, Inc. investigated the application of micro-earthquake and ambient noise seismology methods to imaging and characterizing the structural characteristics and hydrothermal flux of subsurface faults. Significant advances in what could be resolved were enabled by two major developments in seismology: 1) the availability of large-n arrays of low-cost seismometers, and 2) the availability of increased computational power and semi-automated data reduction algorithms. In tandem, these advances may improve the signal-to-noise ratio and spatial precision of the data collected and enable higher-resolution characterization of subsurface fracture systems and their spatio-temporal evolution. These tools supported efforts to reduce dry-hole risk and to improve wellfield productivity for geothermal resource development. In particular, two applications of these advances were evaluated: 1) fracture-seismic imaging, which was used to detect ambient emissions from fluid-filled fractures, and 2) reservoir tomography, which used information about travel paths, source locations, and source parameters of micro-earthquakes to identify areas of enhanced permeability. Integration of these methods provided guidance for siting wells and served as prior constraints for reservoir models, informing forecasts of power potential and production and injection strategies aimed at minimizing temperature decline and improving overall resource productivity.

15 GEOTHERMAL ENERGY↗

A linear method for analyzing lightning field changes

A constrained, least-squares method for analyzing multiple-station measurements of lightning field changes (delta Es) is introduced. Previous methods have attempted to fit the spatial pattern of lightning delta Es using nonlinear models, such as a point charge (Q) or a point dipole (P) model. With the linear method, the delta Es are described not by models but by a general volume charge distribution that is deposited on a large (40 x 40 x 20 cu km) Cartesian grid above the measuring network. A linear system of equations is used to relate the measured delta Es to the charges that are deposited at each grid point. With this approach, the information content of the measurements can be quantified by an eigenanalysis of the covariance matrix of the linear system. Constraints can be used to reduce the infinity of possible solutions to the linear system and also to reduce systematic biases that can be introduced by the method of solution. It is shown that a Landweber iterative method, derived from the general method of steepest descent, can be used to solve the linear system and that the resulting volume charge distributions are generally consistent with computer-simulated charge sources, when these sources are over the measuring network. The Landweber iteration has also provided solutions for natural lightning events that are consistent with Q- and P-model results.

Koshak, William J.↗