Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Uncertainty Propagation”

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.

442 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Precision Measurement of the Nickel Spectrum in Cosmic Rays from 8.8 GeV/n to 240 GeV/n with CALET on the International Space Station

The relative abundance of nickel with respect to iron is by far larger than all other trans-iron elements, therefore it provides a favourable opportunity for a low background measurement of its spectrum. Since nickel, as well as iron, is one of the most stable nuclei, the nickel energy spectrum and its relative abundance with respect to iron provide important information to estimate the abundances at the source and to model the propagation of heavy nuclei. However, only a few direct measurements of cosmic-ray nickel at energy larger than∼3 GeV/n are available at present in the literature and they are affected by strong limitations in both energy reach and statistics. In this paper we present a measurement of the differential energy spectrum of nickel from 8.8 GeV/nto 240 GeV/n, carried out with unprecedented precision by the Calorimetric Electron1 Telescope(CALET) in operation on the International Space Station since 2015. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number Z = 40). The particle’s energy is measured by a homogeneous calorimeter (1.2 proton interaction lengths, 27 radiation lengths) preceded by a thin imaging section (3 radiation lengths) providing tracking and energy sampling. This paper follows our previous measurement of the iron spectrum [1] and it extends our investigation on the energy dependence of the spectral index of heavy elements. It reports the analysis of nickel data collected from November 2015 to May 2021 and a detailed assessment of the systematic uncertainties. In the region from 20 GeV/nto 240 GeV/n our present data are compatible within the errors with a single power law with spectral index−2.51±0.07.

O Adriani↗
NASA NTRS2003

Modeling Candle Flame Behavior In Variable Gravity

The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.

Alsairafi, A.↗
NASA NTRS2015

Advanced Concept Studies for Supersonic Commercial Transports Entering Service in the 2018-2020 Period Phase 2

Lockheed Martin Aeronautics Company (LM), working in conjunction with General Electric Global Research (GE GR) and Stanford University, executed a 19 month program responsive to the NASA sponsored "N+2 Supersonic Validation: Advanced Concept Studies for Supersonic Commercial Transports Entering Service in the 2018-2020 Period" contract. The key technical objective of this effort was to validate integrated airframe and propulsion technologies and design methodologies necessary to realize a supersonic vehicle capable of meeting the N+2 environmental and performance goals. The N+2 program is aligned with NASA's Supersonic Project and is focused on providing system level solutions capable of overcoming the efficiency, environmental, and performance barriers to practical supersonic flight. The N+2 environmental and performance goals are outlined in the technical paper, AIAA-2014-2138 (Ref. 1) along with the validated N+2 Phase 2 results. Our Phase 2 efforts built upon our Phase 1 studies (Ref. 2) and successfully demonstrated the ability to design and test realistic configurations capable of shaped sonic booms over the width of the sonic boom carpet. Developing a shaped boom configuration capable of meeting the N+2 shaped boom targets is a key goal for the N+2 program. During the LM Phase 1 effort, LM successfully designed and tested a shaped boom trijet configuration (1021) capable of achieving 85 PLdB under track (forward and aft shock) and up to 28 deg off-track at Mach 1.6. In Phase 2 we developed a refined configuration (1044-2) that extended the under 85 PLdB sonic boom level over the entire carpet of 52 deg off-track at a cruise Mach number of 1.7. Further, the loudness level of the configuration throughout operational conditions calculates to an average of 79 PLdB. These calculations rely on propagation employing Burger's (sBOOM) rounding methodology, and there are indications that the configuration average loudness would actually be 75 PLdB. We also added significant fidelity to the design of the configuration in this phase by performing a low speed wind tunnel test at our LTWT facility in Palmdale, by more complete modelling of propulsion effects in our sonic boom analysis, and by refining our configuration packaging and performance assessments. Working with General Electric, LM performed an assessment of the impact of inlet and nozzle effects on the sonic boom signature of the LM N+2 configurations. Our results indicate that inlet/exhaust streamtube boundary conditions are adequate for conceptual design studies, but realistic propulsion modeling at similar stream-tube conditions does have a small but measurable impact on the sonic boom signature. Previous supersonic transport studies have identified aeroelastic effects as one of the major challenges associated with the long, slender vehicles particularly common with shaped boom aircraft (Ref. 3). Under the Phase 2 effort, we have developed a detailed structural analysis model to evaluate the impact of flexibility and structural considerations on the feasibility of future quiet supersonic transports. We looked in particular at dynamic structural modes and flutter as a failure that must be avoided. We found that for our N+2 design in particular, adequate flutter margin existed. Our flutter margin is large enough to cover uncertainties like large increases in engine weight and the margin is relatively easy to increase with additional stiffening mass. The lack of major aeroelastic problems probably derives somewhat from an early design bias. While shaped boom aircraft require long length, they are not required to be thin. We intentionally developed our structural depths to avoid major flexibility problems. So at the end of Phase 2, we have validated that aeroelastic problems are not necessarily endemic to shaped boom designs. Experimental validation of sonic boom design and analysis techniques was the primary objective of the N+2 Supersonic Validations contract; and in this Phase, LM participated in four high speed wind tunnel tests. The first so-called Parametric Test in the Ames 9x7 tunnel did an exhaustive look at variation effects of the parameters: humidity, total pressure, sample time, spatial averaging distance and number of measurement locations, and more. From the results we learned to obtain data faster and more accurately, and made test condition tolerances easy to meet (eliminating earlier 60 percent wasted time when condition tolerances could not be held). The next two tests used different tunnels. The Ames 11 ft tunnel was used to test lower Mach numbers of 1.2 and 1.4. There were several difficulties using this tunnel for the first time for sonic boom including having to shift the measurement Mach numbers to 1.15 and 1.3 to avoid flow problems. It is believed that the 11 ft could be used successfully to measure sonic boom but there are likely to be a number of test condition restrictions. The Glenn 8x6 ft tunnel was used next and the tunnel has a number of desirable features for sonic boom measurement. While the Ames 9x7 can only test Mach 1.55 to 2.55 and the 11 ft can only test Mach 1.3 and lower, the Glenn 8x6 can test continuously from Mach 0.3 to 2.0. Unfortunately test measurement accuracy was compromised by a reference pressure drift. Post-test analysis revealed that the drift occurred when Mach number drifted slightly. Test measurements indicated that if Mach number drift is eliminated, results from the 8x6 would be more accurate, especially at longer distances, than results from the 9x7. The fourth test in the 9x7, called LM4, used everything we learned to comprehensively and accurately measure our new 1044-02 configuration with a full-carpet shaped signature design. Productivity was 8 times greater than our Phase 1 LM3 test. Measurement accuracy and repeatability was excellent out to 42 in. However, measurements at greater distances require the rail in the aft position and become substantially less accurate. Further signature processing or measurement improvements are needed for beyond near-field signature validation.

Morgenstern, John↗
NASA NTRS2012

Spectrum and Morphology of the Two Brightest Milagro Sources in the Cygnus Region: MGRO J2019+37 and MGRO J2031+41

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1-100 TeV using the last three years of Milagro data (2005-2008), with the detector in its final configuration. An F-test indicates that MGRO J2019+37 is better fit by a power law with an exponential cutoff than by a simple power law. The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of 7(sup +5 sub −2) × 10(exp −10)/ s /sq m/ TeV, a spectral index of 2.0(sup +0.5 sub -10), and a cutoff energy of 29(sup +50 sub −16) TeV. MGRO J2031+41 shows no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of 2.1(sup +0.6 sub −0.6) × 10(exp −10)/ s/sq m/ TeV and a spectral index of 3.22(sup +0.23 sub −0.18. The overall flux is subject to a approx.. 30% systematic uncertainty. The systematic uncertainty on the power-law indices is approx. 0.1. Both uncertainties have been verified with cosmic-ray data. A comparison with previous results from TeV J2032+4130, MGRO J2031+41, and MGRO J2019+37 is also presented.

BRIGHTEST↗
NASA NTRSDate not supplied

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗
NASA NTRS2006

The Ebb and Flow of Tidal Science, and the Impact of Satellite Altimetry

In the years immediately preceding the launches of Geosat and Topex/Poseidon, tidal science had lapsed into a period of uncertainty and discouragement, brought about by the failure of once-exciting new ideas that eventually proved overly optimistic. A long list of outstanding problems presented themselves, but progress had reached a "low water mark". What was lacking was a high-quality global dataset of tidal measurements, which satellite altimetry -- and especially Topex/Poseidon -- provided. With these data in hand, a "flood tide" of marked progress resulted. In this paper we review some of that progress. An important area of progress, with potentially important implications for other areas of physical oceanography, falls under the topic of "energy dissipation." With precise global constraints provided by altimetry -- combined with precise laser tracking of the altimeter, other geodetic satellites like Lageos, as well as the moon -- the planetary energy budgets of both Earth and ocean tides are now well determined. Moreover, the local energy balances, and thus local estimates of tidal dissipation, have now been mapped, although somewhat coarsely, throughout the ocean. This work has pointed to internal-tide generation in the deep ocean as the once missing sink of tidal energy, and has led to a plethora of new observational and theoretical studies of internal tides, and their role in vertical mixing of the deep ocean. The discovery that internal tides, or some part of them, can be directly mapped with an altimeter opens new lines of research on this topic. Low-mode internal tides have been found, at least in some regions, to propagate several thousand kilometers across open ocean. The study of such waves with altimetry gives us a global view heretofore unattainable, allowing strong observational constraints to be placed on possible ocean mixing processes, such as subharmonic instabilities.

Ray, Richard↗
NASA NTRS1994

Scaling of impact damage in fiber composites

Impact damage in fiber composite structures remains of much concern, and is often the limiting factor in establishing allowable strain levels. The complexity of impact damage formation usually dictates that experiments are required, but scaling of results from small laboratory scale specimens to large structures introduces additional uncertainty into the analysis. This presentation gives the results of an analytical and experimental investigation intended to develop procedures for prediction of damage formation and subsequent strength loss, with particular emphasis on scaling of results with respect to structure size. The experimental investigation involved both drop-weight and airgun impact on carbon/epoxy plates and cylinders. Five sizes of plates ranging from 50 by 50 by 1.072 mm to 250 by 250 by 5.36 mm, and two sizes of cylinders with diameters of 96.5 and 319 mm, were employed in the experimental program. Impact tests were carried out over a range of impact conditions, and specimens were inspected for damage by C-scan and deplying. Analysis procedures were developed for both quasistatic and dynamic impacts for both the plates and cylinders. As has been reported previously, comparison of predicted structural response and measured surface strains was quite good over the entire range of sizes employed in the program. The damage formation and strength loss after impact showed a number of interesting features that are significant with respect to scaling of size. The extent of delamination was observed to increase with specimen size more than would be expected if stresses controlled the delamination extent. This was explained on the basis that delamination is controlled by energy release rates, and thus incorporates the usual dependence on the absolute size characteristic of fracture mechanics. Additionally, the experiments indicated that delamination initiated at matrix cracks and is dependent on the absolute size of the ply group thicknesses. Both the initiation and propagation of delamination are seen to be controlled by fracture mechanics parameters, and thus show specific dependence on size that must be accounted for in extrapolating results from laboratory scale tests to full size structures.

Swanson, Stephen R.↗
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
NASA NTRSDate not supplied

Wall-Modeled Large-Eddy Simulations of Jet Noise in Flight Conditions

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics (LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction (PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology (CST) project. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several operating points both for static (no ambient co-flow), and in-flight (M ͚ =0.3 co-flow) conditions. The operating conditions range in jet exit Mach number from 0.38 to 1.1 with nozzle temperature ratios (NTR) from 0.84 to 2.7. The effect of the flight-stream on far-field noise sound spectra is assessed. Comparison of LES predictions to microphone array measurements demonstrate good agreement within the resolved frequency range. A dip in the predicted low frequency noise spectra for observers between 120° and 145° is observed. This dip is smaller for lower Mach numbers and seems to be correlated to the Mach wave radiation angle. The Mach 1.1 jet shows broadband-shock associated noise. While the onset of BBSN appears to be captured correctly in WMLES, some differences in its magnitude and the prominent frequency at which it occurs persist between the experiment and the simulations. The effect of the outer nozzle boundary layer state created by the co-flow is assessed and shows to be important for accurate comparisons with experiments. A change of 2.5dB between a slip-wall condition and a artificially thickened turbulent boundary layer was observed. In addition, simulations were performed with an alternative nozzle geometry that includes a internal plug and has twice the nozzle exit diameter. These two configurations resulted in very comparable spectra which is consistent with experimental observations. A generally stronger deviation from experimental results is observed for in-flight cases compared to static conditions. The applicability and correct usage of acoustic analogies used for far-field propagation with strong turbulent co-flows needs to be investigated more systematically using canonical problems to improve comparisons between experiments and WMLES. This is especially true for coherent noise sources seen in BBSN.

CST↗
NASA NTRS2021

Seasonal Bias in Global Ocean Color Observations

In this study we identify a seasonal bias in the ocean color satellite derived remote sensing reflectances (Rrs(λ); sr-1) at the ocean color validation site, MOBY (Marine Optical BuoY). The seasonal bias in Rrs(λ) is present to varying degrees in all ocean color satellites examined, including VIIRS (Visible Infrared Imaging Radiometer Suite), SeaWiFS (Sea-viewing Wide Field-of-view Sensor), and MODIS (Moderate Resolution Imaging Spectrometer). The relative bias in Rrs has spectral dependence. Products derived from Rrs(λ) are affected by the bias to varying degrees, with particulate backscattering varying up to 50% over a year, chlorophyll varying up to 25% over a year, and absorption from phytoplankton or dissolved material varying by up to 15%. The propagation of Rrs(λ) bias into derived products is broadly confirmed on regional and global scales using Argo floats and data from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) instrument aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder) satellite. The artifactual seasonality in ocean color is prominent in areas of low biomass (i.e., subtropical gyres) and is not easily discerned in areas of high biomass. While we have eliminated several candidates that could cause the biases in Rrs(λ), there are still outstanding questions regarding potential contributions from atmospheric corrections. Specifically, we provide evidence that the aquatic bidirectional reflectance distribution function may in part cause the observed seasonal bias, but this does not preclude an additional effect of the aerosol estimation. Our investigation highlights the contributions that atmospheric correction schemes can make in introducing biases in Rrs (λ) and we recommend more simulations to discern these influence Rrs (λ) biases. Community efforts are needed to find the root cause of the seasonal bias because all past, present, and future data are or will be affected until a solution is implemented.

satellite ocean color↗