Search NASA⌕ Search

SEARCH · Search NASA

Results for “zone refining”

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

Summary of Data from the Sixth AIAA CFD Drag Prediction Workshop: CRM Cases 2 to 5

Results from the Sixth AIAA CFD Drag Prediction Workshop Common Research Model Cases 2 to 5 are presented. As with past workshops, numerical calculations are performed using industry-relevant geometry, methodology, and test cases. Cases 2 to 5 focused on force/moment and pressure predictions for the NASA Common Research Model wing-body and wing-body-nacelle-pylon configurations, including Case 2 - a grid refinement study and nacelle-pylon drag increment prediction study; Case 3 - an angle-of-attack buffet study; Case 4 - an optional wing-body grid adaption study; and Case 5 - an optional wing-body coupled aero-structural simulation. The Common Research Model geometry differed from previous workshops in that it was deformed to the appropriate static aeroelastic twist and deflection at each specified angle-of-attack. The grid refinement study used a common set of overset and unstructured grids, as well as user created Multiblock structured, unstructured, and Cartesian based grids. For the supplied common grids, six levels of refinement were created resulting in grids ranging from 7x10(exp 6) to 208x10(exp 6) cells. This study (Case 2) showed further reduced scatter from previous workshops, and very good prediction of the nacelle-pylon drag increment. Case 3 studied buffet onset at M=0.85 using the Medium grid (20 to 40x10(exp 6) nodes) from the above described sequence. The prescribed alpha sweep used finely spaced intervals through the zone where wing separation was expected to begin. Although the use of the prescribed aeroelastic twist and deflection at each angle-of-attack greatly improved the wing pressure distribution agreement with test data, many solutions still exhibited premature flow separation. The remaining solutions exhibited a significant spread of lift and pitching moment at each angle-of-attack, much of which can be attributed to excessive aft pressure loading and shock location variation. Four Case 4 grid adaption solutions were submitted. Starting with grids less than 2x10(exp 6) grid points, two solutions showed a rapid convergence to an acceptable solution. Four Case 5 coupled aerostructural solutions were submitted. Both showed good agreement with experimental data. Results from this workshop highlight the continuing need for CFD improvement, particularly for conditions with significant flow separation. These comparisons also suggest the need for improved experimental diagnostics to guide future CFD development.

Tinoco, Edward N.↗
NASA NTRS2021

Reviews and syntheses: Arctic fire regimes and emissions in the 21st century

In recent years, the pan-Arctic region has experienced increasingly extreme fire seasons. Fires in the northern high latitudes are driven by current and future climate change, lightning, fuel conditions, and human activity. In this context, conceptualizing and parameterizing current and future Arctic fire regimes will be important for fire and land management as well as understanding current and predicting future fire emissions. The objectives of this review were driven by policy questions identified by the Arctic Monitoring and Assessment Programme (AMAP) Working Group and posed to its Expert Group on Short-Lived Climate Forcers. This review synthesizes current understanding of the changing Arctic and boreal fire regimes, particularly as fire activity and its response to future climate change in the pan-Arctic have consequences for Arctic Council states aiming to mitigate and adapt to climate change in the north. The conclusions from our synthesis are the following. (1) Current and future Arctic fires, and the adjacent boreal region, are driven by natural (i.e. lightning) and human-caused ignition sources, including fires caused by timber and energy extraction, prescribed burning for landscape management, and tourism activities. Little is published in the scientific literature about cultural burning by Indigenous populations across the pan-Arctic, and questions remain on the source of ignitions above 70_ N in Arctic Russia. (2) Climate change is expected to make Arctic fires more likely by increasing the likelihood of extreme fire weather, increased lightning activity, and drier vegetative and ground fuel conditions. (3) To some extent, shifting agricultural land use and forest transitions from forest–steppe to steppe, tundra to taiga, and coniferous to deciduous in a warmer climate may increase and decrease open biomass burning, depending on land use in addition to climate-driven biome shifts. However, at the country and landscape scales, these relationships are not well established. (4) Current black carbon and PM2:5 emissions from wildfires above 50 and 65_ N are larger than emissions from heanthropogenic sectors of residential combustion, transportation, and flaring. Wildfire emissions have increased from 2010 to 2020, particularly above 60_ N, with 56% of black carbon emissions above 65_ N in 2020 attributed to open biomass burning – indicating how extreme the 2020 wildfire season was and how severe future Arctic wildfire seasons can potentially be. (5) What works in the boreal zones to prevent and fight wildfires may not work in the Arctic. Fire management will need to adapt to a changing climate, economic development, the Indigenous and local communities, and fragile northern ecosystems, including permafrost and peatlands. (6) Factors contributing to the uncertainty of predicting and quantifying future Arctic fire regimes include underestimation of Arctic fires by satellite systems, lack of agreement between Earth observations and official statistics, and still needed refinements of location, conditions, and previous fire return intervals on peat and permafrost landscapes. This review highlights that much research is needed in order to understand the local and regional impacts of the changing Arctic fire regime on emissions and the global climate, ecosystems, and pan-Arctic communities.

Arctic↗
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗
DOE OSTI2021

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , a New Mineral Species from the Sisim Placer Zone, Eastern Sayans, Russia, and the Torryweiserite–Ferrotorryweiserite Series

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , occurs as small grains (≤20 µm) among droplet-like inclusions (up to 50 μm in diameter) of platinum-group minerals (PGM), in association with oberthürite or Rh-bearing pentlandite, laurite, and a Pt-Pd-Fe alloy (likely isoferroplatinum and Fe-Pd-enriched platinum), hosted by placer grains of Os-Ir alloy (≤0.5 mm) in the River Ko deposit. The latter is a part of the Sisim placer zone, which is likely derived from ultramafic units of the Lysanskiy layered complex, southern Krasnoyarskiy kray, Russia. The mineral is opaque, gray to brownish gray in reflected light, very weakly bireflectant, not pleochroic to weakly pleochroic (grayish to light brown tints), and weakly anisotropic. The calculated density is 5.93 g·cm –3 . Mean results (and ranges) of four WDS analyses are: Ir 18.68 (15.55–21.96), Rh 18.34 (16.32–20.32), Pt 0.64 (0.19–1.14), Ru 0.03 (0.00–0.13), Os 0.07 (0.02–0.17), Fe 14.14 (13.63–14.64), Ni 13.63 (12.58–14.66), Cu 4.97 (3.42–6.41), Co 0.09 (0.07–0.11), S 29.06 (28.48–29.44), and total 99.66 wt.%. They correspond to the following formula calculated for a total of 31 atoms per formula unit: (Rh 3.16 Ir 1.72 Pt 0.06 Ru 0.01 Os 0.01 ) Σ 4.95 (Fe 4.48 Ni 4.11 Cu 1.38 Co 0.03 ) Σ10.00 S 16.05 . The results of synchrotron micro-Laue diffraction studies indicate that ferrotorryweiserite is trigonal; its probable space group is $R\bar{3}m$ (#166) based on its Ni-analog, torryweiserite. The unit-cell parameters refined from 177 reflections are a = 7.069 (2) Å, c = 34.286 (11) Å, V = 1484 (1) Å 3 , and Z = 3. The c:a ratio is 4.8502. The strongest eight peaks in the X-ray diffraction pattern derived from results of micro-Laue diffraction study [d in Å(hkil)(I)] are 2.7950 ($20\bar{2}5$) (100); 5.7143 (0006) (60); 1.7671 ($22\bar{4}0$) (44.4); 3.0486 ($20\bar{2}1$) (39.4); 5.7650 ($10\bar{1}2$) (38.6); 2.5956 ($20\bar{2}7$) (37.8); 3.0058 ($11\bar{2}6$) (36.5); and 1.5029 ($4\bar{2}\bar{2}12$) (35.3). Ferrotorryweiserite and the associated PGM crystallized from microvolumes of residual melt at late stages of crystallization of grains of Os- and Ir-dominant alloys occurred in lode zones of chromitites of the Lysanskiy layered complex. In a particular case, the residual melt is disposed peripherally around a core containing a disequilibrium association of magnesian olivine (Fo 72.9–75.6 ) and albite (Ab 81.6–86.4 ), with the development of skeletal crystals of titaniferous augite: Wo 40.8–43.2 En 26.5–29.3 Fs 20.3–22.6 Aeg 6.9–9.5 (2.82–3.12 wt.% TiO 2 ). Ferrotorryweiserite represents the Fe-dominant analog of torryweiserite. We also report occurrences of ferrotorryweiserite in the Marathon deposit, Coldwell Complex, Ontario, Canada, and infer the existence of the torryweiserite–ferrotorryweiserite solid solution in other deposits and complexes.

58 GEOSCIENCES↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2023

Seeing the Disturbed Forest for the Trees: Remote Sensing Is Underutilized to Quantify Critical Zone Response to Unprecedented Disturbance

Understanding the severity and extent of near surface critical zone (CZ) disturbances and their ecosystem response is a pressing concern in the face of increasing human and natural disturbances. Predicting disturbance severity and recovery in a changing climate requires comprehensive understanding of ecosystem feedbacks among vegetation and the surrounding environment, including climate, hydrology, geomorphology, and biogeochemistry. Field surveys and satellite remote sensing have limited ability to effectively capture the spatial and temporal variability of disturbance and CZ properties. Technological advances in remote sensing using new sensors and new platforms have improved observations of changes in vegetation canopy structure and productivity; however, integrating measures of forest disturbance from various sensing platforms is complex. By connecting the potential for remote sensing technologies to observe different CZ disturbance vectors, we show that lower severity disturbance and slower vegetation recovery are more difficult to quantify. Case studies in montane forests from the western United States highlight new opportunities, including evaluating post-disturbance forest recovery at multiple scales, shedding light on understory vegetation regrowth, detecting specific physiological responses, and refining ecohydrological modeling. Learning from regional CZ disturbance case studies, we propose future directions to synthesize fragmented findings with (a) new data analysis using new or existing sensors, (b) data fusion across multiple sensors and platforms, (c) increasing the value of ground-based observations, (d) disturbance modeling, and (e) synthesis to improve understanding of disturbance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An operando synchrotron study on the effect of wire melting state on solidification microstructures of Inconel 718 in wire-laser directed energy deposition

Directed energy deposition (DED) with a coaxial wire-laser configuration has gained significant attention in recent years for the production of large-scale metallic components because of its low directional dependence, fast deposition rate, high feedstock efficiency, and low manufacturing costs. This work studies the coaxial wire-laser DED process of Inconel 718 alloy under a stable deposition condition with a relatively low input volumetric energy density (55.5J/mm 3 ). Post characterization reveals a cluster of refined grains at the center-bottom region of the as-printed track. Operando high-energy synchrotron X-ray experiments and multi-physics modeling are applied innovatively to study the fundamental mechanism responsible for the formation of this microstructure. The X-ray diffraction experiment provides direct evidence, which is supported by the simulation, that the feeding wire can reach the melt pool bottom and release solid particles (primarily carbides) near the mushy zone owing to insufficient melting. Consequently, these sub-micron sized particles suppress the growth of large columnar grains and cause the formation of unique microstructural heterogeneity. In conclusion, this discovery offers new opportunities for tailoring the solidification microstructure by controlling the melting state of the feedstock wire in DED process, in addition to commonly known factors such as the thermal gradient and solidification velocity.

42 ENGINEERING↗
NASA NTRS2003

Reduction of Defects in Germanium-Silicon

Crystals grown without contact with a container have far superior quality to otherwise similar crystals grown in direct contact with a container. In addition to float-zone processing, detached- Bridgman growth is a promising tool to improve crystal quality, without the limitations of float zoning or the defects introduced by normal Bridgman growth. Goals of this project include the development of the detached Bridgman process to be reproducible and well understood and to quantitatively compare the defect and impurity levels in crystals grown by these three methods. Germanium (Ge) and germanium-silicon (Ge-Si) alloys are being used. At MSFC, we are responsible for the detached Bridgman experiments intended to differentiate among proposed mechanisms of detachment, and to confirm or refine our understanding of detachment. Because the contact angle is critical to determining the conditions for detachment, the sessile drop method was used to measure the contact angles as a function of temperature and composition for a large number of substrates made of potential ampoule materials. Growth experiments have used pyrolytic boron nitride (pBN) and fused silica ampoules with the majority of the detached results occurring predictably in the pBN. Etch pit density (EPD) measurements of normal and detached Bridgman-grown Ge samples show a two order of magnitude improvement in the detached-grown samples. The nature and extent of detachment is determined by using profilometry in conjunction with optical and electron microscopy. The stability of detachment has been analyzed, and an empirical model for the conditions necessary to achieve sufficient stability to maintain detached growth for extended periods has been developed. We have investigated the effects on detachment of ampoule material, pressure difference above and below the melt, and Si concentration; samples that are nearly completely detached can be grown repeatedly in pBN. Current work is concentrated on developing a method to make in situ pressure measurements in the growth ampoules.

Source record↗
DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗
DOE OSTI2021

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗
NASA NTRS2010

The Fourth SeaWiFS HPLC Analysis Round-Robin Experiment (SeaHARRE-4)

Ten international laboratories specializing in the determination of marine pigment concentrations using high performance liquid chromatography (HPLC) were intercompared using in situ samples and a mixed pigment sample. Although prior Sea-viewing Wide Field-of-view Sensor (SeaWiFS) High Performance Liquid Chromatography (HPLC) Round-Robin Experiment (SeaHARRE) activities conducted in open-ocean waters covered a wide dynamic range in productivity, and some of the samples were collected in the coastal zone, none of the activities involved exclusively coastal samples. Consequently, SeaHARRE-4 was organized and executed as a strictly coastal activity and the field samples were collected from primarily eutrophic waters within the coastal zone of Denmark. The more restrictive perspective limited the dynamic range in chlorophyll concentration to approximately one and a half orders of magnitude (previous activities covered more than two orders of magnitude). The method intercomparisons were used for the following objectives: a) estimate the uncertainties in quantitating individual pigments and higher-order variables formed from sums and ratios; b) confirm if the chlorophyll a accuracy requirements for ocean color validation activities (approximately 25%, although 15% would allow for algorithm refinement) can be met in coastal waters; c) establish the reduction in uncertainties as a result of applying QA procedures; d) show the importance of establishing a properly defined referencing system in the computation of uncertainties; e) quantify the analytical benefits of performance metrics, and f) demonstrate the utility of a laboratory mix in understanding method performance. In addition, the remote sensing requirements for the in situ determination of total chlorophyll a were investigated to determine whether or not the average uncertainty for this measurement is being satisfied.

Hooker, Stanford B.↗
NASA NTRS2002

New Concepts for Refinement of Cumulus Parameterization in GCM's the Arakawa-Schubert Framework

Several state-of-the-art models including the one employed in this study use the Arakawa-Schubert framework for moist convection, and Sundqvist formulation of stratiform. clouds, for moist physics, in-cloud condensation, and precipitation. Despite a variety of cloud parameterization methodologies developed by several modelers including the authors, most of the parameterized cloud-models have similar deficiencies. These consist of: (a) not enough shallow clouds, (b) too many deep clouds; (c) several layers of clouds in a vertically demoralized model as opposed to only a few levels of observed clouds, and (d) higher than normal incidence of double ITCZ (Inter-tropical Convergence Zone). Even after several upgrades consisting of a sophisticated cloud-microphysics and sub-grid scale orographic precipitation into the Data Assimilation Office (DAO)'s atmospheric model (called GEOS-2 GCM) at two different resolutions, we found that the above deficiencies remained persistent. The two empirical solutions often used to counter the aforestated deficiencies consist of a) diffusion of moisture and heat within the lower troposphere to artificially force the shallow clouds; and b) arbitrarily invoke evaporation of in-cloud water for low-level clouds. Even though helpful, these implementations lack a strong physical rationale. Our research shows that two missing physical conditions can ameliorate the aforestated cloud-parameterization deficiencies. First, requiring an ascending cloud airmass to be saturated at its starting point will not only make the cloud instantly buoyant all through its ascent, but also provide the essential work function (buoyancy energy) that would promote more shallow clouds. Second, we argue that training clouds that are unstable to a finite vertical displacement, even if neutrally buoyant in their ambient environment, must continue to rise and entrain causing evaporation of in-cloud water. These concepts have not been invoked in any of the cloud parameterization schemes so far. We introduced them into the DAO-GEOS-2 GCM with McRAS (Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme).

Sud, Y. C.↗
NASA NTRS2006

Assessing Hurricane Katrina Vegetation Damage at Stennis Space Center using IKONOS Image Classification Techniques

Hurricane Katrina inflicted widespread damage to vegetation in southwestern coastal Mississippi upon landfall on August 29, 2005. Storm damage to surface vegetation types at the NASA John C. Stennis Space Center (SSC) was mapped and quantified using IKONOS data originally acquired on September 2, 2005, and later obtained via a Department of Defense ClearView contract. NASA SSC management required an assessment of the hurricane s impact to the 125,000-acre buffer zone used to mitigate rocket engine testing noise and vibration impacts and to manage forestry and fire risk. This study employed ERDAS IMAGINE software to apply traditional classification techniques to the IKONOS data. Spectral signatures were collected from multiple ISODATA classifications of subset areas across the entire region and then appended to a master file representative of major targeted cover type conditions. The master file was subsequently used with the IKONOS data and with a maximum likelihood algorithm to produce a supervised classification later refined using GIS-based editing. The final results enabled mapped, quantitative areal estimates of hurricane-induced damage according to general surface cover type. The IKONOS classification accuracy was assessed using higher resolution aerial imagery and field survey data. In-situ data and GIS analysis indicate that the results compare well to FEMA maps of flooding extent. The IKONOS classification also mapped open areas with woody storm debris. The detection of such storm damage categories is potentially useful for government officials responsible for hurricane disaster mitigation.

Spruce, Joseph P.↗
DOE OSTI2024

Simulation of Compound Flooding Using River‐Ocean Two‐Way Coupled E3SM Ensemble on Variable‐Resolution Meshes

Abstract Coastal zone compound flooding (CF) can be caused by the interactive fluvial and oceanic processes, particularly when coastal backwater propagates upstream and interacts with high river discharge. The modeling of CF is limited in existing Earth System Models (ESMs) due to coarse mesh resolutions and one‐way coupled river‐ocean components. In this study, we present a novel multi‐scale coupling framework within the Energy Exascale Earth System Model (E3SM), integrating global atmosphere and land with interactively coupled river and ocean models using different meshes with refined resolutions near the coastline. To evaluate this framework, we conducted ensemble simulations of a CF event (Hurricane Irene in 2011) in a Mid‐Atlantic estuary. The results demonstrate that the novel E3SM configuration can reasonably reproduce river discharge and sea surface height variations. The two‐way river‐ocean coupling improves the representation of coastal backwater effects at the terrestrial‐aquatic interface that are caused by the combined actions of tide and storm surge during the CF event, thus providing a valuable modeling tool for better understanding the river‐estuary‐ocean dynamics in extreme events under climate change. Notably, our results show that the most significant CF impacts occur when the highest storm surge generated by a tropical cyclone meets with a moderate river discharge. This study highlights the state‐of‐the‐art advancements developed within E3SM for simulating multi‐scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

[MODIS Investigation]

The objectives of the last six months were: (1) Revise the algorithms for the Fluorescence Line Height (FLH) and Chlorophyll Fluorescence Efficiency (CFE) products, especially the data quality flags; (2) Revise the MOCEAN validation plan; (3) Deploy and recover bio-optical instrumentation at the Hawaii Ocean Time-series (HOT) site as part of the Joint Global Ocean Flux Study (JGOFS); (4) Prepare for field work in the Antarctic Polar Frontal Zone as part of JGOFS; (5) Submit manuscript on bio-optical time scales as estimated from Lagrangian drifters; (6) Conduct chemostat experiments on fluorescence; (7) Interface with the Global Imager (GLI) science team; and (8) Continue development of advanced data system browser. We are responsible for the delivery of two at-launch products for AM-1: Fluorescence line height (FLH) and chlorophyll fluorescence efficiency (CFE). We also considered revising the input chlorophyll, which is used to determine the degree of binning. We have refined the quality flags for the Version 2 algorithms. We have acquired and installed a Silicon Graphics Origin 200. We are working with the University of Miami team to develop documentation that will describe how the MODIS ocean components are linked together.

Abbott, Mark R.↗