Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

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.

1,076 records · Page 60

The Continuous Intercomparison of Radiation Codes (CIRC): Phase I Cases

CIRC aspires to be the successor to ICRCCM (Intercomparison of Radiation Codes in Climate Models). It is envisioned as an evolving and regularly updated reference source for GCM-type radiative transfer (RT) code evaluation with the principle goal to contribute in the improvement of RT parameterizations. CIRC is jointly endorsed by DOE's Atmospheric Radiation Measurement (ARM) program and the GEWEX Radiation Panel (GRP). CIRC's goal is to provide test cases for which GCM RT algorithms should be performing at their best, i.e, well characterized clear-sky and homogeneous, overcast cloudy cases. What distinguishes CIRC from previous intercomparisons is that its pool of cases is based on observed datasets. The bulk of atmospheric and surface input as well as radiative fluxes come from ARM observations as documented in the Broadband Heating Rate Profile (BBHRP) product. BBHRP also provides reference calculations from AER's RRTM RT algorithms that can be used to select the most optimal set of cases and to provide a first-order estimate of our ability to achieve radiative flux closure given the limitations in our knowledge of the atmospheric state.

Oreopoulos, Lazaros↗

Dynamics of three-dimensional plasma clouds with coupling to the background ionosphere

A three-dimensional, time-dependent model with a two-grid system was developed to study the expansion of a plasma cloud in the F region and topside ionosphere. The model maintains an adequate resolution for the released cloud motion and its interaction with the immediate environment, and it includes the effect due to the coupling with the distant part of the ionosphere (i.e., E region). Simulations were performed using realistic background ionospheric density profiles in both the E and F regions. The results show that the cloud coupling to the underlying E region affects the perpendicular cloud motion the most. The distant coupling acts to reduce the perturbation potential and perpendicular velocity and delays or eliminates the striations. These simulation results are consistent with simple analytical approximations. The simulation results also show that the distant coupling has a very small effect on 'localized' phenomena, such as the cloud expansion along the B RIGHT ARROW field and the electrostatic snowplow. The cloud-induced electric potential is attenuated in the lower E region. The electrons flow along the B RIGHT ARROW field, carrying the current to the E region and back to the cloud. The current closure is demonstrated in three dimensions for the first time for such a problem. The perpendicular current flowing through the plasma cloud is closed by the field-aligned electron current and the background perpendicular (mainly Pedersen) current in both the E and F regions. The 'image cloud' formation in t he E region is also clearly demonstrated. The variation of the density change in the 'image cloud' along the B RIGHT ARROW field and the features of the image cloud are shown.

Ma, T.-Z.↗

Aerodynamic Performance and Flow-Field Characteristics of Two Waverider-Derived Hypersonic Cruise Configurations

The component integration of a class of hypersonic high-lift configurations known as waveriders into hypersonic cruise vehicles was evaluated. A wind-tunnel model was developed which integrates realistic vehicle components with two waverider shapes, referred to as the straight-wing and cranked-wing shapes. Both shapes were conical-flow-derived waveriders for a design Mach number of 4.0. Experimental data and limited computational fluid dynamics (CFD) predictions were obtained over a Mach number range of 1.6 to 4.63 at a Reynolds number of 2.0 x 10(exp 6) per foot. The CFD predictions and flow visualization data confirmed the shock attachment characteristics of the baseline waverider shapes and illustrated the waverider flow-field properties. Experimental data showed that no significant performance degradations, in terms of maximum lift-to-drag ratios, occur at off-design Mach numbers for the waverider shapes and the integrated configurations. A comparison of the fully-integrated waverider vehicles to the baseline shapes showed that the performance was significantly degraded when all of the components were added to the waveriders, with the most significant degradation resulting from aftbody closure and the addition of control surfaces. Both fully-integrated configurations were longitudinally unstable over the Mach number range studied with the selected center of gravity location and for unpowered conditions. The cranked-wing configuration provided better lateral-directional stability characteristics than the straight-wing configuration.

Cockrell, Charles E., Jr.↗

A Comparison of Cumulus Parameterizations in Idealized Sea-Breeze Simulations

Four cumulus parameterizations in the Penn State-NCAR model MM5 are compared in idealized sea-breeze simulations, with the aim of discovering why they work as they do. The most realistic results appear to be those using the Kain-Fritsch scheme. Rainfall is significantly delayed with the Betts-Miller-Janjic scheme, due to the method of computing the reference sounding. This method can be corrected, but downdrafts should be added in a physically realistic manner. Even without downdrafts, a corrected version of the BMJ scheme produces nearly the same timing and location of deep convection as the KF scheme, despite the very different physics. In order to simulate the correct timing of the rainfall, a minimum amount of mass is required in the layer that is the source of a parameterized updraft. The Grell parameterization, in the present simulation, always derives the updraft from the top of the mixed layer, where vertical advection predominates over horizontal advection in increasing the moist static energy. This makes the application of the quasi-equilibrium closure more correct than it would be if the updrafts were always derived from the most unstable layer, but it evades the question of whether or not horizontal advection generates instability. Using different physics, the parameterizations produce significantly different cloud-top heights.

Cohen, Charles↗

Safety Culture at the World’s Premier Multi-User Spaceport

NASA’s Agency-wide Safety Culture is implemented at the Kennedy Space Center (KSC) using a unique strategy due an unparalleled approach in making human spaceflight history. Kennedy Space Center, the world’s premier multi-user spaceport, enables U.S. government and commercial space access, while providing the world a resource to allow the exploration of and the ability to work in space. A consistently healthy safety culture at KSC is imperative for mission success: desired achievements, protection of space flight hardware, and ultimately, the preservation of human life requires the support of a healthy safety culture. Emphasis on the NASA Agency-wide development of Safety Culture began with the conception of the NASA Agency Safety Culture Working Group. After the devastating loss of life and mission of the Space Shuttle Colombia, a broken safety culture was identified as an organizational cause by the Columbia Accident Investigation Board Report. Thus, the Agency Safety Culture Working Group was developed in 2009 to assess the status of the Agency’s Safety Culture, while addressing safety culture concerns at the NASA Center-level. A Five-factor model was developed to serve as the guiding principles for Safety Culture: 1) Reporting Culture, 2) Just Culture, 3) Flexible Culture, 4) Learning Culture, and 5) Engaged Culture. These five factors are included in the NASA Safety Culture logo, which was intentionally designed as a DNA double helix to prompt the permeation of safety into day-to-day work. KSC specifically implements the NASA Agency-wide Safety Culture principles in a tailored approach that is relevant to the diversity of work being performed. An emphasis is placed on the implementation to include safety at home, not exclusively at work. This emphasis is a KSC-specific element that has been intentionally added to promote a closed loop Safety Culture. To advertise the safety culture, various safety and health events are held throughout the calendar year, providing innovative speakers and engaging activities, while also promoting a wide range of curated safety initiatives. Development and continuous improvement of the KSC safety tracking database, allows for advanced tracking-to-closure, along with providing data sets used to identify areas of emphasis. Other safety initiatives rely solely on employee participation, such as the photo challenges; participants are encouraged to identify and capture themselves, coworkers, or family members participating in safe or healthy activities to share with others within the Center and at Agency levels. Fabrication of exclusive videos and graphics are utilized to advertise and inform employee of safety initiatives, upcoming safety events, and general dispersion of safety information. In addition, an anonymous Agency-wide Safety Culture Survey is advertised, administered, and analyzed at Kennedy Space Center, with the purpose of receiving basic feedback on Safety Culture perceptions to help prevent future incidents from occurring. Through these briefly identified means, and many other forms of employee engagement, Kennedy Space Center aims to maintain safety in the forefront, while creating an environment where everyone trusts that safety is a priority.

Larrin E. Moody↗

Assessing the Variability of Radiation, Water and Energy in the Deep Tropics Over the Last 3 Decades

Global surface energy closure and its variability depends heavily upon the surface radiative energy budget in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance is observed between energy and water fluxes in terms of energy closure in some atmospheric reanalysis systems, implying that processes are not simulated well. A recent paper by Hsaio et al (2023, pre-published version) found that the transition process has wind shear and longwave cloud radiative feedback signatures. In this presentation, we assess the surface radiative budget in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES Ed4.1, ISSCP FH, etc.) and describing the variability across the deep tropics for the period from 1988 to near present. We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear. We include further analysis contrasting 6 key tropical oceanic regions (Indian Ocean east and west, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean) and the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. We conclude with recommendations for continuing work.

Radiation budget↗

Computational Studies of Magnetic Nozzle Performance

An extensive literature review of magnetic nozzle research has been performed, examining previous work, as well as a review of fundamental principles. This has allow us to catalog all basic physical mechanisms which we believe underlie the thrust generation process. Energy conversion mechanisms include the approximate conservation of the magnetic moment adiabatic invariant, generalized hall and thermoelectric acceleration, swirl acceleration, thermal energy transformation into directed kinetic energy, and Joule heating. Momentum transfer results from the interaction of the applied magnetic field with currents induced in the plasma plume., while plasma detachment mechanisms include resistive diffusion, recombination and charge exchange collisions, magnetic reconnection, loss of adiabaticity, inertial forces, current closure, and self-field detachment. We have performed a preliminary study of Hall effects on magnetic nozzle jets with weak guiding magnetic fields and weak expansions (p(sub jet) approx. = P(sub background)). The conclusion from this study is that the Hall effect creates an azimuthal rotation of the plasma jet and, more generally, creates helical structures in the induced current, velocity field, and magnetic fields. We have studied plasma jet expansion to near vacuum without a guiding magnetic field, and are presently including a guiding magnetic field using a resistive MHD solver. This research is progressing toward the implementation of a full generalized Ohm's law solver. In our paper, we will summarize the basic principle, as well as the literature survey and briefly review our previous results. Our most recent results at the time of submittal will also be included. Efforts are currently underway to construct an experiment at the University of Michigan Plasmadynamics and Electric Propulsion Laboratory (PEPL) to study magnetic nozzle physics for a RF-thruster. Our computational study will work directly with this experiment to validate the numerical model, in order to study magnetic nozzle physics and optimize magnetic nozzle design. Preliminary results from the PEPL experiment will also be presented.

Ebersohn, Frans H.↗

Next-Generation MKIII Lightweight HUT/Hatch Assembly

The MK III (H-1) carbon-graphite/ epoxy Hard Upper Torso (HUT)/Hatch assembly was designed, fabricated, and tested in the early 1990s. The spacesuit represented an 8.3 psi (≈58 kPa) technology demonstrator model of a zero prebreathe suit. The basic torso shell, brief, and hip areas of the suit were composed of a carbon-graphite/epoxy composite lay-up. In its current configuration, the suit weighs approximately 120 lb (≈54 kg). However, since future planetary suits will be designed to operate at 0.26 bar (≈26 kPa), it was felt that the suit's re-designed weight could be reduced to 79 lb (≈35 kg) with the incorporation of lightweight structural materials. Many robust, lightweight structures based on the technologies of advanced honeycomb materials, revolutionary new composite laminates, metal matrix composites, and recent breakthroughs in fullerene fillers and nanotechnology lend themselves well to applications requiring materials that are both light and strong. The major problem involves the reduction in weight of the HUT/ Hatch assembly for use in lunar and/or planetary applications, while at the same time maintaining a robust structural design. The technical objective is to research, design, and develop manufacturing methods that support fa b rica - tion of a lightweight HUT/Hatch assembly using advanced material and geometric redesign as necessary. Additionally, the lightweight HUT/Hatch assembly will interface directly with current MK III hardware. Using the new operating pressure and current MK III (H-1) interfaces as a starting block, it is planned to maximize HUT/Hatch assembly weight reduction through material selection and geometric redesign. A hard upper torso shell structure with rear-entry closure and corresponding hatch will be fabricated. The lightweight HUT/Hatch assembly will retrofit and interface with existing MK III (H-1) hardware elements, providing NASA with immediate "plug-andplay" capability. NASA crewmembers will have a lightweight, robust, life-support system that will minimize fatigue during extraterrestrial surface sojourns. Its unique feature is the utilization of a new and innovative family of materials used by the aerospace industry, which at the time of this reporting has not been used for the proposed application.

McCarthy, Mike↗

A Combined Al-Mg/Pb-Pb Age of the Solar System

Astrophysical models of planet formation and protoplanetary disk evolution demand precise and accurate timing of the sequence of events in the solar nebula, relative to a time t=0, usually taken to be during the short epoch of CAI (Ca-rich, Al-rich inclusion) formation. Most CAIs formed withlive26Al (mean-life τ26= 1.034 Myr [1]), with an abundance 26Al/27Al ≈ (26Al/27Al)SS= 5.23 × 10-5[2]. We adopt this as the widespread level of 26Al in the solar nebula at t=0. Assuming spatial homogeneity of 26Al, an inclusion that had less 26Al, (26Al/27Al)0, formed a time Δt26= τ26ln[(26Al/27Al)SS/ (26Al/27Al)0] after t=0.These ages are typical precise to within ±0.1 Myr. Igneous bulk meteorites and inclusions can be relatively dated by the Al-Mg chronometer, but only ifΔt26<6 Myr. The Pb-Pb system is useful as a longer relative chronometer. It yields absolute ages tPb using 207Pb/206Pb, 206Pb/204Pb, and 238U/235U ratios measured indifferent portions of a sample, assuming certain half-lives [4]. These absolute ages are uncertain to within ±9 Myr due to uncertainties in the 235U half-life[3], but times of formation ΔtPb= tCAI–tPb relative to t=0, are more precise(±0.5Myr),iftCAI can be found. Here, tCAI means the Pb-Pb age that would be measured in CAIs using the half-lives the community typically uses, if they achieved isotopic closure at t=0. Unfortunately, direct Pb-Pb dating of CAIs has not definitively determined tCAI. Based on four CAIs with canonical (26Al/27Al)0,[5,6] found tPb= 4567.30 ± 0.16 Myr. No other CAI ages with measured 238U/235U have been reported in the refereed literature, but there are hints of other CAIs with ages tPb= 4568.0 ± 0.2 Myr [7] and tPb= 4568.3 ± 0.2 Myr [8].It is unclear whether anyof these igneous type B CAIs isotopically closed at t=0 or represents tCAI. Instead of measurements, we advocate finding tCAI by minimizing the discrepancies between the Al-Mg and Pb-Pb chronometers. Assuming Δt26=ΔtPb, we find the implied t’CAI= tPb+Δt26, then define t*CAIas the weighted mean of the t’CAI. t*CAIis the best guess for the Pb-Pb age of t=0; the assumption of homogeneity is justified if the t’CAI cluster within errors around t*CAI. This statistical approach is similar to, but improves on, that of[9]. We find t*CAI= 4568.73 ± 0.16 Myr. Below we discuss our methodology and the implications of this age for CAIs, 1.4 Myr older than the reported and typically used age 4567.30±0.16 Myr. Methods: We base our estimate of t*CAIon five achondrites for which published (26Al/27Al)0and Pb-Pb ages exist: the quenched angrites D’Orbigny, Sahara 99555 (SAH 99555), and Northwest Africa (NWA) 1670; the pseudo-eucrite Asuka 881394; and the inner disk achondrite. All are “NC” (non-carbonaceous) achondrites that likely cooled quickly enough that the Al-Mg and Pb-Pb systems achieved isotopic closure simultaneously. We also considered the “CC” (carbonaceous chondrite-like) achondrites NWA 2796 and NWA 6704, butdo not include them in our fit. Al-Mg and Pb-Pb seem not to have closed simultaneously, possibly because formation in the outer disk from volatile-rich composition led to slower cooling. Of the 8 chondrules from NWA 5697 measured by [20,21], we also consider the 4 for which 238U/235U was measured: 2-C1, 5-C2, 3-C5, 11-C1.Depending on their post-formation thermal histories, the Al-Mg and Pb-Pb systems in chondrules may or may not have closed simultaneously. Table 1: (26Al/27Al)0, Pb-Pb ages of selected samplesSample(26Al/27Al)0/ 10-6RefPb-PbRefD’Orbigny3.98±0.15104563.43±0.19♮10-12SAH 995553.64±0.18104563.88±0.2712NWA16705.92±0.59104564.39±0.24*10Asuka 88139413.1±0.5613-154564.98±0.1715NWA 73253.03±0.14164563.4±2.616NWA 27963.94±0.16174562.89±0.5917NWA 67043.03±0.14184562.76±0.26192-C17.56±1.53204567.57±0.56*215-C27.04±1.51204567.54±0.52*213-C58.85±1.83204566.20±0.63*2111-C15.55±1.84204565.84±0.72*21*regression based on one subset of data points ♮weighted mean of two datasets Pb-Pb ages are proportional to the intercept of the line formed by linear regression of 207Pb/206Pb vs. 206Pb/204Pb data from various washes, leachates and residues of acid dissolution of a sample. Because contamination by terrestrial or primordial Pb is pervasive, some fractions must be excluded from regressions to ensure a fit with acceptable mean squares weighted deviation (MSWD). Usually points are excluded based on low [Pb], or low 206Pb/204Pb ratio(low radiogenic component), with single outliers identified [11,12,15,16,17]. In the starred examples (Table 1)and the case of 3 CAI Pb-Pb ages [5], up to half the points were excluded solely because did not fit a pre-determined line. This approach is vulnerable to confirmation bias and produces fits with low MSWD and too-low Pb-Pb age uncertainty. Regressing the same data points as [10], were produce the Pb-Pb age of NWA 1670 of 4564.39±0.24 Myr. But selecting other combinations of data points, other, equally valid, isochrons yield ages from 4563.77±0.21 Myr to 4564.64±0.23 Myr. Similar arguments apply to the Pb-Pb isochrons built by [21] for chondrules 2-C1 (we find 4567.33±0.44 to 4567.85±0.46 Myr), 5-C2 (4566.84±0.53 to 4567.70±0.44 Myr), 3-C5 (4565.84±0.54to 4567.04±0.54)and 11-C1 (4565.36±0.51 to 4565.74±0.45 Myr). Our adopted ages for these and NWA 1670 are listed in Table 2.Table 2. tCAI estimated from various components, using our regressions for the chondrules & NWA 1670.SampleΔt26(Myr)tPb(Myr)t’CAI(Myr)D’Orbigny5.05±0.044563.43±0.194568.48±0.19SAH 995555.14±0.054563.88±0.274569.02±0.27NWA16704.64±0.104564.21±0.634568.85±0.67Asuka 8813943.81±0.044564.98±0.174568.79±0.17NWA 73255.33±0.054563.4±2.64568.7±2.6NWA 27965.06±0.044562.89±0.594567.95±0.59NWA 67045.29±0.134562.76±0.264568.05±0.292-C12.00±0.214567.59±0.704569.59±0.725-C22.07±0.224567.23±0.914569.30±0.933-C51.84±0.214566.44±1.124568.28±1.1411-C12.32±0.344565.52±0.664567.84±0.73achondrite4568.72±0.16chondrules4568.76±0.58combined4568.73±0.16A weighted average of the five NC achondrites(or just D’Orbigny, SAH 99555 and Asuka 881394)yields t*CAI= 4568.72 ± 0.16Myr. All are consistent with this value to within 1.8σ, and MSWD=1.5. Including the 4 U-corrected chondrules, t*CAI= 4568.73± 0.16Myrwith MSWD=1.66, which is statistically significant. All chondrules and NC achondrites are consistent with this to within 1.8σ, (Figure 1).Figure 1. Al-Mg formation times after t=0 vs. Pb-Pb ages. The five NC achondrites and four chondrules are consistent with a Pb-Pb age of t=0 of 4568.7 Myr. Discussion: The data from achondrites and chondrules are consistent with a single Pb-Pb age at t=0, justifying the assumption of 26Al homogeneity. The age, 4568.7 Myr, is ≈1.4 Myr older than the commonly accepted Pb-Pb age of CAIs that formed with canonical 26Al/27Al at t=0 [3]. Others have interpreted the discrepancy to signify 26Al heterogeneity in the CAI-forming region[5,21]. We suggest instead that CAIs were exposed to transient heating events that reset the Pb-Pb system without disturbing the Al-Mg system. Notably, chondrules typically experienced transient heating at these times in the nebula [22]. If so, direct measurements of CAIs will not yield as reliable a Pb-Pb age of t=0 as statistical approaches like this and that of [9].References:[1] Auer et al. 2009. [2] Jacobsen, B et al. 2008, EPSL 272, 353-364. [3] Tissot, Fet al. 2017, GCA 213, 593-617.[4] Villa, I et al. 2016, GCA 172, 387-392.[5] Amelin, Y et al. 2010, EPSL 300, 343-350.[6] Connelly, Jet al. 2012, Science 338, 651.[7] Bouvier, Aet al. 2011, LPICo 1639, 9054. [8] Bouvier, Aand Wadhwa, M2010, Nat Geosci 3, 637-641.[9] Nyquist, Let al. 2009, GCA 73, 5115-5136.[10] Schiller et al. 2015. [11] Wadhwa & Brennecka 2012. [12] Tissot et al. 2017. [13] Nyquist et al. 2003. [14] Wadhwa et al. 2009. [15] Wimpenny et al. 2019, GCA 244, 478-501.[16] Koefoed et al. 2016, GCA 183, 31-45. [17] Bouvier, A et al. 2011, GCA 75, 5310-5323. [18] Sanborn, Met al. 2019, GCA 245, 577-596. [19] Amelin, Yet al. 2019, GCA 245, 628-642.[20] Bollard, Jet al. 2017, Sci Adv 3 ,e1700407. [21] Bollard, Jet al. 2019, GCA 260, 62-83.[22] Villeneuve, J et al. 2009, Science 325, 985

S. J. Desch↗

Let's Exploit Available Knowledge on Vegetation Fluorescence

The potential to measure vegetation fluorescence from space (1) and to derive from it direct information on the gross primary productivity (GPP) of terrestrial ecosystems is probably the most thrilling development in remote sensing and global ecology of recent years, as it moves Earth observation techniques from the detection of canopy biophysics (e.g., fraction of absorbed radiation) and biochemistry (chlorophyll and nitrogen content) to the realm of ecosystem function. The existence of a functional relationship between fluorescence and photosynthesis has been elucidated over the last decade by several laboratories, notably as part of the preliminary studies of the European Space Agency Fluorescence Explorer (FLEX) Earth Explorer Mission. The empirical observation presented by Guanter et al. (2) of a linear relationship between fluorescence radiance and GPP, however, provides the first experimental confirmation of the feasibility of the approach— already thoroughly tested at leaf level—at the desired scale, despite the confounding effects associated with the satellite detection of such a faint signal. A word of clarification is needed here. The use of fluorescence as a probe of leaf photochemistry has been a staple of plant ecophysiology for decades, rooted in a sound understanding of photosynthetic energy dissipation. However, most past studies had to rely for the interpretation of results on active (pulse-saturated) techniques, making them unsuitable for remote-sensing applications. Over recent years, however, novel process based models have been developed for the interpretation of steady-state, solar-induced fluorescence at the leaf to canopy scale (3). We are therefore in a position to move beyond the mere empirical observation of an association between GPP and fluorescence radiance. In particular, Guanter et al. (2) base their analysis on the assumption of a constant ratio between photosynthetic and fluorescence light use efficiencies (equation 3 in ref. 2).We know, however, that the ratio is not constant, but changes widely in response to light, CO2, stomatal limitations, and extreme stress (4, 5). What’s more, we can make sense of such changes, thus extracting valuable information from the very scatter that is apparent in their data. However, this process will require the availability of more tailored instruments, such as the one planned for the FLEX mission. As already stressed by Guanter et al. (2), the spatial resolution of the Global Ozone Monitoring Experiment-2 sensor (40 × 80 km) makes it difficult to compare meaningfully the fluorescence signal with ground measurements, when only 60–70% of the footprint consists of the desired land-cover type (table S1 in ref. 2), suggesting that this could be largely responsible for the low signals observed in European grasslands. Moreover, the overpass time of the MetOp-A satellite (9:30 AM) implies that fluorescence is generally measured under light-limiting conditions, when fluorescence is only marginally affected by stomatal closure even under stress conditions. This result could explain the seasonal mismatch with daily GPP observed in natural ecosystems in the absence of irrigation (figure 4 in ref. 2). We hope, therefore, that this welcome contribution to this fast-advancing field will help demonstrate the potential of the new technique, and pave the way for more refined studies under both a technological and scientific point of view.

Space↗

Exploration of the Theoretical Physical Capacity of the John F. Kennedy International Airport Runway System

A design study was completed to explore the theoretical physical capacity (TPC) of the John F. Kennedy International Airport (KJFK) runway system for a northflow configuration assuming impedance-free (to throughput) air traffic control functionality. Individual runways were modeled using an agent-based, airspace simulation tool, the Airspace Concept Evaluation System (ACES), with all runways conducting both departures and arrivals on a first-come first-served (FCFS) scheduling basis. A realistic future flight schedule was expanded to 3.5 times the traffic level of a selected baseline day, September 26, 2006, to provide a steady overdemand state for KJFK runways. Rules constraining departure and arrival operations were defined to reflect physical limits beyond which safe operations could no longer be assumed. Safety buffers to account for all sources of operational variability were not included in the TPC estimate. Visual approaches were assumed for all arrivals to minimize inter-arrival spacing. Parallel runway operations were assumed to be independent based on lateral spacing distances. Resulting time intervals between successive airport operations were primarily constrained by same-runway and then by intersecting-runway spacing requirements. The resulting physical runway capacity approximates a theoretical limit that cannot be exceeded without modifying runway interaction assumptions. Comparison with current KJFK operational limits for a north-flow runway configuration indicates a substantial throughput gap of approximately 48%. This gap may be further analyzed to determine which part may be feasibly bridged through the deployment of advanced systems and procedures, and which part cannot, because it is either impossible or not cost-effective to control. Advanced systems for bridging the throughput gap may be conceptualized and simulated using this same experimental setup to estimate the level of gap closure achieved.

Neitzke, Kurt W.↗

Arc electrode interaction study

The project consisted of two parts: (1) the cathode interaction studies which were a continuation of previous work and had the objective of increasing our understanding of the microscopic phenomena controlling cathode erosion in arc jet thrusters, and (2) the studies of the anode attachment in arc jet thrusters. The cathode interaction studies consisted of (1) a continuation of some modeling work in which the previously derived model for the cathode heating was applied to some specific gases and electrode materials, and (2) experimental work in which various diagnostics was applied to the cathode. The specific diagnostics used were observation of the cathode tip during arcing using a Laser Strobe Video system in conjunction with a tele-microscope, a monochromator with an optical multichannel analyzer for the determination of the cathode temperature distribution, and various ex situ materials analysis methods. The emphasis of our effort was shifted to the cathode materials analysis because a parallel project was in place during the second half of 1993 with a visiting scientist pursuing arc electrode materials studies. As a consequence, the diagnostic investigations of the arc in front of the cathode had to be postponed to the first half of 1994, and we are presently preparing these measurements. The results of last year's study showed some unexpected effects influencing the cathode erosion behavior, such as increased erosion away from the cathode tip, and our understanding of these effects should improve our ability to control cathode erosion. The arc jet anode attachment studies concentrated on diagnostics of the instabilities in subsonic anode attachment arc jet thrusters, and were supplemental measurements to work which was performed by one of the authors who spent the summer as an intern at NASA Lewis Research Center. A summary of the results obtained during the internship are included because they formed an integral part of the study. Two tasks for 1994, the diagnostics of the anode closure phenomenon, and the use of arc jet thrusters for the deposition of c-BN, are being prepared.

Zhou, X.↗

Cyber-Threat Assessment for the Air Traffic Management System: A Network Controls Approach

Air transportation networks are being disrupted with increasing frequency by failures in their cyber- (computing, communication, control) systems. Whether these cyber- failures arise due to deliberate attacks or incidental errors, they can have far-reaching impact on the performance of the air traffic control and management systems. For instance, a computer failure in the Washington DC Air Route Traffic Control Center (ZDC) on August 15, 2015, caused nearly complete closure of the Centers airspace for several hours. This closure had a propagative impact across the United States National Airspace System, causing changed congestion patterns and requiring placement of a suite of traffic management initiatives to address the capacity reduction and congestion. A snapshot of traffic on that day clearly shows the closure of the ZDC airspace and the resulting congestion at its boundary, which required augmented traffic management at multiple locations. Cyber- events also have important ramifications for private stakeholders, particularly the airlines. During the last few months, computer-system issues have caused several airlines fleets to be grounded for significant periods of time: these include United Airlines (twice), LOT Polish Airlines, and American Airlines. Delays and regional stoppages due to cyber- events are even more common, and may have myriad causes (e.g., failure of the Department of Homeland Security systems needed for security check of passengers, see [3]). The growing frequency of cyber- disruptions in the air transportation system reflects a much broader trend in the modern society: cyber- failures and threats are becoming increasingly pervasive, varied, and impactful. In consequence, an intense effort is underway to develop secure and resilient cyber- systems that can protect against, detect, and remove threats, see e.g. and its many citations. The outcomes of this wide effort on cyber- security are applicable to the air transportation infrastructure, and indeed security solutions are being implemented in the current system. While these security solutions are important, they only provide a piecemeal solution. Particular computers or communication channels are protected from particular attacks, without a holistic view of the air transportation infrastructure. On the other hand, the above-listed incidents highlight that a holistic approach is needed, for several reasons. First, the air transportation infrastructure is a large scale cyber-physical system with multiple stakeholders and diverse legacy assets. It is impractical to protect every cyber- asset from known and unknown disruptions, and instead a strategic view of security is needed. Second, disruptions to the cyber- system can incur complex propagative impacts across the air transportation network, including its physical and human assets. Also, these implications of cyber- events are exacerbated or modulated by other disruptions and operational specifics, e.g. severe weather, operator fatigue or error, etc. These characteristics motivate a holistic and strategic perspective on protecting the air transportation infrastructure from cyber- events. The analysis of cyber- threats to the air traffic system is also inextricably tied to the integration of new autonomy into the airspace. The replacement of human operators with cyber functions leaves the network open to new cyber threats, which must be modeled and managed. Paradoxically, the mitigation of cyber events in the airspace will also likely require additional autonomy, given the fast time scale and myriad pathways of cyber-attacks which must be managed. The assessment of new vulnerabilities upon integration of new autonomy is also a key motivation for a holistic perspective on cyber threats.

Complex Networks↗

4.32 Billion Year Old Impact Melts at Apollo 14: Dating the Procellarum Basin?

Introduction: Recent work has improved our understanding of lunar crustal structure and basin evolution on the Moon, but its early impact history and implications for solar system dynamics remain poorly established [1,2]. In particular, there is a derth of absolute ages, especially for the period >4.0 Ga, that could be related to basin-forming events. Here we present U-Pb isotopic data for Zr-rich minerals found in impact-melt fragments from Apollo 14 soil sample 14163 that, when combined with previously published data, dates an impact event at 4324±15 Ma. The compositions and ages of these impact-melt fragments provide unique information about the timing of early impact events and the composition of the lunar crust. Results: Fourteen rocklets ranging in size from 1 to 3 mm were extracted from soil 14163. Most are impactmelt rocks with 10-30% clasts of mostly pyroxene and plagioclase (100 to 500 μm) in a crystalline matrix. Textures of the matrix vary from subophitic, formed by intergrowths of 10-20 μm plagioclase and pyroxene crystals, to poikilitic, with plagioclase and pyroxene reaching 20-50 μm size. Some fragments that contain no visible clasts have similar textures are also interpreted as impact-melt rocks. All fragments contain notably large proportions of euhedral to subhedral ilmenite, zircon, apatite/merrillite and less abundant zirconolite and baddeleyite in the melt matrices. These grains often form intricate intergrowths with each other and rock-forming minerals, indicating their crystallization from the melt (Fig. 1). However, some slightly larger (~50 μm) zircon and phosphate grains can be interpreted as relict clasts based on their granular textures and relationships with the surrounding phases (Fig. 1). This implies that zircon and phosphate minerals were present in the target rocks. Some of these grains were profoundly remelted during the impact, which resulted in an oversaturation of the melt in Zr and P and crystallization of new grains of Zrrich minerals and phosphates during melt solidification. The presence of zircon and phosphates indicates that the melts were significantly enriched in KREEP components and ties their provenance to the Procellarum- KREEP Terrane (PKT) [3]. Combined U-Pb zircon data obtained for all fragments form two main clusters on a concordia diagram (Fig. 2), between about 4.3 and 3.9 Ga. Phosphate data mostly concentrate on the younger end of this age range but some analyses are almost as old as the older zircon grains (Fig. 2). Figure 1: Zircon (Zr) and phosphate (Apt; Mer) grains in impact melt fragments from Apollo 14 soil 14161. A-zircon grains crystallized from impact melt; B-phosphate grain inherited from the target; C- granular zircon grains Our interpretation of these data is that all U-rich minerals experienced variable resetting of the U-Pb system, with phosphate, where closure temperature is significantly lower than that in zircon, affected more profoundly by Pb loss than zircon. Our best estimate of the time of formation of zircon and phosphate is based on statistically valid analysis of grains interpreted texturally as grown from the impact melt (Fig. 2). Ten of these analyses define an age of 4324±15 Ma (MSWD=3.0, probability of fit P=0.002). Our best estimate for the time of resetting is obtained by combining data from phosphates that are statistically indistinguishable from 3.9 Ga within the analytical uncertainties. This group is represented by 18 analyses of 15 phosphate grains from different fragments (Fig. 2) and defines an age of 3922±6 Ma (MSWD=1.2, P=0.23). Conversely, a minimum age of the target lithologies, remelted in the impact that produced the rocklets studied here (Fig. 2), can be determined from the five oldest analyses of zircon clasts at 4338±13 Ma (MSWD=1.5, P=0.2), which is indistinguishable from the age of the impact melt within the uncertainties. The obtained ages, combined with textural evidence, imply that the impact melt was formed at 4324±15 Ma, and that it occurred in a zircon-rich target with a minimum age of 4338±13 Ma. Further reworking occurred during a second impact event at 3922±6 Ma. Discussion: The younger age of 3922±6 Ma can be interpreted as the time of the Imbrium impact. The older age of 4324±15 Ma would then be the time of formation of the impact melt, which was then caught in the Imbrium ejecta either at the Apollo 14 landing site or within the target rocks of the Imbrium impact. The presence of abundant accessory phases such as zircons and phosphates is consistent with a substantial KREEP component in the analyzed particles. The current distribution of KREEP on the lunar surface appears to be strongly influenced by Imbrum ejecta [4] so interpretation of the 4.32 Ga age depends in part on assumed structure of the pre-impact crust and distribution of KREEP within the crust. If KREEP was present only in the deep crust at 4.32 Ga, then a basin-scale impact possibly analogous to Imbrium or larger may be necessary to excavate a KREEPy impact melt at this time. However, if KREEPy materials were present closer to the surface perhaps due to redistribution related to Mg-suite magmatism, then smaller impacts might be able to rework KREEPy compositions at shallower depths. The coherence of the data on the particles analysed here suggests a large volume of melt that has been preserved since 4.32 Ga, consistent with a large impact event. Relict zircons and mineral clasts suggests that the igneous crust in the vicinity of this impact was well developed by at least 4.34 Ga, similar to the model age of KREEP and older than the isochron ages of many Mg-suite cumulates [5]. A problem that confronts all lunar sample studies using the current collection is that the pre-Imbrium geology of the PKT (the source of these 4.32 Ga impact melts) is not well constrained. The South Pole-Aitken basin contains regions that are moderately enriched in Th, but its ejecta is Th-poor [6]; therefore these A14 fragments are probably not SPA ejecta. Alternatively, these fragments may represent formation of a hypothesized Procellarum basin [7] although the lack of a clearly defined basin ejecta signature is a potential problem with that interpretation. In any case, they provide a younger limit on the age of lunar differentiation and formation of KREEP within the lunar crust. Figure 2: U-Pb data for zircon and phosphate grains from impact melt fragments. A-all data; B-data used for age calculations. References: [1] Orgel C., Michael G., Fassett C. I., van der Bogert C. H., Riedel C., Kneissl T., and Hiesinger H. (2018) J. Geophys. Res. Planets 123, 748- 762. [2] Evans A. J., Andrews-Hanna J. C., Head J. W., Soderblom J. M., Solomon S. C., and Zuber M. T. (2018) J. Geophys. Planets. 123, 1596-1617. [3] Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L. and Wieczorek, M.A. (2000) J. Geophys. Res: Planets 105, 4197-4216. [4] Haskin L. A. (1998) J. Geophys. Res. Planets 103, 1679-1689. [5] Borg L.E., Gaffney A.M., and Shearer C.K. (2015) MAPS 50, 715-732. [6] Moriarty, D.P., Watkins, R.N., Valencia, S.N., Kendall, J.D., Evans, A.J., Dygert, N. and Petro, N.E. (2021) J. Geophys. Res. Planets 126. [7] Zhu, M.H., Wünnemann, K., Potter, R.W., Kleine, T. and Morbidelli, A. (2019) J. Geophys. Res: Planets 124, 2117-2140

M D Norman↗