90.00.00 GEOPHYSICS, ASTRONOMY, AND ASTROPHYSICS (for more detailed headings, see the Geophysics Appendix)
Refine
Has Fulltext
- yes (2)
Is part of the Bibliography
- yes (2)
Year of publication
- 2022 (2) (remove)
Document Type
- Doctoral Thesis (2)
Keywords
- Bodenparameter (1)
- Datenanalyse (1)
- Klimamodell (1)
- Klimamodellierung (1)
- Modellierung (1)
- REMO (1)
- Topographie (1)
- Topographiedaten (1)
- Vergleich verschiedener Modellparameterisierungen (1)
- blazar (1)
Institute
Sonstige beteiligte Institutionen
Context. In active galaxies, matter is accreted onto super massive black holes (SMBH). This accretion process causes a region roughly the size of our solar system to outshine the entire host galaxy, forming an active galactic nucleus (AGN). In some of these active galaxies, highly relativistic particle jets are formed parallel to the rotation axis of the super massive black hole. A fraction of these sources is observed under a small inclination angle between the pointing direction of the jet and the observing line of sight. These sources are called blazars. Due to the small inclination angle and the highly relativistic speeds of the particles in the jet, beaming effects occur in the radiation of these particles. Blazars can be subdivided into the high luminosity flat spectrum radio quasars (FSRQs) and the low luminosity BL Lacertae objects (BL Lacs). As all AGN, blazars are broadband emitters and therefore observable from the longest wavelengths in the radio regime to the shortest wavelengths in the gamma-ray regime. In this thesis I will analyze blazars at these two extremes with respect to their parsec-scale properties in the radio and their time evolution properties in gamma-ray flux.
Method. In the radio regime the technique of very long baseline interferometry (VLBI) can be used in order to spatially resolve the synchrotron radiation coming from those objects down to sub-parsec scales. This information can be used to observe the time evolution of the structure of such sources. This is done in large monitoring programs such as the MOJAVE (15 GHz) and the Boston University blazar monitoring program (43 GHz). In this thesis I utilize data of 28 sources from these monitoring programs spanning 10 years of observation from 2003 to 2013, resulting in over 1800 observed epochs, to study the brightness temperature and diameter gradients of these jets. I conduct a search for systematic geometry transitions in the radio jets. The synchrotron cooling time in the radio core of the jets is used to determine the magnetic field strength in the radio core. Considering the jet geometry, these magnetic field strengths are scaled to the ergosphere of the SMBH in order to obtain the distance of the radio core to the SMBH.
In the gamma-regime these blazars cannot be spatially resolved. Due to this, it is hard to put strong constrains onto where the gamma-ray emitting region is. Blazars have shown to be variable at high energies on time scales down to minutes. The nature of this variability can be studied in order to put constrains on the particle acceleration mechanism and possibly the region and size of the gamma-ray emitting region. The variability of blazars in the energy range between 0.1 GeV and 1 GeV can for example be observed with the pair-conversion telescope on board the Fermi satellite. I use 10 years of data from the Fermi-LAT (LAT: Large Area Telescope) satellite in order to study the variability of a large sample of blazars (300-800, depending on the used significance filters for data points). I quantify this variability with the Ornstein-Uhlenbeck (OU) parameters and the power spectral density (PSD) slopes. The same procedure is applied to 20 light curves available for the radio sample.
Results. The diameter evolution along the jet axis of the radio sources suggests, that FSRQs feature flatter gradients than BL Lacs. Fitting these gradients, it is revealed that BL Lacs are systematically better described by a simple single power law than FSRQs. I found 9 sources with a strongly constrained geometry transition. The sources are 0219+421, 0336-019, 0415+379, 0528+134, 0836+710, 1101+384, 1156+295, 1253-055 and 2200+420. In all of these sources, the geometry transition regions are further out in the jet than the Bondi sphere. The magnetic field strengths of BL Lacs is systematically larger than that of FSRQs. However the scaling of these fields suggest that the radio cores of BL Lac objects are closer to the SMBHs than the radio cores of FSRQs. Analyzing the variability of Fermi-LAT light curves yields consistent results for all samples. FSRQs show systematically steeper PSD slopes and feature OU parameters more favorable to strong variability than BL Lacs. The Fermi-LAT light curves of the sub-sample of radio jets, suggest an anticorrelation between the jet complexity from the radio observations and the OU-parameters as well as the PSD slopes from the gamma-ray observations.
Conclusion.
The flatter diameter gradients of FSRQs suggest that these sources are more collimated further down the jet than BL Lacs. The systematically better description of the diameter and brightness temperature gradient by a single power law of BL Lacs, suggest that FSRQs are more complex with respect to the diameter evolution along the jet and the surface brightness distribution than BL Lac objects. FSRQs often feature regions where recollimation can occur in distinct knots within the jets. For the sources where a geometry transition could be constrained, the Bondi radius, being systematically smaller than the position of the transition region along the jet axis, suggest that changing pressure gradients are not the sole cause for these systematic geometry transitions. Nevertheless they may be responsible for recollimation regions, found typically downstream the jet, beyond the Bondi radius and the transition zone. The difference in the distance of the radio cores between FSRQs and BL Lacs is most likely due to the combination of differences in SMBH masses and systematically smaller jet powers in BL Lacs. The variability in energy ranges above 100 MeV and above 1 GeV-regime suggest that many light curves of BL Lac objects are more likely to be white noise while the PSD slopes and the OU parameters of FSRQ gamma-ray light curves favor stronger variability on larger time scales with respect to the time binning of the analyzed light curve. Although the anticorrelation of the jet complexity acquired from the radio observations and the PSD slopes and OU parameters from the gamma-observations suggest that more complex sources favor OU parameters and PSD slopes resulting in more variability (not white noise) it is beyond the scope of this thesis to pinpoint whether this correlation results from causation. The question whether a complex jet causes more gamma-ray variability or more gamma-ray variability causes more complex jets cannot be answered at this point. Nevertheless the computed correlation measures suggest that this dependence is most likely not linear and therefore an indication that these effects might even interact.
Das Ziel dieser Arbeit war neue Eingangsdaten für die Landoberflächenbeschreibung des regionalen Klimamodells REMO zu finden und ins Modell zu integrieren, um die Vorhersagequalität des Modells zu verbessern. Die neuen Daten wurden so in das Modell eingebaut, dass die bisherigen Daten weiterhin als Option verfügbar sind. Dadurch kann überprüft werden, ob und in welchem Umfang sich die von jedem Klimamodell benötigten Rahmendaten auf Modellergebnisse auswirken. Im Zuge der Arbeit wurden viele unterschiedliche Daten und Methoden zur Generierung neuer Parameter miteinander verglichen, denn neben dem Ersetzen der konstanten Eingangswerte für verschiedene Oberflächenparameter und den damit verbundenen Änderungen wurden als zusätzliche Verbesserung auch Veränderungen an der Parametrisierung des Bodens speziell in Hinblick auf die Bodentemperaturen in REMO vorgenommen. Im Rahmen dieser Arbeit wurden die durch die verschiedenen Änderungen ausgelösten Auswirkungen für das CORDEX-Gebiet EUR-44 mit einer Auflösung von ca. 50km und für das in dem darin eingebetteten neu definierten Deutschlandgebiet GER-11 mit einer Auflösung von ca. 12km getestet sowie alle Änderungen anhand von verschiedenen Beobachtungsdatensätzen validiert.
Die vorgenommenen Arbeiten gliederten sich in drei Hauptteile. Der erste Teil bestand in dem vom eigentlichen Klimamodell unabhängigen Vergleich der verschiedenen Eingangsdaten auf unterschiedlichen Auflösungen und deren Performanz in allen Teilen der Erde, wobei ein besonderer Fokus auf der Qualität in den späteren Modellgebieten lag. Unter Berücksichtigung der Faktoren, wie einer globalen Verfügbarkeit der Daten, einer verbesserten räumlichen Auflösung und einer kostenlosen Nutzung der Daten sowie verschiedener Validationsergebnissen von anderen Studien, wurden in dieser Arbeit vier neue Topographiedatensätze (SRTM, ALOS, TANDEM und ASTER) und drei neue Bodendatensätze (FAOn, Soilgrid und HWSD) für die Verwendung im Präprozess von REMO aufbereitet und miteinander sowie mit den bisher in REMO verwendeten Daten verglichen. Auf Grundlage dieser Vergleichsstudien schieden bei den Topographiedaten die verwendeten Datensatz-Versionen von SRTM, ALOS und TANDEM für die in dieser Arbeit durchgeführten REMO-Läufe aus. Bei den neuen Bodendatensätzen wurde ausgenutzt, dass diese verschiedenen Bodeneigenschaften für unterschiedliche Tiefen als Karten zur Verfügung stellen. In REMO wurden bisher alle benötigten Bodenparameter abhängig von fünf verschiedenen Bodentexturklassen und einer zusätzlichen Torfklasse ausgewiesen und als konstant über die gesamte Modellbodensäule (bis ca. 10m) angenommen. Im zweiten Teil wurden auf Basis der im ersten Teil ausgewählten neuen Datensätze und den neu verfügbaren Bodenvariablen verschiedene Sensitivitätsstudien über das Beispieljahr 2000 durchgeführt. Dabei wurden verschiedene neue Parametrisierungen für die bisher aus der Textur abgeleiteten Bodenvariablen und die Parametrisierung von weiteren hydrologischen und thermalen Bodeneigenschaften verglichen. Ferner wurde aufgrund der neuen nicht über die Tiefe konstanten Bodeneigenschaften eine neue numerische Methode zur Berechnung der Bodentemperaturen der fünf Schichten in REMO getestet, welche wiederum andere Anpassungen erforderte. Der Test und die Auswahl der verschiedenen Datensatz- und Parametrisierungsversionen auf die Modellperformanz wurde in drei Experimentpläne unterteilt. Im ersten Plan wurden die Auswirkungen der ausgewählten Topographie- und Bodendatensätze überprüft. Der zweite Plan behandelte die Unterschiede der verschiedenen Parametrisierungsarten der Bodenvariablen hinsichtlich der verwendeten Variablen zur Berechnung der Bodeneigenschaften, der über die Tiefe variablen oder konstanten Eigenschaften und der verwendeten Berechnungsmethode der Bodentemperaturänderungen. Durch die Erkenntnisse aus diesen beiden Experimentplänen, die für beide Untersuchungsgebiete durchgeführt wurden, ergaben sich im dritten Plan weitere Parametrisierungsänderungen. Alle Änderungen dieses dritten Experimentplans wurden sukzessiv getestet, sodass der paarweise Vergleich von zwei aufeinanderfolgenden Modellläufen die Auswirkungen der Neuerung im jeweils zweiten Lauf widerspiegelt. Der letzte Teil der Arbeit bestand aus der Analyse von fünf längeren Modellläufen (2000-2018), die zur Überprüfung der Ergebnisse aus den Sensitivitätsstudien sowie zur Einschätzung der Performanz in weiteren teilweise extremen atmosphärischen Bedingungen durchgeführt wurden. Hierfür wurden die bisherige Modellversion von REMO (id01) für die beiden Untersuchungsgebiete EUR-44 und GER-11 als Referenzläufe, zwei aufgrund der Vergleichsergebnisse von Experimentplan 3 selektierte Modellversionen (id06 und id15a für GER-11) sowie die finale Version (id18a für GER-11), die alle vorgenommenen Änderungen dieser Arbeit enthält, ausgewählt.
Es stellte sich heraus, dass sowohl die neuen Topographiedaten als auch die neuen Bodendaten große Differenzen zu den bisherigen Daten in REMO haben. Zudem änderten sich die von diesen konstanten Eingangsdaten abgeleiteten Hilfsvariablen je nach verwendeter Parametrisierung sehr deutlich. Dies war besonders gut anhand der Bodenparameter zu erkennen. Sowohl die räumliche Verteilung als auch der Wertebereich der verschiedenen Modellversionen unterschieden sich stark. Eine Einschätzung der Qualität der resultierenden Parameter wurde jedoch dadurch erschwert, dass auch die verschiedenen zur Validierung herangezogenen Bodendatensätze für diese Parameter deutlich voneinander abweichen. Die finale Modellversion id18a ähnelte trotz der umfassenden Änderungen in den meisten Variablen den Ergebnissen der bisherigen REMO-Version. Je nach zeitlicher und räumlicher Aggregation sowie unterschiedlichen Regionen und Jahreszeiten wurden leichte Verbesserungen, aber auch leichte Verschlechterungen im Vergleich zu den klimatologischen Validationsdaten festgestellt. Größere Veränderungen im Vergleich zur bisherigen Modellversion konnten in den tieferen Bodenschichten aufgezeigt werden, welche allerdings aufgrund von fehlenden Validationsdaten nicht beurteilt werden konnten. Für alle 2m-Temperaturen konnte eine tendenzielle leichte Erwärmung im Vergleich zum bisherigen Modelllauf beobachtet werden, was sich einerseits negativ auf die ohnehin durchschnittlich zu hohe Minimumtemperatur, aber andererseits positiv auf die bisher zu niedrige Maximumtemperatur des Modells in den betrachteten Gebieten auswirkte. Im Niederschlagssignal und in den 10m-Windvariablen konnten keine signifikanten Änderungen nachgewiesen werden, obwohl die neue Topographie an manchen Stellen im Modellgebiet deutlich von der bisherigen abweicht. Des Weiteren variierte das Ranking der verschiedenen Modellversionen jeweils nach dem angewendeten Qualitätsindex.
Um diese Ergebnisse besser einordnen zu können, muss berücksichtigt werden, dass die neuen Daten für Modellgebiete mit 50 bzw. 12km räumlicher Auflösung und der damit verbundenen hydrostatischen Modellversion getestet wurden. Zudem sind vor allem in Fall der Topographie die bisher enthaltenen GTOPO-Daten (1km Auflösung) für die Aggregation auf diese gröbere Modellauflösung geeignet. Die bisherigen Bodendaten stoßen jedoch mit 50km Auflösung bereits an ihre Grenzen. Zusätzlich ist zu beachten, dass nicht nur die Mittelwerte dieser Daten, sondern auch deren Subgrid-Variabilität als Variablen im Modell für verschiedene Parametrisierungen verwendet werden. Daher ist es essentiell, dass die Eingangsdaten eine deutlich höhere Auflösung bereitstellen als die zur Modellierung definierte Auflösung. Für lokale Klimasimulationen mit Auflösungen im niedrigen Kilometerbereich spielen auch die Vertikalbewegungen (nicht-hydrostatische Modellversion) eine wichtige Rolle, die stark von der Topographie sowie deren horizontaler und vertikaler Änderungsrate beeinflusst werden, was die in dieser Arbeit eingebauten wesentlich höher aufgelösten Daten für die zukünftige Weiterentwicklung von REMO wertvoll machen kann.