Observations of the solar wind-magnetosphere-ionosphere coupling

In this thesis, the solar wind-magnetosphere-ionosphere coupling is studied observationally, with the main focus on the ionospheric currents in the auroral region. The thesis consists of five research articles and an introductory part that summarises the most important results reached in the articles and places them in a wider context within the field of space physics. Ionospheric measurements are provided by the International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network, by the low-orbit CHAllenging Minisatellite Payload (CHAMP) satellite, by the European Incoherent SCATter (EISCAT) radar, and by the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) satellite. Magnetospheric observations, on the other hand, are acquired from the four spacecraft of the Cluster mission, and solar wind observations from the Advanced Composition Explorer (ACE) and Wind spacecraft. Within the framework of this study, a new method for determining the ionospheric currents from low-orbit satellite-based magnetic field data is developed. In contrast to previous techniques, all three current density components can be determined on a matching spatial scale, and the validity of the necessary one-dimensionality approximation, and thus, the quality of the results, can be estimated directly from the data. The new method is applied to derive an empirical model for estimating the Hall-to-Pedersen conductance ratio from ground-based magnetic field data, and to investigate the statistical dependence of the large-scale ionospheric currents on solar wind and geomagnetic parameters. Equations describing the amount of field-aligned current in the auroral region, as well as the location of the auroral electrojets, as a function of these parameters are derived. Moreover, the mesoscale (10-1000 km) ionospheric equivalent currents related to two magnetotail plasma sheet phenomena, bursty bulk flows and flux ropes, are studied. Based on the analysis of 22 events, the typical equivalent current pattern related to bursty bulk flows is established. For the flux ropes, on the other hand, only two conjugate events are found. As the equivalent current patterns during these two events are not similar, it is suggested that the ionospheric signatures of a flux rope depend on the orientation and the length of the structure, but analysis of additional events is required to determine the possible ionospheric connection of flux ropes.

Influence of Spatial Scales on the Performance of Saturated Subsurface Flow and Transport Models

Numerical modeling of saturated subsurface flow and transport has been widely used in the past using different numerical schemes such as finite difference and finite element methods. Such modeling often involves discretization of the problem in spatial and temporal scales. The choice of the spatial and temporal scales for a modeling scenario is often not straightforward. For example, a basin-scale saturated flow and transport analysis demands larger spatial and temporal scales than a meso-scale study, which in turn has larger scales compared to a pore-scale study. The choice of spatial-scale is often dictated by the computational capabilities of the modeler as well as the availability of fine-scale data. In this study, we analyze the impact of different spatial scales and scaling procedures on saturated subsurface flow and transport simulations.

Processes of 30-90 days sea surface temperature variability in the northern Indian Ocean during boreal summer

During summer, the northern Indian Ocean exhibits significant atmospheric intraseasonal variability associated with active and break phases of the monsoon in the 30-90 days band. In this paper, we investigate mechanisms of the Sea Surface Temperature (SST) signature of this atmospheric variability, using a combination of observational datasets and Ocean General Circulation Model sensitivity experiments. In addition to the previously-reported intraseasonal SST signature in the Bay of Bengal, observations show clear SST signals in the Arabian Sea related to the active/break cycle of the monsoon. As the atmospheric intraseasonal oscillation moves northward, SST variations appear first at the southern tip of India (day 0), then in the Somali upwelling region (day 10), northern Bay of Bengal (day 19) and finally in the Oman upwelling region (day 23). The Bay of Bengal and Oman signals are most clearly associated with the monsoon active/break index, whereas the relationship with signals near Somali upwelling and the southern tip of India is weaker. In agreement with previous studies, we find that heat flux variations drive most of the intraseasonal SST variability in the Bay of Bengal, both in our model (regression coefficient, 0.9, against similar to 0.25 for wind stress) and in observations (0.8 regression coefficient); similar to 60% of the heat flux variation is due do shortwave radiation and similar to 40% due to latent heat flux. On the other hand, both observations and model results indicate a prominent role of dynamical oceanic processes in the Arabian Sea. Wind-stress variations force about 70-100% of SST intraseasonal variations in the Arabian Sea, through modulation of oceanic processes (entrainment, mixing, Ekman pumping, lateral advection). Our similar to 100 km resolution model suggests that internal oceanic variability (i.e. eddies) contributes substantially to intraseasonal variability at small-scale in the Somali upwelling region, but does not contribute to large-scale intraseasonal SST variability due to its small spatial scale and random phase relation to the active-break monsoon cycle. The effect of oceanic eddies; however, remains to be explored at a higher spatial resolution.

Space-Time Scales of Northward Propagation of Convection during Boreal Summer

This study uses precipitation estimates from the Tropical Rainfall Measuring Mission to quantify the spatial and temporal scales of northward propagation of convection over the Indian monsoon region during boreal summer. Propagating modes of convective systems in the intraseasonal time scales such as the Madden-Julian oscillation can interact with the intertropical convergence zone and bring active and break spells of the Indian summer monsoon. Wavelet analysis was used to quantify the spatial extent (scale) and center of these propagating convective bands, as well as the time period associated with different spatial scales. Results presented here suggest that during a good monsoon year the spatial scale of this oscillation is about 30 degrees centered around 10 degrees N. During weak monsoon years, the scale of propagation decreases and the center shifts farther south closer to the equator. A strong linear relationship is obtained between the center/scale of convective wave bands and intensity of monsoon precipitation over Indian land on the interannual time scale. Moreover, the spatial scale and its center during the break monsoon were found to be similar to an overall weak monsoon year. Based on this analysis, a new index is proposed to quantify the spatial scales associated with propagating convective bands. This automated wavelet-based technique developed here can be used to study meridional propagation of convection in a large volume of datasets from observations and model simulations. The information so obtained can be related to the interannual and intraseasonal variation of Indian monsoon precipitation.

An exploration framework to identify and track movement of cloud systems

We describe a framework to explore and visualize the movement of cloud systems. Using techniques from computational topology and computer vision, our framework allows the user to study this movement at various scales in space and time. Such movements could have large temporal and spatial scales such as the Madden Julian Oscillation (MJO), which has a spatial scale ranging from 1000 km to 10000 km and time of oscillation of around 40 days. Embedded within these larger scale oscillations are a hierarchy of cloud clusters which could have smaller spatial and temporal scales such as the Nakazawa cloud clusters. These smaller cloud clusters, while being part of the equatorial MJO, sometimes move at speeds different from the larger scale and in a direction opposite to that of the MJO envelope. Hitherto, one could only speculate about such movements by selectively analysing data and a priori knowledge of such systems. Our framework automatically delineates such cloud clusters and does not depend on the prior experience of the user to define cloud clusters. Analysis using our framework also shows that most tropical systems such as cyclones also contain multi-scale interactions between clouds and cloud systems. We show the effectiveness of our framework to track organized cloud system during one such rainfall event which happened at Mumbai, India in July 2005 and for cyclone Aila which occurred in Bay of Bengal during May 2009.

In a hot bubble: why does superbubble feedback work, but isolated supernovae do not?

Using idealized one-dimensional Eulerian hydrodynamic simulations, we contrast the behaviour of isolated supernovae with the superbubbles driven by multiple, collocated supernovae. Continuous energy injection via successive supernovae exploding within the hot/dilute bubble maintains a strong termination shock. This strong shock keeps the superbubble over-pressured and drives the outer shock well after it becomes radiative. Isolated supernovae, in contrast, with no further energy injection, become radiative quite early (less than or similar to 0.1Myr, tens of pc), and stall at scales less than or similar to 100 pc. We show that isolated supernovae lose almost all of their mechanical energy by 1 Myr, but superbubbles can retain up to similar to 40 per cent of the input energy in the form of mechanical energy over the lifetime of the star cluster (a few tens of Myr). These conclusions hold even in the presence of realistic magnetic fields and thermal conduction. We also compare various methods for implementing supernova feedback in numerical simulations. For various feedback prescriptions, we derive the spatial scale below which the energy needs to be deposited in order for it to couple to the interstellar medium. We show that a steady thermal wind within the superbubble appears only for a large number (greater than or similar to 10(4)) of supernovae. For smaller clusters, we expect multiple internal shocks instead of a smooth, dense thermalized wind.

Critical Sensitivity and Trans-Scale Fluctuations in Catastrophic Rupture

Rupture in the heterogeneous crust appears to be a catastrophe transition. Catastrophic rupture sensitively depends on the details of heterogeneity and stress transfer on multiple scales. These are difficult to identify and deal with. As a result, the threshold of earthquake-like rupture presents uncertainty. This may be the root of the difficulty of earthquake prediction. Based on a coupled pattern mapping model, we represent critical sensitivity and trans-scale fluctuations associated with catastrophic rupture. Critical sensitivity means that a system may become significantly sensitive near catastrophe transition. Trans-scale fluctuations mean that the level of stress fluctuations increases strongly and the spatial scale of stress and damage fluctuations evolves from the mesoscopic heterogeneity scale to the macroscopic scale as the catastrophe regime is approached. The underlying mechanism behind critical sensitivity and trans-scale fluctuations is the coupling effect between heterogeneity and dynamical nonlinearity. Such features may provide clues for prediction of catastrophic rupture, like material failure and great earthquakes. Critical sensitivity may be the physical mechanism underlying a promising earthquake forecasting method, the load-unload response ratio (LURR).

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

Scale-Dependent Habitat Selection of Nesting Great Egrets and Snowy Egrets

Foraging habitat selection of nesting Great Egrets ( Ardea alba ) and Snowy Egrets ( Egretta thula ) was investigated within an estuary with extensive impounded salt marsh habitat. Using a geographic information system, available habitat was partitioned into concentric bands at five, ten, and 15 km radius from nesting colonies to assess the relative effects of habitat composition and distance on habitat selection. Snowy Egrets were more likely than Great Egrets to depart colonies and travel to foraging sites in groups, but both species usually arrived at sites that were occupied by other wading birds. Mean flight distances were 6.2 km (SE = 0.4, N = 28, range 1.8-10.7 km) for Great Egrets and 4.7 km (SE = 0.48, N = 31, range 0.7-12.5 km) for Snowy Egrets. At the broadest spatial scale both species used impounded (mostly salt marsh) and estuarine edge habitat more than expected based on availability while avoiding unimpounded (mostly fresh water wetland) habitat. At more local scales habitat use matched availability. Interpretation of habitat preference differed with the types of habitat that were included and the maximum distance that habitat was considered available. These results illustrate that caution is needed when interpreting the results of habitat preference studies when individuals are constrained in their choice of habitats, such as for central place foragers.

Multiscale analysis of factors that affect the distribution of sharks throughout the northern Gulf of Mexico

Identification of the spatial scale at which marine communities are organized is critical to proper management, yet this is particularly difficult to determine for highly migratory species like sharks. We used shark catch data collected during 2006–09 from fishery-independent bottom-longline surveys, as well as biotic and abiotic explanatory data to identify the factors that affect the distribution of coastal sharks at 2 spatial scales in the northern Gulf of Mexico. Centered principal component analyses (PCAs) were used to visualize the patterns that characterize shark distributions at small (Alabama and Mississippi coast) and large (northern Gulf of Mexico) spatial scales. Environmental data on temperature, salinity, dissolved oxygen (DO), depth, fish and crustacean biomass, and chlorophyll-a (chl-a) concentration were analyzed with normed PCAs at both spatial scales. The relationships between values of shark catch per unit of effort (CPUE) and environmental factors were then analyzed at each scale with co-inertia analysis (COIA). Results from COIA indicated that the degree of agreement between the structure of the environmental and shark data sets was relatively higher at the small spatial scale than at the large one. CPUE of Blacktip Shark (Carcharhinus limbatus) was related positively with crustacean biomass at both spatial scales. Similarly, CPUE of Atlantic Sharpnose Shark (Rhizoprionodon terraenovae) was related positively with chl-a concentration and negatively with DO at both spatial scales. Conversely, distribution of Blacknose Shark (C. acronotus) displayed a contrasting relationship with depth at the 2 scales considered. Our results indicate that the factors influencing the distribution of sharks in the northern Gulf of Mexico are species specific but generally transcend the spatial boundaries used in our analyses.

Integrating remote sensing products and GIS tools to support marine spatial management in West Hawai'i

Marine protected areas (MPAs) represent a form of spatial management, and geospatial information on living marine resources and associated habitat is extremely important to support best management practices in a spatially discrete MPA. Benthic habitat maps provide georeferenced information on the geomorphic structure and biological cover types in the marine environment. This information supports an enhanced understanding of ecosystem function and species habitat utilization patterns. Benthic habitat maps are most useful for marine management and spatial planning purposes when they are created at a scale that is relevant to management actions. We sought to improve the resolution of existing benthic habitat maps created during a regional mapping effort in Hawai`i. Our results complemented these existing regional maps and provided more detailed, finer-scale habitat maps for a network of MPAs in West Hawai`i. The map products created during this study allow local planners and managers to extract information at a spatial scale relevant to the discrete management units, and appropriate for local marine management efforts on the Kona Coast. The resultant benthic habitat maps were integrated in a geographic information system (GIS) that also included aerial imagery, underwater video, MPA regulations, summarized ecological data and other relevant and spatially explicit information. The integration of the benthic habitat maps with additional “value added” geospatial information into a dynamic GIS provide a decision support tool with pertinent marine resource information available in one central location and support the application of a spatial approach to the management of marine resources. Further, this work can serve as a case study to demonstrate the integration of remote sensing products and GIS tools at a fine spatial scale relevant to local-level marine spatial planning and management efforts.

Interannual variability in the PFEG coastal upwelling indices

Numerous studies examine decadal-scale variability in basin-scale parameters in the Northern Pacific. Characterizing such interannual-to-interdecadal variability is essential to identifying long-term climate changes. The Pacific Fisheries Environmental Group (PFEG) coastal upwelling indices display variability on these time scales and may help explain the mechanisms responsible for such climate variability. ... In this study, examination of 49-year time series of monthly mean upwelling indices at the 15 PFEG-standard positions along the west coast of North America revealed variability on large spatial scales as well as temporal scales.

Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

中国沙棘(Hippophae rhamnoides subsp. sinensis)天然群体的遗传变异分析

近十年，植物群体遗传学的研究飞速发展，然而与海拔相关的植物群体遗传结构和遗传变异研究却相对较少。到目前为止，还不清楚遗传变异与海拔之间是否有一个通用的格局。在山区，各种生态因子，如温度、降水、降雪、紫外线辐射强度以及土壤成分都随海拔梯度急剧变化，造成了即使在一个小的空间区域，植被类型变化显著，这种高山环境的异质性和复杂性为我们研究植物群体遗传结构和分化提供了方便。沙棘(Hippophea)属于胡颓子科(Elaeagnaceae)为多年生落叶灌木或乔木，雌雄异株，天然种群分布极为广泛。中国沙棘(H. rhamnoides subsp. sinensis)是沙棘属植物中分布较广的一个亚种，种内形态变异非常丰富，加之其具有独特的繁育系统和广泛的生态地理分布，是研究沙棘属植物遗传变异和系统分化的理想材料。本文从1,800 m 到3,400 m 分5 个海拔梯度进行取样，用RAPD 和cpSSR 分子标记研究了卧龙自然保护区中国沙棘天然群体的遗传结构和遗传变异。5 个取样群体依次标记为A、B、C、D 和E，它们分别代表分布在海拔1,800，2,200，2,600，3,000 和3,400 m 的5 个天然群体。RAPD实验用11 条寡核苷酸引物，扩增得到151 个重复性好的位点，其中143 个多态位点，多态率达94.7%。在5 个沙棘群体中，总遗传多样性值(HT)为0.289，B群体内的遗传多样性值为0.315，这完全符合沙棘这种多年生、远交的木本植物具有高遗传变异的特性。5 个群体内遗传多样性随海拔升高呈低－高－低变异趋势，在2,200 m海拔处的B群体遗传多样性达最大值0.315，3,400 m海拔处的E群体则表现最小仅0.098。5 个群体间的遗传分化值GST＝0.406，也即是说有40.6%的遗传变异存在于群体间，1,800 m海拔处的A群体与其它群体的明显分离是造成群体间遗传分化大的原因。UPGMA聚类图和PCoA散点图进一步确证了5 个群体间的关系和所有个体间的关系。最后，经过Mantel检测，遗传距离与海拔表现了明显的相关性(r = 0.646, P = 0.011)。cpSSR 实验中，经过对24 对cpSSR 通用引物筛选，11 对引物能扩增出特异性条带，只有2 对引物(ccmp2 和ARCP4)呈现多态性。4 个等位基因共组合出4 种单倍型，单倍型Ⅰ出现在A 群体的所有个体和B 群体的8 个个体中，C、D、E 三个群体均不含有，而单倍型Ⅱ出现在C、D、E 三个群体的所有个体及B 群体的18 个个体中，A 群体不含有。另外两种单倍型Ⅲ和Ⅳ为稀有类型，仅B 群体中的4 个个体拥有。这种单倍型分布模式和TFPGA 群体聚类图揭示了，C、D、E 群体可能来源于同一祖先种，而A 群体却是由另一祖先种发展起来的，B 群体则兼具了这两种起源种的信息，这可能是因为在历史上的某一时期，在中国沙棘群体高山分化的过程中，B 群体处某个或者某些个体发生了基因突变，具备了适应高海拔环境的能力，产生了高海拔沙棘群体的祖先种。 In recent ten years, studies about population genetics of plants developed rapidly,whereas their genetic structure and genetic variation along altitudinal gradients have beenstudied relatively little. So far, it is uncleared whether there is a common pattern betweengenetic variation and altitudinal gradients. In the mountain environments, importantecological factors, e.g., temperature, rainfall, snowfall, ultraviolet radiation and soil substratesetc., change rapidly with altitudes, which cause the vegetation distribution varying typically,even on a small spatial scale. The mountain environments, which are heterogeneous andcomplex, facilitate and offer a good opportunity to characterize population genetic structureand population differentiation.The species of the genus Hippophae L. (Elaeagnaceae) are perennial deciduous shrubs ortrees, which are dioecious, wind-pollinated pioneer plants. The natural genus has a widedistribution extending from Northern Europe through Central Europe and Central Asia toChina. According to the latest taxonomy, the genus Hippophae is divided into six species and12 subspecies. The subspecies H. rhamnoides ssp. sinensis shows significant morphologicalvariations, large geographic range and dominantly outcrossing mating system. Thesecharacteristics of the subspecies are favourable to elucidate genetic variation and systemevolution. To estimate genetic variation and genetic structure of H. rhamnoides ssp. sinensisat different altitudes, we surveyed five natural populations in the Wolong Natural Reserve at altitudes ranging from 1,800 to 3,400 m above sea level (a.s.l.) using random amplifiedpolymorphic DNA markers (RAPDs) and cpSSR molecular methods. The five populations A,B, C, D, and E correspond to the altitudes 1,800, 2,200, 2,600, 3,000 and 3,400 m,respectively.Based on 11 decamer primers, a total of 151 reproducible DNA loci were yielded, ofwhich 143 were polymorphic and the percentage of polymorphic loci equaled 94.7%. Amongthe five populations investigated, the total gene diversity (HT) and gene diversity within population B equaled 0.289 and 0.315, respectively, which are modest for a subspecies of H.rhamnoides, which is an outcrossing, long-lived, woody plant. The amount of geneticvariation within populations varied from 0.098 within population E (3,400 m a.s.l.) to 0.315within population B (2,200 m a.s.l.). The coefficient of gene differentiation (GST) amongpopulations equaled 0.406 and revealed that 40.6% of the genetic variance existed amongpopulations and 59.4% within populations. The population A (1,800 m a.s.l.) differed greatlyfrom the other four populations, which contributes to high genetic differentiation. A UPGMAcluster analysis and principal coordinate analyses based on Nei's genetic distances furthercorroborated the relationships among the five populations and all the sampling individuals,respectively. Mantel tests detected a significant correlation between genetic distances andaltitudinal gradients (r = 0.646, P = 0.011).Eleven of the original 24 cpSSR primer pairs tested produced good PCR products, onlytwo (ccmp2 and ARCP4) of which were polymorphic. Four total length variants (alleles) werecombined resulting in 4 haplotypes. The haplotype was present in all individuals of Ⅰpopulation A and 8 individuals of populations B, the other three populations (C, D and Epopulations) did not share. The haplotype was present in all individuals of populations C, D Ⅱand E and 18 individuals of populations B, population A did not share. The other twohaplotypes and were rare haplotypes, which were only shared in 4 individuals of Ⅲ Ⅳpopulation B. The distribution of haplotypes and TFPGA population clustering map showedthat the populations C, D and E might be origined from one ancestor seed and population Amight be from another, whereas population B owned information of the two ancestor seeds. Itwas because that gene mutation within some individual or seed in the location of population Bwas likely to happen in the history of H. rhamnoides, which was the original ancestor of thehigh-altitude populations.

Land-cover classification of China: integrated analysis of AVHRR imagery and geophysical data

Over last two decades, numerous studies have used remotely sensed data from the Advanced Very High Resolution Radiometer (AVHRR) sensors to map land use and land cover at large spatial scales, but achieved only limited success. In this paper, we employed an approach that combines both AVHRR images and geophysical datasets (e.g. climate, elevation). Three geophysical datasets are used in this study: annual mean temperature, annual precipitation, and elevation. We first divide China into nine bio-climatic regions, using the long-term mean climate data. For each of nine regions, the three geophysical data layers are stacked together with AVHRR data and AVHRR-derived vegetation index (Normalized Difference Vegetation Index) data, and the resultant multi-source datasets were then analysed to generate land-cover maps for individual regions, using supervised classification algorithms. The nine land-cover maps for individual regions were assembled together for China. The existing land-cover dataset derived from Landsat Thematic Mapper (TM) images was used to assess the accuracy of the classification that is based on AVHRR and geophysical data. Accuracy of individual regions varies from 73% to 89%, with an overall accuracy of 81% for China. The results showed that the methodology used in this study is, in general, feasible for large-scale land-cover mapping in China.