Accepted for the Graduate School Abu S.M. Masud, Interim Dean iii DEDICATION To my parents iv ACKNOWLEDGEMENTS I would first like to acknowledge God and then, secondly, all my family who have supported me and given me strength to come to this point. This dissertation would not have been possible without the guidance and help of several individuals who, in one way or another, have contributed and extended their valuable assistance in the preparation and completion of this study. My first debt of gratitude goes to my advisor, Professor Chunseng Ma, who has truly been an inspiration. Without his invaluable guidance, this dissertation would not have been possible. I would also like to express my gratitude to all the members of my committee for the contribution of their valuable time.v ABSTRACT There is a great demand for analyzing multivariate measurements observed across space and over time, due to an increasing wealth of multivariate spatial, temporal, or spatiotemporal data, which may be treated as the realizations of vector (multivariate) random fields. As one of the most important random fields in theory and application, Gaussian random field has been extensively investigated in the literature. Non-Gaussian models and random fields are often encountered in many natural and applied science areas, with specific reasons for assuming particular non-Gaussian finite-dimensional distributions in practice.One of the objectives of this dissertation is to introduce a new non-Gaussian vector random field, which belongs to the family of elliptically contoured vector random fields.This new field is named the K-differenced vector random fields, because its finite-dimensional densities are the difference of two Bessel K functions. A K-differenced vector random field is of second-order and allows for any possible correlation structure, just as a Gaussian one does. It includes a Laplace vector random field as a limiting case. This dissertation studies the properties of the K-differenced vector random field and proposes some covariance matrix structures for not only a K-differenced vector random field but also a second-order elliptically contoured one.