The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Effects of large-scale human land use on Capercaillie (Tetrao urogallus L.) populations in Finland

The Capercaillie (Tetrao urogallus L.) is often used as a focal species for landscape ecological studies: the minimum size for its lekking area is 300 ha, and the annual home range for an individual may cover 30 80 km2. In Finland, Capercaillie populations have decreased by approximately 40 85%, with the declines likely to have started in the 1940s. Although the declines have partly stabilized from the 1990s onwards, it is obvious that the negative population trend was at least partly caused by changes in human land use. The aim of this thesis was to study the connections between human land use and Capercaillie populations in Finland, using several spatial and temporal scales. First, the effect of forest age structure on Capercaillie population trends was studied in 18 forestry board districts in Finland, during 1965 1988. Second, the abundances of Capercaillie and Moose (Alces alces L.) were compared in terms of several land-use variables on a scale of 50 × 50 km grids and in five regions in Finland. Third, the effects of forest cover and fine-grain forest fragmentation on Capercaillie lekking area persistence were studied in three study locations in Finland, on 1000 and 3000 m spatial scales surrounding the leks. The analyses considering lekking areas were performed with two definitions for forest: > 60 and > 152 m3ha 1 of timber volume. The results show that patterns and processes at large spatial scales strongly influence Capercaillie in Finland. In particular, in southwestern and eastern Finland, high forest cover and low human impact were found to be beneficial for this species. Forest cover (> 60 m3ha 1 of timber) surrounding the lekking sites positively affected lekking area persistence only at the larger landscape scale (3000 m radius). The effects of older forest classes were hard to assess due to scarcity of older forests in several study areas. Young and middle-aged forest classes were common in the vicinity of areas with high Capercaillie abundances especially in northern Finland. The increase in the amount of younger forest classes did not provide a good explanation for Capercaillie population decline in 1965 1988. In addition, there was no significant connection between mature forests (> 152 m3ha 1 of timber) and lekking area persistence in Finland. It seems that in present-day Finnish landscapes, area covered with old forest is either too scarce to efficiently explain the abundance of Capercaillie and the persistence of the lekking areas, or the effect of forest age is only important when considering smaller spatial scales than the ones studied in this thesis. In conclusion, larger spatial scales should be considered for assessing the future Capercaillie management. According to the proposed multi-level planning, the first priority should be to secure the large, regional-scale forest cover, and the second priority should be to maintain fine-grained, heterogeneous structure within the separate forest patches. A management unit covering hundreds of hectares, or even tens or hundreds of square kilometers, should be covered, which requires regional-level land-use planning and co-operation between forest owners.

Risk factors for open angle glaucoma; a clinical and molecular genetic study

Glaucoma, optic neuropathy with excavation in the optic nerve head and corresponding visual field defect, is one of the leading causes for blindness worldwide. However, visual disability can often be avoided or delayed if the disease is diagnosed at an early stage. Therefore, recognising the risk factors for development and progression of glaucoma may prevent further damage. The purpose of the present study was to evaluate factors associated with visual disability caused by glaucoma and the genetic features of two risk factors, exfoliation syndrome (ES) and a positive family history of glaucoma. The present study material consisted of three study groups 1) deceased glaucoma patients from the Ekenäs practice 2) glaucoma families from the Ekenäs region and 3) population based families with and without exfoliation syndrome from Kökar Island. For the retrospective study, 106 patients with open angle glaucoma (OAG) were identified. At the last visit, 17 patients were visually impaired. Blindness induced by glaucoma was found in one or both eyes in 16 patients and in both eyes in six patients. The cumulative incidence of glaucoma caused blindness for one eye was 6% at 5 years, 9% at 10 years, and 15% at 15 years from initialising the treatment. The factors associated with blindness caused by glaucoma were an advanced stage of glaucoma at diagnosis, fluctuation in intraocular pressure during treatment, the presence of exfoliation syndrome, and poor patient compliance. A cross-sectional population based study performed in 1960-1962 on Kökar Island and the same population was followed until 2002. In total 965 subjects (530 over 50 years) have been examined at least once. The prevalence of exfoliation syndrome (ES) was 18% among subjects older than 50 years. Seventy-five of all 78 ES-positives belonged to the same extended pedigree. According to the segregation and family analysis, exfoliation syndrome seemed to be inherited as an autosomal dominant trait with reduced penetrance. The penetrance was more reduced for males, but the risk for glaucoma was higher in males than in females. To find the gene or genes associated with exfoliation syndrome, a genome wide scan was performed for 64 members (28 ES affected and 36 controls) of the Kökar pedigree. A promising result was found: the highest two-point LOD score of 3.45 (θ=0.04) in chromosome18q12.1-21.33. The presence of mutations in glaucoma genes TIGR/MYOC (myocilin) and OPTN (optineurin) was analysed in eight glaucoma families from the Ekenäs region. An inheritance pattern resembling autosomal dominant mode was detected in all these families. Primary open angle glaucoma or exfoliation glaucoma was found in 35% of 136 family members and 28% were suspected to have glaucoma. No mutations were detected in these families.

Prognostic molecular factors and algorithms in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common of the non-Hodgkin lymphomas. As DLBCL is characterized by heterogeneous clinical and biological features, its prognosis varies. To date, the International Prognostic Index has been the strongest predictor of outcome for DLBCL patients. However, no biological characters of the disease are taken into account. Gene expression profiling studies have identified two major cell-of-origin phenotypes in DLBCL with different prognoses, the favourable germinal centre B-cell-like (GCB) and the unfavourable activated B-cell-like (ABC) phenotypes. However, results of the prognostic impact of the immunohistochemically defined GCB and non-GCB distinction are controversial. Furthermore, since the addition of the CD20 antibody rituximab to chemotherapy has been established as the standard treatment of DLBCL, all molecular markers need to be evaluated in the post-rituximab era. In this study, we aimed to evaluate the predictive value of immunohistochemically defined cell-of-origin classification in DLBCL patients. The GCB and non-GCB phenotypes were defined according to the Hans algorithm (CD10, BCL6 and MUM1/IRF4) among 90 immunochemotherapy- and 104 chemotherapy-treated DLBCL patients. In the chemotherapy group, we observed a significant difference in survival between GCB and non-GCB patients, with a good and a poor prognosis, respectively. However, in the rituximab group, no prognostic value of the GCB phenotype was observed. Likewise, among 29 high-risk de novo DLBCL patients receiving high-dose chemotherapy and autologous stem cell transplantation, the survival of non-GCB patients was improved, but no difference in outcome was seen between GCB and non-GCB subgroups. Since the results suggested that the Hans algorithm was not applicable in immunochemotherapy-treated DLBCL patients, we aimed to further focus on algorithms based on ABC markers. We examined the modified activated B-cell-like algorithm based (MUM1/IRF4 and FOXP1), as well as a previously reported Muris algorithm (BCL2, CD10 and MUM1/IRF4) among 88 DLBCL patients uniformly treated with immunochemotherapy. Both algorithms distinguished the unfavourable ABC-like subgroup with a significantly inferior failure-free survival relative to the GCB-like DLBCL patients. Similarly, the results of the individual predictive molecular markers transcription factor FOXP1 and anti-apoptotic protein BCL2 have been inconsistent and should be assessed in immunochemotherapy-treated DLBCL patients. The markers were evaluated in a cohort of 117 patients treated with rituximab and chemotherapy. FOXP1 expression could not distinguish between patients, with favourable and those with poor outcomes. In contrast, BCL2-negative DLBCL patients had significantly superior survival relative to BCL2-positive patients. Our results indicate that the immunohistochemically defined cell-of-origin classification in DLBCL has a prognostic impact in the immunochemotherapy era, when the identifying algorithms are based on ABC-associated markers. We also propose that BCL2 negativity is predictive of a favourable outcome. Further investigational efforts are, however, warranted to identify the molecular features of DLBCL that could enable individualized cancer therapy in routine patient care.

Iodine Brachytherapy for Large Uveal Melanomas

Uveal melanoma is the most common primary intraocular malignancy in adults. Vision in the affected eye is threatened by both the tumor and side-effects from the treatments currently available. Poor prognosis for saving vision increases with tumor size and, consequently, enucleation has been the treatment of choice for large uveal melanomas in most centers. However, increasing evidence suggests that no survival benefit is gained (nor lost) by enucleation as compared to eye-conserving methods. The Helsinki University Eye Hospital has since 1990 offered episcleral iodine-125 plaque brachytherapy (IBT) for all patients unwilling to undergo enucleation for a large uveal melanoma. The primary aim of this study was to assess survival, local tumor recurrence and preservation of the eye and vision after IBT in a population-based series of 97 patients with uveal melanomas classified as large by the Collaborative Ocular Melanoma Study (COMS) criteria. Further aims included reporting the incidence of side-effects and assessing the role of intraocular dose distribution and clinical risk factors in their development. Finally, means to improve the current treatment were investigated by using computer models to compare existing plaques with collimating ones and by comparing the outcome of a subgroup of 54 IBT patients with very thick tumors with 33 patients with similarly-sized tumors managed with transscleral local resection (TSR) in Liverpool, United Kingdom. Kaplan-Meier estimates of all-cause and melanoma-specific survival at 5 years after IBT were 62% and 65%, respectively, and visually comparable with the survival experience of patients reported after enucleation by the COMS. Local recurrence developed in 6% of eyes and 84% of eyes were conserved at 5 years. Visual prognosis was guarded with 11% avoiding loss of 20/70 vision and 26% avoiding loss of 20/400 vision in the tumor eye at 2 years. Large tumor height and short distance from the posterior pole were independently associated with loss of vision. Using cumulative incidence analysis to account for competing risks, such as enucleation and metastatic death, the 5-year incidence of cataract after IBT was 79%, glaucoma 60%, optic neuropathy 46%, maculopathy 52%, persistent or recurring retinal detachment (RD) 25%, and vitreous hemorrhage 36%. In multivariate competing risks regression models, increasing tumor height was associated with cataract, iris neovascularization and RD. Maculopathy and optic neuropathy were associated with distance from the tumor to the respective structure. Median doses to the tumor apex, macula and optic disc were 81 Gy (range, 40-158), 79 Gy (range, 12-632), and 83 Gy (range, 10-377), respectively. Dose to the optic disc was independently associated with optic neuropathy, and both dose to the optic disc and dose to the macula predicted vision loss after IBT. Simulated treatment using collimating plaques resulted in clinically meaningful reduction in both optic disc (median reduction, 30 Gy) and macular (median reduction, 36 Gy) doses as compared to the actual treatment with standard plaques. In the subgroup of patients with uveal melanomas classified as large because of tumor height, cumulative incidence analysis revealed that while long-term preservation of 20/70 vision was rare after both IBT and TSR, preservation of 20/400 vision was better after TSR (32% vs. 5% at 5 years). In multivariate logistic regression models, TSR was independently associated with better preservation of 20/400 vision (OR 0.03 at 2 years, P=0.005) No cases of secondary glaucoma were observed after TSR and optic neuropathy was rare. However, local tumor recurrence was more common after TSR than it was after IBT (Cumulative incidence 41% vs. 7% at 5 years, respectively). In terms of survival, IBT seems to be a safe alternative to enucleation in managing large uveal melanomas. Local tumor control is no worse than with medium-sized tumors and the chances of avoiding secondary enucleation are good. Unfortunately, side-effects from radiotherapy are frequent, especially in thick tumors, and long-term prognosis of saving vision is consequently guarded. Some complications can be limited by using collimating plaques and by managing uveal melanomas that are large because of tumor height with TSR instead of IBT. However, the patient must be willing to accept a substantial risk of local tumor recurrence after TSR and it is best suited for cases in which the preservation of vision in the tumor eye is critical.

On the mechanical behavior of compacted pack ice : A theoretical and numerical investigation

Pack ice is an aggregate of ice floes drifting on the sea surface. The forces controlling the motion and deformation of pack ice are air and water drag forces, sea surface tilt, Coriolis force and the internal force due to the interaction between ice floes. In this thesis, the mechanical behavior of compacted pack ice is investigated using theoretical and numerical methods, focusing on the three basic material properties: compressive strength, yield curve and flow rule. A high-resolution three-category sea ice model is applied to investigate the sea ice dynamics in two small basins, the whole Gulf Riga and the inside Pärnu Bay, focusing on the calibration of the compressive strength for thin ice. These two basins are on the scales of 100 km and 20 km, respectively, with typical ice thickness of 10-30 cm. The model is found capable of capturing the main characteristics of the ice dynamics. The compressive strength is calibrated to be about 30 kPa, consistent with the values from most large-scale sea ice dynamic studies. In addition, the numerical study in Pärnu Bay suggests that the shear strength drops significantly when the ice-floe size markedly decreases. A characteristic inversion method is developed to probe the yield curve of compacted pack ice. The basis of this method is the relationship between the intersection angle of linear kinematic features (LKFs) in sea ice and the slope of the yield curve. A summary of the observed LKFs shows that they can be basically divided into three groups: intersecting leads, uniaxial opening leads and uniaxial pressure ridges. Based on the available observed angles, the yield curve is determined to be a curved diamond. Comparisons of this yield curve with those from other methods show that it possesses almost all the advantages identified by the other methods. A new constitutive law is proposed, where the yield curve is a diamond and the flow rule is a combination of the normal and co-axial flow rule. The non-normal co-axial flow rule is necessary for the Coulombic yield constraint. This constitutive law not only captures the main features of forming LKFs but also takes the advantage of avoiding overestimating divergence during shear deformation. Moreover, this study provides a method for observing the flow rule for pack ice during deformation.

Scattering and Refraction as Contrast Mechanisms in X-Ray Imaging

Differentiation of various types of soft tissues is of high importance in medical imaging, because changes in soft tissue structure are often associated with pathologies, such as cancer. However, the densities of different soft tissues may be very similar, making it difficult to distinguish them in absorption images. This is especially true when the consideration of patient dose limits the available signal-to-noise ratio. Refraction is more sensitive than absorption to changes in the density, and small angle x-ray scattering on the other hand contains information about the macromolecular structure of the tissues. Both of these can be used as potential sources of contrast when soft tissues are imaged, but little is known about the visibility of the signals in realistic imaging situations. In this work the visibility of small-angle scattering and refraction in the context of medical imaging has been studied using computational methods. The work focuses on the study of analyzer based imaging, where the information about the sample is recorded in the rocking curve of the analyzer crystal. Computational phantoms based on simple geometrical shapes with differing material properties are used. The objects have realistic dimensions and attenuation properties that could be encountered in real imaging situations. The scattering properties mimic various features of measured small-angle scattering curves. Ray-tracing methods are used to calculate the refraction and attenuation of the beam, and a scattering halo is accumulated, including the effect of multiple scattering. The changes in the shape of the rocking curve are analyzed with different methods, including diffraction enhanced imaging (DEI), extended DEI (E-DEI) and multiple image radiography (MIR). A wide angle DEI, called W-DEI, is introduced and its performance is compared with that of the established methods. The results indicate that the differences in scattered intensities from healthy and malignant breast tissues are distinguishable to some extent with reasonable dose. Especially the fraction of total scattering has large enough differences that it can serve as a useful source of contrast. The peaks related to the macromolecular structure come to angles that are rather large, and have intensities that are only a small fraction of the total scattered intensity. It is found that such peaks seem to have only limited usefulness in medical imaging. It is also found that W-DEI performs rather well when most of the intensity remains in the direct beam, indicating that dark field imaging methods may produce the best results when scattering is weak. Altogether, it is found that the analysis of scattered intensity is a viable option even in medical imaging where the patient dose is the limiting factor.

On the homogeneous and heterogeneous nucleation of some organic compounds

The conversion of a metastable phase into a thermodynamically stable phase takes place via the formation of clusters. Clusters of different sizes are formed spontaneously within the metastable mother phase, but only those larger than a certain size, called the critical size, will end up growing into a new phase. There are two types of nucleation: homogeneous, where the clusters appear in a uniform phase, and heterogeneous, when pre-existing surfaces are available and clusters form on them. The nucleation of aerosol particles from gas-phase molecules is connected not only with inorganic compounds, but also with nonvolatile organic substances found in atmosphere. The question is which ones of the myriad of organic species have the right properties and are able to participate in nucleation phenomena. This thesis discusses both homogeneous and heterogeneous nucleation, having as theoretical tool the classical nucleation theory (CNT) based on thermodynamics. Different classes of organics are investigated. The members of the first class are four dicarboxylic acids (succinic, glutaric, malonic and adipic). They can be found in both the gas and particulate phases, and represent good candidates for the aerosol formation due to their low vapor pressure and solubility. Their influence on the nucleation process has not been largely investigated in the literature and it is not fully established. The accuracy of the CNT predictions for binary water-dicarboxylic acid systems depends significantly on the good knowledge of the thermophysical properties of the organics and their aqueous solutions. A large part of the thesis is dedicated to this issue. We have shown that homogeneous and heterogeneous nucleation of succinic, glutaric and malonic acids in combination with water is unlikely to happen in atmospheric conditions. However, it seems that adipic acid could participate in the nucleation process in conditions occurring in the upper troposphere. The second class of organics is represented by n-nonane and n-propanol. Their thermophysical properties are well established, and experiments on these substances have been performed. The experimental data of binary homogeneous and heterogeneous nucleation have been compared with the theoretical predictions. Although the n-nonane - n-propanol mixture is far from being ideal, CNT seems to behave fairly well, especially when calculating the cluster composition. In the case of heterogeneous nucleation, it has been found that better characterization of the substrate - liquid interaction by means of line tension and microscopic contact angle leads to a significant improvement of the CNT prediction. Unfortunately, this can not be achieved without well defined experimental data.

Lastensuojelun sosiaalityö asiakirjoina : Dokumentoinnin ja tiedonmuodostuksen dynamiikka

Tutkimuksen keskeinen tehtävä on selvittää, mikä on dokumentoinnin merkitys lastensuojelun sosiaalityön tiedonmuodostuksessa ja ammattikäytännöissä. Asiakirjateksteistä koostuvaa tutkimusaineistoa tarkastellaan kolmesta eri suunnasta kysymällä: 1)Miten asiakirjoja kirjoitetaan? 2) Mitä asiakirjoihin kirjoitetaan? 3) Miksi asiakirjoja kirjoitetaan niin kuin kirjoitetaan? Tutkimusaineisto muodostuu lastensuojelun sosiaalityöntekijöiden laatimista asiakastietojärjestelmään tallennetuista muistiinpanoista ja huostaanottopäätöksistä. Tutkimukseen on valittu 20 huostaanotetun eri-ikäisen lapsen ja heidän perheensä asiakirjat yhteensä 1613 asiakirjatulostussivua. Tekstit ajoittuvat vuodesta 1989 vuoteen 2000. Tutkimusmenetelmä on diskurssianalyyttinen ja tukeutuu Fairclough`n (1997)esittämään kolmiulotteiseen malliin, jossa diskurssi määritellään tekstin, käytäntöjen ja sosiokulttuurisen ympäristön suhteeksi. Diskurssianalyysi on näiden rakenteiden ja niiden välisten suhteiden kuvaamista, tulkintaa ja selittämistä. Fairclough’n mallia mukaillen tutkimuksen analyysi koostuu retoriikan ja tematiikan analyyseistä sekä pragmatiikan näkökulman sisältävästä tarkastelusta. Asiakirjatekstien pilkkominen puhujakategorioihin osoitti tekstien olevan moniäänisiä, useiden henkilöiden näkemyksiä ja mielipiteitä sisältäviä tekstipintoja. Retoriikan analyysi näytti, että lastensuojelun sosiaalityön asiakirjat sisältävät paljon dynaamisia kuvauksia työstä. Asiakirjojen kirjoittaminen moniäänisiksi tuo tekstiin uskottavuutta, ja se on myös yksi retorinen vaikuttamiskeino. Tematiikan tarkastelu osoitti,että asiakirjojen sisällölliset teemat (lapsen hoiva, arjen hallinta, yhteistyö ja päihteiden käyttö) ja kokemukselliset teemat (huoli, vastuu, yhteys ja moraali) toistuvat sisäkkäisinä ja päällekkäisinä säikeinä dynaamisesti vaihdellen. Sosiaalityöntekijät kirjaavat teksteihin monia yhtäaikaisia teemoja, joiden avulla rakentavat ammatillista ymmärrystä kyseessä olevasta tilanteesta. Asiakirjojen tutkiminen pragmatiikan suunnasta toi esiin, kirjoittamisen ja lukemisen kontekstiulottuvuudet sekä tiedonmuodostusprosessin. Asiakirjojen laatiminen on osa sosiaalityön käytäntöjä. Se on myös keskeinen alue ammattikunnan yhteisen ammatillisen ymmärryksen luomisessa ja ylläpitämisessä. Muistiinpanot, huostaanottopäätökset ja lakitekstit ovat intertekstuaalisia. Lastensuojelun sosiaalityön asiakirjojen tutkiminen on avannut uusia mahdollisuuksia ymmärtää sosiaalityön dokumentointiprosessia, merkitystä ja roolia sekä tiedonmuodostuksen dynamiikkaa. Tekstien kirjoittaminen, niiden lukeminen, tietojen siirtäminen ja asiakkaan kuuleminen samoin kuin kuulemisen kirjaaminen ovat sosiaalityön dokumentoinnin keskeisiä haasteita. Tutkimus pyrkii avaamaan ymmärrystä asiakirjatekstien monivivahteiseen ja dynaamiseen maailmaan ja siten myös sosiaalityön dokumentoinnin arkeen. Tarkastelut mahdollistavat työn kehittämisen erityisesti sosiaalityön asiakasvaikuttavuuden mittaamisen ja parantamisen suuntaan. Asiakirjoissa ilmenevä tiedonmuodostuksen dynamiikka syntyy kirjoittamiskäytäntöjen, kirjoittamisen ja lukemisen sekä toimintakäytäntöjen yhteisessä alueessa. Avainsanat: sosiaalityö, lastensuojelu, dokumentointi, asiakirja, diskurssianalyysi, tiedonmuodostus.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.