Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kappa”

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 397 records · Page 22

Variation in Global Chemical Composition of PM2.5: Emerging Results from SPARTAN

The Surface PARTiculate mAtter Network (SPARTAN) is a long-term project that includes characterization of chemical and physical attributes of aerosols from filter samples collected worldwide. This paper discusses the ongoing efforts of SPARTAN to define and quantify major ions and trace metals found in fine particulate matter (PM (sub 2.5). Our methods infer the spatial and temporal variability of PM (sub 2.5) in a cost-effective manner. Gravimetrically weighed filters represent multi-day averages of PM (sub 2.5), with a collocated nephelometer sampling air continuously. SPARTAN instruments are paired with AErosol RObotic NETwork (AERONET) sun photometers to better understand the relationship between ground-level PM (sub 2.5) and columnar aerosol optical depth (AOD). We have examined the chemical composition of PM (sub 2.5) at 12 globally dispersed, densely populated urban locations and a site at Mammoth Cave (US) National Park used as a background comparison. So far, each SPARTAN location has been active between the years 2013 and 2016 over periods of 2-26 months, with an average period of 12 months per site. These sites have collectively gathered over 10 years of quality aerosol data. The major PM (sub 2.5) constituents across all sites (relative contribution plus or minus Standard Deviation) are ammoniated sulfate (20 percent plus or minus 11 percent), crustal material (13.4 percent plus or minus 9.9 percent), equivalent black carbon (11.9 percent plus or minus 8.4 percent), ammonium nitrate (4.7 percent plus or minus 3.0 percent), sea salt (2.3 percent plus or minus 1.6 percent), trace element oxides (1.0 percent plus or minus 1.1 percent), water (7.2 percent plus or minus 3.3 percent) at 35 percent relative humidity, and residual matter (40 percent plus or minus 24 percent). Analysis of filter samples reveals that several PM (sub 2.5) chemical components varied by more than an order of magnitude between sites. Ammoniated sulfate ranges from 1.1 microns per cubic meter (Buenos Aires, Argentina) to 17 microns per cubic meter (Kanpur, India in the dry season). Ammonium nitrate ranged from 0.2 microns per cubic meter (Mammoth Cave, in summer) to 6.8 microns per cubic meter (Kanpur, dry season). Equivalent black carbon ranged from 0.7 microns per cubic meter (Mammoth Cave) to over 8 microns per cubic meter (Dhaka, Bangladesh and Kanpur, India). Comparison of SPARTAN vs. coincident measurements from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network at Mammoth Cave yielded a high degree of consistency for daily PM (sub 2.5) (r squared equals 0.76, slope equals 1.12), daily sulfate (r squared equals 0.86, slope equals 1.03), and mean fractions of all major PM (sub 2.5) components (within 6 percent). Major ions generally agree well with previous studies at the same urban locations (e.g. sulfate fractions agree within 4 percent for 8 out of 11 collocation comparisons). Enhanced anthropogenic dust fractions in large urban areas (e.g. Singapore, Kanpur, Hanoi, and Dhaka) are apparent from high Zn to Al ratios. The expected water contribution to aerosols is calculated via the hygroscopicity parameter kappa (sub v (volume)) for each filter. Mean aggregate values ranged from 0.15 (Ilorin) to 0.28 (Rehovot). The all-site parameter mean is 0.20 plus or minus 0.04. Chemical composition and water retention in each filter measurement allows inference of hourly PM (sub 2.5) at 35 percent relative humidity by merging with nephelometer measurements. These hourly PM (sub 2.5) estimates compare favourably with a beta attenuation monitor (MetOne) at the nearby US embassy in Beijing, with a coefficient of variation r squared equals 0.67 (number equals 3167), compared to r squared equals 0.62 when v (volume) was not considered. SPARTAN continues to provide an open-access database of PM (sub 2.5) compositional filter information and hourly mass collected from a global federation of instruments.

Snider, Graydon↗

HOTEX: An Approach for Global Mapping of Human Built-Up and Settlement Extent

Understanding the impacts of urbanization requires accurate and updatable urban extent maps. Here we present an algorithm for mapping urban extent at global scale using Landsat data. An innovative hierarchical object-based texture (HOTex) classification approach was designed to overcome spectral confusion between urban and nonurban land cover types. VIIRS nightlights data and MODIS vegetation index datasets are integrated as high-level features under an object-based framework. We applied the HOTex method to the GLS-2010 Landsat images to produce a global map of human built-up and settlement extent. As shown by visual assessments, our method could effectively map urban extent and generate consistent results using images with inconsistent acquisition time and vegetation phenology. Using scene-level cross validation for results in Europe, we assessed the performance of HOTex and achieved a kappa coefficient of 0.91, compared to 0.74 from a baseline method using per-pixel classification using spectral information.

Wang, Panshi↗

Formation and Past Evolution of the Showers of 96P/Machholz Complex

In this work we model the dynamical evolution of meteoroid streams of comet 96P/Machholz, and the largest member of the Marsden sunskirters, comet P/1999 J6. We simultaneously fit the characteristics of eight meteor showers which have been proposed to be linked to the complex, using observations from a range of techniques - visual, video, TV and radar. The aim is to obtain a self-consistent scenario of past capture of a large comet into a short-period orbit, and its subsequent fragmentation history. Moreover, we also aim to constrain the dominant parent of these showers. The fit of our simulated shower characteristics to observations is consistent with the scenario of a capture of a proto-comet 96P/Machholz by Jupiter circa 20000 BCE (Before Christian Era), and a subsequent major breakup around 100-950 CE (Christian Era) which resulted in the formation of the Marsden group of comets. We find that the Marsden group of comets are not the immediate parents of the daytime Arietids and Northern and Southern delta-Aquariids, as previously suggested. In fact, the hypothesis that the Northern delta-Aquariids are related to the Marsden group of comets is not supported by this study. The bulk of the observational characteristics of all eight showers can be explained by meteoroid ejection primarily from comet 96P/Machholz between 10000 BCE and 20000 BCE. Assuming the Marsden group of comets originated between 100 CE-950 CE, we conclude that sunskirting comets contribute mainly to the meteoroid stream near the time of the peak of the daytime Arietids, Southern delta-Aquariids, kappa-Velids. Finally, we find that the meteor showers identified by Babadzhanov and Obrubov (1992) as the alpha-Cetids, the Ursids and Carinids correspond to the daytime lambda-Taurids, the November iota-Draconids or December alpha-Draconids and the theta-Carinids.

Abedin, Abedin↗

Mapping Mangrove Extent and Change: A Globally Applicable Approach

This study demonstrates a globally applicable method for monitoring mangrove forest extent at high spatial resolution. A 2010 mangrove baseline was classified for 16 study areas using a combination of ALOS PALSAR and Landsat composite imagery within a random forests classifier. A novel map-to-image change method was used to detect annual and decadal changes in extent using ALOS PALSAR/JERS-1 imagery. The map-to-image method presented makes fewer assumptions of the data than existing methods, is less sensitive to variation between scenes due to environmental factors (e.g., tide or soil moisture) and is able to automatically identify a change threshold. Change maps were derived from the 2010 baseline to 1996 using JERS-1 SAR and to 2007, 2008 and 2009 using ALOS PALSAR. This study demonstrated results for 16 known hotspots of mangrove change distributed globally, with a total mangrove area of 2,529,760 ha. The method was demonstrated to have accuracies consistently in excess of 90% (overall accuracy: 92.2–93.3%, kappa: 0.86) for mapping baseline extent. The accuracies of the change maps were more variable and were dependent upon the time period between images and number of change features. Total change from 1996 to 2010 was 204,850 ha (127,990 ha gain, 76,860 ha loss), with the highest gains observed in French Guiana (15,570 ha) and the highest losses observed in East Kalimantan, Indonesia (23,003 ha). Changes in mangrove extent were the consequence of both natural and anthropogenic drivers, yielding net increases or decreases in extent dependent upon the study site. These updated maps are of importance to the mangrove research community, particularly as the continual updating of the baseline with currently available and anticipated spaceborne sensors. It is recommended that mangrove baselines are updated on at least a 5-year interval to suit the requirements of policy makers.

Thomas, Nathan↗

Mark Twain National Forest Ecological Forecasting - Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest

The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.

Kaitlyn Bretz↗

Mark Twain National Forest Ecological Forecasting: Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest

The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.

DEVELOP Tech Paper↗

Remote monitoring of agricultural systems using NDVI time series and machine learning methods: a tool for an adaptive agricultural policy

This study aims to provide accurate information about changes in agricultural systems (AS) using phenological metrics derived from the NDVI time series. Use of such information could help land managers optimize land use choices and monitor the status of agricultural lands, under a variety of environmental and socioeconomic conditions. For this purpose, the Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI data were used to derive phenological metrics over the Oum Er-Rbia basin (central Morocco). Random forest (RF), support vector machine (SVM), and K-nearest neighbor (KNN) classifiers were explored and compared on their ability to classify AS classes over the study area. Four main AS classes have been considered: (1) irrigated annual crop (IAC), (2) irrigated perennial crop (IPC), (3) rainfed area (RA), and (4) fallow (FA). By comparing the accuracy of the three classifiers, the RF method showed the best performance with an overall accuracy of 0.97 and kappa coefficient of 0.96.The RF method was then chosen to examine time variations in AS over a 16-year period (2000–2016). The AS main variations were detected and evaluated for the four AS classes. These variations have been found to be linked well with other indicators of local agricultural land management, as well as the historical agricultural drought changes over the study area. Overall, the results present a tool for decision makers to improve agricultural management and provide a different perspective in understanding the spatiotemporal dynamics of agricultural systems.

Youssef Lebrini↗

Improved assessment of mangrove forests in Sundarbans East Wildlife Sanctuary using WorldView 2 and TanDEM-X high resolution imagery

Recent developments of remote sensing techniques which can capture both the structure and function of the ecosystem provide a more representative view of the landscape. These unique Earth observations were used to help improve traditional forestry surveys by providing species-specific land cover classes for mangrove forests in the Sundarbans East Wildlife Sanctuary. By combining optical data from WorldView2 (WV2; 2 m pixel) and a canopy height model derived using radar data from TanDEM-X (TDX; 12 m pixel), we identified nine mangrove and five non-mangrove classes by following an Iterative Self-Organizing Data Analysis Algorithm. Three dominant mangrove species accounted for nearly 50% of the sanctuary. Heritieria fomes disproportionately covered the largest area at 43%, overturning previous field-based estimates of Excoecaria agallocha dominance. E. agallocha and Sonneratia apetala, covered 3% and 1.47% of the sanctuary, respectively. Four mixed species classes were also identified with clear vegetation zonation patterns that trended toward species homogeneity with increasing distance from shore. The overall land cover accuracy (WV2: 89.33%; WV2-TDX: 89.89%), the Kappa Coefficient (WV2:0.88; WV2-TDX: 0.89) and change statistics between WV2 and WV2-TDX landcover classifications indicate that the WV2 imagery can separate mangrove community types without structural data. The combination of the land cover classifications and the canopy height model indicated that H. fomes were not only the most dominant forest but also, on average, the tallest (12.3 m) among the other eight mangrove types. Our large-scale mapping with high resolution optical and radar platforms can capture subtle changes in mangrove vegetation and canopy structural gradients more accurately and be used to monitor biodiversity changes and Aichi Biodiversity Targets and Indicators, which would contribute to biodiversity policy updating.

Md Mizanur Rahman↗

Assessing Cyanobacterial Frequency and Abundance at Surface Waters Near Drinking Water Intakes Across the United States

This study presents the first large-scale assessment of cyanobacterial frequency and abundance of surface water near drinking water intakes across the United States. Public water systems serve drinking water to nearly 90% of the United States population. Cyanobacteria and their toxins may degrade the quality of finished drinking water and can lead to negative health consequences. Satellite imagery can serve as a cost-effective and consistent monitoring technique for surface cyanobacterial blooms in source waters and can provide drinking water treatment operators information for managing their systems. This study uses satellite imagery from the European Space Agency’s Ocean and Land Colour Instrument (OLCI) spanning June 2016 through April 2020. At 300-m spatial resolution, OLCI imagery can be used to monitor cyanobacteria in 685 drinking water sources across 285 lakes in 44 states, referred to here as resolvable drinking water sources. First, a subset of satellite data was compared to a subset of responses (n = 84) submitted as part of the U.S. Environmental Protection Agency’s fourth Unregulated Contaminant Monitoring Rule (UCMR 4). These UCMR 4 qualitative responses included visual observations of algal bloom presence and absence near drinking water intakes from March 2018 through November 2019. Overall agreement between satellite imagery and UCMR 4 qualitative responses was 94% with a Kappa coefficient of 0.70. Next, temporal frequency of cyanobacterial blooms at all resolvable drinking water sources was assessed. In 2019, bloom frequency averaged 2% and peaked at 100%, where 100% indicated a bloom was always present at the source waters when satellite imagery was available. Monthly cyanobacterial abundances were used to assess short-term trends across all resolvable drinking water sources and effect size was computed to provide insight on the number of years of data that must be obtained to increase confidence in an observed change. Generally, 2016-2020 was an insufficient time period for confidently observing changes at these source waters; on average, a decade of satellite imagery would be required for observed environmental trends to outweigh variability in the data. However, five source waters did demonstrate a sustained short-term trend, with one increasing in cyanobacterial abundance from June 2016 to April 2020 and four decreasing.

remote sensing↗

Remote Sensing of Lineage Functional Types for Modeling and Monitoring Biodiversity

Hyperspectral remote sensing has the potential to continuously scale plant function and plant diversity information from landscape to global extents. Numerous studies have indicated that VSWIR (400-2500 nm) reflectance properties of vegetation capture evolutionarily conserved biochemical, structural, and other functional attributes of plant species. Spectral properties conserved in plants provide the opportunity to both 1) aggregate species into lineages with improved classification accuracy and 2) link those lineages directly to plant traits. Full realization of this goal will enable parameterization of Land Surface Models (LSMs) with remotely sensed information, e.g., canopy nitrogen, and better representations of biodiversity and functional diversity in biogeographic studies. In this study, we use hyperspectral AVIRIS data from the 2013 HyspIRI campaign over the Southern Sierra Nevada, California flight box to investigate the potential for incorporating evolutionary thinking into landcover classification. We link the airborne hyperspectral data with vegetation plot data from roughly 1372 surveys and a phylogeny representing 1361 species. We aggregate species into lineages ranging from species level groups down to similar number of Plant Functional Types as often used in LSMs. We assessed the ability of Random Forest and Partial Least Squares Discriminant Analysis to discriminate across these different phylogenetic scales and determine the optimal number of lineages to classify. Although there are some temporal and spatial differences in our training data, our best approaches achieved moderate classification accuracy (Kappa > 0.65). Given an optimal number of lineages, we explored approaches to improve classifications including machine learning and unmixing approaches. This work suggests that lineage-based methods may be a promising way to leverage the huge amounts of data that will come from high resolution and high return interval hyperspectral data planned for the Surface Biology and Geology mission with sparsely sampled existing ground-based ecological data.

Hyperspectral↗

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Solar Flare Heating with Turbulent Suppression of Thermal Conduction

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppmission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function.

Joel C. Allred↗

Solar Flare Heating with Turbulent Suppression of Thermal Conduction

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function

Solar flares↗

Bhutan Agriculture III: Monitoring Cropland Changes in Bhutan using Remote Sensing to Bolster Food Security and Support Crop Monitoring

The Bhutan Agriculture III team aimed to improve agricultural efficiency in Bhutan. Bhutan is a nation heavily reliant on agriculture, but it faces challenges such as geophysical limitations and lack of scientific agricultural practice. The team partnered with a primary end user, Bhutan’s Department of Agriculture (DoA), and with collaborators; the Bhutan Foundation, National Plant Protection Centre (NPPC), Agricultural Research Department Centre (ARDC), National Statistics Bureau (NSB), and the Ugyen Wangchuck Institute for Conservation and Environment Research (UWICER). Advised by NASA SERVIR, the team developed crop masks and monitored rice distribution from 2015 to 2022 utilizing Earth observations such as Landsat 8 Operational Land Imager (OLI), Landsat 9 OLI-2, Sentinel-1 C-Band Synthetic Aperture Radar (C-SAR), Sentinel-2 MultiSpectral Instrument (MSI) and Shuttle Radar Topography Mission (SRTM). The team gathered 5,000 points from the five dzongkhags that yield the most rice in Bhutan (Paro, Punakha, Samtse, Sarpang and Wangue Phodrang) using Collect Earth Online (CEO). With the data collected, the team split the data into training and validation data on Google Earth Engine (GEE) for a random forest (RF) classifier for rice and non-rice classification. After running the data on the Random Forest (RF) model, the team got an accuracy score of 81.48%, a kappa score of 55.75% and an F1 score of 86.11%. This data supports better agricultural decision-making for the governing body of Bhutan, helps enhance farming efficiency and foster sustainable practices, assists in overcoming data inaccuracy and bolsters food security in the country.

Sonam Seldon Tshering↗

Jefferson County Ecological Conservation: Quantifying the Effects of Hydrologic Restoration in the Camas National Wildlife Refuge and Mud Lake Wildlife Management Area

Wetlands in the Camas National Wildlife Refuge and Mud Lake Wildlife Management Area, a rare landscape feature in the Intermountain West, are crucial for migratory birds along the Pacific Flyway, yet these wetlands have experienced a noticeable decline in extent and inundation over the last 40+ years. Partners at the U.S. Fish and Wildlife Service and Idaho Department of Fish and Game have ongoing restoration efforts that are yet to be quantified. To assist these partners, we used NASA Earth observations to quantify change in wetland extent and determine if the restoration projects met the intended impact. The QA_PIXEL band from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) was used to identify water. Landsat 8 OLI imagery was also used for classifying landcover in ArcGIS Pro for 2016 and 2020 and to forecast landcover in 2060 using TerrSet’s Land Change Modeler. We also compared lidar data from 2011 and 2019 to detect changes in vegetation height, ground elevation, and surface water levels, and Sentinel-2 Multispectral Instrument (MSI) imagery to detect changes in vegetation health. Classification results indicated overall accuracies over 84% and Kappa indices above 0.80, well surpassing random classification performance. From 2016 to 2020 our classifications showed an 825 acre increase in wetland extent, and our forecasted 2060 model showed a 2795 acre decrease. Combining our partners’ knowledge of the area with our analyses, we created a remote sensing workflow to enhance future monitoring and decision-making. This study demonstrated the feasibility of NASA Earth observations in quantifying wetlands, informing restoration projects, and supporting avian conservation.

Cassidy Bromka↗

Equivalent Fundamental-Mode Source Simulations for Spherical Uranium and Plutonium Systems

It is important in neutron noise measurements and reactor startup to understand how the source distribution in a system compares to the equivalent fundamental-mode (EFM) source. Deviation from the fundamental mode source distribution increases as the system $\kappa$ eff increases. A correction factor, g *, is a measure of this deviation, and is the focus of this work. Two methods were utilized to simulate g *. These methods are compared to each other and (where applicable) previous results. The effects associated with nuclear material inner radius, materials of reflectors, and reflector thickness are also investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗