El trabajo realizado en esta tesis doctoral se basa en caracterizar la variabilidad espacial y temporal que presenta el plasma que se encuentra en la alta atmósfera, región conocida como la ionosfera, mediante el modelado de las leyes físicas que rigen su comportamiento. El parámetro más relevante a la hora de estudiar a la ionosfera es la densidad de electrones, cuya presencia se debe a la ionización de varias de las componentes neutras de la atmósfera al absorber la radiación solar ultravioleta, dando como resultado pares de electrón-ion libre. Luego, como por cada electrón existente en la ionosfera se encuentra un ion, la densidad de electrones libre es simplemente igual a la suma sobre las densidades de las distintas especies ionizadas. Dado que la distribución de los iones es inhomogénea respecto de la altura, la densidad de los electrones está controlada por distintas especies en función de la altura. Por esta razón, en lugar de tratar de modelar de manera directa la densidad de electrones, se la modela de manera indirecta, al calcular la densidad de las especies ionizadas. Este procedimiento permite un mejor entendimiento de los distintos procesos físicos y químicos que suceden en las distintas regiones de la ionosfera. Para este trabajo se eligío modelar 7 especies ionizadas que abarcan las diferentes regiones de la ionosfera: NO+, N2 + y O2 + para modelar la región E y la parte baja de la región F, el O+ y el N+ para modelar la región F y el He+ y el H+ para modelar la región por encima de la capa F. La ecuación que regula el comportamiento de la densidad es la llamada “ecuación de continuidad”, en la que se tiene en cuenta varios procesos: i) la producción de iones por fotoionización y reacciones de intercambio, ii) la pérdida de iones debido a interacciones químicas con los electrones y con los elementos neutros de la atmósfera, iii) efectos de transporte por movimientos convectivos originados en la interacción entre el campo magnético terrestre, los campos eléctricos presentes en la ionosfera y la acción de los vientos neutros. Se elige un sistema coordenado asociado con el campo magnético de la Tierra lo que permite llevar el problema de 3 dimensiones a 2 dimensiones y presenta la ventaja de desacoplar los movimientos en las direcciones paralela y perpendicular a las líneas de campo magnético. Como las especies ionizadas interactúan entre si a través de reacciones químicas y de colisiones, las ecuaciones para cada especie están acopladas y es necesario resolver el sistema en forma simultánea. De este manera, el sistema a resolver queda conformado por ecuaciones diferenciales parciales en su forma conservativa. Se utiliza el Método de Volúmenes Finitos con un esquema de <i>upwind</i> para su resolución. El resultado final es un modelo que permite estimar el comportamiento que presenta la ionosfera a latitudes bajas y medias a través de la evolución temporal y espacial de las densidades para 7 especies ionizadas y para los electrones sin la necesidad de alimentarlas con mediciones.