Atmospheric turbulence as an agency affecting the propagation of electromagnetic (EM) waves in different regions of the earth relative to the ground plane has been studied extensively over the past several decades. Mathematical models describing turbulence itself relative to EM waves have been developed by a variety of investigators in the last 50 or more years. It turns out that the majority of these models are essentially in the spatial domain, involving transverse spatial coordinates and their spatial frequency counterparts in the spectral domain. Most turbulence models start out by assuming a random dependence of the medium permittivity on the turbulence. This leads to a random model describing what is commonly referred to as the refractive index power density spectrum. It is well known that propagation through standard atmospheric turbulence creates ripples, random distortions, phase variations and also for monochromatic cases scintillations in the recovered signals. One idea that was proposed to the investigators of this research was that perhaps pre-packaging the EM signal inside a trackable chaos waveform might offer some measure of shielding for the signal even as the overall EM wave passes through turbulence. With this objective in mind, this work began by first establishing standard numerical simulations of EM propagation through homogeneous regions upon passage through a variety of apertures. This standard application involved the use of the Fresnel-Kirchhoff diffraction integral implemented in two ways: (a) as a direct propagation from an object to an image plane, and (b) segmented propagation over uniform incremental layers of the medium in the longitudinal direction. The latter approach was put into place in anticipation of the later introduction of a turbulent layer in the system. Following successful implementation of this technique, turbulence was inserted once again in two different ways: (a) assuming a relatively narrow region of turbulence, modeled as a planar random phase screen derived from the use of the well-known modified von Karman spectrum (MVKS) for refractive index; and (b) the case of an extended random region which is modeled by inserting multiple planar random screens along the propagation path. These initial approaches led to the determination of the resulting scalar output fields numerically derived as complex entities. In the first half of this work, time statistics of the scalar fields were obtained by repeating the simulation multiple times on the basis of an assumed (relatively low) frequency of variation of the turbulence phenomenon (of the order of, say, 20-100 Hz). These time statistics were then incorporated into a transfer function model involving two random processes: (a) the MVKS phase turbulence for which the time statistics were derived as mentioned; and (b) a purely time-dependent chaos wave generated via an acousto-optic (A-O) Bragg cell under feedback, whose first-order optical output is thereby encrypted by an input signal waveform. Use the transfer function approach, cross spectral densities and corresponding cross-correlation functions between the two random phenomena were numerically derived with the final cross correlation product containing the vital message information. Retrieving the message signal from the turbulence-chaos cross correlation product became a prohibitive task, and therefore, even though further investigations are needed, a new approach was developed to complete the intended work. In this approach, a modulated carrier wave which is both time and space-dependent, is propagated through a region of homogeneous space, and upon diffraction through the region, is picked up at the receiver and the embedded message is recovered using appropriate electronics. Thereafter, the same process is repeated in the presence of spatial turbulence, and the recovered signal waveforms are averaged over multiple runs of the simulation representing the time statistics of the turbulence. It is demonstrated that signals recovered under varying degrees of turbulence indeed suffer moderate to severe phase and amplitude distortion, as expected. It must be noted that all numerical simulations reported here are based on strictly near-isoplanatic and paraxial or low propagation angle basis, such that the essential turbulence parameter, C_n^(2 ) is h-independent for all practical purposes. In the final application of this strategy, an encrypted chaos wave riding on an optical carrier is propagated through narrow turbulence of varying strengths, and recovered using a chaos-based heterodyne detection technique. It is shown that indeed encapsulation of the message inside the chaos reduces the distortions in the recovered signal which occur when chaos is not used.