Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Biomass Gasification - Process analysis and dimensioning aspects for downdraft units and gas cleaning lines

In such territories where food production is mostly scattered in several small / medium size or even domestic farms, a lot of heterogeneous residues are produced yearly, since farmers usually carry out different activities in their properties. The amount and composition of farm residues, therefore, widely change during year, according to the single production process periodically achieved. Coupling high efficiency micro-cogeneration energy units with easy handling biomass conversion equipments, suitable to treat different materials, would provide many important advantages to the farmers and to the community as well, so that the increase in feedstock flexibility of gasification units is nowadays seen as a further paramount step towards their wide spreading in rural areas and as a real necessity for their utilization at small scale. Two main research topics were thought to be of main concern at this purpose, and they were therefore discussed in this work: the investigation of fuels properties impact on gasification process development and the technical feasibility of small scale gasification units integration with cogeneration systems. According to these two main aspects, the present work was thus divided in two main parts. The first one is focused on the biomass gasification process, that was investigated in its theoretical aspects and then analytically modelled in order to simulate thermo-chemical conversion of different biomass fuels, such as wood (park waste wood and softwood), wheat straw, sewage sludge and refuse derived fuels. The main idea is to correlate the results of reactor design procedures with the physical properties of biomasses and the corresponding working conditions of gasifiers (temperature profile, above all), in order to point out the main differences which prevent the use of the same conversion unit for different materials. At this scope, a gasification kinetic free model was initially developed in Excel sheets, considering different values of air to biomass ratio and the downdraft gasification technology as particular examined application. The differences in syngas production and working conditions (process temperatures, above all) among the considered fuels were tried to be connected to some biomass properties, such elementary composition, ash and water contents. The novelty of this analytical approach was the use of kinetic constants ratio in order to determine oxygen distribution among the different oxidation reactions (regarding volatile matter only) while equilibrium of water gas shift reaction was considered in gasification zone, by which the energy and mass balances involved in the process algorithm were linked together, as well. Moreover, the main advantage of this analytical tool is the easiness by which the input data corresponding to the particular biomass materials can be inserted into the model, so that a rapid evaluation on their own thermo-chemical conversion properties is possible to be obtained, mainly based on their chemical composition A good conformity of the model results with the other literature and experimental data was detected for almost all the considered materials (except for refuse derived fuels, because of their unfitting chemical composition with the model assumptions). Successively, a dimensioning procedure for open core downdraft gasifiers was set up, by the analysis on the fundamental thermo-physical and thermo-chemical mechanisms which are supposed to regulate the main solid conversion steps involved in the gasification process. Gasification units were schematically subdivided in four reaction zones, respectively corresponding to biomass heating, solids drying, pyrolysis and char gasification processes, and the time required for the full development of each of these steps was correlated to the kinetics rates (for pyrolysis and char gasification processes only) and to the heat and mass transfer phenomena from gas to solid phase. On the basis of this analysis and according to the kinetic free model results and biomass physical properties (particles size, above all) it was achieved that for all the considered materials char gasification step is kinetically limited and therefore temperature is the main working parameter controlling this step. Solids drying is mainly regulated by heat transfer from bulk gas to the inner layers of particles and the corresponding time especially depends on particle size. Biomass heating is almost totally achieved by the radiative heat transfer from the hot walls of reactor to the bed of material. For pyrolysis, instead, working temperature, particles size and the same nature of biomass (through its own pyrolysis heat) have all comparable weights on the process development, so that the corresponding time can be differently depending on one of these factors according to the particular fuel is gasified and the particular conditions are established inside the gasifier. The same analysis also led to the estimation of reaction zone volumes for each biomass fuel, so as a comparison among the dimensions of the differently fed gasification units was finally accomplished. Each biomass material showed a different volumes distribution, so that any dimensioned gasification unit does not seem to be suitable for more than one biomass species. Nevertheless, since reactors diameters were found out quite similar for all the examined materials, it could be envisaged to design a single units for all of them by adopting the largest diameter and by combining together the maximum heights of each reaction zone, as they were calculated for the different biomasses. A total height of gasifier as around 2400mm would be obtained in this case. Besides, by arranging air injecting nozzles at different levels along the reactor, gasification zone could be properly set up according to the particular material is in turn gasified. Finally, since gasification and pyrolysis times were found to considerably change according to even short temperature variations, it could be also envisaged to regulate air feeding rate for each gasified material (which process temperatures depend on), so as the available reactor volumes would be suitable for the complete development of solid conversion in each case, without even changing fluid dynamics behaviour of the unit as well as air/biomass ratio in noticeable measure. The second part of this work dealt with the gas cleaning systems to be adopted downstream the gasifiers in order to run high efficiency CHP units (i.e. internal engines and micro-turbines). Especially in the case multi–fuel gasifiers are assumed to be used, weightier gas cleaning lines need to be envisaged in order to reach the standard gas quality degree required to fuel cogeneration units. Indeed, as the more heterogeneous feed to the gasification unit, several contaminant species can simultaneously be present in the exit gas stream and, as a consequence, suitable gas cleaning systems have to be designed. In this work, an overall study on gas cleaning lines assessment is carried out. Differently from the other research efforts carried out in the same field, the main scope is to define general arrangements for gas cleaning lines suitable to remove several contaminants from the gas stream, independently on the feedstock material and the energy plant size The gas contaminant species taken into account in this analysis were: particulate, tars, sulphur (in H2S form), alkali metals, nitrogen (in NH3 form) and acid gases (in HCl form). For each of these species, alternative cleaning devices were designed according to three different plant sizes, respectively corresponding with 8Nm3/h, 125Nm3/h and 350Nm3/h gas flows. Their performances were examined on the basis of their optimal working conditions (efficiency, temperature and pressure drops, above all) and their own consumption of energy and materials. Successively, the designed units were combined together in different overall gas cleaning line arrangements, paths, by following some technical constraints which were mainly determined from the same performance analysis on the cleaning units and from the presumable synergic effects by contaminants on the right working of some of them (filters clogging, catalysts deactivation, etc.). One of the main issues to be stated in paths design accomplishment was the tars removal from the gas stream, preventing filters plugging and/or line pipes clogging At this scope, a catalytic tars cracking unit was envisaged as the only solution to be adopted, and, therefore, a catalytic material which is able to work at relatively low temperatures was chosen. Nevertheless, a rapid drop in tars cracking efficiency was also estimated for this same material, so that an high frequency of catalysts regeneration and a consequent relevant air consumption for this operation were calculated in all of the cases. Other difficulties had to be overcome in the abatement of alkali metals, which condense at temperatures lower than tars, but they also need to be removed in the first sections of gas cleaning line in order to avoid corrosion of materials. In this case a dry scrubber technology was envisaged, by using the same fine particles filter units and by choosing for them corrosion resistant materials, like ceramic ones. Besides these two solutions which seem to be unavoidable in gas cleaning line design, high temperature gas cleaning lines were not possible to be achieved for the two larger plant sizes, as well. Indeed, as the use of temperature control devices was precluded in the adopted design procedure, ammonia partial oxidation units (as the only considered methods for the abatement of ammonia at high temperature) were not suitable for the large scale units, because of the high increase of reactors temperature by the exothermic reactions involved in the process. In spite of these limitations, yet, overall arrangements for each considered plant size were finally designed, so that the possibility to clean the gas up to the required standard degree was technically demonstrated, even in the case several contaminants are simultaneously present in the gas stream. Moreover, all the possible paths defined for the different plant sizes were compared each others on the basis of some defined operational parameters, among which total pressure drops, total energy losses, number of units and secondary materials consumption. On the basis of this analysis, dry gas cleaning methods proved preferable to the ones including water scrubber technology in al of the cases, especially because of the high water consumption provided by water scrubber units in ammonia adsorption process. This result is yet connected to the possibility to use activated carbon units for ammonia removal and Nahcolite adsorber for chloride acid. The very high efficiency of this latter material is also remarkable. Finally, as an estimation of the overall energy loss pertaining the gas cleaning process, the total enthalpy losses estimated for the three plant sizes were compared with the respective gas streams energy contents, these latter obtained on the basis of low heating value of gas only. This overall study on gas cleaning systems is thus proposed as an analytical tool by which different gas cleaning line configurations can be evaluated, according to the particular practical application they are adopted for and the size of cogeneration unit they are connected to.

Influence of nitrogen and soil physical characteristics on belowground carbon flux dynamics of woody plants

At ecosystem level soil respiration (Rs) represents the largest carbon (C) flux after gross primary productivity, being mainly generated by root respiration (autotrophic respiration, Ra) and soil microbial respiration (heterotrophic respiration, Rh). In the case of terrestrial ecosystems, soils contain the largest C-pool, storing twice the amount of C contained in plant biomass. Soil organic matter (SOM), representing the main C storage in soil, is decomposed by soil microbial community. This process produces CO2 which is mainly released as Rh. It is thus relevant to understand how microbial activity is influenced by environmental factors like soil temperature, soil moisture and nutrient availability, since part of the CO2 produced by Rh, directly increases atmospheric CO2 concentration and therefore affects the phenomenon of climate change. Among terrestrial ecosystems, agricultural fields have traditionally been considered as sources of atmospheric CO2. In agricultural ecosystems, in particular apple orchards, I identified the role of root density, soil temperature, soil moisture and nitrogen (N) availability on Rs and on its two components, Ra and Rh. To do so I applied different techniques to separate Rs in its two components, the ”regression technique” and the “trenching technique”. I also studied the response of Ra to different levels of N availability, distributed either in a uniform or localized way, in the case of Populus tremuloides trees. The results showed that Rs is mainly driven by soil temperature, to which it is positively correlated, that high levels of soil moisture have inhibiting effects, and that N has a negligible influence on total Rs, as well as on Ra. Further I found a negative response of Rh to high N availability, suggesting that microbial decomposition processes in the soil are inhibited by the presence of N. The contribution of Ra to Rs was of 37% on average.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

Transcriptional dynamics of the human DMD locus

The human DMD locus encodes dystrophin protein. Absence or reduced levels of dystrophin (DMD or BMD phenotype, respectively) lead to progressive muscle wasting. Little is known about the complex coordination of dystrophin expression and its transcriptional regulation is a field of intense interest. In this work we found that DMD locus harbours multiple long non coding RNAs which orchestrate and control transcription of muscle dystrophin mRNA isoforms. These lncRNAs are tissue-specific and highly expressed during myogenesis, suggesting a possible role in tissue-specific expression of DMD gene isoforms. Their forced ectopic expression in human muscle and neuronal cells leads to a specific and negative regulation of endogenous dystrophin full lenght isoforms. An intriguing aspect regarding the transcription of the DMD locus is the gene size (2.4Mb). The mechanism that ensures the complete synthesis of the primary transcript and the coordinated splicing of 79 exons is still completely unknown. By ChIP-on-chip analyses, we discovered novel regions never been involved before in the transcription regulation of the DMD locus. Specifically, we observed enrichments for Pol II, P-Ser2, P-Ser5, Ac-H3 and 2Me-H3K4 in an intronic region of 3Kb (approximately 21Kb) downstream of the end of DMD exon 52 and in a region of 4Kb spanning the DMD exon 62. Interestingly, this latter region and the TSS of Dp71 are strongly marked by 3Me-H3K36, an histone modification associated with the regulation of splicing process. Furthermore, we also observed strong presence of open chromatin marks (Ac-H3 and 2Me-H3K4) around intron 34 and the exon 45 without presence of RNA pol II. We speculate that these two regions may exert an enhancer-like function on Dp427m promoter, although further investigations are necessary. Finally, we investigated the nuclear-cytoplasmic compartmentalization of the muscular dystrophin mRNA and, specifically, we verified whether the exon skipping therapy could influence its cellular distribution.

Identity dynamics and the emergence of new organizational arrangements: A multi-level study

Organizational and institutional scholars have advocated the need to examine how processes originating at an individual level can change organizations or even create new organizational arrangements able to affect institutional dynamics (Chreim et al., 2007; Powell & Colyvas, 2008; Smets et al., 2012). Conversely, research on identity work has mainly investigated the different ways individuals can modify the boundaries of their work in actual occupations, thus paying particular attention to ‘internal’ self-crafting (e.g. Wrzesniewski & Dutton, 2001). Drawing from literatures on possible and alternative self and on positive organizational scholarship (e.g., Obodaru, 2012; Roberts & Dutton, 2009), my argument is that individuals’ identity work can go well beyond the boundaries of internal self-crafting to the creation of new organizational arrangements. In this contribution I analyze, through multiple case studies, healthcare professionals who spontaneously participated in the creation of new organizational arrangements, namely health structures called Community Hospitals. The contribution develops this form of identity work by building a grounded model. My findings disclose the process that leads from the search for the enactment of different self-concepts to positive identities, through the creation of a new organizational arrangement. I contend that this is a particularly complex form of collective identity work because it requires, to be successful, concerted actions of several internal, external and institutional actors, and it also requires balanced tensions that – at the same time - enable individuals’ aspirations and organizational equilibrium. I name this process organizational collective crafting. Moreover I inquire the role of context in supporting the triggering power of those unrealized selves. I contribute to the comprehension of the consequences of self-comparisons, organizational identity variance, and positive identity. The study bears important insights on how identity work originating from individuals can influence organizational outcomes and larger social systems.

Simulation of photoinduced processes in organic chromophores: mapping photochemically relevant states for accurate dynamics

The simulation of ultrafast photoinduced processes is a fundamental step towards the understanding of the underlying molecular mechanism and interpretation/prediction of experimental data. Performing a computer simulation of a complex photoinduced process is only possible introducing some approximations but, in order to obtain reliable results, the need to reduce the complexity must balance with the accuracy of the model, which should include all the relevant degrees of freedom and a quantitatively correct description of the electronic states involved in the process. This work presents new computational protocols and strategies for the parameterisation of accurate models for photochemical/photophysical processes based on state-of-the-art multiconfigurational wavefunction-based methods. The required ingredients for a dynamics simulation include potential energy surfaces (PESs) as well as electronic state couplings, which must be mapped across the wide range of geometries visited during the wavepacket/trajectory propagation. The developed procedures allow to obtain solid and extended databases reducing as much as possible the computational cost, thanks to, e.g., specific tuning of the level of theory for different PES regions and/or direct calculation of only the needed components of vectorial quantities (like gradients or nonadiabatic couplings). The presented approaches were applied to three case studies (azobenzene, pyrene, visual rhodopsin), all requiring an accurate parameterisation but for different reasons. The resulting models and simulations allowed to elucidate the mechanism and time scale of the internal conversion, reproducing or even predicting new transient experiments. The general applicability of the developed protocols to systems with different peculiarities and the possibility to parameterise different types of dynamics on an equal footing (classical vs purely quantum) prove that the developed procedures are flexible enough to be tailored for each specific system, and pave the way for exact quantum dynamics with multiple degrees of freedom.

Smart Farming in Italian agriculture: essays on adoption and diffusion dynamics shaping the agricultural digital transition

Smart Farming Technologies (SFT) is a term used to define the set of digital technologies able not only to control and manage the farm system, but also to connect it to the many disruptive digital applications posed at multiple links along the value chain. The adoption of SFT has been so far limited, with significant differences at country-levels and among different types of farms and farmers. The objective of this thesis is to analyze what factors contributes to shape the agricultural digital transition and to assess its potential impacts in the Italian agri-food system. Specifically, this overall research objective is approached under three different perspectives. Firstly, we carry out a review of the literature that focuses on the determinants of adoption of farm-level Management Information Systems (MIS), namely the most adopted smart farming solutions in Italy. Secondly, we run an empirical analysis on what factors are currently shaping the adoption of SFT in Italy. In doing so, we focus on the multi-process and multi-faceted aspects of the adoption, by overcoming the one-off binary approach often used to study adoption decisions. Finally, we adopt a forward-looking perspective to investigate what the socio-ethical implications of a diffused use of SFT might be. On the one hand, our results indicate that bigger, more structured farms with higher levels of commercial integration along the agri-food supply chain are those more likely to be early adopters. On the other hand, they highlight the need for the institutional and organizational environment around farms to more effectively support farmers in the digital transition. Moreover, the role of several other actors and actions are discussed and analyzed, by highlighting the key role of specific agri-food stakeholders and ad-hoc policies, with the aim to propose a clearer path towards an efficient, fair and inclusive digitalization of the agrifood sector.

Chemical looping reforming from methane to syngas: investigation on CeO2 oxygen carrier lifetime and optimal process times

Emissions of CO2 are constantly growing since the beginning of industrial era. Interruption of the production of major emitters sectors (energy and agriculture) is not a viable way and reducing all the emission through carbon capture and storage (CCS) is not economically viable and little publicly accepted, therefore, it becomes fundamentals to take actions like retrofitting already developed infrastructure employing cleanest resources, modify the actual processes limiting the emissions, and reduce the emissions already present through direct air capture. The present thesis will deeply discuss the aspects mentioned in regard to syngas and hydrogen production since they have a central role in the market of energy and chemicals. Among the strategies discussed, greater emphasis is given to the application of looping technologies and to direct air capture processes, as they have been the main point of this work. Particularly, chemical looping methane reforming to syngas was studied with Aspen Plus thermodynamic simulations, thermogravimetric analysis characterization (TGA) and testing in a fixed bed reactor. The process was studied cyclically exploiting the redox properties of a Ce-based oxide oxygen carrier synthetized with a simple forming procedure. The two steps of the looping cycles were studied isothermally at 900 °C and 950° C with a mixture of 10 %CH4 in N2 and of 3% O2 in N2, for carrier reduction and oxidation, respectively. During the stay abroad, in collaboration with the EHT of Zurich, a CO2 capture process in presence of amine solid sorbents was investigated, studying the difference in the performance achievable with the use of contactors of different geometry. The process was studied at two concentrations (382 ppm CO2 in N2 and 5.62% CO2 in N2) and at different flow rates, to understand the dynamics of the adsorption process and to define the mass transfer limiting step.

Navigation design and flight dynamics operations of the ArgoMoon cis-lunar CubeSat.

On November 16, 2022, the NASA’s Space Launch System (SLS) has been launched for the first time in the context of Artemis-1 mission where, together with the Orion Multi-Purpose Crew Vehicle, a set of 10 CubeSats have been delivered into a translunar trajectory. Among the small satellites deployed during Artemis-1 there is ArgoMoon, a 6U CubeSat built by the Italian company Argotec and coordinated by Italian Space Agency (ASI). The primary goal of ArgoMoon is to capture images of the Interim Cryogenic Propulsion Stage. The ArgoMoon trajectory has been designed as a highly elliptical geocentric orbit, with several encounters with the Moon. In order to successfully fly ArgoMoon along the designed cis-lunar trajectory, a ground-based navigation system has been developed exploiting the guidance techniques also used for regular deep space missions. The navigation process is subdivided into Orbit Determi- nation (OD) and a Flight Path Control (FPC), and it is designed to follow the reference trajectory, prevent impacts with the Earth and the Moon, intensively test the navigation techniques, and guarantee the spacecraft disposal at the end of the mission. The work done in this thesis has accomplished the navigation of ArgoMoon, covering all aspects of the project life, from pre-launch design and analysis to actual operations. Firstly, the designed navigation process and the pre-mission assessment of its performance will be presented. Then, the results of the ArgoMoon navigation operations performed after the launch in November 2022 will be described in detail by discussing the main encountered challenges and the adopted solutions. The results of the operations confirmed the robustness of the designed navigation which allowed to accurately estimate the trajectory of ArgoMoon despite a series of complex events.