Modulation of growth by aromatase inhibitor treatment in boys: Efficacy and safety

Without estrogen action, the fusion of the growth plates is postponed and statural growth continues for an exceptionally long time. Aromatase inhibitors, blockers of estrogen biosynthesis, have therefore emerged as a new potential option for the treatment of children with short stature. We investigated the efficacy of the aromatase inhibitor letrozole in the treatment of boys with idiopathic short stature (ISS) using a randomised, placebo-controlled, double-blind research setting. A total of 30 boys completed the two-year treatment. By decreasing estrogen-mediated central negative feedback, letrozole increased gonadotrophin and testosterone secretion in pubertal boys, whereas the pubertal increase in IGF-I was inhibited. Treatment with letrozole effectively delayed bone maturation and increased predicted adult height by 5.9 cm (P0.001), while placebo had no effect on either parameter. The effect of letrozole treatment on near-final height was studied in another population, in boys with constitutional delay of puberty, who received letrozole (n=9) or placebo (n=8) for one year, in combination with low-dose testosterone for six months during adolescence. The mean near-final height of boys randomised to receive testosterone and letrozole was significantly greater than that of boys who received testosterone and placebo (175.8 vs. 169.1 cm, P=0.04). As regards safety, treatment effects on bone health, lipid metabolism, insulin sensitivity, and body composition were monitored in boys with ISS. During treatment, no differences in bone mass accrual were evident between the treatment groups, as evaluated by dual-energy x-ray absorptiometry measurements of the lumbar spine and femoral neck. Bone turnover and cortical bone growth, however, were affected by letrozole treatment. As indicated by differences in markers of bone resorption (U-INTP) and formation (S-PINP and S-ALP), the long-term rate of bone turnover was lower in letrozole-treated boys, despite their more rapid advancement in puberty. Letrozole stimulated cortical bone growth in those who progressed in puberty: the metacarpal index (MCI), a measure of cortical bone thickness, increased more in letrozole-treated pubertal boys than in placebo-treated pubertal boys (25% vs. 9%, P=0.007). The change in MCI correlated positively with the mean testosterone-to-estradiol ratio. In post-treatment radiographic evaluation of the spine, a high rate of vertebral deformities - mild anterior wedging and mild compression deformities - were found in both placebo and letrozole groups. In pubertal boys with ISS treated with letrozole, stimulated testosterone secretion was associated with a decrease in the percentage of fat mass and in HDL-cholesterol, while LDL-cholesterol and triglycerides remained unchanged. Insulin sensitivity, as evaluated by HOMA-IR, was not significantly affected by the treatment. In summary, treatment with the aromatase inhibitor letrozole effectively delayed bone maturation and increased predicted adult height in boys with ISS. Long-term follow-up data of boys with constitutional delay of puberty, treated with letrozole for one year during adolescence, suggest that the achieved gain in predicted adult height also results in increased adult height. However, until the safety of aromatase inhibitor treatment in children and adolescents is confirmed, such treatment should be considered experimental.

Intrahepatic cholestasis of pregnancy : genetic background, epidemiology and hepatobiliary consequences

Intrahepatic cholestasis of pregnancy (ICP) is the most common cholestatic liver disease during pregnancy. The reported incidence varies from 0.4 to 15% of full-term pregnancies. The etiology is heterogeneous but familial clustering is known to occur. Here we have studied the genetic background, epidemiology, and long-term hepatobiliary consequences of ICP. In a register-based nation-wide study (n=1 080 310) the incidence of ICP was 0.94% during 1987-2004. A slightly higher incidence, 1.3%, was found in a hospital-based series (n=5304) among women attending the University Hospital of Helsinki in 1992-1993. Of these 16% (11/69) were familial and showed a higher (92%) recurrence rate than the sporadic (40%) cases. In the register-based epidemiological study, advanced maternal age and, to a lesser degree, parity were identified as new risk factors for ICP. The risk was 3-fold higher in women >39 years of age compared to women <30 years. Multiple pregnancy also associated with an elevated risk. In a genetic study we found no association of ICP with the genes regulating bile salt transport (ABCB4, ABCB11 and ATP8B1). The livers of postmenopausal women with a history of ICP tolerated well the short-term exposure to oral and transdermal estradiol, although the doses used were higher than those in routine clinical use. The response of serum levels of sex hormone-binding globulin (SHBG) to oral estradiol was slightly reduced in the ICP group. Transdermal estradiol had no effect on C-reactive protein (CRP) or SHBG. A number of liver and biliary diseases were found to be associated with ICP. Women with a history of ICP showed elevated risks for non-alcoholic liver cirrhosis (8.2 CI 1.9-36), cholelithiasis and cholecystitis (3.7 CI 3.2-4.2), hepatitis C (3.5 CI 1.6-7.6) and non-alcoholic pancreatitis (3.2 CI 1.7-5.7). In conclusion, ICP complicates around 1% of all full-term pregnancies in Finland and its incidence has remained unchanged since 1987. It is familial in 16% of cases with a higher recurrence rate. Although the cause remains unknown, several risk factors, namely advanced maternal age, parity and multiple pregnancies, can be identified. Both oral and transdermal regimens of postmenopausal hormone therapy (HT) are safe for women with a history of ICP when liver function is considered. Some ICP patients are at risk of other liver and biliary diseases and, contrary to what has been thought, a follow-up is warranted.

Acute Lymphoblastic Leukemia in Adolescents and Young Adults in Finland

In recent reports, adolescents and young adults (AYA) with acute lymphoblastic leukemia (ALL) have had a better outcome with pediatric treatment than with adult protocols. ALL can be classified into biologic subgroups according to immunophenotype and cytogenetics, with different clinical characteristics and outcome. The proportions of the subgroups are different in children and adults. ALL subtypes in AYA patients are less well characterized. In this study, the treatment and outcome of ALL in AYA patients aged 10-25 years in Finland on pediatric and adult protocols was retrospectively analyzed. In total, 245 patients were included. The proportions of biologic subgroups in different age groups were determined. Patients with initially normal or failed karyotype were examined with oligonucleotide microarray-based comparative genomic hybridization (aCGH). Also deletions and instability of chromosome 9p were screened in ALL patients. In addition, patients with other hematologic malignancies were screened for 9p instability. aCGH data were also used to determine a gene set that classifies AYA patients at diagnosis according to their risk of relapse. Receiver operating characteristic analysis was used to assess the value of the set of genes as prognostic classifiers. The 5-year event-free survival of AYA patients treated with pediatric or adult protocols was 67% and 60% (p=0.30), respectively. White blood cell count larger than 100x109/l was associated with poor prognosis. Patients treated with pediatric protocols and assigned to an intermediate-risk group fared significantly better than those of the pediatric high-risk or adult treatment groups. Deletions of 9p were detected in 46% of AYA ALL patients. The chromosomal region 9p21.3 was always affected, and the CDKN2A gene was always deleted. In about 15% of AYA patients, the 9p21.3 deletion was smaller than 200 kb in size, and therefore, probably undetectable with conventional methods. Deletion of 9p was the most common aberration of AYA ALL patients with initially normal karyotype. Instability of 9p, defined as multiple separate areas of copy number loss or homozygous loss within a larger heterozygous area in 9p, was detected in 19% (n=27) of ALL patients. This abnormality was restricted to ALL; none of the patients with other hematologic malignancies had the aberration. The prognostic model identification procedure resulted in a model of four genes: BAK1, CDKN2B, GSTM1, and MT1F. The copy number profile combinations of these genes differentiated between AYA ALL patients at diagnosis depending on their risk of relapse. Deletions of CDKN2B and BAK1 in combination with amplification of GSTM1 and MT1F were associated with a higher probability of relapse. Unlike all previous studies, we found that the outcome of AYA patients with ALL treated using pediatric or adult therapeutic protocols was comparable. The success of adult ALL therapy emphasizes the benefit of referral of patients to academic centers and adherence to research protocols. 9p deletions and instability are common features of ALL and may act together with oncogene-activating translocations in leukemogenesis. New and more sensitive methods of molecular cytogenetics can reveal previously cryptic genetic aberrations with an important role in leukemic development and prognosis and that may be potential targets of therapy. aCGH also provides a viable approach for model design aiming at evaluation of risk of relapse in ALL.

Biology and Molecular Markers of Malignant Gonadal Germ Cell Tumors

Germ cell tumors occur both in the gonads of both sexes and in extra-gonadal sites during adoles-cence and early adulthood. Malignant ovarian germ cell tumors are rare neoplasms accounting for less than 5% of all cases of ovarian malignancy. In contrast, testicular cancer is the most common malignancy among young males. Most of patients survive the disease. Prognostic factors of gonadal germ cell tumors include histology, clinical stage, size of the primary tumor and residua, and levels of tumor markers. Germ cell tumors include heterogeneous histological subgroups. The most common subgroup includes germinomas (ovarian dysgerminoma and testicular seminoma); other subgroups are yolk sac tumors, embryonal carcinomas, immature teratomas and mixed tumors. The origin of germ cell tumors is most likely primordial germ cells. Factors behind germ cell tumor development and differentiation are still poorly known. The purpose of this study was to define novel diagnostic and prognostic factors for malignant gonadal germ cell tumors. In addition, the aim was to shed further light into the molecular mechanisms regulating gonadal germ cell tumorigenesis and differentiation by studying the roles of GATA transcription factors, pluripotent factors Oct-3/4 and AP-2γ, and estrogen receptors. This study revealed the prognostic value of CA-125 in malignant ovarian germ cell tumors. In addition advanced age and residual tumor had more adverse outcome. Several novel markers for histological diagnosis were defined. In the fetal development transcription factor GATA-4 was expressed in early fetal gonocytes and in testicular carcinoma precursor cells. In addition, GATA-4 was expressed in both gonadal germinomas, thus it may play a role in the development and differentiation of the germinoma tumor subtype. Pluripotent factors Oct-3/4 and AP-2γ were expressed in dysgerminomas, thus they could be used in the differential diagnosis of the germ cell tumors. Malignant ovarian germ cell tumors expressed estrogen receptors and their co-regulator SNURF. In addition, estrogen receptor expression was up-regulated by estradiol stimulation. Thus, gonadal steroid hormone burst in puberty may play a role in germ cell tumor development in the ovary. This study shed further light in to the molecular pathology of malignant gonadal germ cell tumors. In addition, some novel diagnostic and prognostic factors were defined. This data may be used in the differential diagnosis of germ cell tumor patients.

Imaging of knee injuries with special focus on tibial plateau fractures

Acute knee injury is a common event throughout life, and it is usually the result of a traffic accident, simple fall, or twisting injury. Over 90% of patients with acute knee injury undergo radiography. An overlooked fracture or delayed diagnosis can lead to poor patient outcome. The major aim of this thesis was retrospectively to study imaging of knee injury with a special focus on tibial plateau fractures in patients referred to a level-one trauma center. Multi-detector computed tomography (MDCT) findings of acute knee trauma were studied and compared to radiography, as well as whether non-contrast MDCT can detect cruciate ligaments with reasonable accuracy. The prevalence, type, and location of meniscal injuries in magnetic resonance imaging (MRI) were evaluated, particularly in order to assess the prevalence of unstable meniscal tears in acute knee trauma with tibial plateau fractures. The possibility to analyze with conventional MRI the signal appearance of menisci repaired with bioabsorbable arrows was also studied. The postoperative use of MDCT was studied in surgically treated tibial plateau fractures: to establish the frequency and indications of MDCT and to assess the common findings and their clinical impact in a level-one trauma hospital. This thesis focused on MDCT and MRI of knee injuries, and radiographs were analyzed when applica-ble. Radiography constitutes the basis for imaging acute knee injury, but MDCT can yield information beyond the capabilities of radiography. Especially in severely injured patients , sufficient radiographs are often difficult to obtain, and in those patients, radiography is unreliable to rule out fractures. MDCT detected intact cruciate ligaments with good specificity, accuracy, and negative predictive value, but the assessment of torn ligaments was unreliable. A total of 36% (14/39) patients with tibial plateau fracture had an unstable meniscal tear in MRI. When a meniscal tear is properly detected preoperatively, treatment can be combined with primary fracture fixation, thus avoiding another operation. The number of meniscal contusions was high. Awareness of the imaging features of this meniscal abnormality can help radiologists increase specificity by avoiding false-positive findings in meniscal tears. Postoperative menisci treated with bioabsorbable arrows showed no difference, among different signal intensities in MRI, among menisci between patients with operated or intact ACL. The highest incidence of menisci with an increased signal intensity extending to the meniscal surface was in patients whose surgery was within the previous 18 months. The results may indicate that a rather long time is necessary for menisci to heal completely after arrow repair. Whether the menisci with an increased signal intensity extending to the meniscal surface represent improper healing or re-tear, or whether this is just the earlier healing feature in the natural process remains unclear, and further prospective studies are needed to clarify this. Postoperative use of MDCT in tibial plateau fractures was rather infrequent even in this large trauma center, but when performed, it revealed clinically significant information, thus benefitting patients in regard to treatment.

Cyclosporine population pharmacokinetics in pediatric renal transplant recipients

Cyclosporine is an immunosuppressant drug with a narrow therapeutic index and large variability in pharmacokinetics. To improve cyclosporine dose individualization in children, we used population pharmacokinetic modeling to study the effects of developmental, clinical, and genetic factors on cyclosporine pharmacokinetics in altogether 176 subjects (age range: 0.36–20.2 years) before and up to 16 years after renal transplantation. Pre-transplantation test doses of cyclosporine were given intravenously (3 mg/kg) and orally (10 mg/kg), on separate occasions, followed by blood sampling for 24 hours (n=175). After transplantation, in a total of 137 patients, cyclosporine concentration was quantified at trough, two hours post-dose, or with dose-interval curves. One-hundred-four of the studied patients were genotyped for 17 putatively functionally significant sequence variations in the ABCB1, SLCO1B1, ABCC2, CYP3A4, CYP3A5, and NR1I2 genes. Pharmacokinetic modeling was performed with the nonlinear mixed effects modeling computer program, NONMEM. A 3-compartment population pharmacokinetic model with first order absorption without lag-time was used to describe the data. The most important covariate affecting systemic clearance and distribution volume was allometrically scaled body weight i.e. body weight**3/4 for clearance and absolute body weight for volume of distribution. The clearance adjusted by absolute body weight declined with age and pre-pubertal children (< 8 years) had an approximately 25% higher clearance/body weight (L/h/kg) than did older children. Adjustment of clearance for allometric body weight removed its relationship to age after the first year of life. This finding is consistent with a gradual reduction in relative liver size towards adult values, and a relatively constant CYP3A content in the liver from about 6–12 months of age to adulthood. The other significant covariates affecting cyclosporine clearance and volume of distribution were hematocrit, plasma cholesterol, and serum creatinine, explaining up to 20%–30% of inter-individual differences before transplantation. After transplantation, their predictive role was smaller, as the variations in hematocrit, plasma cholesterol, and serum creatinine were also smaller. Before transplantation, no clinical or demographic covariates were found to affect oral bioavailability, and no systematic age-related changes in oral bioavailability were observed. After transplantation, older children receiving cyclosporine twice daily as the gelatine capsule microemulsion formulation had an about 1.25–1.3 times higher bioavailability than did the younger children receiving the liquid microemulsion formulation thrice daily. Moreover, cyclosporine oral bioavailability increased over 1.5-fold in the first month after transplantation, returning thereafter gradually to its initial value in 1–1.5 years. The largest cyclosporine doses were administered in the first 3–6 months after transplantation, and thereafter the single doses of cyclosporine were often smaller than 3 mg/kg. Thus, the results suggest that cyclosporine displays dose-dependent, saturable pre-systemic metabolism even at low single doses, whereas complete saturation of CYP3A4 and MDR1 (P-glycoprotein) renders cyclosporine pharmacokinetics dose-linear at higher doses. No significant associations were found between genetic polymorphisms and cyclosporine pharmacokinetics before transplantation in the whole population for which genetic data was available (n=104). However, in children older than eight years (n=22), heterozygous and homozygous carriers of the ABCB1 c.2677T or c.1236T alleles had an about 1.3 times or 1.6 times higher oral bioavailability, respectively, than did non-carriers. After transplantation, none of the ABCB1 SNPs or any other SNPs were found to be associated with cyclosporine clearance or oral bioavailability in the whole population, in the patients older than eight years, or in the patients younger than eight years. In the whole population, in those patients carrying the NR1I2 g.-25385C–g.-24381A–g.-205_-200GAGAAG–g.7635G–g.8055C haplotype, however, the bioavailability of cyclosporine was about one tenth lower, per allele, than in non-carriers. This effect was significant also in a subgroup of patients older than eight years. Furthermore, in patients carrying the NR1I2 g.-25385C–g.-24381A–g.-205_-200GAGAAG–g.7635G–g.8055T haplotype, the bioavailability was almost one fifth higher, per allele, than in non-carriers. It may be possible to improve individualization of cyclosporine dosing in children by accounting for the effects of developmental factors (body weight, liver size), time after transplantation, and cyclosporine dosing frequency/formulation. Further studies are required on the predictive value of genotyping for individualization of cyclosporine dosing in children.

The effect of blood glucose on the vasculature in young patients with type 1 diabetes

Background. Patients with type 1 diabetes are at markedly increased risk of vascular complications. In this respect it is noteworthy that hyperglycaemia that is shown to cause endothelial dysfunction, has clearly been shown to be a risk factor for diabetic microvascular disease. However, the role of hyperglycaemia as a predictor of macrovascular disease is not as clear as for microvascular disease, although type 1 diabetes itself increases the risk of cardiovascular disease substantially. Furthermore, it is not known whether it is the short-term or the long-term hyperglycaemia that confers possible risk. In addition, the role of glucose variability as a predictor of complications is to a large extent unexplored. Interestingly, although hyperglycaemia increases the risk of pre-eclampsia in women with type 1 diabetes, it is unclear whether pre-eclampsia, a condition characterized by endothelial dysfunction, is also a risk factor for microvascular complication, diabetic nephropathy. Aims. This doctoral thesis investigated the role of acute hyperglycaemia and glucose variability on arterial stiffness and cardiac ventricular repolarisation in male patients with type 1 diabetes as well as in healthy male volunteers. The thesis also explored whether acute hyperglycaemia leads to an inflammatory response, endothelial dysfunction and oxidative stress. Finally, the role of pre-eclampsia, as a predictor of diabetic nephropathy in type 1 diabetes was examined. Subjects and methods. In order to study glucose variability and the daily glycaemic control, 22 male patients with type 1 diabetes, without any diabetic complications, were monitored for 72-h with a continuous glucose monitoring system. At the end of the 72-h glucose monitoring period a 2-h hyperglycaemic clamp was performed both in the patients with type 1 diabetes and in the 13 healthy age-matched male volunteers. Blood pressure, arterial stiffness and QT time were measured to detect vascular changes during acute hyperglycaemia. Blood samples were drawn at baseline (normoglycaemia) and during acute hyperglycaemia. In another patient sample, women with type 1 diabetes were followed during their pregnancy and restudied eleven years later to elucidate the role of pre-eclampsia and pregnancy-induced hypertension as potential risk factors for diabetic nephropathy. Results and conclusions. Acute hyperglycaemia increased arterial stiffness as well as caused a disturbance in the myocardial ventricular repolarisation, emphasizing the importance of a strict daily glycaemic control in male patients with type 1 diabetes. An inflammatory response was also observed during acute hyperglycaemia. Furthermore, a high mean daily blood glucose but not glucose variability per se is associated with arterial stiffness. While glucose variability in turn correlated with central blood pressure, the results suggest that the glucose metabolism is closely linked to the haemodynamic changes in male patients with uncomplicated type 1 diabetes. Notably, the results are not directly applicable to females. Finally, a history of a pre-eclamptic pregnancy, but not pregnancy-induced hypertension was associated with increased risk of diabetic nephropathy.

Towards the Use of Radar Winds in Numerical Weather Prediction

Numerical weather prediction (NWP) models provide the basis for weather forecasting by simulating the evolution of the atmospheric state. A good forecast requires that the initial state of the atmosphere is known accurately, and that the NWP model is a realistic representation of the atmosphere. Data assimilation methods are used to produce initial conditions for NWP models. The NWP model background field, typically a short-range forecast, is updated with observations in a statistically optimal way. The objective in this thesis has been to develope methods in order to allow data assimilation of Doppler radar radial wind observations. The work has been carried out in the High Resolution Limited Area Model (HIRLAM) 3-dimensional variational data assimilation framework. Observation modelling is a key element in exploiting indirect observations of the model variables. In the radar radial wind observation modelling, the vertical model wind profile is interpolated to the observation location, and the projection of the model wind vector on the radar pulse path is calculated. The vertical broadening of the radar pulse volume, and the bending of the radar pulse path due to atmospheric conditions are taken into account. Radar radial wind observations are modelled within observation errors which consist of instrumental, modelling, and representativeness errors. Systematic and random modelling errors can be minimized by accurate observation modelling. The impact of the random part of the instrumental and representativeness errors can be decreased by calculating spatial averages from the raw observations. Model experiments indicate that the spatial averaging clearly improves the fit of the radial wind observations to the model in terms of observation minus model background (OmB) standard deviation. Monitoring the quality of the observations is an important aspect, especially when a new observation type is introduced into a data assimilation system. Calculating the bias for radial wind observations in a conventional way can result in zero even in case there are systematic differences in the wind speed and/or direction. A bias estimation method designed for this observation type is introduced in the thesis. Doppler radar radial wind observation modelling, together with the bias estimation method, enables the exploitation of the radial wind observations also for NWP model validation. The one-month model experiments performed with the HIRLAM model versions differing only in a surface stress parameterization detail indicate that the use of radar wind observations in NWP model validation is very beneficial.

X-ray scattering and diffraction enhanced imaging studies of in vitro breast tissues

This thesis is a study of the x-ray scattering properties of tissues and tumours of the breast. Clinical radiography is based on the absorption of the x-rays when passing right through the human body and gives information about the densities of the tissues. Besides being absorbed, x-rays may change their direction within the tissues due to elastic scattering or even to refraction. The phenomenon of scattering is a nuisance to radiography in general, and to mammography in particular, because it reduces the quality of the images. However, scattered x-rays bear very useful information about the structure of the tissues at the supra-molecular level. Some pathologies, like breast cancer, produce alterations to the structures of the tissues, being especially evident in collagen-rich tissues. On the other hand, the change of direction due to refraction of the x-rays on the tissue boundaries can be mapped. The diffraction enhanced imaging (DEI) technique uses a perfect crystal to convert the angular deviations of the x-rays into intensity variations, which can be recorded as images. This technique is of especial interest in the cases were the densities of the tissues are very similar (like in mammography) and the absorption images do not offer enough contrast. This thesis explores the structural differences existing in healthy and pathological collagen in breast tissue samples by the small-angle x-ray scattering (SAXS) technique and compares these differences with the morphological information found in the DEI images and the histo-pathology of the same samples. Several breast tissue samples were studied by SAXS technique in the European Synchrotron Radiation Facility (ESRF) in Grenoble, France. Scattering patterns of the different tissues of the breast were acquired and compared with the histology of the samples. The scattering signals from adipose tissue (fat), connective tissue (collagen) and necrotic tissue were identified. Moreover, a clear distinction could be done between the scattering signals from healthy collagen and from collagen from an invasive tumour. Scattering from collagen is very characteristic. It includes several scattering peaks and scattering features that carry information about the size and the spacing of the collagen fibrils in the tissues. It was found that the collagen fibrils in invaded tumours were thinner and had a d-spacing length 0,7% longer that fibrils from healthy tumours. The scattering signals from the breast tissues were compared with the histology by building colour-coded maps across the samples. They were also imaged with the DEI technique. There was a total agreement between the scattering maps, the morphological features seen in the images and the information of the histo- pathological examination. The thesis demonstrates that the x-ray scattering signal can be used to characterize tissues and that it carries important information about the pathological state of the breast tissues, thus showing the potential of the SAXS technique as a possible diagnostic tool for breast cancer.

Rekonnektio Maan magnetosfäärissä

Tässä väitöskirjassa perehdytään magneettisen rekonnektion ilmenemismuotoihin ja vaikutuksiin Maan magnetosfäärissä. Keskeisenä tutkimusvälineenä käytetään magnetohydrodynaamista (MHD) Gumics-magnetosfäärisimulaatiota. Työssä kehitetään myös uusia menetelmiä simulaatiossa ilmenevän rekonnektion tunnistamiseksi ja mittaamiseksi. MHD-simulaatio sopii suuren mittakaavan ilmiöiden tarkasteluun, joten kuvaa rekonnektiosta täydennetään pienen mittakaavan piirteiden osalta Cluster-satelliittien avulla. Tärkein tutkimuksen tuoma edistysaskel menetelmien saralla on rekonnektioviivan paikallistaminen topologisesti erityyppisten magneettikenttäviivojen alueiden liitoskohdassa olevana erottajaviivana neljän kentän tienoon menetelmää käyttäen. Tämä topologinen lähestymistapa on hyödyllinen erityisesti magnetopausilla, jonka monimutkainen geometria tekee magneettikentän paikalliseen käyttäytymiseen perustuvien rekonnektioviivan etsintätapojen soveltamisen hankalaksi. Topologisesti määritelty rekonnektioviiva on myös helppo tunnistaa magnetosfäärin globaalin konvektion solmukohdaksi. Magnetopausin rekonnektioviivan käyttäytyminen Gumicsissa noudattaa komponenttirekonnektio-olettamaan pohjautuvia teoreettisia ennusteita. Rekonnektion kvantitatiivinen tarkastelu Gumics-simulaatiossa perustuu energian muuntumiseen, joka lasketaan Poyntingin vektorin divergenssinä tai Poyntingin vuona valitun umpinaisen pinnan läpi. Rekonnektioon liittyvän energian muuntumisen jakautumista magnetopausilla tarkastellaan energian muuntumisen pintatiheyden avulla ja rekonnektion kokonaismäärää rekonnektiotehon avulla. Magnetopausin ja pyrstön rekonnektiotehot ovat simulaatiossa samaa suuruusluokkaa. Tärkeimmät magnetopausin rekonnektiotehoa säätelevät parametrit ovat aurinkotuulen nopeus ja aurinkotuulen magneettikentän suunta. Magnetopausin rekonnektio puolestaan säätelee energian ja aineen pääsyä magnetosfääriin, joskaan magnetopausin läpäisevät vuot eivät ole aivan suoraan verrannollisia rekonnektiotehoon. Pyrstön rekonnektioteho sen sijaan on suoraan verrannollinen magnetopausilta tulevaan energiavuohon; pyrstörekonnektio Gumicsissa on siis ulkoista pakotetta seuraava passiivinen energian käsittelijä. Simulaation tuottama rekonnektio on realistinen magnetosfäärin globaalissa mittakaavassa tarkasteltuna, mutta satelliittihavainnot paljastavat rekonnektiosta simulaation erottelukykyä pienimittakaavaisempia piirteitä. Havaintopuolella tämän väitöstutkimuksen tärkein löytö on protonien diffuusioalueen rakenteeseen kuuluvien Hallin kenttien kääntyminen pyrstön virtalevyn aaltoilun mukana.

Electromagnetic signatures of lightning near the HF frequency band

The dissertation deals with remote narrowband measurements of the electromagnetic radiation emitted by lightning flashes. A lightning flash consists of a number of sub-processes. The return stroke, which transfers electrical charge from the thundercloud to to the ground, is electromagnetically an impulsive wideband process; that is, it emits radiation at most frequencies in the electromagnetic spectrum, but its duration is only some tens of microseconds. Before and after the return stroke, multiple sub-processes redistribute electrical charges within the thundercloud. These sub-processes can last for tens to hundreds of milliseconds, many orders of magnitude longer than the return stroke. Each sub-process causes radiation with specific time-domain characteristics, having maxima at different frequencies. Thus, if the radiation is measured at a single narrow frequency band, it is difficult to identify the sub-processes, and some sub-processes can be missed altogether. However, narrowband detectors are simple to design and miniaturize. In particular, near the High Frequency band (High Frequency, 3 MHz to 30 MHz), ordinary shortwave radios can, in principle, be used as detectors. This dissertation utilizes a prototype detector which is essentially a handheld AM radio receiver. Measurements were made in Scandinavia, and several independent data sources were used to identify lightning sub-processes, as well as the distance to each individual flash. It is shown that multiple sub-processes radiate strongly near the HF band. The return stroke usually radiates intensely, but it cannot be reliably identified from the time-domain signal alone. This means that a narrowband measurement is best used to characterize the energy of the radiation integrated over the whole flash, without attempting to identify individual processes. The dissertation analyzes the conditions under which this integrated energy can be used to estimate the distance to the flash. It is shown that flash-by-flash variations are large, but the integrated energy is very sensitive to changes in the distance, dropping as approximately the inverse cube root of the distance. Flashes can, in principle, be detected at distances of more than 100 km, but since the ground conductivity can vary, ranging accuracy drops dramatically at distances larger than 20 km. These limitations mean that individual flashes cannot be ranged accurately using a single narrowband detector, and the useful range is limited to 30 kilometers at the most. Nevertheless, simple statistical corrections are developed, which enable an accurate estimate of the distance to the closest edge of an active storm cell, as well as the approach speed. The results of the dissertation could therefore have practical applications in real-time short-range lightning detection and warning systems.

Features of the geomagnetic field during the Holocene and Proterozoic

The geomagnetic field is one of the most fundamental geophysical properties of the Earth and has significantly contributed to our understanding of the internal structure of the Earth and its evolution. Paleomagnetic and paleointensity data have been crucial in shaping concepts like continental drift, magnetic reversals, as well as estimating the time when the Earth's core and associated geodynamo processes begun. The work of this dissertation is based on reliable Proterozoic and Holocene geomagnetic field intensity data obtained from rocks and archeological artifacts. New archeomagnetic field intensity results are presented for Finland, Estonia, Bulgaria, Italy and Switzerland. The data were obtained using sophisticated laboratory setups as well as various reliability checks and corrections. Inter-laboratory comparisons between three laboratories (Helsinki, Sofia and Liverpool) were performed in order to check the reliability of different paleointensity methods. The new intensity results fill up considerable gaps in the master curves for each region investigated. In order to interpret the paleointensity data of the Holocene period, a novel and user-friendly database (GEOMAGIA50) was constructed. This provided a new tool to independently test the reliability of various techniques and materials used in paleointensity determinations. The results show that archeological artifacts, if well fired, are the most suitable materials. Also lavas yield reliable paleointensity results, although they appear more scattered. This study also shows that reliable estimates are obtained using the Thellier methodology (and its modifications) with reliability checks. Global paleointensity curves during Paleozoic and Proterozoic have several time gaps with few or no intensity data. To define the global intensity behavior of the Earth's magnetic field during these times new rock types (meteorite impact rocks) were investigated. Two case histories are presented. The Ilyinets (Ukraine) impact melt rocks yielded a reliable paleointensity value at 440 Ma (Silurian), whereas the results from Jänisjärvi impact melts (Russian Karelia, ca. 700 Ma) might be biased towards high intensity values because of non-ideal magnetic mineralogy. The features of the geomagnetic field at 1.1 Ga are not well defined due to problems related to reversal asymmetries observed in Keweenawan data of the Lake Superior region. In this work new paleomagnetic, paleosecular variation and paleointensity results are reported from coeval diabases from Central Arizona and help understanding the asymmetry. The results confirm the earlier preliminary observations that the asymmetry is larger in Arizona than in Lake Superior area. Two of the mechanisms proposed to explain the asymmetry remain plausible: the plate motion and the non-dipole influence.

Aerosol size distribution and its connection to cloud droplet number concentration

In order to predict the current state and future development of Earth s climate, detailed information on atmospheric aerosols and aerosol-cloud-interactions is required. Furthermore, these interactions need to be expressed in such a way that they can be represented in large-scale climate models. The largest uncertainties in the estimate of radiative forcing on the present day climate are related to the direct and indirect effects of aerosol. In this work aerosol properties were studied at Pallas and Utö in Finland, and at Mount Waliguan in Western China. Approximately two years of data from each site were analyzed. In addition to this, data from two intensive measurement campaigns at Pallas were used. The measurements at Mount Waliguan were the first long term aerosol particle number concentration and size distribution measurements conducted in this region. They revealed that the number concentration of aerosol particles at Mount Waliguan were much higher than those measured at similar altitudes in other parts of the world. The particles were concentrated in the Aitken size range indicating that they were produced within a couple of days prior to reaching the site, rather than being transported over thousands of kilometers. Aerosol partitioning between cloud droplets and cloud interstitial particles was studied at Pallas during the two measurement campaigns, First Pallas Cloud Experiment (First PaCE) and Second Pallas Cloud Experiment (Second PaCE). The method of using two differential mobility particle sizers (DMPS) to calculate the number concentration of activated particles was found to agree well with direct measurements of cloud droplet. Several parameters important in cloud droplet activation were found to depend strongly on the air mass history. The effects of these parameters partially cancelled out each other. Aerosol number-to-volume concentration ratio was studied at all three sites using data sets with long time-series. The ratio was found to vary more than in earlier studies, but less than either aerosol particle number concentration or volume concentration alone. Both air mass dependency and seasonal pattern were found at Pallas and Utö, but only seasonal pattern at Mount Waliguan. The number-to-volume concentration ratio was found to follow the seasonal temperature pattern well at all three sites. A new parameterization for partitioning between cloud droplets and cloud interstitial particles was developed. The parameterization uses aerosol particle number-to-volume concentration ratio and aerosol particle volume concentration as the only information on the aerosol number and size distribution. The new parameterization is computationally more efficient than the more detailed parameterizations currently in use, but the accuracy of the new parameterization was slightly lower. The new parameterization was also compared to directly observed cloud droplet number concentration data, and a good agreement was found.

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.

Two views on interstellar dust: near-infrared scattering and polarized thermal dust emission

Interstellar clouds are not featureless, but show quite complex internal structures of filaments and clumps when observed with high enough resolution. These structures have been generated by 1) turbulent motions driven mainly by supernovae, 2) magnetic fields working on the ions and, through neutral-ion collisions, on neutral gas as well, and 3) self-gravity pulling a dense clump together to form a new star. The study of the cloud structure gives us information on the relative importance of each of these mechanisms, and helps us to gain a better understanding of the details of the star formation process. Interstellar dust is often used as a tracer for the interstellar gas which forms the bulk of the interstellar matter. Some of the methods that are used to derive the column density are summarized in this thesis. A new method, which uses the scattered light to map the column density in large fields with high spatial resolution, is introduced. This thesis also takes a look at the grain alignment with respect to the magnetic fields. The aligned grains give rise to the polarization of starlight and dust emission, thus revealing the magnetic field. The alignment mechanisms have been debated for the last half century. The strongest candidate at present is the radiative torques mechanism. In the first four papers included in this thesis, the scattered light method of column density estimation is formulated, tested in simulations, and finally used to obtain a column density map from observations. They demonstrate that the scattered light method is a very useful and reliable tool in column density estimation, and is able to provide higher resolution than the near-infrared color excess method. These two methods are complementary. The derived column density maps are also used to gain information on the dust emissivity within the observed cloud. The two final papers present simulations of polarized thermal dust emission assuming that the alignment happens by the radiative torques mechanism. We show that the radiative torques can explain the observed decline of the polarization degree towards dense cores. Furthermore, the results indicate that the dense cores themselves might not contribute significantly to the polarized signal, and hence one needs to be careful when interpreting the observations and deriving the magnetic field.