O modelo numérico Delft3D-Flow foi utilizado nos estudos da distribuição e variabilidade de propriedades termohalinas, e da circulação tridimensional do sistema estuarino dos rios Caravelas e Peruípe-SERCP. Os resultados do modelo foram validados quantitativamente usando dados experimentais de marés de quadratura e de sizígia de correntes e salinidades em quatro estações fixas (duas em cada um dos estuários de Caravelas e Nova Viçosa). No processamento foram aplicadas inicialmente condições iniciais homogêneas da salinidade, densidade e dos coeficientes cinemáticos de viscosidade e difusividade vertical, e após quatro semanas de simulação foram extraídos resultados de diferentes condições. Esses resultados foram utilizados como condições iniciais mais realistas, tendo o campo espacial variado das propriedades termohalinas, e dos parâmetros físicos. Essas novas condições iniciais permitiram pular o tempo transiente otimizando assim as simulações subsequentes. A validação das simulações foi feita com base nas medições da campanha de verão, por causa do atraso do período chuvoso a vazão do rio Peruípe foi bem representativa para condições de inverno e, portanto, seria redundante usar as forçantes dos rios com valores tão próximos na validação, além da necessidade da demanda de um tempo muito maior nas analise comparativas dos resultados. Durante o verão (janeiro de 2008), os dados de vazão do rio Peruípe obtidos na Agência Nacional de Águas (ANA) levaram a estimativa da descarga fluvial de ~ 20 m3.s-1, com a extrapolação de ~ 4 m3.s-1 para vazão conjunta dos dois principais afluentes do canal estuarino do rio Caravelas (rios Cúpido e Jaburuna). As oscilações da maré foram simuladas satisfatoriamente para as quatro estações de controle utilizadas na validação, com valores médios do parâmetro Skill superiores a 0,97. As amplitudes da maré calculadas numericamente foram de aproximadamente 1,3 m e 2,5 m na quadratura e sizígia, respectivamente. Os resultados numéricos das simulações das velocidades foram melhores em condições de sizígia com valores médios do parâmetro Skill entre 0,77 e 0,93, enquanto que na quadratura esse parâmetro variou entre 0,38 a 0,65. Para a estrutura de salinidade durante as condições de maré de sizígia, melhores resultados no SERCP do skill e, portanto abrangendo as quatro estações de controle, foram confirmados pelos parâmetros Skill médios superiores a 0,83. Na maré de quadratura os correspondentes valores médios também foram relativamente altos, variando entre 0,73-0,85. Entretanto, houve dificuldades em simular adequadamente a alta estratificação (vertical e longitudinal) do estuário de Nova Vicosa, podendo-se atribuir esse fato à maior vazão do rio Peruípe, que é muito maior do que a do rio de Caravelas. Com os resultados da estrutura de salinidade a intrusão da massa de Água Tropical (AT) foi adequadamente representada nos quilômetros iniciais do canal estuarino do rio Caravelas, e na parte interna da desembocadura na região de Nova Viçosa. Movimentos bidirecionais foram observados nos resultados das simulações na preamar e baixamar, na região de desembocadura dos estuários da Caravelas e Nova Viçosa, movimentos estes forçados pela componente baroclínica da força de gradiente de pressão. Vale ressaltar que o efeito baroclínico ficou mais visível nos instantes da baixamar. Na preamar para a maior parte destas duas regiões verificou-se apenas um pequeno desvio no sentido das correntes entre superfície (Z = -0,1) e fundo (Z = -0,9) e em geral inferior a 30 graus. Usando os traçadores lagrangeanos virtuais lançados ao longo dos canais estuarinos das entradas norte (Caravelas) e sul (Nova Viçosa) do SERCP foi obtido como resultado um tempo de residência relativamente pequeno e comparável ao estuário tropical do rio Curimataú. No trecho analisado no canal estuarino do rio Caravelas, a ~ 3 km e ~ 12 km distantes da Boca do Tomba o tempo de residência médio foi ~ 4,2 dias e ~10,3 dias, respectivamente. Já no estuário de Nova Viçosa, ou canal estuarino do rio Peruípe, apenas o trecho inicial de 5 km foi considerado para o lançamento dos traçadores, e os tempos de residência de ~ 1,5 dias e 2,5 dias foram estimados para as posições de ~2,5 km e 5,0 km distantes da desembocadura. O modelo analítico proposto e usado no cálculo de tempo de residência teve resultados comparáveis aos obtidos pelas simulações pelo método lagrangeano do Delft-3D Flow. Diferentemente, o modelo LOICZ apresentou valores bem diferentes do tempo de residência para os seis estuários analisados (Caravelas, Nova Viçosa, Curimataú, Hudson, Conwy and Mersey), em geral bem inferiores aos do modelo proposto, indicando assim uma possível estimativa maior do fluxo difusivo na formulação do modelo LOICZ. / The numerical model Delft-Flow was used to study the spatial distribution and variability of the termohaline properties, and the tridimensional circulation in the estuarine system of the rivers Caravelas and Peruípe. The model results were validated quantitatively using field data of tidal oscilation and currents, and salinity measurements during neap and spring tides at four mooring stations, two in each estuary. In the data processing homogeneous conditions were used initially for the fields of salinity, density and the vertical kinematic coeficients of viscosity and diffusivity, then after four weeks of running simulation those fields were saved with spatially varied conditions. Those results were used for a more realistic initial condition, thus having a varied field of the termohaline properties, and of the physical parameters. The new initial conditions allowed avoiding the transient time and thus optimizing the following simulations. The model evaluation was done based on measurements undertaken during the summer season, because of the delay in the raining season, the Peruípe river discharge was low and more representative of dry season condition, for those reasons we just used the measurement taken in the January 2008. In the summer the measurements of the Peruíbe river discharge taken from the Agência Nacional de Águas (ANA), and the river flow of was estimated at ~ 20 m3.s-1, with the extrapolation of ~ 4 m3.s-1 for the unified flow from the rivers Cúpido and Jaburuna, with contribute to the flow in the Caravelas estuary. Tidal oscillation at neap and spring tides were well simulated for the 4 controling sites used to the validation, and with the mean skill values over 0.97. The tidal amplitudes were numerical calculated and nearly between 1.3 and 2.5 meters at neap and spring tides, respectively. The numerical results of the velocities were better simulated in spring tides with the mean skill values in the range of 0.77 and 0.93, while at neap tides this parameter varied between 0.38 and 0.65. For the salinity structure at spring tides good results in the estuarine system were achieved, hence comprising all the controling station for validation, with all mean skill values over 0.83. At neap tides the corresponding mean skill values were relatively high, varing in the range of 0.73-0.85. In addition, there were difficulties in simulation adequately the highly vertical and longitudinal stratification of the Nova Vicosa estuary. That was caused by the stronger river inflow of the Peruípe River, which is much higher than the inflow by Caravelas River. The intrusion of the Tropical Water mass (TW) was properly simulated given the right distribution of the salinity along the first quilometers of the estuarine channel of Caravelas, and the internal part of the Nova Viçosa estuarine mouth. Bidirectional movements were observed from the model outputs at the high and low tides, at the region near the mouth of both, Caravelas and Nova Viçosa estuary. Those moviments are forced by the baroclinic component of the pressure gradient. Is worth to mention the the baroclinic effect was better observed at the low tides of these regions. At high tide most of areas showed only a small change in the current direction between surface (Z = -0.1) and bottom (Z = -0.9) and usually smaller than 30 degrees. The virtual lagrangian drifters released along the estuarine channels of Caravelas (north) and Nova Viçosa (south) have shown a residence time relatively short, and comparable to the Curimataú tropical estuary. In the released location between 3 and 12 km of distance from the Caravelas mouth, the residence time was in the range of 4.2 and 10.3 days. Differently, at the Nova Viçosa estuary, only the first 5 km were studied and the residence time of 1.5 and 2.5 days was estimated for the respectives positions of 2.5 and 5.0 km away from the mouth. The analytical model proposed and used to calculate the residence time leaded to results comparable to the results from the numerical simulation using the Lagrangian method of the Delft-3D Flow. Differently, the LOICZ model has shown results really different of the residence time for all the six estuaries (Caravelas, Nova Viçosa, Curimataú, Hudson, Conwy and Mersey), and in general much smaller then the results from the proposed formula, thus indicating a possible overestimation of the diffusive flux assumed in the LOICZ formulation.