Exoplanet discovery has grown more quickly in recent years. However, the nature of their discovery leaves many unanswered in questions regarding exoplanetary habitability. Proxima Centauri b, an exoplanet which orbits the Sun's closest stellar neighbour, Proxima Centauri, was recently discovered with a subzero equilibrium temperature. Although not considered habitable based on the classical definition of the liquid water range, there may be fractions of Proxima Centauri b which are habitable. A prior study simulated the climate conditions of Proxima Centauri b until equilibrium was reached, using a variety of initial conditions. In this project, various metrics for calculating the fractional habitability of Proxima Centauri b are presented and applied to the results of the prior study's simulations. Colormaps are used to show the ice and temperature distributions that produce the calculated values of fractional habitability. The fractional habitabilities calculated show that while the value is both case and metric dependent, for the vast majority of all cases and metrics the value is nonzero implying that Proxima Centauri b is likely to have habitable regions. / Upptäckandet av exoplaneter har ökat i takt över de senaste åren. Samtidigt, på grund av sättet som de upptäcks finns många obesvarade frågor angående planeternas beboelighet. Proxima Centauri b är en exoplanet som kretsar kring solens närmsta granne, Proxima Centauri. Exoplaneten upptäcktes nyligen med en jämviktstemperatur under $0\degree$C. Trots att exoplaneten inte anses beboelig enligt klassisk definition kan det finnas delar av Proxima Centauri b som är beboeliga. En tidigare studie simulerade klimatförhållandena av Proxima Centarui b till jämvikt nåddes, med varierade begynnelsetillstånd. I detta projekt beräknas andelen av Proxima Centauri b som är beboelig genom flera olika mått för "fractional habitability". Måtten jämförs med den tidigare studien och dess simuleringar. Grafiskt åsikdligörs resultaten via färgkartor över planeten för istjocklek och yttemperatur. De beräknade värdena på Proxima Centauri b's "fractional habitability" påvisar beroende på mått och begynnelsetillstånd. Däremot, för en majoritet av både fall och mått är värdet nollskilt vilket antyder att Proxima Centauri b är delvist beboelig.

The influence of anthropogenic emissions on aerosol distributions and the hydrological cycle are examined with a focus on monsoon precipitation over the Indian subcontinent, during January 2001 to December 2005, using the European Centre for Medium-Range Weather Forecasts-Hamburg (ECHAM5.5) general circulation model extended by the Hamburg Aerosol Module (HAM). The seasonal variability of aerosol optical depth (AOD) retrieved from the MODerate Resolution Imaging Spectroradiometer (MODIS) on board the Terra and Aqua satellite is broadly well simulated (R 0.6–0.85) by the model. The spatial distribution and seasonal cycle of the precipitation observed over the Indian region are reasonably well simulated (R 0.5 to 0.8) by the model, while in terms of absolute magnitude, the model underestimates precipitation, in particular in the south-west (SW) monsoon season. The model simulates significant anthropogenic aerosol-induced changes in clear-sky net surface solar radiation (dimming greater than -7 W m-2), which agrees well with the observed trends over the Indian region. A statistically significant decreasing precipitation trend is simulated only for the SWmonsoon season over the central-north Indian region, which is consistent with the observed seasonal trend over the Indian region. In the model, this decrease results from a reduction in convective precipitation, where there is an increase in stratiform cloud droplet number concentration (CDNC) and solar dimming that resulted from increased stability and reduced evaporation. Similarities in spatial patterns suggest that surface cooling, mainly by the aerosol indirect effect, is responsible for this reduction in convective activity. When changes in large-scale dynamics are allowed by slightly disturbing the initial state of the atmosphere, aerosol absorption in addition leads to a further stabilization of the lower troposphere, further reducing convective precipitation.

The scale dependency of variance of total water mixing ratio is explored by analyzing data from a general circulation model (GCM), a numerical weather prediction model (NWP), and large-eddy simulations (LESs). For clarification, direct numerical simulation (DNS) data are additionally included, but the focus is placed on defining a general scaling behavior for scales ranging from global down to cloud resolving. For this, appropriate power-law exponents are determined by calculating and approximating the power density spectrum. The large-scale models (GCM and NWP) show a consistent scaling with a power-law exponent of approximately 22. For the high-resolution LESs, the slope of the power density spectrum shows evidence of being somewhat steeper, although the estimates are more uncertain. Also the transition between resolved and parameterized scales in a current GCM is investigated. Neither a spectral gap nor a strong scale break is found, but a weak scale break at high wavenumbers cannot be excluded. The evaluation of the parameterized total water variance of a state-of-the-art statistical scheme shows that the scale dependency is underestimated by this parameterization. This study and the discovered general scaling behavior emphasize the need for a development of scale-dependent parameterizations.

This study explores the viability of parameter estimation in the comprehensive general circulation model ECHAM6 using ensemble Kalman filter data assimilation techniques. Four closure parameters of the cumulus-convection scheme are estimated using increasingly less idealized scenarios ranging from perfect-model experiments to the assimilation of conventional observations. Updated parameter values from experiments with real observations are used to assess the error of the model state on short 6 h forecasts and on climatological timescales. All parameters converge to their default values in single parameter perfect-model experiments. Estimating parameters simultaneously has a neutral effect on the success of the parameter estimation, but applying an imperfect model deteriorates the assimilation performance. With real observations, single parameter estimation generates the default parameter value in one case, converges to different parameter values in two cases, and diverges in the fourth case. The implementation of the two converging parameters influences the model state: Although the estimated parameter values lead to an overall error reduction on short timescales, the error of the model state increases on climatological timescales.

Documenter la variabilité climatique récente est nécessaire à la compréhension des mécanismes en jeu, associés au rôle du bilan de masse de l’Antarctique pour l’élévation du niveau des mers globale. Les enregistrements issus des carottes peu profondes d’Antarctique sont des données précieuses, complémentaires aux observations instrumentales et satellitaires, pour couvrir en continu l’ensemble du continent. Mesurés le long de ces carottes de glace, les isotopes stables de l’eau sont traditionnellement utilisés pour quantifier les changements passés de la température locale.Cette thèse doctorale a été initiée dans le cadre du programme de l’Agence Nationale de la Recherche ASUMA (“Improving the Accurancy of SUrface Mass balance of Antarctica”), ayant pour objectif de reconstruire et identifier les processus contrôlant la variabilité spatio-temporelle du bilan de masse de surface (BMS) de la Terre Adélie. J’ai utilisé des données d’isotopes stables de l’eau enregistrées dans des carottes de névé, des simulations atmosphériques produites par le modèle atmosphérique de circulation générale de haute résolution ECHAM5-wiso équipé des isotopes stables de l’eau, des réanalyses atmosphériques, des rétro-trajectoires, ainsi que des observations instrumentales satellitaires et de surface.Dans une première partie, j’ai évalué les capacités du modèle ECHAM5-wiso à simuler les températures de l’Antarctique, le BMS, le δ18O et le d-excess (ci-après, d-excess), comme prérequis à l’exploitation du modèle pour interpréter les compositions isotopiques. J’ai développé des diagnostics pour les relations δ18O-température et d-excess- δ18O sur l’ensemble du continent de l’Antarctique, en montrant que les différences issues des pentes des relations δ18O-température spatiales, inter-annuelles et saisonnières. Au sein du groupe de travail international de PAGES (Past Global Changes) Antarctica2k, j’ai utilisé des calibrations établies issues du modèle ECHAM5-wiso pour reconstruire la température de 7 régions d’Antarctique à partir d’une synthèse d’enregistrements de δ18O issus de carottes de glace couvrant les 2 000 dernières années.Dans une seconde partie, de nouveaux enregistrements issus de deux carottes de névé extraites en Terre Adélie, la S1C1 et la TA192A, ont été exploités, couvrant respectivement les périodes 1947-2007 et 1998-2014. Les BMS reconstruits décrivent une grande variabilité spatiale (74,11 ± 14,1 cm w.e. y-1 et 21,8 ± 6,9 cm w.e. y-1 pour la TA192A et la S1C1 respectivement), cohérente avec les données de balise disponibles. En utilisant une base de données mise à jour des isotopes stables de l’eau de l’Antarctique, j’ai montré que les valeurs moyennes isotopiques de Terre Adélie appartiennent à l’intervalle des valeurs côtières de l’Antarctique. Des analyses statistiques montrent une absence de relation entre nos enregistrements avec la température de surface locale à l’échelle inter-annuelle, mais des relations significatives avec des rétro-trajectoires atmosphériques et des simulations isotopiques suggérant que les isotopes de l’eau de la Terre Adélie fournissent des indications de la variabilité de la dynamique atmosphérique et du transport d’humidité, aux échelles saisonnière et inter-annuelle.Les analyses de cette thèse ont été limitées par la quantité d’enregistrements isotopiques disponibles pour la Terre Adélie, ainsi que par le manque de compréhension des effets de dépôt et de post-dépôt. Il est donc nécessaire d’exploiter les nouvelles carottes de névé extraites au cours du programme ASUMA, et d’effectuer en continu des mesures de la composition isotopique des précipitations, de la vapeur d’eau et de la neige de surface de Terre Adélie, en combinaison avec des outils de simulations atmosphériques, tels que des rétro-trajectoires associées à un diagnostic des sources d’humidité, et des modèles atmosphériques de circulation générale et régionaux équipés des isotopes stables de l’eau. / Documenting recent Antarctic climate variability is needed in order to understand the mechanisms at play, associated with the role of Antarctic mass balance for global sea level rise. Proxy records from Antarctic shallow firn cores are precious data, which complement instrumental and remote sensing observations to continuously cover the whole continent. Within these ice cores, water stable isotopes are commonly used to quantify past changes in local temperature.This PhD thesis was initiated within the French Agence Nationale de la Recherche “Improving the Accurancy of SUrface Mass balance of Antarctica” (ASUMA) project, which aims to reconstruct and to identify the processes controlling the spatio-temporal variability of the surface mass balance (SMB) in Adélie Land. I used water stable isotopes records from recently drilled shallow firn cores, as well as atmospheric simulations performed with the high resolution atmospheric general circulation model ECHAM5-wiso model, equipped with water stable isotopes, atmospheric reanalyses and back-trajectories, instrumental and remote sensing climate observations.In a first part, I assessed the skills of the ECHAM5-wiso with respect to Antarctic temperature, SMB, δ18O and deuterium excess (hereafter d-excess), as a prerequisite for the exploitation of the model to interpret isotope compositions. I developed Antarctic-wide diagnostics of the δ18O-temperature and d-excess- δ18O relationships, showing differences in the spatial, seasonal and interannual δ18O-temperature slopes. Within the international working group of PAGES (Past Global Changes) Antarctica 2k, I used the calibrations inferred from ECHAM5-wiso to reconstruct temperatures over 7 Antarctic regions from a synthesis of ice core δ18O records spanning the past 2,000 years.In a second part, new water stable isotope records from two firn core drilled in Adélie Land, the S1C1 and the TA192A, were investigated, covering the periods 1947-2007 and 1998-2014 respectively. The reconstructed SMB display a high spatial variability (74.1 ± 14.1 cm w.e. y-1 and 21.8 ± 6.9 cm w.e. y-1 for the TA192A and S1C1 respectively), consistent with Adélie Land stake data. Using an updated database of Antarctic water stable isotope datasets, I showed that the mean isotopic values (δ18O and d-excess) in Adélie Land are in line with the range of Antarctic coastal values. Statistical analyses show no relationship between our records and local surface air temperature, at the inter-annual scale, but significant relationships with atmospheric back-trajectories and isotopic simulations, suggesting that water stable isotopes in Adélie Land provide fingerprints of the variability of atmospheric dynamics and moisture transport, at the seasonal and inter-annual scales.The analyses performed during this PhD thesis have been limited by the few available Adélie Land water stable isotope records, and by the lack of understanding of deposition and post-deposition processes. Further work is thus needed to exploit the new firn cores drilled within the ASUMA project, and to monitor continuously Adélie Land water stable isotopes in precipitation, surface water vapour and surface snow, in combination with tools of atmospheric simulations such as back-trajectory simulations provided with moisture sources diagnostics, as well as water stable isotopes-enabled atmospheric general and regional circulation models.

Water resource planning and management practices in the southeastern United States may be vulnerable to climate change. This vulnerability has not been quantified, and decision makers, although generally concerned, are unable to appreciate the extent of the possible impact of climate change nor formulate and adopt mitigating management strategies. Thus, this dissertation aims to fulfill this need by generating decision worthy data and information using an integrated climate change assessment framework. To begin this work, we develop a new joint variable spatial downscaling technique for statistically downscaling gridded climatic variables to generate high-resolution, gridded datasets for regional watershed modeling and assessment. The approach differs from previous statistical downscaling methods in that multiple climatic variables are downscaled simultaneously and consistently to produce realistic climate projections. In the bias correction step, JVSD uses a differencing process to create stationary joint cumulative frequency statistics of the variables being downscaled. The functional relationship between these statistics and those of the historical observation period is subsequently used to remove GCM bias. The original variables are recovered through summation of bias corrected differenced sequences. In the spatial disaggregation step, JVSD uses a historical analogue approach, with historical analogues identified simultaneously for all atmospheric fields and over all areas of the basin under study. In the second component of the integrated assessment framework, we develop a data-driven, downward hydrological watershed model for transforming the climate variables obtained from the downscaling procedures to hydrological variables. The watershed model includes several water balance elements with nonlinear storage-release functions. The release functions and parameters are data driven and estimated using a recursive identification methodology suitable for multiple, inter-linked modeling components. The model evolves from larger spatial/temporal scales down to smaller spatial/temporal scales with increasing model structure complexity. For ungauged or poorly-gauged watersheds, we developed and applied regionalization hydrologic models based on stepwise regressions to relate the parameters of the hydrological models to observed watershed responses at specific scales. Finally, we present the climate change assessment results for six river basins in the southeastern United States. The historical (baseline) assessment is based on climatic data for the period 1901 through 2009. The future assessment consists of running the assessment models under all IPCC A1B and A2 climate scenarios for the period from 2000 through 2099. The climate assessment includes temperature, precipitation, and potential evapotranspiration; the hydrology assessment includes primary hydrologic variables (i.e., soil moisture, evapotranspiration, and runoff) for each watershed.

Climatological winds over the equatorial Indian Ocean (EqlO) are westerly most of the year. Twice a year, in April-May ("spring") and October-December ("fall"), strong, sustained westerly winds generate eastward equatorial jets in the ocean. There are several unresolved issues related to the equatorial jets. They accelerate rapidly to speeds over lms"1 when westerly wind stress increases to about 0.7 dyne cm"2 in spring and fall, but decelerate while the wind stress continues to be westerly; each jet is followed by westward flow in the upper ocean lasting a month or longer. In addition to the semi-annual cycle, the equatorial winds and currents have strong in-traseasonal fluctuations. Observations show strong 30-60 day variability of zonal flow, and suggest that there might be variability with periods shorter than 20 days in the central EqlO. Observations from moored current meter arrays along 80.5°E south of Sri Lanka showed a distinct 15 day oscillation of equatorial meridional velocity (v) and off-equatorial zonal velocity (u). Recent observations from current meter moorings at the equator in the eastern EqlO show continuous 10-20 day, or biweekly, oscillations of v. The main motivation for the present study is to understand the dynamics of intraseasonal variability in the Indian Ocean that has been documented in the observational literature. What physical processes are responsible for the peculiar behavior of the equatorial jets? What are the relative roles of wind stress and large scale ocean dynamics? Does intraseasonal variability of wind stress force intraseasonal jets? What is the structure and origin of the biweekly variability? The intraseasonal and longer timescale variability of the equatorial Indian Ocean circulation is studied using an ocean general circulation model (OGCM) and recent in Abstract ii situ observations. The OGCM simulations are validated against other available observations. In this thesis, we document the space-time structure of the variability of equatorial Indian Ocean circulation, and attempt to find answers to some of the questions raised above. The main results are based on OGCM simulations forced by high frequency reanalysis and satellite scatterometer (QuikSCAT) winds. Several model experiments with idealized winds are used to interpret the results of the simulations. In addition to the OGCM simulations, the origin of observed intraseasonal anomalies of sea surface temperature (SST) in the eastern EqlO and Bay of Bengal, and related air-sea interaction, are investigated using validated satellite data. The main findings of the thesis can be summarized as: • High frequency accurate winds are required for accurate simulation of equatorial Indian Ocean currents, which have strong variability on intraseasonal to interannual time scales. • The variability in the equatorial waveguide is mainly driven by variability of the winds; there is some intraseasonal variability near the western boundary and in the equatorial waveguide due to dynamic instability of seasonal "mean" flows. • The fall equatorial jet is generally stronger and longer lived than the spring jet; the fall jet is modulated on intraseasonal time scales. Westerly wind bursts can drive strong intraseasonal equatorial jets in the eastern EqlO during the summer monsoon. • Eastward equatorial jets create a westward zonal pressure gradient force by raising sea level, and deepening the thermocline, in the east relative to the west. The zonal pressure force relaxes via Rossby wave radiation from the eastern boundary. • The zonal pressure force exerts strong control on the evolution of zonal flow; the decel eration of the eastward jets, and the subsequent westward flow in the upper ocean in the presence of westerly wind stress, is due to the zonal pressure force. • Neither westward currents in the upper ocean nor subsurface eastward flow (the ob served spring and summer "undercurrent") requires easterly winds; they can be gener ated by equatorial adjustment due to Kelvin (Rossby) waves generated at the western (eastern) boundary. • The biweekly variability in the EqlO is associated with forced mixed Rossby-gravity (MRG) waves generated by intraseasonal variability of winds. The biweekly MRG wave in has westward and upward phase propagation, zonal wavelength of 3000-4500 km and phase speed of 4 m s"1; it is associated with deep off equatorial upwelling/downwelling. Intraseasonal SST anomalies are forced mainly by net heat flux anomalies in the central and eastern EqlO; the large northward propagating SST anomalies in summer in the Bay of Bengal are due to net heat flux anomalies associated with the monsoon active-break cycle. Coherent variability in the atmosphere and ocean suggests air-sea interaction.

Today, ocean modeling is fast developing as a versatile tool for the study of earth’s climate, local marine ecosystems and coastal engineering applications. Though the field of ocean modeling began in the early 1950s along with the development of climate models and primitive computers, even today, the state-of-the-art ocean models have their own limitations. Many issues still remain such as the uncertainity in the parameterisation of essential processes that occur on spatial and temporal scales smaller than that can be resolved in model calculations, atmospheric forcing of the ocean and the boundary and initial conditions. The advent of data assimilation into ocean modeling has heralded a new era in the field of ocean modeling and oceanic sciences. “Data assimilation” is a methodology in which observations are used to improve the forecasting skill of operational meteorological models. The study in the present thesis mainly focuses on obtaining a four dimensional realization (the spatial description coupled with the time evolution) of the oceanic flow that is simultaneously consistent with the observational evidence and with the dynamical equations of motion and to provide initial conditions for predictions of oceanic circulation and tracer distribution. A good implementation of data assimilation can be achieved with the availability of large number of good quality observations of the oceanic fields as both synoptic and in-situ data. With the technology in satellite oceanography and insitu measurements advancing by leaps over the past two decades, good synoptic and insitu observations of oceanic fields have been achieved. The current and expected explosion in remotely sensed and insitu measured oceanographic data is ushering a new age of ocean modeling and data assimilation. The thesis presents results of analysis of the impact of data assimilation in an ocean general circulation model of the North Indian Ocean. In this thesis we have studied the impact of assimilation of temperature and salinity profiles from Argo floats and Sea Surface height anomalies from satellite altimeters in a Sigma-coordinate Indian Ocean model. An ocean data assimilation system based on the Regional Ocean Modeling System (ROMS) for the Indian Ocean is used. This model is implemented, validated and applied in a climatological simulation experiment to study the circulation in the Indian Ocean. The validated model is then used for the implementation of the data assimilation system for the Indian Ocean region. This dissertation presents the qualitative and quantitative comparisons of the model simulations with and without subsurface temperature and salinity profiles and sea surface height anamoly data assimilation for the Indian Ocean region. This is the first ever reported data assimilation studies of the Argo subsurface temperature and salinity profile data with ROMS in the Indian Ocean region.

Implementing a nonlinear gravity wave (GW) parameterization into a mechanistic middle and upper atmosphere model, which extends to the lower thermosphere (160 km), we study the response of the atmosphere in terms of the circulation patterns, temperature distribution, and migrating terdiurnal solar tide activity to the upward propagating smallscale internal GWs originating in the lower atmosphere. We perform three test simulations for the Northern Hemisphere winter conditions in order to assess the effects of variations in the initial GWspectrum on the climatology and tidal patterns of the mesosphere and lower thermosphere. We find that the overall strength of the source level momentum flux has a relatively small impact on the zonal mean climatology. The tails of the GW source level spectrum, however, are crucial for the lower thermosphere climatology. With respect to the terdiurnal tide, we find a strong dependence of tidal amplitude on the induced GW drag, generally being larger when GW drag is increased.

I present herein an ensemble of ROCKE-3D aqua-planet simulations which I compare with the simulations presented in the work of Yang et al. (2019) and other similar works. The focus was on contrasting differences in the greenhouse effect between the models. In contrast to their work, I examined simulations with a dynamic ocean instead of a slab ocean, as well as the inclusion of dynamic sea ice for 2 out of 4 of them. A subset of the simulations examined prevented the formation of sea ice to make them more comparable to Yang et al. (2019), but they never reached radiative equilibrium and this made it difficult to utilize their results. When contrasting the sea ice simulations of ROCKE-3D with the CAM4_Wolf/ExoCAM simulation of Komacek & Abbot(2019), I found that the inclusion of ocean heat transport through a dynamic ocean increases the ice-free region around the sub-tropics for the rapidly rotating aqua-planet around a G-star, thus, resulting in a lower Bond albedo and more surface warming. Supporting previous intercomparisons (Sergeevet al. 2021), ROCKE-3D produces less low- to midlevel clouds toward the equator/substellar point, compared to other models. Consequently, this leads to less cooling through the shortwave cloud radiative forcing. Lastly, I looked at the specific humidity. ROCKE-3D produced the highest stratospheric water vapor content in the M-star scenario, which suggests that ROCKE-3D is closer to the moist greenhouse limit of Kasting et al. (1993); although, the model is still far off. / I detta arbete så presenterar jag vattenplanet simulationer producerade med hjälp av ROCKE-3D, som jag sedan jämför med simulationerna som presenteras i Yang et al. (2019). Fokuset för jämförelsen låg på att jämföra skillnader gällande den producerade växthuseffekten. Alla simulationer utnyttjade ett dynamiskt hav i stället för ett enklare "platt hav", varav 2 utav 4 av dom simulationer som presenteras tillåter havsis att formas. De simulationer som inte tillät is uppnådde aldrig termisk jämvikt vilket gjorde det svårt att antyda något utifrån dom. Vid jämförelse av is-simulationerna som producerades av ROCKE-3D med de is-simulationer producerade med hjälp av CAM4/ExoCAM i Komacek & Abbot (2019), så visades det sig att den is-fria regionen runtomkring de subtropiska områderna vart större för den snabbt roterande vattenplaneten runt en G-klassad stjärna om man inkluderar ett dynamiskt hav i stället för ett "platt hav". Vidare, så stödjer detta arbete dom resultat presenterade i Sergeev et al. (2021), där ROCKE-3D producerar mindre låg- och medelhöga molnformationer runtom ekvatorn/substellära regionen, jämfört med andra modeller. Vidare, så leder detta till en mindre kylningseffekt genom molnens reflektion av kortvågsstrålning. Sist, så undersökte jag den specifika fuktigheten, där ROCKE-3D visar på ett högre värde av stratosfärisk vattenånga i fallet av en tidsvattenlåst havsplanet runt en M-stjärna. Detta tyder på att ROCKE-3D är närmare den fuktiga växthusgränsen som presenteras i Kasting et al. (1993). Dock, så är den fortfarande långt ifrån att uppnås.