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

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

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.

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.

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.

This dissertation presents three different aspects of the incorporation of directionality into the Monte Carlo ray-trace (MCRT) environment: (1) the development of a methodology for using directional surface optical data, (2) the measurement of the bi-directional reflectivity functions for two different surfaces, and (3) MCRT simulations performed using these directional data sets. The methodology presented is based upon a rigorous analytical formulation and is capable of performing simulations of radiation exchange involving directional emission, absorption and reflection given the bi-directional reflectivity functions (BDRF) of the participating surfaces. A wavelet compression technique is presented for the management of extremely large directional data sets. The BDRFs of two different surfaces were acquired using a Surface Optics Corporation model SOC-250 bi-directional reflectometer. These data were processed according to the methodology presented and an MCRT code was used to simulate the action of the SOC-250 in measuring radiant energy reflected from the surfaces of the two samples when illuminated by the source of the SOC-250. Another MCRT code was used to simulate the radiant energy reflected into a plane at the exit of an open-ended rectangular box when the entrance to the box is illuminated by source of the SOC-250. The RMS error between the MCRT simulations of sampling using the SOC-250 and the measured data were determined and then divided by the mean BDRF level of the measured data (RMS/mean[rho]) to provide an estimate of convergence. The RMS/mean[rho] was observed to fall from as much as 138 to 0.84 for the aluminum substrate coated with Krylon Shortcuts Hunter Green Satin aerosol paint as the number of energy bundles emitted in the MCRT simulation went from 103 to 106 at an incident zenith angle of 40 deg. The RMS/mean[rho] was observed to fall from as much as 2.2 to 0.2 for the Norton (150 Fine grit) all-purpose sandpaper coated with Krylon Shortcuts Hunter Green Satin aerosol paint as the number of energy bundles emitted in the MCRT simulation went from 103 to 106 at an incident zenith angle of 40 deg. / Ph. D.

Two aspects of the study of the Earth radiation budget and the effects of clouds on our climate system are considered in this dissertation : instrumentation and data interpretation. Numerical models have been developed to characterize the optical/thermal-radiative behavior, the dynamic electrothermal response and the structural thermal transients of radiometric channels. These models, applied to a satellite-borne scanning radiometer, are used to determine the instrument point spread function and the potential for optical and thermal-radiative contamination of the signal due to out-of-field radiation and emission from the radiometer structure. The capabilities of the model are demonstrated by scanning realistic Earth scenes. In addition, the optical/thermal-radiative model is used for the development of an infrared field radiometer to interpret results from the experimental characterization of the instrument. The model allowed the sensitivity of the instrument response to assembly uncertainties to be determined. Data processing consists of converting radiometric data into estimates of the flux at the top of the atmosphere. Primary error sources are associated with the procedures used to compensate for unsampled data. The time interpolation algorithm applied to a limited number of observations can produce significantly biased estimates of monthly mean fluxes. A diurnal interpolation protocol using correlative ISCCP cloudiness data is developed to compensate for sparse temporal sampling of Earth radiation budget data. The bias is shown to be significantly reduced in regions where the variability of the cloud cover is well accounted for by ISCCP data. / Ph. D.