Dependence of climate forcing and response on the altitude of black carbon aerosols

Black carbon aerosols absorb solar radiation and decrease planetary albedo, and thus can contribute to climate warming. In this paper, the dependence of equilibrium climate response on the altitude of black carbon is explored using an atmospheric general circulation model coupled to a mixed layer ocean model. The simulations model aerosol direct and semi-direct effects, but not indirect effects. Aerosol concentrations are prescribed and not interactive. It is shown that climate response of black carbon is highly dependent on the altitude of the aerosol. As the altitude of black carbon increases, surface temperatures decrease; black carbon near the surface causes surface warming, whereas black carbon near the tropopause and in the stratosphere causes surface cooling. This cooling occurs despite increasing planetary absorption of sunlight (i.e. decreasing planetary albedo). We find that the trend in surface air temperature response versus the altitude of black carbon is consistent with our calculations of radiative forcing after the troposphere, stratosphere, and land surface have undergone rapid adjustment, calculated as ``regressed'' radiative forcing. The variation in climate response from black carbon at different altitudes occurs largely from different fast climate responses; temperature dependent feedbacks are not statistically distinguishable. Impacts of black carbon at various altitudes on the hydrological cycle are also discussed; black carbon in the lowest atmospheric layer increases precipitation despite reductions in solar radiation reaching the surface, whereas black carbon at higher altitudes decreases precipitation.

Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Year-to-year changes of vertical temperature distribution in the California Current region: 1954 to 1986

Studies by Enfield and Allen (1980), McLain et al (1985), and others have shown that anomalously warm years in the northern coastal California Current correspond to El Niño conditions in the equatorial Pacific Ocean. Ocean model studies suggest a mechanical link between the northern coastal California Current and the equatorial ocean through long waves that propagate cyclonically along the ocean boundary (McCreary 1976; Clarke 1983; Shriver et al 1991). However, distinct observational evidence of such an oceanic connection is not extensive. Much of the supposed El Niño variation in temperature and sea level data from the coastal California Current region can be associated with the effects of anomalously intense north Pacific atmospheric cyclogenesis, which is frequently augmented during El Niño years (Wallace and Gutzler 1981; Simpson 1983; Emery and Hamilton 1984). This study uses time series of ocean temperature data to distinguish between locally forced effects, initiated by north Pacific atmospheric changes, and remotely forced effects, initiated by equatorial Pacific atmospheric changes related to El Niño events.

Dynamics of intraseasonal sea level and thermocline variability in the equatorial Atlantic during 2002-03

Satellite and in situ observations in the equatorial Atlantic Ocean during 2002-03 show dominant spectral peaks at 40-60 days and secondary peaks at 10-40 days in sea level and thermocline within the intraseasonal period band (10-80 days). A detailed investigation of the dynamics of the intraseasonal variations is carried out using an ocean general circulation model, namely, the Hybrid Coordinate Ocean Model (HYCOM). Two parallel experiments are performed in the tropical Atlantic Ocean basin for the period 2000-03: one is forced by daily scatterometer winds from the Quick Scatterometer (QuikSCAT) satellite together with other forcing fields, and the other is forced by the low-passed 80-day version of the above fields. To help in understanding the role played by the wind-driven equatorial waves, a linear continuously stratified ocean model is also used. Within 3 degrees S-3 degrees N of the equatorial region, the strong 40-60-day sea surface height anomaly (SSHA) and thermocline variability result mainly from the first and second baroclinic modes equatorial Kelvin waves that are forced by intraseasonal zonal winds, with the second baroclinic mode playing a more important role. Sharp 40-50-day peaks of zonal and meridional winds appear in both the QuikSCAT and Pilot Research Moored Array in the Tropical Atlantic (PIRATA) data for the period 2002-03, and they are especially strong in 2002. Zonal wind anomaly in the central-western equatorial basin for the period 2000-06 is significantly correlated with SSHA across the equatorial basin, with simultaneous/ lag correlation ranging from-0.62 to 0.74 above 95% significance. Away from the equator (3 degrees-5 degrees N), however, sea level and thermocline variations in the 40-60-day band are caused largely by tropical instability waves (TIWs). On 10-40-day time scales and west of 10 degrees W, the spectral power of sea level and thermocline appears to be dominated by TIWs within 5 degrees S-5 degrees N of the equatorial region. The wind-driven circulation, however, also provides a significant contribution. Interestingly, east of 10 W, SSHA and thermocline variations at 10 40- day periods result almost entirely from wind-driven equatorial waves. During the boreal spring of 2002 when TIWs are weak, Kelvin waves dominate the SSHA across the equatorial basin (2 degrees S-2 degrees N). The observed quasi-biweekly Yanai waves are excited mainly by the quasi-biweekly meridional winds, and they contribute significantly to the SSHA and thermocline variations in 1 degrees-5 degrees N and 1 degrees-5 degrees S regions.

Water transports through the four main straits around the South China Sea

A quasi-global high-resolution HYbrid Coordinate Ocean Model (HYCOM) is used to investigate seasonal variations of water transports through the four main straits in the South China Sea. The results show that the annual transports through the four straits Luzon Strait, Taiwan Strait, Sunda Shelf and Mindoro Strait are -4.5, 2.3, 0.5 and 1.7 Sv (1 Sv=10(6) m(3)s(-1)), respectively. The Mindoro Strait has an important outflow that accounts for over one third of the total inflow through the Luzon Strait. Furthermore, it indicates that there are strong seasonal variations of water transport in the four straits. The water transport through the Luzon Strait (Taiwan Strait, Sunda Shelf, Mindoro Strait) has a maximum value of -7.6 Sv in December (3.1 Sv in July, 2.1S v in January, 4.5Sv in November), a minimum value of -2.1 Sv in June (1.5 Sv in October, -1.0 Sv in June, -0.2 Sv in May), respectively.

Factors influencing the climatological mixed layer depth in the South China Sea: numerical simulations

The mixed layer depth (MLD) in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth (CMLD), and established a numerical simulation in the South China Sea (SCS) using the Regional Ocean Model System (ROMS) with a high-resolution (1/12A degrees x1/12A degrees) grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress (SSWS), sea surface net heat flux (SSNHF), and the difference between evaporation and precipitation (DEP). The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF (13%) for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third (only 2%). Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

Ecohydrodynamics of cold-water coral reefs: A case study of the Mingulay Reef complex (western Scotland)

Ecohydrodynamics investigates the hydrodynamic constraints on ecosystems across different temporal and spatial scales. Ecohydrodynamics play a pivotal role in the structure and functioning of marine ecosystems, however the lack of integrated complex flow models for deep-water ecosystems beyond the coastal zone prevents further synthesis in these settings. We present a hydrodynamic model for one of Earth's most biologically diverse deep-water ecosystems, cold-water coral reefs. The Mingulay Reef Complex (western Scotland) is an inshore seascape of cold-water coral reefs formed by the scleractinian coral Lophelia pertusa. We applied single-image edge detection and composite front maps using satellite remote sensing, to detect oceanographic fronts and peaks of chlorophyll a values that likely affect food supply to corals and other suspension-feeding fauna. We also present a high resolution 3D ocean model to incorporate salient aspects of the regional and local oceanography. Model validation using in situ current speed, direction and sea elevation data confirmed the model's realistic representation of spatial and temporal aspects of circulation at the reef complex including a tidally driven current regime, eddies, and downwelling phenomena. This novel combination of 3D hydrodynamic modelling and remote sensing in deep-water ecosystems improves our understanding of the temporal and spatial scales of ecological processes occurring in marine systems. The modelled information has been integrated into a 3D GIS, providing a user interface for visualization and interrogation of results that allows wider ecological application of the model and that can provide valuable input for marine biodiversity and conservation applications.

Southern Annular Mode and westerly-wind-driven changes in Indian-Atlantic exchange mechanisms

The dynamical link between the Indian Ocean and Atlantic Meridional Overturning Circulation (AMOC) remains poorly understood. This partly arises from the complex Agulhas leakage, which occurs via rings, cyclones, and non-eddy flux. Hindcast simulations suggest that leakage has recently increased but have not decomposed this signal into its constituent mechanisms. Here these are isolated in a realistic ocean model. Increases in simulated leakage are attributed to stronger eddy and non-eddy-driven transports, and a strong warming and salinification, especially within Agulhas rings. Variability in both regimes is associated with strengthening Indian Ocean westerly winds, reflecting an increasingly positive Southern Annular Mode. While eddy and non-eddy flux signals are tied through turbulent eddy dissipation, the ratio between the two varies decadally. Consequently, while altimetry suggests a recent increase in retroflection turbulence and implied leakage, non-eddy flux may also play a significant role in modulating the leakage AMOC connection.

Quantifying pursuit diving seabirds' associations with fine-scale physical features in a high tidal energy environment

The rapid increase in the number of tidal stream turbine arrays will create novel and unprecedented levels of anthropogenic activity within habitats characterized by horizontal current speeds exceeding 2 ms−1. However, the potential impacts on pursuit‐diving seabirds exploiting these tidal stream environments remain largely unknown. Identifying similarities between the fine‐scale physical features (100s of metres) suitable for array installations, and those associated with foraging pursuit‐diving seabirds, could identify which species are most vulnerable to either collisions with moving components, or displacement from these installations. A combination of vessel‐based observational surveys, Finite Volume Community Ocean Model outputs and hydroacoustic seabed surveys provided concurrent measures of foraging distributions and physical characteristics at a fine temporal (15 min) and spatial (500 m) resolution across a tidal stream environment suitable for array installations, during both breeding and non‐breeding seasons. These data sets were then used to test for associations between foraging pursuit‐diving seabirds (Atlantic puffins Fratercula arctica, black guillemots Cepphus grylle, common guillemots Uria aalge, European shags Phalacrocorax aristotelis) and physical features. These species were associated with areas of fast horizontal currents, slow horizontal currents, high turbulence, downward vertical currents and also hard–rough seabeds. The identity and strength of associations differed among species, and also within species between seasons, indicative of interspecific and intraspecific variations in habitat use. However, Atlantic puffins were associated particularly strongly with areas of fast horizontal currents during breeding seasons, and European shags with areas of rough–hard seabeds and downward vertical currents during non‐breeding seasons. Synthesis and applications. Atlantic puffins’ strong association with fast horizontal current speeds indicates that they are particularly likely to interact with installations during breeding seasons. Any post‐installation monitoring and mitigation measures should therefore focus on this species and season. The multi‐species associations with high turbulence and downward vertical currents, which often coincide with fast horizontal current speeds, also highlight useful pre‐installation mitigation measures via the omission of devices from these areas, reducing the overall likelihood of interactions. Environmental impact assessments (EIA) generally involve once‐a‐month surveys across 2‐year periods. However, the approaches used in this study show that more focussed surveys can greatly benefit management strategies aiming to reduce the likelihood of negative impacts by facilitating the development of targeted mitigation measures. It is therefore recommended that these approaches contribute towards EIA within development sites.

Modelling the cycling of persistent organic pollutants (POPs) in the North Sea system: Fluxes, loading, seasonality, trends

The fate and cycling of two selected POPs is investigated for the North Sea system with an improved version of a fate and transport ocean model (FANTOM). The model uses atmospheric data from the EMEP MSC East POP model (Gusev et al., 2009), giving reasonable concentrations and seasonal distributions for the entire region, as opposed to the three observation stations that Ilyina et al. (2006) were limited to. Other model improvements include changes in the calculation of POP exchange between the water column and sediment.

We chose to simulate the fate of two POPs with very different properties, ?-HCH and PCB 153. Since the fate and cycling of POPs are strongly affected by hydrodynamic processes, a high resolution version of the Hamburg Shelf Ocean Model (HAMSOM) was developed and utilised. Simulations were made for the period 1996–2005. Both models were validated by comparing results with available data, which showed that the simulations were of very satisfactory quality.

Model results show that the North Sea is a net sink for ?-HCH and a net source to the atmosphere of PCB 153. Total masses of ?-HCH and PCB 153 in 2005 are reduced to 30% and 50%, respectively, of 1996 values.

Storms resuspending bottom sediments into the water column mobilise POPs into the atmosphere and have the potential to deliver substantial loads of these POPs into Europe.

The effects of global climate change on the cycling and processes of persistent organic pollutants (POPs) in the North Sea

The fate and cycling of two selected legacy persistent organic pollutants (POPs), PCB 153 and gamma-HCH, in the North Sea in the 21st century have been modelled with combined hydrodynamic and fate and transport ocean models
(HAMSOM and FANTOM, respectively). To investigate the impact of climate variability on POPs in the North Sea in the 21st century, future scenario model runs for three 10-year periods to the year 2100 using plausible levels of both in
situ concentrations and atmospheric, river and open boundary inputs are performed. This slice mode under a moderate scenario (A1B) is sufficient to provide a basis for further analysis. For the HAMSOM and atmospheric forcing, results of the IPCC A1B (SRES) 21st century scenario are utilized, where surface forcing is provided by the REMO downscaling of the ECHAM5 global atmospheric model, and open boundary conditions are provided by the MPIOM global ocean model.
Dry gas deposition and volatilization of gamma-HCH increase in the future relative to the present by up to 20% (in the spring and summer months for deposition and in summer for volatilization). In the water column, total mass of
gamma-HCH and PCB 153 remain fairly steady in all three runs. In sediment,
gamma-HCH increases in the future runs, relative to the present, while PCB 153 in sediment decreases exponentially in all three runs, but even faster in the future, due to the increased number of storms, increased duration of gale wind conditions and increased water and air temperatures, all of which are the result of climate change. Annual net sinks exceed sources at the ends of all periods.
Overall, the model results indicate that the climate change scenarios considered here generally have a negligible influence on the simulated fate and transport of the two POPs in the North Sea, although the increased number and magnitude of storms in the 21st century will result in POP resuspension and ensuing revolatilization events. Trends in emissions from primary and secondary sources will remain the key driver of levels of these contaminants over time.

Modelling the fate of persistent organic pollutants (POPs) in the North Sea system

The environmental fate of selected persistent organic pollutants (POPs) in the North Sea system is modelled with a high resolution Fate and Transport Ocean Model (FANTOM) that uses hydrodynamic model output from the Hamburg Shelf Ocean Model (HAMSOM). Large amounts of POPs enter the North Sea from the surrounding highly populated, industrialised and agricultural countries. Major pathways to the North Sea are atmospheric deposition and river inputs, with additional contributions coming from bottom sediments and adjacent seas. The model domain covers the entire North Sea region, extending northward as far as the Shetland Islands, and includes adjacent basins such as the Skagerrak, Kattegat, and the westernmost part of the Baltic Sea. Model resolution (for both models) is 1.5’ latitude x 2.5’ longitude (approximately 3 km horizontal resolution) with 30 vertical levels. The POP model also has 20 sediment layers. Important model processes controlling the fate of POPs in the North Sea system are discussed. Results focus on Lindane gamma- HCH or gamma-hexachlorocyclohexane) and PCB 153.

Modeling the Western Iberian margin circulation: present and future

The purpose of this work is to carry out a comprehensive study on the Western Iberian Margin (WIM) circulation my means of numerical modeling, and to postulate what this circulation will be in the future. The adopted approach was the development of a regional ocean model configuration with high resolution, capable of reproducing the largeand small-scale dynamics of the coastal transition zone. Four numerical experiences were carried out according to these objectives: (1) a climatological run, in order to study the system’s seasonal behavior and its mean state; (2) a run forced with real winds and fluxes for period 2001-2011 in order to study the interannual variability of the system; (3) a run forced with mean fields from Global Climate Models (GCMs) for the present, in order to validate GCMs as adequate forcing for regional ocean modeling; (4) a similar run (3) for period 2071-2100, in order to assess possible consequences of a future climate scenario on the hydrography and dynamics of the WIM. Furthermore, two Lagrangian particle studies were carried out: one in order to trace the origin of the upwelled waters along the WIM; the other in order to portrait the patterns of larval dispersal, accumulation and connectivity. The numerical configuration proved to be adequate in the reproduction of the system’s mean state, seasonal characterization and an interannual variability study. There is prevalence of poleward flow at the slope, which coexists with the upwelling jet during summer, although there is evidence of its shifting offshore, and which is associated with the Mediterranean Water flow at deeper levels, suggesting a barotropic character. From the future climate scenario essay, the following conclusions were drawn: there is general warming and freshening of upper level waters; there is still poleward tendency, and despite the upwellingfavorable winds strengthening in summer the respective coastal band becomes more restricted in width and depth. In what concerns larval connectivity and dispersion along the WIM, diel vertical migration was observed to increase recruitment throughout the domain, and while smooth coastlines are better suppliers, there is higher accumulation where the topography is rougher.

Acoustic-oceanographic data fusion for prediction of oceanic acoustic elds

Tese de doutoramento, Engenharia Electrónica e Telecomunicações (Processamento de Sinal), Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2014

Running Climate Models On The NERC Cluster Grid Using G-Rex

Compute grids are used widely in many areas of environmental science, but there has been limited uptake of grid computing by the climate modelling community, partly because the characteristics of many climate models make them difficult to use with popular grid middleware systems. In particular, climate models usually produce large volumes of output data, and running them usually involves complicated workflows implemented as shell scripts. For example, NEMO (Smith et al. 2008) is a state-of-the-art ocean model that is used currently for operational ocean forecasting in France, and will soon be used in the UK for both ocean forecasting and climate modelling. On a typical modern cluster, a particular one year global ocean simulation at 1-degree resolution takes about three hours when running on 40 processors, and produces roughly 20 GB of output as 50000 separate files. 50-year simulations are common, during which the model is resubmitted as a new job after each year. Running NEMO relies on a set of complicated shell scripts and command utilities for data pre-processing and post-processing prior to job resubmission. Grid Remote Execution (G-Rex) is a pure Java grid middleware system that allows scientific applications to be deployed as Web services on remote computer systems, and then launched and controlled as if they are running on the user's own computer. Although G-Rex is general purpose middleware it has two key features that make it particularly suitable for remote execution of climate models: (1) Output from the model is transferred back to the user while the run is in progress to prevent it from accumulating on the remote system and to allow the user to monitor the model; (2) The client component is a command-line program that can easily be incorporated into existing model work-flow scripts. G-Rex has a REST (Fielding, 2000) architectural style, which allows client programs to be very simple and lightweight and allows users to interact with model runs using only a basic HTTP client (such as a Web browser or the curl utility) if they wish. This design also allows for new client interfaces to be developed in other programming languages with relatively little effort. The G-Rex server is a standard Web application that runs inside a servlet container such as Apache Tomcat and is therefore easy to install and maintain by system administrators. G-Rex is employed as the middleware for the NERC1 Cluster Grid, a small grid of HPC2 clusters belonging to collaborating NERC research institutes. Currently the NEMO (Smith et al. 2008) and POLCOMS (Holt et al, 2008) ocean models are installed, and there are plans to install the Hadley Centre’s HadCM3 model for use in the decadal climate prediction project GCEP (Haines et al., 2008). The science projects involving NEMO on the Grid have a particular focus on data assimilation (Smith et al. 2008), a technique that involves constraining model simulations with observations. The POLCOMS model will play an important part in the GCOMS project (Holt et al, 2008), which aims to simulate the world’s coastal oceans. A typical use of G-Rex by a scientist to run a climate model on the NERC Cluster Grid proceeds as follows :(1) The scientist prepares input files on his or her local machine. (2) Using information provided by the Grid’s Ganglia3 monitoring system, the scientist selects an appropriate compute resource. (3) The scientist runs the relevant workflow script on his or her local machine. This is unmodified except that calls to run the model (e.g. with “mpirun”) are simply replaced with calls to "GRexRun" (4) The G-Rex middleware automatically handles the uploading of input files to the remote resource, and the downloading of output files back to the user, including their deletion from the remote system, during the run. (5) The scientist monitors the output files, using familiar analysis and visualization tools on his or her own local machine. G-Rex is well suited to climate modelling because it addresses many of the middleware usability issues that have led to limited uptake of grid computing by climate scientists. It is a lightweight, low-impact and easy-to-install solution that is currently designed for use in relatively small grids such as the NERC Cluster Grid. A current topic of research is the use of G-Rex as an easy-to-use front-end to larger-scale Grid resources such as the UK National Grid service.