The total solar eclipse over the continental United States on August 21, 2017 offered a unique opportunity to study the dependence of the ionospheric density and morphology on incident solar radiation. Unique responses may be witnessed during eclipses, including changes in radio frequency (RF) propagation at high frequency (HF). Such changes in RF propagation were observed by the Super Dual Auroral Radar Network (SuperDARN) radars in Christmas Valley, Oregon and in Fort Hays, Kansas during the 2017 eclipse. At each site, the westward looking radar observed an increase in slant range of the backscattered signal during the eclipse onset followed by a decrease after totality. In order to investigate the underlying processes governing the ionospheric response to the eclipse, we employ the HF propagation toolbox (PHaRLAP), created by Dr. Manuel Cervera, to simulate SuperDARN data for different models of the eclipsed ionosphere. Thus, by invoking different hypotheses and comparing simulated results to SuperDARN measurements, we can study the underlying processes governing the ionosphere and improve our model of the ionospheric responses to an eclipse. This thesis presents three studies using this method: identification of the cause of the increase in slant range observed by SuperDARN during the eclipse; evaluation of different eclipse obscuration models; and quantification of the effect of the neutral wind velocity on the simulated eclipse data. / Master of Science / The ionosphere is the charged layer of the upper atmosphere, which is generated and sustained by sunlight ionizing neutral particles to form a plasma. In the absence of sunlight, ions and electrons can recombine into neutral particles. The total solar eclipse over the continental United States on August 21, 2017 offered a unique opportunity to study the dependence of the ionospheric density and plasma motion on sunlight as the period of the eclipse is much shorter than night. Observations of the ionosphere during past eclipses indicate that unique ionospheric behavior may be witnessed during eclipses, including changes in radio wave propagation for radio waves in the high frequency (HF) regime. Such changes in radio propagation were observed by the Super Dual Auroral Radar Network (SuperDARN) ionospheric HF radars in Christmas Valley, Oregon and in Fort Hays, Kansas during the 2017 eclipse. At each site, the westward looking radar observed an increase in distance that the radio waves traveled before they were reflected back to the radar during the eclipse onset followed by a decrease in this distance after totality. In order to investigate the mechanisms that produce these observed effects, we employed the HF propagation toolbox (PHaRLAP), created by Dr. Manuel Cervera, to simulate radio propagation and generate simulated SuperDARN data for different models of the eclipsed ionosphere. Thus, different models can be tested by comparing simulated data to measured data. Hence, we can study the underlying processes governing the ionosphere and improve our model of the ionospheric responses to an eclipse. This thesis presents three studies using this method to: identify the cause of the increase in the distance radio waves traveled during the eclipse; evaluate different models of change in eclipse magnitude over time; and investigate the effect of the neutral wind velocity on the simulated eclipse data.

The current effort addresses two issues important to the research conducted by the Thermal Radiation Group at Virginia Tech. The first research topic involves the development of a method which can properly model the diffraction of radiation as it enters an instrument aperture. The second topic involves the study of a potential next-generation space-borne radiometric instrument concept. Presented are multiple modeling efforts to describe the diffraction of monochromatic radiant energy passing through an aperture for use in the Monte-Carlo ray-trace environment. Described in detail is a deterministic model based upon Heisenberg's uncertainty principle and the particle theory of light. This method is applicable to either Fraunhofer or Fresnel diffraction situations, but is incapable of predicting the secondary fringes in a diffraction pattern. Also presented is a second diffraction model, based on the Huygens-Fresnel principle with a correcting obliquity factor. This model is useful for predicting Fraunhofer diffraction, and can predict the secondary fringes because it keeps track of phase. NASA is planning for the next-generation of instruments to follow CERES (Clouds and the Earth's Radiant Energy System), an instrument which measures components of the Earth's radiant energy budget in three spectral bands. A potential next-generation concept involves modification of the current CERES instrument to measure in a larger number of wavelength bands. This increased spectral partitioning would be achieved by the addition of filters and detectors to the current CERES geometry. The capacity of the CERES telescope to serve for this purpose is addressed in this thesis. / Master of Science

The survivability of U.S. warships has become a higher priority than ever before. Two ways to improve survivability are to either avoid damage, or to continue to operate after damage has been incurred. This thesis concentrates on the first line of defense, which involves the first of these two approaches. Specifically, this thesis evaluates the extent of threat due to optical contrast with the ocean background. As part of this effort, an MCRT tool was created that allows the user to vary the shape and surface properties of a ship. A reverse MCRT was performed in order to reduce the processing time required to get accurate results. Using this MCRT tool, the user can determine the theoretical contrast with the ocean surface that would be seen at any viewing angle with and without a polarization filter. The contrast due to differential polarization and a change in viewing angle is estimated to determine the extent of threat. These results can be determined for both daytime and nighttime conditions by specifying if the ray trace is in the infrared or visible light range. The location of the sun for daytime conditions, and the temperature of the surfaces for nighttime conditions, can all be adjusted by the user. In order to get an accurate estimation of the signal power coming from the ocean surface, a great deal of time and effort was spent modeling the ocean surface. Many studies have been done concerning the slope statistics of an ocean surface, some more informative than others. This thesis takes two of the most complete studies and brings them together to get accurate slope statistics in both along-wind and crosswind directions. An original idea by the author was used to give a typical shape to the waves of the simulated ocean surface. The surface properties of the ship were determined using Fresnel's equations and the complex index of refraction of water at the particular wavelengths of interest. / Master of Science

Since the early 70's chemistry and transport models (ChTMs) have been proposed and improved. Tropospheric ChTMs for trace species are detailed numerical formulations intended to represent the atmospheric system as a whole, accounting for all the individual processes and phenomena that influence climate changes. The development of computer resources and the retrieval of emission inventories and observational data of the species of interest have enhanced the model evolution towards three-dimensional global models that account for more complicated chemical mechanisms, wet and dry deposition phenomena, and interactions and feedback mechanisms between meteorology and atmospheric chemistry. The purpose of this study is to ascertain the sensitivity of the solar radiative field in the atmosphere to absorption and scattering by aerosols. This effort is preliminary to the study of feedback mechanisms between photolytic processes that create and destroy aerosols and the radiation field itself. In this study, a cloud of water-soluble aerosols, randomly distributed in space within hypothetical 1-cm cubes of atmosphere, is generated. A random radius is assigned to each aerosol according to a lognormal size distribution function. The radiative field characterization is analyzed using a Mie scattering code to determine the scattering phase function and the absorption and scattering coefficients of sulfate aerosols, and a Monte Carlo ray-trace code is used to evaluate the radiative exchange. The ultimate goal of the effort is to create a tool to analyze the vertical distribution of absorption by aerosols in order to determine whether or not feedback between photolytic processes and the radiation field needs to be included in a Third Generation Chemistry and Transport model. / Master of Science

As the planetary climate continues to evolve, it is important to build an accurate long-term climate record. State-of-the-art atmospheric science requires a variety of approaches to the measurement of the atmospheric structure and composition. This thesis supports the possibility of inferring cloud properties from sun photometer observations of the cloud solar aureole using an artificial neural network (ANN). Training of an ANN requires a large number of input and output parameter sets. A cloud radiance model is derived that takes into consideration the cloud depth, the mean size of the cloud water particles, and the cloud liquid water content. The cloud radiance model derived here is capable of considering the wavelength of the incident sunlight and the cloud lateral dimensions as parameters; however, here we consider only one wavelength—550 nm—and one lateral dimension—500 m—to demonstrate its performance. The cloud radiance model is then used to generate solar aureole profiles corresponding to the cloud parameters as they would be observed using a sun photometer. Coefficients representative of the solar aureole profiles may then be used as inputs to a trained ANN to infer the parameters used to generate the profile. This process is demonstrated through examples. A manuscript submitted for possible publication based on an early version of the cloud radiance model was deemed naïve by reviewers, ultimately leading to improvements documented here. / Master of Science / The Earth's climate is driven by heat from the sun and the exchange of heat between the Earth and space. The role of clouds is paramount in this process. One aspect of "cloud forcing" is cloud structure and composition. Required measures may be obtained by satellite or surface-based observations. Described here is the creation of a numerical model that calculates the disposition of individual bundles of light within water clouds. The clouds created in the model are all described by the mean size of the cloud water droplets, the amount of water in the cloud, and cloud depth. Changing these factors relative to each other changes the amount of light that traverses the cloud and the angle at which the individual bundles of light leave the cloud as measured using a device called a sun photometer. The measured amount and angle of bundles of light leaving the cloud are used to recover the parameters that characterize the cloud; i.e., the size of the cloud water droplets, the amount of water in the cloud, and the cloud depth. Two versions of the cloud radiance model are described.

A primary objective of the effort reported here is to develop a radiometric instrument modeling environment to provide complete end-to-end numerical models of radiometric instruments, integrating the optical, electro-thermal, and electronic systems. The modeling environment consists of a Monte Carlo ray-trace (MCRT) model of the optical system coupled to a transient, three-dimensional finite-difference electrothermal model of the detector assembly with an analytic model of the signal-conditioning circuitry. The environment provides a complete simulation of the dynamic optical and electrothermal behavior of the instrument. The modeling environment is used to create an end-to-end model of the CERES scanning radiometer, and its performance is compared to the performance of an operational CERES total channel as a benchmark. A further objective of this effort is to formulate an efficient design environment for radiometric instruments. To this end, the modeling environment is then combined with evolutionary search algorithms known as genetic algorithms (GA's) to develop a methodology for optimal instrument design using high-level radiometric instrument models. GA's are applied to the design of the optical system and detector system separately and to both as an aggregate function with positive results. / Ph. D.

The primary objective of the work reported here is to develop a methodology to predict the uncertainty associated with radiation heat transfer problems solved using the Monte Carlo ray-trace method (MCRT). Four equations are developed to predict the uncertainty of the distribution factor from one surface to another, the global uncertainty of all the distribution factors in an enclosure, the uncertainty of the net heat flux from a surface, and the global uncertainty of the net heat flux from all the surfaces in an enclosure, respectively. Numerical experiments are performed to successfully validate these equations and to study the impact of various parameters such as the number of surfaces in an enclosure, the number of energy bundles traced in the MCRT model, the fractional uncertainty of emissivity and temperature, and the temperature distribution in the enclosure. Finally, the methodology is successfully applied to a detailed MCRT model of a CERES-like radiometer. / Ph. D.

This study was conducted to examine the effect of flow turbulence on sound waves propagating across a velocity field. The resulting information can be used to determine the potential for increasing the accuracy of an ultrasonic flowmeter, and understand the data scatter typically seen when using an ultrasonic flowmeter. A modification of the Ray Trace Method was employed which enabled the use of multiple rays in a very fine grid through a flow field. This technique allowed for the computation of the statistical variation of the propagation times for sound pulses traversing a flow field. The statistical variation was studied using two flow fields: 1) a uniform flow field with a superimposed vortex street and 2) an experimentally measured channel flow. The uniform flow field with a superimposed vortex street allowed for the examination of the effects of a large-scale flow structure on sound wave propagation, and for the verification of the analysis technique. Next by using the measured turbulent channel flow, as an example, the statistical variation of sound pulse propagation time was computed for flow likely to be encountered in actual flow measurement situations. Analysis was also conducted to determine the maximum allowable repetition rate of measurements with regard to the optimal time of flight measurements. Both the propagation time of a sound pulse moving across a uniform flow field with superimposed vortex street, and the resultant computed flow were observed to vary at the same frequency of the vortex street. Further, the magnitude of the variations was proportional with the strength of the individual vortices in the vortex street. A sound pulse propagating back and forth across a measured turbulent channel flow, afforded individual time difference variation from the mean propagation time of up to 5%. It was shown that a minimum variation occurred when the sound pulses were transmitted at a 75 degree angle to the flow axis. It was also determined that the average speed of sound in a flow field affected the final flow measurements by decreasing the measured delta time difference between the upstream and downstream propagating sound waves, and therefore the measured flow. The width of the sound path also contributed to decreasing the variation of the individual measurements by integrating over a larger sound path. These findings suggest that turbulence in a flow field affects ultrasonic flowmeter measurements by creating differences in the propagation times of individual sound pulses. Thus, turbulence and large-scale flow structures can result in variations in volumetric flow rate determination made by an ultrasonic flowmeter system.

The study whose results are presented here was carried out in support of an ongoing larger effort to investigate and understand the impact of coherence and polarization on the performance of instruments intended to monitor the Earth's radiant energy budget. The visibility of fringes produced by a Michelson interferometer is known to be sensitive to the degree to which the incident light beam is monochromatic. Therefore, the Michelson interferometer has significant potential as a tool for quantifying the degree of temporal coherence of a quasi-monochromatic light beam. Simulation of the performance of an optical instrument using the Monte-Carlo ray-trace (MCRT) method has been shown to be an efficient method for transferring knowledge of the coherence state of a beam of light from one instrument to another. The goal of the effort reported here is to create and experimentally validate an MCRT model for the optical performance of a Michelson interferometer. The effort is motivated by the need to consolidate the knowledge and skills of the investigator in the realm of physical optics, and by the need to make a useful analytical tool available to other investigators in the larger effort. / Master of Science

Black coatings essentially absorb incident light at all wavelengths from all directions. They are used when minimal reflection or maximum absorption is desired and therefore are effective in applications that require control of stray light. Our motivation stems from the use of black coating Lord Aeroglaze® Z302 in aerospace and remote sensing applications and the desire to support the development of bidirectional spectral models that can be used successfully to predict the performance of optical instruments such as telescopes. The bidirectional reflectance distribution function (BRDF) is an indispensable parameter in the optical characterization of such coatings. The current effort involves investigation of the BRDF of the commercial black coating Aeroglaze® Z302. An automated goniometer reflectometer has been designed, fabricated and successfully used for performing the BRDF measurements of Z302 at visible and ultraviolet wavelengths and at both polarizations. The current contribution involves study of Z302 samples prepared at different thicknesses and by different methods, which provides insight about influence of surface roughness on BRDF of Z302. / Master of Science / When light falls on different materials it undergoes various phenomenon such as reflection, refraction, absorption and scattering. The amount of each phenomenon varies with the optical nature of a material as well as the wavelength and direction of the light. Therefore, understanding the optical properties of materials at various wavelengths of light is necessary for effectively using those materialsin specific applications which require light to be efficiently reflected or absorbed. This research studies an optical property known as Bidirectional Reflectance Distribution Function (BRDF) of a black coating called Lord Aeroglaze Z302. Black coatings are materials that ideally absorb almost all light that falls on them irrespective of the light’s direction and wavelength. They are used in applications where maximum absorption of light is required. One such application which relates to the motivation for this research is absorbing unwanted light in instruments used in space such as telescopes and radiometers. Z302 is used in the Clouds and the Earth’s Radiant Energy System (CERES) instruments developed by NASA. BRDF is an important parameter which gives information about all other optical properties of a surface and can be used to know optical performance of that surface. The current work describes the experiments and an automated device developed, called reflectometer, to measure the BRDF of Z302 at different angles and wavelengths of light. The results are reported for different thickness samples of Z302 coating, and two different wavelengths of light that belong to the visible and ultraviolet spectrum of light.