Optical Characterization of DebriSat Fragments in Support of Orbital Debris Environmental Models

The NASA Orbital Debris Program Office (ODPO) develops, maintains, and updates orbital debris environmental models, such as the NASA Orbital Debris Engineering Model (ORDEM), to support satellite designers and operators by estimating the risk from orbital debris impacts on their vehicles in orbit. Updates to ORDEM utilize the most recent validated datasets from radar, optical, and in situ sources to provide estimates of the debris flux as a function of size, material density, impact speed, and direction along a mission orbit. On-going efforts within the NASA ODPO to update the next version of ORDEM include a new parameter that highly affects the damage risk – shape. Shape can be binned by material density and size to better understand the damage assessments on spacecraft. The in situ and laboratory research activities at the NASA ODPO are focused on cataloging and characterizing fragments from a laboratory hypervelocity-impact test using a high-fidelity, mock-up satellite, DebriSat, in controlled and instrumented laboratory conditions. DebriSat is representative of present-day, low Earth orbit satellites, having been constructed with modern spacecraft materials and techniques. The DebriSat fragment ensemble provides a variety of shapes, bulk densities, and dimensions. Fragments down to 2 mm in size are being characterized by their physical and derived properties. A subset of fragments is being analyzed further in NASA’s Optical Measurement Center (OMC) using broadband, bidirectional reflectance measurements to provide insight into the optical-based NASA Size Estimation Model. Additionally, pre-impact spectral measurements on a subset of DebriSat materials were acquired for baseline material characterization. This paper provides an overview of DebriSat, the status of the project, and ongoing fragment characterization efforts within the OMC.


Introduction
The DebriSat project is a collaboration of the NASA Orbital Debris Program Office (ODPO), the Air Force Space and Missile Systems Center (SMC), The Aerospace Corporation (Aerospace), the University of Florida (UF), and the Air Force Arnold Engineering Development Complex (AEDC). The project has four primary goals: 1) design and fabricate a 56-kg class spacecraft ("DebriSat") representative of modern spacecraft in the low Earth orbit (LEO) environment; 2) conduct a hypervelocity laboratory impact test to simulate a catastrophic fragmentation event of DebriSat; 3) collect, measure, and characterize all fragments down to 2 mm in size; and 4) use the data to improve space situational awareness applications and satellite breakup models for better orbital debris environment definition [1].
The motivation for the DebriSat project was based on a key impact test series, Satellite Orbital Debris Characterization Impact Test (SOCIT), which was conducted by the Department of Defense (DOD) and NASA at AEDC in 1992 to support the development of satellite breakup models. The main target for SOCIT was a fully functional U.S. Navy Transit 1960's era satellite. The DOD and NASA breakup models based on the SOCIT data have supported many applications and matched on-orbit events reasonably well over the years [1].
New laboratory-based tests to acquire data to improve the existing DOD and NASA breakup models are needed as new materials and construction techniques are developed for modern satellites. The need for such tests is supported also by discrepancies between model predictions and observations of fragments generated from the breakup of modern satellites, including the Iridium 33 and Fengyun 1-C [13].
The DebriSat design was based on a survey of modern satellites in LEO, as discussed in two 2013 University of Florida dissertations: M. Werremeyer's "Design of Sub-systems for a Representative Modern LEO Satellite" and S. Clark's "Design of a Representative LEO Satellite and Hypervelocity Impact Test to Improve the NASA Standard Breakup Model." All major design decisions, including the selection of components, subsystems, mass fractions, structure, and construction methods, were reviewed and approved by Aerospace subject matter experts. In addition, the DebriSat body was covered with multi-layer insulation (MLI) and three solar panels were attached to one side of the main body (Fig. 1).
To reduce the project cost, a decision was made to emulate the majority of components. The emulated components were based on existing designs of flight hardware, including structure, dimensions, materials, and connection mechanisms. At the end of the assembly, DebriSat was subjected to a NASA General Environmental Verification Specification qualification vibration test to ensure the integrity of the structure.
To increase the project's benefits further, Aerospace designed and built a target resembling a launch vehicle upper stage ("DebrisLV") for the pre-test shot. The pre-test shot was conducted two weeks prior to the DebriSat impact test with identical conditions to DebriSat (facility set-up, projectile, and impact speed). The resulting fragments were collected in boxes and will be characterized after DebriSat characterization is complete. Fig. 2 shows the mounting of DebrisLV inside the target chamber. Fig. 1 Illustration of the DebriSat satellite from two different perspectives. View a) shows hexagonal satellite with the composite body panels each measuring 50 cm × 30 cm. The body of the satellite is wrapped in multilayered insulation. View b) shows the satellite with the 50 cm × 90 cm deployable solar panels. In both views various components are identified including the sun sensor, spectrometer, S-band antenna, divert thruster, optical imager, UHF/VHF antenna, and x-band antenna (Credit: University of Florida) Fig. 2 Digital image of DebrisLV mounted inside the impact chamber at AEDC. DebrisLV was used as a pretest shot conducted two-weeks prior to the DebriSat test and built to resemble a vehicle upper stage, named "DebrisLV." The cylindrical target is composed of various metals and is shown hung inside the foam enclosed chamber prior to impact (Credit: Arnold Engineering Development Complex/Air Force) On 1 April and 15 April 2014, respectively, the DebrisLV and DebriSat impacts were successfully carried out at AEDC Range G. To maximize the projectile mass at the 7 km/sec impact speed without a sabot, the AEDC team developed a special projectile design featuring a hollow aluminum cylinder embedded in a nylon cap. The nylon cap served as a bore rider for the aluminum cylinder to prevent hydrogen leakage and to protect the two-stage light gas gun barrel [14]. Table 1 shows a comparison of some SOCIT, DebriSat, and DebrisLV test conditions. After the impacts of DebrisLV and DebriSat, all soft catch panels used inside the chamber to minimize secondary damage, loose fragments, and dust were carefully collected, processed, documented, and placed in bags or plastic containers for shipping to a storage facility. The initial estimates using the NASA Standard Satellite Breakup Model (SSBM) indicated the number of 2 mm (and larger) fragments from DebriSat and DebrisLV were approximately 85,000 and 35,000, respectively. Since August 2014, the UF team has been working to extract fragments from foam panels, characterize each fragment down to 2 mm in size, and upload the data into a database. Figure 3 shows the number of fragments recorded in the database and the estimated numbers of fragments collected to be analyzed in comparison to the SSBM prediction since the year the experiment was conducted through fiscal year (FY) 2019. Note, the need to update the SSBM is evident utilizing laboratory break-up experiments that encompass more modern construction methods and materials, as well as classical represented satellites, such as SOCIT, to improve breakup models used by NASA and DOD.

DebriSat Characterization & Analysis
All fragments with at least one dimension ≥ 2 mm are carefully collected and/or extracted from the foam panels/pieces and assigned unique identification numbers before their physical characteristics are determined. During the characterization process, each  Table 2, associated metadata (e.g., location the fragment was found within the chamber, images of the fragments, etc.) are also archived in the DCS. The initial step involves assessing each uniquely identified fragment in terms of material, shape, and color. Mass measurements are acquired using a suite of scales appropriate for the individual fragment's mass. To optimize the characterization process and implement methodologies to minimize biases and/or errors associated with "human-inthe-loop" activities during measurements, 2D and 3D imagers are utilized to measure the physical size based on the solid body (SB) dimensions (XSB, YSB, and ZSB). From these initial measurements, the characteristic length (L c ), the average of the object's three maximum orthogonal projected dimensions (XDIM, YDIM, and ZDIM), average crosssectional area, volume, and bulk density are calculated and stored in the DCS. The L c is a common parameter used to represent size as it can be translated from observable parameters, such as radar cross section (RCS) via the NASA Size Estimation Model (SEM). This geometrical mean accounts for removing any aspect angle dependencies during observations and provides an averaged size dimension for modeling purposes. Note, the initial set of fragments investigated concentrated on carbon fiber reinforced polymers (CFRP) using the 2D imagers where the third dimension (i.e., ZDIM) could be extracted from the material's density. As the project has progressed from these initial homogenous fragments, mirrors have been implemented in the 2D imaging system to calculate the third dimension. More details on the characterization can be found in [5,17] ( Table 2). The current NASA SSBM was formulated using laboratory tests and ground-based remote measurements of on-orbit fragmentation events to provide an average breakup ensemble for spacecraft and upper stage collisions and explosions. Significant validation efforts against observable data were conducted, including the development of a radar Size Estimation Model (SEM) based on SOCIT fragments. To assess how to best update or improve future environmental models, such as the Orbital Debris Engineering Model (ORDEM) and the SSBM, a statistical assessment must be done to determine the distribution of key parameters (shape, density, and size) using laboratory data.
The SSBM uses L c as the fundamental independent variable in lieu of mass. This choice of L c in preference to mass was driven by the on-orbit observables, namely Two-Line Element (TLE) orbital data sets and RCS time series. More details on the SSBM can be found in [5,17]. In Fig. 4 the cumulative number as a function of L c is shown binned by primary material type. Due to the project focus on CFRP during the initial DebriSat fragment analysis, the number of CFRP fragments may be biased in the overall material assessment. Fig. 5 shows cumulative number as a function of mass, also binned by material type. Currently, CFRP dominates the small fragment material type in both size and mass categories, but as these parameters increase the metals dominate the population. As the characterization process continues, the trend assessment is subject to change. The transition from predominantly CFRP to metal is currently around 2 cm in L c (see Fig. 4) and 10 milligrams (see Fig. 5). Fig. 6 shows the cumulative number distribution of fragments as a function of L c categorized by shape categories. Below L c of approximately 20 mm, the CFRP dominate the material category, which happen to be mostly flat plates and "needlelike" straight rod fragments. As the distribution increases in size, the larger fragments Average cross-sectional area Weighted average of projected areas visible in multiple 2-D images; pixel-to-length (and area) relationship calculated from hardware characteristics and scene geometry (mm 2 ).

Mass Fragment mass (g)
Area-to-mass ratio (AMR) Calculated average cross-sectional area divided by measured mass (mm 2 /g).

Volume
Calculated from point cloud using convex hull and alpha-shape subtraction (mm 3 ).
Bulk Density Measured mass divided by calculated volume (g/mm 3 ).
are classified as parallelepiped/nugget/spheroid. As the fragments continue to be analysed, these classifications will continue to evolve. Further distributions by density and AMR can be viewed in Ref. [17].

Optical Characterization
The NASA ODPO has utilized optical observations of Earth-orbiting objects to better characterize the orbital debris environment from low Earth orbit, to geosynchronous orbit, and to geostationary transfer orbit. Ground-based measurements provide time-dependent orbital parameters and brightness (i.e., magnitudes). The magnitude data is converted into size assuming a range, phase function, and albedo. To better model the optical size to target size, as RCS is equated to L c using the SEM size conversation as a function of radar wavelength, data acquired in the OMC allow optical data products to be measured using space-based illumination conditions, equipment, and techniques that parallel telescopic observations and source-target-sensor orientations (Fig. 7). To improve upon the current optical size estimation, the OMC is investigating new phase functions based on spacecraft material samples and fragments from laboratory hypervelocity impact tests, such as DebriSat, with plans to generate a distribution of phase functions.

Photometric Measurements/BRDF
To improve the sampled brightness of laboratory targets, the aspect-angle dependencies in object measurements are removed by fixing the object to the end effector of a six degree-offreedom robotic arm and sampling a hemisphere of object view orientations. The totality of this sampled hemisphere is analogous to the bidirectional reflectance distribution function (BRDF), defined by Nicodemus, et al. [18], as a radiometric concept that identifies the reflective characteristics of a surface by the ratio of reflected radiance in the direction ( r ,ψ r ) to the radiation flux incident on a surface from direction ( i ,ψ i ) within the solid angle element dω i [10]. In the OMC, a BRDF is constructed as a hemispheric surface interpolated from laboratory measurements made at 21 discrete view orientations (Fig. 8). Intensity curves (i.e., Fig. 6 Cumulative number of fragments versus characteristics length (L c ) as a function of six shape categories as of 16 January 2020. For objects less than 3 mm, the flat plates represent the majority of fragments, followed by straight rods, nuggets, bent rods, bent plates, and flexible shapes, in that order. For fragments larger than 2 cm the nuggets dominate the shape category, with the remaining shapes following with low count numbers linear dataset) represent instrumental data products, whereas magnitude is normalized to a Lambertian calibration target to compute absolute magnitudes for each target's reflectance distribution. By normalizing the BRDF to that of a Lambertian sphere, the BRDF represents all possible values of the object's phase function at each phase angle in the solar range. Previous work studied the relation between an object's surface optical properties on subsequent behaviors exhibited in their phase functions. By imaging several objects with identical sizes and shapes but varying optical properties, it was discovered that the dominance of diffuse versus specular surfaces resulted in unique characteristics in their associated phase function. These preliminary findings suggest that it may be possible to discern an object's material composition from these photometric functions. With the material composition confirmed, the overall geometric albedo, α g , may be subsequently approximated, and the phase function f(ϕ) may also be directly constructed using photometric data collected when sampling the BRDF [11]. Therefore, the OMC intends to image a multitude of fragments generated via high velocity impact testing and cataloging their experimentally derived BRDFs. By imaging these objects, which are representative of objects present in the orbital debris environment, a more general understanding of an object's material composition impact on phase function behavior is sought, with the ultimate goal of providing a means for identifying these characteristics for data collected via telescopic measurements.

Phase Function Investigation
The accepted model for approximating target size from calibrated magnitudes is the optical-based SEM developed by NASA's Orbital Debris Program Office, which is a photometric determination utilizing both observable and assumed properties of the target. While many of the quantities utilized in the SEM are readily accessible via photometric analysis (e.g., range, phase angle, etc.), the quantities α g (geometric albedo) and f(ϕ) (phase function) are largely determined by the target's shape and material composition. These properties cannot be derived solely from telescopic optical measurements, necessitating that assumptions be made to overall optical data products to derive size estimates. Conventionally, a geometric albedo of 0.175 is assumed for all objects, a value found to minimize the fractional error in resulting size estimations [2,16], and the phase function of orbital debris is assumed to be either Lambertian or specular, although it is well understood that for rocket bodies, spacecraft, and complex shapes the phase function is a complex combination of specular and Lambertian phase functions. The geometric albedo is a byproduct of the bond albedo (ratio of incident to scattered power incorporating all wavelengths), where the bond albedo is defined as the product of the geometrical albedo of the surface and the phase function, which defines how the sunlight is scattered by the surface in the direction facing the observer. Of note, the investigation into why the global geometric albedo of 0.175 is used was based on a study of observational data products focused on fragmentation debris acquired from the Liquid Mirror Telescope with an applied Lambertian phase function correction, as well as, SEM-corrected RCS values used to calculate the global albedo value [16]. Several investigations have been initiated to supplement these assumptions in an effort to  [11] improve the accuracy of size estimates. Hejduk, et al., have conducted investigations into correlating material composition with phase function characteristics [9], and the OMC has also sought to provide an experimental approach for estimating an object's phase function [10,11].
The objective of the OMC's phase function investigation was to replace the Lambertian assumption with an experimentally derived phase function approximation generated via imaging. After initially validating the experimental phase function process with a Lambertian sphere, four objects that were representative of orbital debris with varying optical surface properties were imaged. These objects included a small section of CFRP, an intact section of MLI, a telescope frame-cylindrical object resulting from the DebriSat experiment, and a fragment from SOCIT. Imagery and corresponding data results from the telescope frame are shown in Fig. 9 and Table 3. The experimentally determined phase functions resulted in the most accurate size estimations for the CFRP and MLI, and the shapes of their phase functions demonstrated how certain aspects of the object's geometry may be inferred. Using the known size and experimentally determined phase function, the geometric albedo for each object was determined, and this value for α g was then implemented in a second specular/Lambertian phase function size estimation. The experimentally determined values for geometric albedo were much larger than the assumed value of 0.175, likely inflated, in part, by the short distance between the object and the CCD (i.e., range), which created a much larger solid angle subtending the target. The comparative errors using both the assumed and experimentally determined values for α g indicate that the experimental phase function is more sensitive to the geometric albedo than its specular and Lambertian counterparts. Provided these data and results, to date, the path forward will utilize OMC measurements to assess a broader, more statistically relevant sample size that represents orbital debris of various materials, geometries, and sizes. The goal is to utilize the phase function investigations and telescopic data to ascertain appropriate updates and improvements to the current optical SEM. The complexities involved with assessing a "one size fits all" geometric albedo are discussed further in the Section 3.3 Spectral Measurements (Fig. 9).
Further formulations exist which include the assumed quantities incorporated into the SEM, which may prove useful to supplement experimentally derived phase function approximations. The space illuminance ratio [9], for instance, may provide an avenue for further improving phase function approximations by understanding the impact of an object's shape on its phase function characteristics. This formulation includes both the phase function and a constant, k, which is specific to the observed shape, thereby providing an avenue of investigation into this relatively under-studied aspect of phase function behavior. Therefore, the OMC intends to experimentally approximate the shape factors of various basic shapes by imaging a suite of Lambertian objects with varying geometries. By isolating the shape factor as the only unknown term in the space illuminance ratio, it is possible to identify the impact of a target's shape on phase function, which will ultimately aide in the identification of unknown targets, in addition to improving the accuracy of phase function approximations utilized in subsequent size estimations.
While the effort of experimentally approximating the shape factors of basic geometries is a useful demonstration of the impact that a body's shape has on its photometric properties, the utility in applying these approximations directly to size estimations made for unknown orbital debris is complicated. Orbital debris generated via impact, break-up, or any such kinetic mechanism is unlikely to result in bodies that have a simplistic geometric shape, as has been shown with the DebriSat fragments. Therefore, the OMC intends to utilize the multitude of DebriSat fragments, which are categorized into six shape categories. These fragments have been imaged to create precise point clouds used to calculate the L c among other parameters and to validate the laboratory experimentally derived-phase functions against 3D model ray-tracing software. The ray-tracing software will allow for rapid parametric analyses of target characteristics on  subsequent behaviors exhibited in the phase function. While these analyses would initially be confined to the influence of material composition, the OMC intends to pursue a shape factor study in parallel by implementing a triangular meshing approach for object modeling. It is hoped that through utilizing this approach, it is possible to approximate BRDFs of objects that have increasingly complex shapes and thereby increase the accuracy of experimentally derived-phase functions for objects that are more characteristic of orbital debris.

Spectral Measurements
One of the variables is remote optical measurements albedo variations due to space weathering and/or other external events (i.e., collision or explosion) that can affect the optical size calculations. Provided one albedo value is used for optical size calculations, beyond studying phase function and BRDFs, the following provides analysis of the spectral data that was collected on a range of materials pre-laboratory hypervelocity impact of DebriSat [6]. The goal is to provide an overview of typical space craft materials used in the construction of a modern-day satellite to establish if there are groupings or taxonomy-based relationships that could be applied in material assessments using remote observations. Fig. 10 provides an example of a typical set-up in the OMC using an Analytical Spectral Device (ASD) field spectrometer with a range from 350-2500 nm and a resolving power of approximately 200 (corresponding to a spectral resolution of 10 nm at 2000 nm), and 717 channels that serve to obtain baseline spectral data on various material types. The system only needs 210 channels to obtain the desired bandwidth/ spectral resolution, so using 717 channels is over-sampling. Measurements are acquired by placing the target under the quartz lamp illuminator and orienting the spectrometer's fiber feed (mounted in a pistol grip) approximately perpendicular to the target surface. Depending on the reflectance properties of the target/material, the angle of incidence and detector direction are modified to best sample the material with maximum signal and without saturating the system. The data are output to a laptop computer and reduced using in-house developed software to provide the absolute spectral reflectance. The processed data are uploaded to the NASA Johnson Space Center (JSC) Spacecraft Materials Spectral Database and available upon request and approval to U. S. citizens. A subset of spectral reflectance data on typical spacecraft materials used in DebriSat will be presented in Figs. 11, 12, 14 and 15. All figures below are plotted as a function of absolute reflectance on the y-axis, where 1.0 would be a perfect diffuse reflector as calibrated against a Spectralon calibration plate, as a function of wavelength. For values that extend beyond 1.0, the reflectance may have saturated the detector as observed with highly specular material or components, illustrated in Figs. 12 and 15. The Ultra Triple Junction (UTJ) solar cells used for constructing the DebriSat mockup satellite are composed of a GaInP2/GaAs/Ge structure above a germanium substrate and coated with anti-reflective cover glass for maximum absorption [20]. The reflectance spectrum of the DebriSat UTJ solar cell, shown in Fig. 11, was measured and compared to various additional solar cell samples of similar composition for optical data comparison. Germanium-based photovoltaic cells have consistently produced a reflectance spectrum with an apparent peak present between 850 nm to 900 nm, a characteristic seen in all four data samples shown in the figure. The dual junction GaAs/Ge structured solar cell doped above a germanium substrate, and Spectrolab UTJ and Improved Triple Junction (ITJ) solar cells all show a typical carbon-hydrogen (C-H) bond feature near 2300 nm. It is not clear why the DebriSat UTJ solar cell does not show a strong indication of this feature.
Optical reflectance data of DebriSat Aluminum 6061 alloy samples coated with diverse colors are displayed in Fig. 12. Metals are optically opaque, making their distinguished pigmentation discernable by reflectance within the visible region of the spectrum [4]. This is possible because light impinging on the material surface will excite electrons to unoccupied sites of higher energy, resulting in the absorption and reflectance of different light colors [4]. Hence, the blue anodized Aluminum 6061  Reflectance spectrum for DebriSat UTJ solar cell and various solar cell samples of similar material composition. Overall, the DebriSat solar cell had a higher spectral reflectance in comparison to other solar cells, but also did not exhibit a typical increase in reflectivity (slope) near 1600 nm sample produced initial prominent reflectance at~400 nm, whereas the gold-and redcolored alloy samples exhibit heightened reflectance near 550-600 nm (Fig. 12). The magenta-coated aluminum alloy sample resulted in a modest amount of reflectance at 380 nm coupled with substantial reflectance at 600 nm, symbolizing purple and red color components in the magenta pigment. The black anodized Aluminum 6061 sample showed a featureless, low reflectance through the entire visible region (350-700 nm), typical for materials that are black. The clear anodized aluminum alloy produced a higher reflectance response between 350-700 nm, differentiating itself from all other anodized samples, due to the high reflectivity seen in bare metals [4]. A spectral feature common amongst all Aluminum 6061 alloy samples belongs to the absorption feature present at 850 nm; an optical property typical of most aluminum metals [8]. When analyzing optical behavior of the measured alloy samples in the near infrared region beyond 1000 nm, the colored-anodized aluminum samples produced common absorption features at 1400 nm and 1900 nm, associated with water, and a wide absorption feature centered at 2200 nm, indicating the presence of organic content likely resulting from carbon-hydrogen (C-H) bonds or a variation of oxygen-hydrogen (X-O-H) bonds, an unknown element that bonds with O-H [12,15]. Therefore, the clear anodized alloy produced spectrum with no organic features and only the clear indication of aluminum. Using this data to differentiate between types of anodized Al 60601 may be best suited focusing on the visible and near-infrared regime and binning the reflectance spectra into broadband filters to produce color-color index groupings (i.e., B-R vs R-I). The clear anodized/unpainted spectra had a higher reflectance than the colored anodized samples from 350 nm to approximately 600 nm. For wavelengths outside the visible regime, magenta dominated the spectral, followed closely by black and clear anodized aluminum, followed by blue, red, and gold aluminum in that order MLI (i.e., thermal blankets) components have been considered an essential material incorporated on the exterior of countless spacecraft articles. MLI functions as a substance employed to protect components, such as solid rocket motors and propellant tanks, from elevated temperatures due to external environmental sources while hindering the component from experiencing its own loss of heat [7]. A Kapton® polyimide film, manufactured by E.I. du Pont de Nemours and Co., was used in the DebriSat MLI. Trade names and trademarks are used in this report for identification only. Their usage does not constitute an official endorsement, either expressed or implied, by the National Aeronautics and Space Administration.
The DebriSat MLI Kapton polyimide shows optical similarities to that of two other polyimide samples studied previously (Fig. 13). An in-depth study of MLI components, such as the Space Facing Kapton and Spacecraft Facing Kapton data are provided in Rodriguez 2007. In an MLI blanket, the Spacecraft Facing Kapton sample resembles a woven quilt-like structure and it would be used closest to the spacecraft, while the Space Facing Kapton sample refers to the visibly external side of the polyimide sheet covering the spacecraft article and is perforated for thermal regulation [19]. The previously evaluated Space Facing and Spacecraft Facing polyimide data are compared alongside the DebriSat MLI Kapton polyimide for characteristic similarities. Coppercolored Kapton has a known absorption feature near 450 nm, which is seen in all three data samples in Fig. 13. The DebriSat MLI that was measured resembled the space facing structure, but it is interesting to note that the overall reflectance for the DebriSat MLI polyimide sample remains at an intermediate value between Space Facing and Spacecraft Facing samples over the majority of the spectral response. Fig. 13 Reflectance spectrum for DebriSat MLI wrap Kapton polyimide sample compared with two other Kapton polyimide copper-colored samples. The DebriSat sample has an overall higher reflectance throughout the spectral instrument bandpass and exhibits similar absorption features typical of copper-colored polyimide film A Composite Overwrapped Pressure Vessel (COPV) used in DebriSat and visual "gold shimmer" DebriSat composite materials were also compared with similar composite materials in terms of measured absolute spectral reflectance. The DebriSat COPV and gold shimmer composite materials had equal configurations consisting of an Aluminum-6061-T6 honeycombed core embedded within a Toray T1000 composite. Previous studies performed in the realm of material optical properties claim that composites made of carbon fiber produce reflectance independent of surface roughness, are Lambertian rather than specular, and are dependable absorbers comparable to material absorbers designated as good standards for calibration purposes [3]. This is likely the reason we see relatively low absolute reflectance signatures less than 0.45 produced by the DebriSat composites and sample composites within 350-2500 nm in Fig 14. With a tensile strength of 924 ksi and low density of 1.8 g/cm 3 , Toray carbon fiber composites are commonly utilized in space, military, and commercial based applications [21,22]. Therefore, it is essential to discern the optical signature produced by composites to aid in remote observations telescopically. The "DebriSat Composite: Gold Shimmer" appeared dominantly black in color, though had a gold shimmer tint to its surface upon light impinging on it. This gold shimmer was also present on composite samples 1, 2, and 3, plotted in Fig. 14 for comparison, but this was not the case for the black COPV DebriSat composite. Composites are well known to incorporate epoxy resin in their structure for enhanced physical properties and generate a yellow amber hue that is seen in resins having low-molecular-weight or medium viscosity [1,23]. In general composites, predominately black in color, share a Fig. 14 Reflectance spectrum for a set of DebriSat composite material specimen compared with spectra belonging to previously measured spacecraft composite samples. The DebriSat samples (black and gold) were both more reflective over the 350-2500 nm than the comparative samples and showed an overall featureless data curve with a slight increase in slope approaching the near infrared regime featureless spectral signature, although the two DebriSat samples and the flat composite sample show a slight increasing trend through the infrared wavelengths.
Several printed circuit boards of different colors were incorporated into the DebriSat experiment and have had their optical reflectance properties measured. Between all seven circuit boards tested, a total of four different companies supplied the samples, therefore not all circuit boards were provided by the same manufacturer. Yet, similarities in spectral characteristics are identified in most of the circuit board samples, specifically the organic features in the infrared section beyond 800 nm. There are common absorption bands in the different colored circuit boards observed at 1450 nm, 1700 nm, 1900 nm, and~2300 nm in Fig. 15. The absorption band present at 1700 nm, and 2300 nm are likely due to C-H bands, while the features at~1400 nm and 1900 nm suggest the typical H 2 O (water band) presence [15]. Regarding color spectrums, each colored circuit board will have a bandgap due to its respective color; blue has absorption features between~350 nm (below the signal of the ASD spectrometer) and~400 nm, green between 450 nm and 550 nm and red near 600 nm to 650 nm. Like black, brown also has a featureless, low reflectance, but shows a slight, increasing reflectance through the visible and near-infrared wavelengths.

Summary and Path Forward
This paper reviewed the current status of DebriSat, hypervelocity impact experiment on a representative of modern spacecraft in the LEO environment. The data provided the Fig. 15 Reflectance spectrum for multiple DebriSat circuit board components of various colors (blue, red, green, and brown). Circuit boards are nonhomogeneous materials in general, thus the presence, or lack thereof, spectral features may be a result of where the spectrometer field of view was placed on the target's surface. The variation in spectral plots are attributed to the solar of the circuit board measured results to date of the collected materials, sizes, and shapes that represent the fragmentation event.
Photometric-based phase functions and the best-fit for geometric albedo using a target of known size were presented with path forwards for the development of a more sophisticated phase function. These laboratory results would be used in parallel with optical telescopic data to improve the current optical SEM.
Characterizing the complexities using spectral measurements of spacecraft materials was also presented in finding potential groupings/taxonomies. Future research would incorporate color-color index data utilizing the same optical filters (i.e., bandpasses) as optical telescopes from the spectral data to determine if subsets of spacecraft materials could be isolated in this taxonomy research.
The work presented in this paper will be used not only to improve the optical SEM, but to improve the SSBM and incorporate shape parameterization into the next version of ORDEM.
Funding This work was funded by National Aeronautics and Space Administration.

Declarations
Conflict of Interest The authors declare no conflict of, nor competing interest in the creation of this work.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.