Cluster observations and simulations of He+ EMIC triggered emissions

NASA Astrophysics Data System (ADS)

Grison, B.; Shoji, M.; Santolik, O.; Omura, Y.

2012-12-01

EMIC triggered emissions have been reported in the inner magnetosphere at the edge of the plasmapause nightside [Pickett et al., 2010]. The generation mechanism proposed by Omura et al. [2010] is very similar to the one of the whistler chorus emissions and simulation results agree with observations and theory [Shoji et Omura, 2011]. The main characteristics of these emissions generated in the magnetic equatorial plane region are a frequency with time dispersion and a high level of coherence. The start frequency of previously mentioned observations is above half of the proton gyrofrequency. It means that the emissions are generated on the proton branch. On the He+ branch, generation of triggered emissions, in the same region, requests more energetic protons and the triggering process starts below the He+ gyrofrequency. It makes their identification in Cluster data rather difficult. Recent simulation results confirm the possibility of EMIC triggered emission on the He+ branch. In the present contribution we propose to compare a Cluster event to simulation results in order to investigate the possibility to identify observations to a He+ triggered emission. The impact of the observed waves on particle precipitation is also investigated.

A NEO population generation and observation simulation software tool

NASA Astrophysics Data System (ADS)

Müller, Sven; Gelhaus, Johannes; Hahn, Gerhard; Franco, Raffaella

One of the main targets of ESA's Space Situational Awareness (SSA) program is to build a wide knowledge base about objects that can potentially harm Earth (Near-Earth Objects, NEOs). An important part of this effort is to create the Small Bodies Data Centre (SBDC) which is going to aggregate measurement data from a fully-integrated NEO observation sensor network. Until this network is developed, artificial NEO measurement data is needed in order to validate SBDC algorithms. Moreover, to establish a functioning NEO observation sensor network, it has to be determined where to place sensors, what technical requirements have to be met in order to be able to detect NEOs and which observation strategies work the best. Because of this, a sensor simulation software was needed. This paper presents a software tool which allows users to create and analyse NEO populations and to simulate and analyse population observations. It is a console program written in Fortran and comes with a Graphical User Interface (GUI) written in Java and C. The tool can be distinguished into the components ``Population Generator'' and ``Observation Simulator''. The Population Generator component is responsible for generating and analysing a NEO population. Users can choose between creating fictitious (random) and synthetic populations. The latter are based on one of two models describing the orbital and size distribution of observed NEOs: The existing socalled ``Bottke Model'' (Bottke et al. 2000, 2002) and the new ``Granvik Model'' (Granvik et al. 2014, in preparation) which has been developed in parallel to the tool. Generated populations can be analysed by defining 2D, 3D and scatter plots using various NEO attributes. As a result, the tool creates the appropiate files for the plotting tool ``gnuplot''. The tool's Observation Simulator component yields the Observation Simulation and Observation Analysis functions. Users can define sensor systems using ground- or space-based locations as well as

Observing System Simulation Experiments

NASA Technical Reports Server (NTRS)

Prive, Nikki

2015-01-01

This presentation gives an overview of Observing System Simulation Experiments (OSSEs). The components of an OSSE are described, along with discussion of the process for validating, calibrating, and performing experiments. a.

Pit formation observed in a multilayer dielectric coating as a result of simulated space environmental exposure

NASA Astrophysics Data System (ADS)

Fuqua, Peter D.; Presser, Nathan; Barrie, James D.; Meshishnek, Michael J.; Coleman, Dianne J.

2002-06-01

Certain spaceborne telescope designs require that dielectric-coated lenses be exposed to the energetic electrons and protons associated with the space environment. Test coupons that were exposed to a simulated space environment showed extensive pitting as a result of dielectric breakdown. A typical pit was 50-100 mum at the surface and extended to the substrate material, in which a 10-mum-diameter melt region was found. Pitting was not observed on similar samples that had also been overcoated with a transparent conductive thin film. Measurement of the bidirectional reflectance distribution transfer function showed that pitting caused a fivefold to tenfold increase in the scattering of visible light.

Ionospheric Trend Over Wuhan During 1947-2017: Comparison Between Simulation and Observation

NASA Astrophysics Data System (ADS)

Yue, Xinan; Hu, Lianhuan; Wei, Yong; Wan, Weixing; Ning, Baiqi

2018-02-01

Since Roble and Dickinson (1989), who drew the community's attention about the greenhouse gas effect on the ionosphere, huge efforts have been implemented on ionospheric climate study. However, direct comparison between observations and simulations is still rare. Recently, the Wuhan ionosonde observations were digitized and standardized through unified method back to 1947. In this study, the NCAR-TIEGCM was driven by Mauna Loa Observatory observed CO2 level and International Geomagnetic Reference Field (IGRF) geomagnetic field to simulate their effects on ionospheric long-term trend over Wuhan. Only March equinox was considered in both data analysis and simulation. Simulation results show that the CO2 and geomagnetic field have comparable effect on hmF2 trend, while geomagnetic field effect is stronger than CO2 on foF2 trend over Wuhan. Both factors result in obvious but different diurnal variations of foF2/hmF2 long-term trends. The geomagnetic field effect is nonlinear versus years since the long-term variation of geomagnetic field intensity and orientation is complex. Mean value of foF2 and hmF2 trend is (-0.0021 MHz/yr, -0.106 km/yr) and (-0.0022 MHz/yr, -0.0763 km/yr) for observation and simulation, respectively. Regarding the diurnal variation of the trend, the simulation accords well with that of observation except hmF2 results around 12 UT. Overall, good agreement between observation and simulation illustrates the good quality of Wuhan ionosonde long-term data and the validity of ancient ionosphere reconstruction based on realistic indices driving simulation.

Cassini radar : system concept and simulation results

NASA Astrophysics Data System (ADS)

Melacci, P. T.; Orosei, R.; Picardi, G.; Seu, R.

1998-10-01

The Cassini mission is an international venture, involving NASA, the European Space Agency (ESA) and the Italian Space Agency (ASI), for the investigation of the Saturn system and, in particular, Titan. The Cassini radar will be able to see through Titan's thick, optically opaque atmosphere, allowing us to better understand the composition and the morphology of its surface, but the interpretation of the results, due to the complex interplay of many different factors determining the radar echo, will not be possible without an extensive modellization of the radar system functioning and of the surface reflectivity. In this paper, a simulator of the multimode Cassini radar will be described, after a brief review of our current knowledge of Titan and a discussion of the contribution of the Cassini radar in answering to currently open questions. Finally, the results of the simulator will be discussed. The simulator has been implemented on a RISC 6000 computer by considering only the active modes of operation, that is altimeter and synthetic aperture radar. In the instrument simulation, strict reference has been made to the present planned sequence of observations and to the radar settings, including burst and single pulse duration, pulse bandwidth, pulse repetition frequency and all other parameters which may be changed, and possibly optimized, according to the operative mode. The observed surfaces are simulated by a facet model, allowing the generation of surfaces with Gaussian or non-Gaussian roughness statistic, together with the possibility of assigning to the surface an average behaviour which can represent, for instance, a flat surface or a crater. The results of the simulation will be discussed, in order to check the analytical evaluations of the models of the average received echoes and of the attainable performances. In conclusion, the simulation results should allow the validation of the theoretical evaluations of the capabilities of microwave instruments, when

Ocean Carbon States: Data Mining in Observations and Numerical Simulations Results

NASA Astrophysics Data System (ADS)

Latto, R.; Romanou, A.

2017-12-01

Advanced data mining techniques are rapidly becoming widely used in Climate and Earth Sciences with the purpose of extracting new meaningful information from increasingly larger and more complex datasets. This is particularly important in studies of the global carbon cycle, where any lack of understanding of its combined physical and biogeochemical drivers is detrimental to our ability to accurately describe, understand, and predict CO2 concentrations and their changes in the major carbon reservoirs. The analysis presented here evaluates the use of cluster analysis as a means of identifying and comparing spatial and temporal patterns extracted from observational and model datasets. As the observational data is organized into various regimes, which we will call "ocean carbon states", we gain insight into the physical and/or biogeochemical processes controlling the ocean carbon cycle as well as how well these processes are simulated by a state-of-the-art climate model. We find that cluster analysis effectively produces realistic, dynamic regimes that can be associated with specific processes at different temporal scales for both observations and the model. In addition, we show how these regimes can be used to illustrate and characterize the model biases in the model air-sea flux of CO2. These biases are attributed to biases in salinity, sea surface temperature, wind speed, and nitrate, which are then used to identify the physical processes that are inaccurately reproduced by the model. In this presentation, we provide a proof-of-concept application using simple datasets, and we expand to more complex ones, using several physical and biogeochemical variable pairs, thus providing considerable insight into the mechanisms and phases of the ocean carbon cycle over different temporal and spatial scales.

Atmospheric icing of structures: Observations and simulations

NASA Astrophysics Data System (ADS)

Ágústsson, H.; Elíasson, Á. J.; Thorsteins, E.; Rögnvaldsson, Ó.; Ólafsson, H.

2012-04-01

This study compares observed icing in a test span in complex orography at Hallormsstaðaháls (575 m) in East-Iceland with parameterized icing based on an icing model and dynamically downscaled weather at high horizontal resolution. Four icing events have been selected from an extensive dataset of observed atmospheric icing in Iceland. A total of 86 test-spans have been erected since 1972 at 56 locations in complex terrain with more than 1000 icing events documented. The events used here have peak observed ice load between 4 and 36 kg/m. Most of the ice accretion is in-cloud icing but it may partly be mixed with freezing drizzle and wet snow icing. The calculation of atmospheric icing is made in two steps. First the atmospheric data is created by dynamically downscaling the ECMWF-analysis to high resolution using the non-hydrostatic mesoscale Advanced Research WRF-model. The horizontal resolution of 9, 3, 1 and 0.33 km is necessary to allow the atmospheric model to reproduce correctly local weather in the complex terrain of Iceland. Secondly, the Makkonen-model is used to calculate the ice accretion rate on the conductors based on the simulated temperature, wind, cloud and precipitation variables from the atmospheric data. In general, the atmospheric model correctly simulates the atmospheric variables and icing calculations based on the atmospheric variables correctly identify the observed icing events, but underestimate the load due to too slow ice accretion. This is most obvious when the temperature is slightly below 0°C and the observed icing is most intense. The model results improve significantly when additional observations of weather from an upstream weather station are used to nudge the atmospheric model. However, the large variability in the simulated atmospheric variables results in high temporal and spatial variability in the calculated ice accretion. Furthermore, there is high sensitivity of the icing model to the droplet size and the possibility that

Tag team simulation: An innovative approach for promoting active engagement of participants and observers during group simulations.

PubMed

Levett-Jones, Tracy; Andersen, Patrea; Reid-Searl, Kerry; Guinea, Stephen; McAllister, Margaret; Lapkin, Samuel; Palmer, Lorinda; Niddrie, Marian

2015-09-01

Active participation in immersive simulation experiences can result in technical and non-technical skill enhancement. However, when simulations are conducted in large groups, maintaining the interest of observers so that they do not disengage from the learning experience can be challenging. We implemented Tag Team Simulation with the aim of ensuring that both participants and observers had active and integral roles in the simulation. In this paper we outline the features of this innovative approach and provide an example of its application to a pain simulation. Evaluation was conducted using the Satisfaction with Simulation Experience Scale. A total of 444 year nursing students participated from a population of 536 (response rate 83%). Cronbach's alpha for the Scale was .94 indicating high internal consistency. The mean satisfaction score for participants was 4.63 compared to 4.56 for observers. An independent sample t test revealed no significant difference between these scores (t (300) = -1.414, p = 0.16). Tag team simulation is an effective approach for ensuring observers' and participants' active involvement during group-based simulations and one that is highly regarded by students. It has the potential for broad applicability across a range of leaning domains both within and beyond nursing. Copyright © 2015 Elsevier Ltd. All rights reserved.

LES ARM Symbiotic Simulation and Observation (LASSO) Implementation Strategy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gustafson Jr., WI; Vogelmann, AM

2015-09-01

This document illustrates the design of the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) workflow to provide a routine, high-resolution modeling capability to augment the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility’s high-density observations. LASSO will create a powerful new capability for furthering ARM’s mission to advance understanding of cloud, radiation, aerosol, and land-surface processes. The combined observational and modeling elements will enable a new level of scientific inquiry by connecting processes and context to observations and providing needed statistics for details that cannot be measured. The result will be improved process understandingmore » that facilitates concomitant improvements in climate model parameterizations. The initial LASSO implementation will be for ARM’s Southern Great Plains site in Oklahoma and will focus on shallow convection, which is poorly simulated by climate models due in part to clouds’ typically small spatial scale compared to model grid spacing, and because the convection involves complicated interactions of microphysical and boundary layer processes.« less

Biases and systematics in the observational derivation of galaxy properties: comparing different techniques on synthetic observations of simulated galaxies

NASA Astrophysics Data System (ADS)

Guidi, Giovanni; Scannapieco, Cecilia; Walcher, C. Jakob

2015-12-01

We study the sources of biases and systematics in the derivation of galaxy properties from observational studies, focusing on stellar masses, star formation rates, gas and stellar metallicities, stellar ages, magnitudes and colours. We use hydrodynamical cosmological simulations of galaxy formation, for which the real quantities are known, and apply observational techniques to derive the observables. We also analyse biases that are relevant for a proper comparison between simulations and observations. For our study, we post-process the simulation outputs to calculate the galaxies' spectral energy distributions (SEDs) using stellar population synthesis models and also generate the fully consistent far-UV-submillimetre wavelength SEDs with the radiative transfer code SUNRISE. We compared the direct results of simulations with the observationally derived quantities obtained in various ways, and found that systematic differences in all studied galaxy properties appear, which are caused by: (1) purely observational biases, (2) the use of mass-weighted and luminosity-weighted quantities, with preferential sampling of more massive and luminous regions, (3) the different ways of constructing the template of models when a fit to the spectra is performed, and (4) variations due to different calibrations, most notably for gas metallicities and star formation rates. Our results show that large differences can appear depending on the technique used to derive galaxy properties. Understanding these differences is of primary importance both for simulators, to allow a better judgement of similarities and differences with observations, and for observers, to allow a proper interpretation of the data.

Results from Binary Black Hole Simulations in Astrophysics Applications

NASA Technical Reports Server (NTRS)

Baker, John G.

2007-01-01

Present and planned gravitational wave observatories are opening a new astronomical window to the sky. A key source of gravitational waves is the merger of two black holes. The Laser Interferometer Space Antenna (LISA), in particular, is expected to observe these events with signal-to-noise ratio's in the thousands. To fully reap the scientific benefits of these observations requires a detailed understanding, based on numerical simulations, of the predictions of General Relativity for the waveform signals. New techniques for simulating binary black hole mergers, introduced two years ago, have led to dramatic advances in applied numerical simulation work. Over the last two years, numerical relativity researchers have made tremendous strides in understanding the late stages of binary black hole mergers. Simulations have been applied to test much of the basic physics of binary black hole interactions, showing robust results for merger waveform predictions, and illuminating such phenomena as spin-precession. Calculations have shown that merging systems can be kicked at up to 2500 km/s by the thrust from asymmetric emission. Recently, long lasting simulations of ten or more orbits allow tests of post-Newtonian (PN) approximation results for radiation from the last orbits of the binary's inspiral. Already, analytic waveform models based PN techniques with incorporated information from numerical simulations may be adequate for observations with current ground based observatories. As new advances in simulations continue to rapidly improve our theoretical understanding of the systems, it seems certain that high-precision predictions will be available in time for LISA and other advanced ground-based instruments. Future gravitational wave observatories are expected to make precision.

Experimental and simulational result multipactors in 112 MHz QWR injector

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xin, T.; Ben-Zvi, I.; Belomestnykh, S.

2015-05-03

The first RF commissioning of 112 MHz QWR superconducting electron gun was done in late 2014. The coaxial Fundamental Power Coupler (FPC) and Cathode Stalk (stalk) were installed and tested for the first time. During this experiment, we observed several multipacting barriers at different gun voltage levels. The simulation work was done within the same range. The comparison between the experimental observation and the simulation results are presented in this paper. The observations during the test are consisted with the simulation predictions. We were able to overcome most of the multipacting barriers and reach 1.8 MV gun voltage under pulsedmore » mode after several round of conditioning processes.« less

Simulations of star-forming molecular clouds: observational predictions

NASA Astrophysics Data System (ADS)

Zhang, Shangjia; Hartmann, Lee; Kuznetsova, Aleksandra; Abelardo Zamora, Manuel

2018-01-01

Observations of protostellar molecular cloud cores can be used to test theories of star formation. However, observational results can be biased because of limited information: (a) only two spatial dimensions and one velocity dimension can be measured, (b) and cores generally are not spherically symmetric. We use numerical simulations of the formation and collapse of molecular gas with sink particles to make observational predictions. We use the radiative transfer code LIME to predict CO and NH3 channel maps. We find reasonable agreement with observed velocity structures and gradients but occasional large differences depending on viewing angle.

Observing System Simulation Experiments for Fun and Profit

NASA Technical Reports Server (NTRS)

Prive, Nikki C.

2015-01-01

Observing System Simulation Experiments can be powerful tools for evaluating and exploring both the behavior of data assimilation systems and the potential impacts of future observing systems. With great power comes great responsibility - given a pure modeling framework, how can we be sure our results are meaningful? The challenges and pitfalls of OSSE calibration and validation will be addressed, as well as issues of incestuousness, selection of appropriate metrics, and experiment design. The use of idealized observational networks to investigate theoretical ideas in a fully complex modeling framework will also be discussed

Laboratory observations and simulations of phase reddening

NASA Astrophysics Data System (ADS)

Schröder, S. E.; Grynko, Ye.; Pommerol, A.; Keller, H. U.; Thomas, N.; Roush, T. L.

2014-09-01

The visible reflectance spectrum of many Solar System bodies changes with changing viewing geometry for reasons not fully understood. It is often observed to redden (increasing spectral slope) with increasing solar phase angle, an effect known as phase reddening. Only once, in an observation of the martian surface by the Viking 1 lander, was reddening observed up to a certain phase angle with bluing beyond, making the reflectance ratio as a function of phase angle shaped like an arch. However, in laboratory experiments this arch-shape is frequently encountered. To investigate this, we measured the bidirectional reflectance of particulate samples of several common rock types in the 400-1000 nm wavelength range and performed ray-tracing simulations. We confirm the occurrence of the arch for surfaces that are forward scattering, i.e. are composed of semi-transparent particles and are smooth on the scale of the particles, and for which the reflectance increases from the lower to the higher wavelength in the reflectance ratio. The arch shape is reproduced by the simulations, which assume a smooth surface. However, surface roughness on the scale of the particles, such as the Hapke and van Horn (Hapke, B., van Horn, H. [1963]. J. Geophys. Res. 68, 4545-4570) fairy castles that can spontaneously form when sprinkling a fine powder, leads to monotonic reddening. A further consequence of this form of microscopic roughness (being indistinct without the use of a microscope) is a flattening of the disk function at visible wavelengths, i.e. Lommel-Seeliger-type scattering. The experiments further reveal monotonic reddening for reflectance ratios at near-IR wavelengths. The simulations fail to reproduce this particular reddening, and we suspect that it results from roughness on the surface of the particles. Given that the regolith of atmosphereless Solar System bodies is composed of small particles, our results indicate that the prevalence of monotonic reddening and Lommel

Observation and Simulation of Microseisms Offshore Ireland

NASA Astrophysics Data System (ADS)

Le Pape, Florian; Bean, Chris; Craig, David; Jousset, Philippe; Donne, Sarah; Möllhoff, Martin

2017-04-01

Although more and more used in seismic imagery, ocean induced ambient seismic noise is still not so well understood, particularly how the signal propagates from ocean to land. Between January and September 2016, 10 broadband Ocean Bottom Seismometers (OBSs) stations, including acoustic sensors (hydrophone), were deployed across the shelf offshore Donegal and out into the Rockall Trough. The preliminary results show spatial and temporal variability in the ocean generated seismic noise which holds information about changes in the generation source process, including meteorological information, but also in the geological structure. In addition to the collected OBS data, numerical simulations of acoustic/seismic wave propagation are also considered in order to study the spatio-temporal variation of the broadband acoustic wavefield and its connection with the measured seismic wavefield in the region. Combination of observations and simulations appears significant to better understand what control the acoustic/seismic coupling at the sea floor as well as the effect of the water column and sediments thickness on signal propagation. Ocean generated seismic ambient noise recorded at the seafloor appears to behave differently in deep and shallow water and 3D simulations of acoustic/seismic wave propagation look particularly promising for reconciling deep ocean, shelf and land seismic observations.

A Statistical Comparison of PSC Model Simulations and POAM Observations

NASA Technical Reports Server (NTRS)

Strawa, A. W.; Drdla, K.; Fromm, M.; Bokarius, K.; Gore, Warren J. (Technical Monitor)

2002-01-01

A better knowledge of PSC composition and formation mechanisms is important to better understand and predict stratospheric ozone depletion. Several past studies have attempted to compare modeling results with satellite observations. These comparisons have concentrated on case studies. In this paper we adopt a statistical approach. POAM PSC observations from several Arctic winters are categorized into Type Ia and Ib PSCs using a technique based on Strawa et al. The discrimination technique has been modified to employ the wavelengths dependence of the extinction signal at all wavelengths rather than only at 603 and 10 18 nm. Winter-long simulations for the 1999-2000 Arctic winter have been made using the IMPACT model. These simulations have been constrained by aircraft observations made during the SOLVE/THESEO 2000 campaign. A complete set of winter-long simulations was run for several different microphysical and PSC formation scenarios. The simulations give us perfect knowledge of PSC type (Ia, Ib, or II), composition, especially condensed phase HNO3 which is important for denitrification, and condensed phase H2O. Comparisons are made between the simulation and observation of PSC extinction at 1018 rim versus wavelength dependence, winter-long percentages of Ia and Ib occurrence, and temporal and altitude trends of the PSCs. These comparisons allow us to comment on how realistic some modeling scenarios are.

Hemispheric Coupling: Comparing Dynamo Simulations and Observations

NASA Astrophysics Data System (ADS)

Norton, A. A.; Charbonneau, P.; Passos, D.

2014-12-01

Numerical simulations that reproduce solar-like magnetic cycles can be used to generate long-term statistics. The variations in north-south hemispheric solar cycle synchronicity and amplitude produced in simulations has not been widely compared to observations. The observed limits on solar cycle amplitude and phase asymmetry show that hemispheric sunspot area production is no more than 20 % asymmetric for cycles 17-23 and that phase lags do not exceed 20 % (or two years) of the total cycle period, as determined from Royal Greenwich Observatory sunspot data. Several independent studies have found a long-term trend in phase values as one hemisphere leads the other for, on average, four cycles. Such persistence in phase is not indicative of a stochastic phenomenon. We compare these observational findings to the magnetic cycle found in a numerical simulation of solar convection recently produced with the EULAG-MHD model. This long "millennium simulation" spans more than 1600 years and generated 40 regular, sunspot-like cycles. While the simulated cycle length is too long (˜40 yrs) and the toroidal bands remain at too high of latitudes (>30°), some solar-like aspects of hemispheric asymmetry are reproduced. The model is successful at reproducing the synchrony of polarity inversions and onset of cycle as the simulated phase lags do not exceed 20 % of the cycle period. The simulated amplitude variations between the north and south hemispheres are larger than those observed in the Sun, some up to 40 %. An interesting note is that the simulations also show that one hemisphere can persistently lead the other for several successive cycles, placing an upper bound on the efficiency of transequatorial magnetic coupling mechanisms. These include magnetic diffusion, cross-equatorial mixing within latitudinally-elongated convective rolls (a.k.a. "banana cells") and transequatorial meridional flow cells. One or more of these processes may lead to magnetic flux cancellation whereby

Observing the baryon cycle in hydrodynamic cosmological simulations

NASA Astrophysics Data System (ADS)

Vander Vliet, Jacob Richard

An understanding of galaxy evolution requires an understanding of the flow of baryons in and out of a galaxy. The accretion of baryons is required for galaxies to form stars, while stars eject baryons out of the galaxy through stellar feedback mechanisms such as supernovae, stellar winds, and radiation pressure. The interplay between outfiowing and infalling material form the circumgalactic medium (CGM). Hydrodynamic simulations provide understanding of the connection between stellar feedback and the distribution and kinematics of baryons in the CGM. To compare simulations and observations properly the simulated CGI must be observed in the same manner as the real CGM. I have developed the Mockspec code to generate synthetic quasar absorption line observations of the CGM in cosmological hydrodynamic simulations. Mockspec generates synthetic spectra based on the phase; lnetallicity, and kinematics of CGM gas and mimics instrumental effects. Mockspec includes automated analysis of the spectra and identifies the gas responsible for the absorption. Mockspec was applied to simulations of dwarf galaxies at low redshift to examine the observable effect different feedback models have on the CGM. While the different feedback models had strong effects on the galaxy, they all produced a similar CGM that failed match observations. Mockspec was applied to the VELA simulation suite of high redshift, high mass galaxies to examine the variance of the CGM across different galaxies in different environments. The observed CGM showed little variation between the different galaxies and almost no evolution from z=4 to z=1. The VELAs were not able to generate a CGM to match the observations. The properties of cells responsible for the absorption were compared to the derived properties from Voigt Profile decomposition. VP modeling was found to accurately describe the HI and MgII absorbing gas but failed for high ionization species such as CIV and OVI, which do not arise in the coherent

Evaluation of simulated tropical convective updraft hydrometeor properties using aircraft observations

NASA Astrophysics Data System (ADS)

Stanford, McKenna W.

The High Altitude Ice Crystals - High Ice Water Content (HAIC-HIWC) field campaign produced aircraft retrievals of total condensed water content (TWC), hydrometeor particle size distributions, and vertical velocity (w) in high ice water content regions of tropical mesoscale convective systems (MCSs). These observations are used to evaluate deep convective updraft properties in high-resolution nested Weather Research and Forecasting (WRF) simulations of observed MCSs. Because simulated hydrometeor properties are highly sensitive to the parameterization of microphysics, three commonly used microphysical parameterizations are tested, including two bulk schemes (Thompson and Morrison) and one bin scheme (Fast Spectral Bin Microphysics). A commonly documented bias in cloud-resolving simulations is the exaggeration of simulated radar reflectivities aloft in tropical MCSs. This may result from overly strong convective updrafts that loft excessive condensate mass and from simplified approximations of hydrometeor size distributions, properties, species separation, and microphysical processes. The degree to which the reflectivity bias is a separate function of convective dynamics, condensate mass, and hydrometeor size has yet to be addressed. This research untangles these components by comparing simulated and observed relationships between w, TWC, and hydrometer size as a function of temperature. All microphysics schemes produce median mass diameters that are generally larger than observed for temperatures between -10 °C and -40 °C and TWC > 1 g m-3. Observations produce a prominent mode in the composite mass size distribution around 300 microm, but under most conditions, all schemes shift the distribution mode to larger sizes. Despite a much greater number of samples, all simulations fail to reproduce observed high TWC or high w conditions between -20 °C and -40 °C in which only a small fraction of condensate mass is found in relatively large particle sizes. Increasing

Magneto-acoustic wave energy in sunspots: observations and numerical simulations

NASA Astrophysics Data System (ADS)

Felipe, T.; Khomenko, E.; Collados, M.; Beck, C.

2011-11-01

We have reproduced some sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numericalsimulations. The results of the simulations arecompared with the oscillations observed simultaneously at different heights from the SiI lambda10827Å line, HeI lambda10830Å line, the CaII H core and the FeI blends at the wings of the CaII H line. The simulations show a remarkable agreement with the observations, and we have used them to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is 5-10 times lower than the amount needed to balance the chromospheric radiative losses.

Observations and Simulations of Electron Dynamics Near an Active Neutral Line

NASA Technical Reports Server (NTRS)

Goldstein, M. L.; Hwang, Kyoung-Joo; Ashour-Abdalla, Maha; El-Aloui, Mostafa; Schriver, David; Richard, Robert; Zhou, Meng; Walker, Ray

2010-01-01

Recent observations in the Earth's magnetotail have shown rapid increases in the fluxes of energetic electrons with energies up to 100's of keV associated with dipolarization fronts that propagate into the inner magnetosphere. On August 15, 2001 the four Cluster spacecraft located slightly dawnward of midnight (yGSM approx. -5.4RE) at xGSM approx. -18RE observed a series of earthward propagating dipolarization fronts [Hwang et al., 2010]. At least 6 dipolarization fronts were observed in a 20m interval. Unlike previously reported cases the fluxes of electrons up to 95keV decreased during the passage of the first three fronts over the spacecraft. The energetic electron fluxes increased during the passage of the last three fronts. We have performed a global magnetohydrodynamic simulation of this event using solar wind observations from the ACE satellite to drive the simulation. In the simulation a very complex reconnection system in the near-Earth tail at XGSM approx. -20RE launched a series of earthward propagating dipolarization fronts that are similar to those observed on Cluster. The simulation results indicate that the Cluster spacecraft were just earthward of the reconnection site. In this paper we will present a study of the dynamics of electrons associated with these events by using the large-scale kinetic simulation approach in which we launch a large number of electrons into the electric and magnetic fields from this simulation.

Simulating AIA observations of a flux rope ejection

NASA Astrophysics Data System (ADS)

Pagano, P.; Mackay, D. H.; Poedts, S.

2014-08-01

Context. Coronal mass ejections (CMEs) are the most violent phenomena observed on the Sun. Currently, extreme ultraviolet (EUV) images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) are providing new insights into the early phase of CME evolution. In particular, observations now show the ejection of magnetic flux ropes from the solar corona and how they evolve into CMEs. While this is the case, these observations are difficult to interpret in terms of basic physical mechanisms and quantities. To fully understand CMEs we need to compare equivalent quantities derived from both observations and theoretical models. This will aid in bridging the gap between observations and models. Aims: To this end, we aim to produce synthesised AIA observations from simulations of a flux rope ejection. To carry this out we include the role of thermal conduction and radiative losses, both of which are important for determining the temperature distribution of the solar corona during a CME. Methods: We perform a simulation where a flux rope is ejected from the solar corona. From the density and temperature of the plasma in the simulation we synthesise AIA observations. The emission is then integrated along the line of sight using the instrumental response function of AIA. Results: We sythesise observations of AIA in the channels at 304 Å, 171 Å, 335 Å, and 94 Å. The synthesised observations show a number of features similar to actual observations and in particular reproduce the general development of CMEs in the low corona as observed by AIA. In particular we reproduce an erupting and expanding arcade in the 304 Å and 171 Å channels with a high density core. Conclusions: The ejection of a flux rope reproduces many of the features found in the AIA observations. This work is therefore a step forward in bridging the gap between observations and models, and can lead to more direct interpretations of EUV observations in terms of flux rope

Simulations of extragalactic magnetic fields and of their observables

NASA Astrophysics Data System (ADS)

Vazza, F.; Brüggen, M.; Gheller, C.; Hackstein, S.; Wittor, D.; Hinz, P. M.

2017-12-01

The origin of extragalactic magnetic fields is still poorly understood. Based on a dedicated suite of cosmological magneto-hydrodynamical simulations with the ENZO code we have performed a survey of different models that may have caused present-day magnetic fields in galaxies and galaxy clusters. The outcomes of these models differ in cluster outskirts, filaments, sheets and voids and we use these simulations to find observational signatures of magnetogenesis. With these simulations, we predict the signal of extragalactic magnetic fields in radio observations of synchrotron emission from the cosmic web, in Faraday rotation, in the propagation of ultra high energy cosmic rays, in the polarized signal from fast radio bursts at cosmological distance and in spectra of distant blazars. In general, primordial scenarios in which present-day magnetic fields originate from the amplification of weak (⩽nG ) uniform seed fields result in more homogeneous and relatively easier to observe magnetic fields than astrophysical scenarios, in which present-day fields are the product of feedback processes triggered by stars and active galaxies. In the near future the best evidence for the origin of cosmic magnetic fields will most likely come from a combination of synchrotron emission and Faraday rotation observed at the periphery of large-scale structures.

Validation: Codes to compare simulation data to various observations

NASA Astrophysics Data System (ADS)

Cohn, J. D.

2017-02-01

Validation provides codes to compare several observations to simulated data with stellar mass and star formation rate, simulated data stellar mass function with observed stellar mass function from PRIMUS or SDSS-GALEX in several redshift bins from 0.01-1.0, and simulated data B band luminosity function with observed stellar mass function, and to create plots for various attributes, including stellar mass functions, and stellar mass to halo mass. These codes can model predictions (in some cases alongside observational data) to test other mock catalogs.

Reconciling Simulated and Observed Views of Clouds: MODIS, ISCCP, and the Limits or Instrument Simulators

NASA Technical Reports Server (NTRS)

Ackerman, Steven A.; Hemler, Richard S.; Hofman, Robert J. Patrick; Pincus, Robert; Platnick, Steven

2011-01-01

The properties of clouds that may be observed by satellite instruments, such as optical depth and cloud top pressure, are only loosely related to the way clouds m-e represented in models of the atmosphere. One way to bridge this gap is through "instrument simulators," diagnostic tools that map the model representation to synthetic observations so that differences between simulator output and observations can be interpreted unambiguously as model error. But simulators may themselves be restricted by limited information available from the host model or by internal assumptions. This paper considers the extent to which instrument simulators are able to capture essential differences between MODIS and ISCCP, two similar but independent estimates of cloud properties. The authors review the measurements and algorithms underlying these two cloud climatologies, introduce a MODIS simulator, and detail data sets developed for comparison with global models using ISCCP and MODIS simulators, In nature MODIS observes less mid-level doudines!> than ISCCP, consistent with the different methods used to determine cloud top pressure; aspects of this difference are reproduced by the simulators running in a climate modeL But stark differences between MODIS and ISCCP observations of total cloudiness and the distribution of cloud optical thickness can be traced to different approaches to marginal pixels, which MODIS excludes and ISCCP treats as homogeneous. These pixels, which likely contain broken clouds, cover about 15 k of the planet and contain almost all of the optically thinnest clouds observed by either instrument. Instrument simulators can not reproduce these differences because the host model does not consider unresolved spatial scales and so can not produce broken pixels. Nonetheless, MODIS and ISCCP observation are consistent for all but the optically-thinnest clouds, and models can be robustly evaluated using instrument simulators by excluding ambiguous observations.

Observational Tracers of Hot and Cold Gas in Isolated Galaxy Simulations

NASA Astrophysics Data System (ADS)

Brzycki, Bryan; Silvia, Devin

2018-01-01

We present results from an analysis comparing simulations of isolated spiral galaxies with recent observations of the circumgalactic medium (CGM). As the interface containing inflows and outflows between the interstellar and intergalactic media, the CGM plays an important role in the composition and evolution of galaxies. Using a set of isolated galaxy simulations over different initial conditions and star formation and feedback parameters, we investigate the evolution of CGM gas. Specifically, in light of recent observational studies, we compute the radial column density profiles and covering fractions of various observable ion species (H I, C IV, O VI, Mg II, Si III) for each simulated galaxy. Taking uniformly random sightlines through the CGM of each simulated galaxy, we find the abundance of gas absorbers and analyze their contribution to the overall column density along each sightline. By identifying the prevalence of high column density absorbers, we seek to characterize the distribution and evolution of observable ion species in the CGM. We also highlight a subset of our isolated galaxy simulations that produce and maintain a stable precipitating CGM that fuels high rates of sustained star formation. This project was supported in part by the NSF REU grant AST-1358980 and by the Nantucket Maria Mitchell Association.

Observations and statistical simulations of a proposed solar cycle/QBO/weather relationship

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baldwin, M.P.; Dunkerton, T.J.

1989-08-01

The 10.7 cm solar flux is observed to be highly correlated with north pole stratospheric temperatures when partitioned according to the phase of the equatorial stratospheric winds (the quasi-biennial oscillation, or QBO). The authors supplement observations with calculations showing that temperatures over most of the northern hemisphere are highly correlated or anticorrelated with north pole temperatures. The observed spatial pattern of solar cycle correlations at high latitudes is shown to be not unique to the solar cycle. The authors present results, similar to the observed solar cycle correlations, with simulated harmonics of various periods replacing the solar cycle. These calculationsmore » demonstrate the correlations at least as high as those for the solar cycle results may be obtained using simulated harmonics.« less

Science Results from a Mars Drilling Simulation (Río Tinto, Spain) and Ground Truth for Remote Science Observations

NASA Astrophysics Data System (ADS)

Bonaccorsi, Rosalba; Stoker, Carol R.

2008-10-01

Science results from a field-simulated lander payload and post-mission laboratory investigations provided "ground truth" to interpret remote science observations made as part of the 2005 Mars Astrobiology Research and Technology Experiment (MARTE) drilling mission simulation. The experiment was successful in detecting evidence for life, habitability, and preservation potential of organics in a relevant astrobiological analogue of Mars. Science results. Borehole 7 was drilled near the Río Tinto headwaters at Peña de Hierro (Spain) in the upper oxidized remnant of an acid rock drainage system. Analysis of 29 cores (215 cm of core was recovered from 606 cm penetrated depth) revealed a matrix of goethite- (42-94%) and hematite-rich (47-87%) rocks with pockets of phyllosilicates (47-74%) and fine- to coarse-grained loose material. Post-mission X-ray diffraction (XRD) analysis confirmed the range of hematite:goethite mixtures that were visually recognizable (˜1:1, ˜1:2, and ˜1:3 mixtures displayed a yellowish-red color whereas 3:1 mixtures displayed a dark reddish-brown color). Organic carbon was poorly preserved in hematite/goethite-rich materials (Corg <0.05 wt %) beneath the biologically active organic-rich soil horizon (Corg ˜3-11 wt %) in contrast to the phyllosilicate-rich zones (Corg ˜0.23 wt %). Ground truth vs. remote science analysis. Laboratory-based analytical results were compared to the analyses obtained by a Remote Science Team (RST) using a blind protocol. Ferric iron phases, lithostratigraphy, and inferred geologic history were correctly identified by the RST with the exception of phyllosilicate-rich materials that were misinterpreted as weathered igneous rock. Adenosine 5‧-triphosphate (ATP) luminometry, a tool available to the RST, revealed ATP amounts above background noise, i.e., 278-876 Relative Luminosity Units (RLUs) in only 6 cores, whereas organic carbon was detected in all cores. Our manned vs. remote observations based on automated

Science results from a Mars drilling simulation (Río Tinto, Spain) and ground truth for remote science observations.

PubMed

Bonaccorsi, Rosalba; Stoker, Carol R

2008-10-01

Science results from a field-simulated lander payload and post-mission laboratory investigations provided "ground truth" to interpret remote science observations made as part of the 2005 Mars Astrobiology Research and Technology Experiment (MARTE) drilling mission simulation. The experiment was successful in detecting evidence for life, habitability, and preservation potential of organics in a relevant astrobiological analogue of Mars. SCIENCE RESULTS: Borehole 7 was drilled near the Río Tinto headwaters at Peña de Hierro (Spain) in the upper oxidized remnant of an acid rock drainage system. Analysis of 29 cores (215 cm of core was recovered from 606 cm penetrated depth) revealed a matrix of goethite- (42-94%) and hematite-rich (47-87%) rocks with pockets of phyllosilicates (47-74%) and fine- to coarse-grained loose material. Post-mission X-ray diffraction (XRD) analysis confirmed the range of hematite:goethite mixtures that were visually recognizable (approximately 1:1, approximately 1:2, and approximately 1:3 mixtures displayed a yellowish-red color whereas 3:1 mixtures displayed a dark reddish-brown color). Organic carbon was poorly preserved in hematite/goethite-rich materials (C(org) <0.05 wt %) beneath the biologically active organic-rich soil horizon (C(org) approximately 3-11 wt %) in contrast to the phyllosilicate-rich zones (C(org) approximately 0.23 wt %). GROUND TRUTH VS. REMOTE SCIENCE ANALYSIS: Laboratory-based analytical results were compared to the analyses obtained by a Remote Science Team (RST) using a blind protocol. Ferric iron phases, lithostratigraphy, and inferred geologic history were correctly identified by the RST with the exception of phyllosilicate-rich materials that were misinterpreted as weathered igneous rock. Adenosine 5'-triphosphate (ATP) luminometry, a tool available to the RST, revealed ATP amounts above background noise, i.e., 278-876 Relative Luminosity Units (RLUs) in only 6 cores, whereas organic carbon was detected in all

Understanding observed and simulated historical temperature trends in California

NASA Astrophysics Data System (ADS)

Bonfils, C. J.; Duffy, P. B.; Santer, B. D.; Lobell, D. B.; Wigley, T. M.

2006-12-01

In our study, we attempt 1) to improve our understanding of observed historical temperature trends and their underlying causes in the context of regional detection of climate change and 2) to identify possible neglected forcings and errors in the model response to imposed forcings at the origin of inconsistencies between models and observations. From eight different observational datasets, we estimate California-average temperature trends over 1950- 1999 and compare them to trends from a suite of IPCC control simulations of natural internal climate variability. We find that the substantial night-time warming occurring from January to September is inconsistent with model-based estimates of natural internal climate variability, and thus requires one or more external forcing agents to be explained. In contrast, we find that a significant day-time warming occurs only from January to March. Our confidence in these findings is increased because there is no evidence that the models systematically underestimate noise on interannual and decadal timescales. However, we also find that IPCC simulations of the 20th century that include combined anthropogenic and natural forcings are not able to reproduce such a pronounced seasonality of the trends. Our first hypothesis is that the warming of Californian winters over the second half of the twentieth century is associated with changes in large-scale atmospheric circulation that are likely to be human-induced. This circulation change is underestimated in the historical simulations, which may explain why the simulated warming of Californian winters is too weak. We also hypothesize that the lack of a detectable observed increase in summertime maximum temperature arises from a cooling associated with large-scale irrigation. This cooling may have, until now, counteracted the warming induced by increasing greenhouse gases and urbanization effects. Omitting to include this forcing in the simulations can result in overestimating the

Observations and simulations of the ionospheric lunar tide: Seasonal variability

NASA Astrophysics Data System (ADS)

Pedatella, N. M.

2014-07-01

The seasonal variability of the ionospheric lunar tide is investigated using a combination of Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) observations and thermosphere-ionosphere-mesosphere electrodynamics general circulation model (TIME-GCM) simulations. The present study focuses on the seasonal variability of the lunar tide in the ionosphere and its potential connection to the occurrence of stratosphere sudden warmings (SSWs). COSMIC maximum F region electron density (NmF2) and total electron content observations reveal a primarily annual variation of the ionospheric lunar tide, with maximum amplitudes occurring at low latitudes during December-February. Simulations of the lunar tide climatology in TIME-GCM display a similar annual variability as the COSMIC observations. This leads to the conclusion that the annual variability of the lunar tide in the ionosphere is not solely due to the occurrence of SSWs. Rather, the annual variability of the lunar tide in the ionosphere is generated by the seasonal variability of the lunar tide at E region altitudes. However, compared to the observations, the ionospheric lunar tide annual variability is weaker in the climatological simulations which is attributed to the occurrence of SSWs during the majority of the years included in the observations. Introducing a SSW into the TIME-GCM simulation leads to an additional enhancement of the lunar tide during Northern Hemisphere winter, increasing the lunar tide annual variability and resulting in an annual variability that is more consistent with the observations. The occurrence of SSWs can therefore potentially bias lunar tide climatologies, and it is important to consider these effects in studies of the lunar tide in the atmosphere and ionosphere.

Towards Observational Astronomy of Jets in Active Galaxies from General Relativistic Magnetohydrodynamic Simulations

NASA Astrophysics Data System (ADS)

Anantua, Richard; Roger Blandford, Jonathan McKinney and Alexander Tchekhovskoy

2016-01-01

We carry out the process of "observing" simulations of active galactic nuclei (AGN) with relativistic jets (hereafter called jet/accretion disk/black hole (JAB) systems) from ray tracing between image plane and source to convolving the resulting images with a point spread function. Images are generated at arbitrary observer angle relative to the black hole spin axis by implementing spatial and temporal interpolation of conserved magnetohydrodynamic flow quantities from a time series of output datablocks from fully general relativistic 3D simulations. We also describe the evolution of simulations of JAB systems' dynamical and kinematic variables, e.g., velocity shear and momentum density, respectively, and the variation of these variables with respect to observer polar and azimuthal angles. We produce, at frequencies from radio to optical, fixed observer time intensity and polarization maps using various plasma physics motivated prescriptions for the emissivity function of physical quantities from the simulation output, and analyze the corresponding light curves. Our hypothesis is that this approach reproduces observed features of JAB systems such as superluminal bulk flow projections and quasi-periodic oscillations in the light curves more closely than extant stylized analytical models, e.g., cannonball bulk flows. Moreover, our development of user-friendly, versatile C++ routines for processing images of state-of-the-art simulations of JAB systems may afford greater flexibility for observing a wide range of sources from high power BL-Lacs to low power quasars (possibly with the same simulation) without requiring years of observation using multiple telescopes. Advantages of observing simulations instead of observing astrophysical sources directly include: the absence of a diffraction limit, panoramic views of the same object and the ability to freely track features. Light travel time effects become significant for high Lorentz factor and small angles between

The SELGIFS data challenge: generating synthetic observations of CALIFA galaxies from hydrodynamical simulations

NASA Astrophysics Data System (ADS)

Guidi, G.; Casado, J.; Ascasibar, Y.; García-Benito, R.; Galbany, L.; Sánchez-Blázquez, P.; Sánchez, S. F.; Rosales-Ortega, F. F.; Scannapieco, C.

2018-06-01

In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and scattering. The results of the radiative transfer simulations have been post-processed to reproduce the main properties of the CALIFA V500 and V1200 observational setups. The data has been further formatted to mimic the CALIFA survey in terms of field of view size, spectral range and sampling. We have included the effect of the spatial and spectral Point Spread Functions affecting CALIFA observations, and added detector noise after characterizing it on a sample of 367 galaxies. The simulated datacubes are suited to be analysed by the same algorithms used on real IFS data. In order to provide a benchmark to compare the results obtained applying IFS observational techniques to our synthetic datacubes, and test the calibration and accuracy of the analysis tools, we have computed the spatially-resolved properties of the simulations. Hence, we provide maps derived directly from the hydrodynamical snapshots or the noiseless spectra, in a way that is consistent with the values recovered by the observational analysis algorithms. Both the synthetic observations and the product datacubes are public and can be found in the collaboration website http://astro.ft.uam.es/selgifs/data_challenge/.

Charaterizing the O+ ion plume from Hybrid simulations: comparison to MAVEN observations

NASA Astrophysics Data System (ADS)

Modolo, R.; Leblanc, F.; Chaufray, J. Y.; Leclercq, L.; Esteban-Hernandez, R.; Curry, S.; Dong, Y.; Brain, D. A.; Bowers, C.; Luhmann, J. G.; McFadden, J. P.; Halekas, J. S.; Espley, J. R.; Connerney, J. E. P.; Jakosky, B. M.

2015-12-01

MAVEN observations show a substantial plume-like distribution of escaping ions from the Martian atmosphere. It represents an important ion escape channel with large fluxes (Brain et al, 2015; Dong et al, 2015, Curry et al, 2015). Such structure is organized by the solar wind convection electric field and it is located in the MSE northward hemisphere. Global hybrid models (eg Modolo et al, 2005, 2012; Kallio et al, 2006; Brecht et al, 2006) reproduce nicely this plume. To further characterize this population, hybrid simulations have been performed with upstream solar wind conditions observed by MAVEN. Simulation results along the spacecraft track present signatures of high energetic O+ ions similar to MAVEN measurements. Comparison of simulated 3D distribution functions of this population are compared to STATIC and SWIA observations. Moreover a comparison of hybrid results with statistical ion fluxes maps derived from MAVEN (Dong et al, 2015; Brain et al, 2015) have been conducted and a reasonable agreement is found .

Properties of the circumgalactic medium in simulations compared to observations

NASA Astrophysics Data System (ADS)

Machado, R. E. G.; Tissera, P. B.; Lima Neto, G. B.; Sodré, L.

2018-01-01

Context. Galaxies are surrounded by extended gaseous halos that store significant fractions of chemical elements. These are syntethized by the stellar populations and later ejected into the circumgalactic medium (CGM) by different mechanism, of which supernova feedback is considered one of the most relevant. Aims: We aim to explore the properties of this metal reservoir surrounding star-forming galaxies in a cosmological context aiming to investigate the chemical loop between galaxies and their CGM, and the ability of the subgrid models to reproduce observational results. Methods: Using cosmological hydrodynamical simulations, we have analysed the gas-phase chemical contents of galaxies with stellar masses in the range 109-1011 M⊙. We estimated the fractions of metals stored in the different CGM phases, and the predicted O VI and Si III column densities within the virial radius. Results: We find roughly 107 M⊙ of oxygen in the CGM of simulated galaxies having M⋆ 1010 M⊙, in fair agreement with the lower limits imposed by observations. The Moxy is found to correlate with M⋆, at odds with current observational trends but in agreement with other numerical results. The estimated profiles of O VI column density reveal a substantial shortage of that ion, whereas Si III, which probes the cool phase, is overpredicted. Nevertheless, the radial dependences of both ions follow the respective observed profiles. The analysis of the relative contributions of both ions from the hot, warm and cool phases suggests that the warm gas (105 K < T < 106 K) should be more abundant in order to bridge the mismatch with the observations, or alternatively, that more metals should be stored in this gas-phase. These discrepancies provide important information to improve the subgrid physics models. Our findings show clearly the importance of tracking more than one chemical element and the difficulty of simultaneously satisfying the observables that trace the circumgalactic gas at

Percolation analyses of observed and simulated galaxy clustering

NASA Astrophysics Data System (ADS)

Bhavsar, S. P.; Barrow, J. D.

1983-11-01

A percolation cluster analysis is performed on equivalent regions of the CFA redshift survey of galaxies and the 4000 body simulations of gravitational clustering made by Aarseth, Gott and Turner (1979). The observed and simulated percolation properties are compared and, unlike correlation and multiplicity function analyses, favour high density (Omega = 1) models with n = - 1 initial data. The present results show that the three-dimensional data are consistent with the degree of filamentary structure present in isothermal models of galaxy formation at the level of percolation analysis. It is also found that the percolation structure of the CFA data is a function of depth. Percolation structure does not appear to be a sensitive probe of intrinsic filamentary structure.

Magnetoacoustic Wave Energy from Numerical Simulations of an Observed Sunspot Umbra

NASA Astrophysics Data System (ADS)

Felipe, T.; Khomenko, E.; Collados, M.

2011-07-01

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through three-dimensional MHD numerical simulations. A magnetostatic sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere, introducing the fluctuations measured with the Si I λ10827 line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the He I λ10830 line, the Ca II H core, and the Fe I blends in the wings of the Ca II H line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pairs of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magnetoacoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost Si I λ10827 line in the form of slow magnetoacoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.

Ionospheric Simulation System for Satellite Observations and Global Assimilative Modeling Experiments (ISOGAME)

NASA Technical Reports Server (NTRS)

Pi, Xiaoqing; Mannucci, Anthony J.; Verkhoglyadova, Olga P.; Stephens, Philip; Wilson, Brian D.; Akopian, Vardan; Komjathy, Attila; Lijima, Byron A.

2013-01-01

ISOGAME is designed and developed to assess quantitatively the impact of new observation systems on the capability of imaging and modeling the ionosphere. With ISOGAME, one can perform observation system simulation experiments (OSSEs). A typical OSSE using ISOGAME would involve: (1) simulating various ionospheric conditions on global scales; (2) simulating ionospheric measurements made from a constellation of low-Earth-orbiters (LEOs), particularly Global Navigation Satellite System (GNSS) radio occultation data, and from ground-based global GNSS networks; (3) conducting ionospheric data assimilation experiments with the Global Assimilative Ionospheric Model (GAIM); and (4) analyzing modeling results with visualization tools. ISOGAME can provide quantitative assessment of the accuracy of assimilative modeling with the interested observation system. Other observation systems besides those based on GNSS are also possible to analyze. The system is composed of a suite of software that combines the GAIM, including a 4D first-principles ionospheric model and data assimilation modules, an Internal Reference Ionosphere (IRI) model that has been developed by international ionospheric research communities, observation simulator, visualization software, and orbit design, simulation, and optimization software. The core GAIM model used in ISOGAME is based on the GAIM++ code (written in C++) that includes a new high-fidelity geomagnetic field representation (multi-dipole). New visualization tools and analysis algorithms for the OSSEs are now part of ISOGAME.

Observer roles that optimise learning in healthcare simulation education: a systematic review.

PubMed

O'Regan, Stephanie; Molloy, Elizabeth; Watterson, Leonie; Nestel, Debra

2016-01-01

Simulation is widely used in health professional education. The convention that learners are actively involved may limit access to this educational method. The aim of this paper is to review the evidence for learning methods that employ directed observation as an alternative to hands-on participation in scenario-based simulation training. We sought studies that included either direct comparison of the learning outcomes of observers with those of active participants or identified factors important for the engagement of observers in simulation. We systematically searched health and education databases and reviewed journals and bibliographies for studies investigating or referring to observer roles in simulation using mannequins, simulated patients or role play simulations. A quality framework was used to rate the studies. We sought studies that included either direct comparison of the learning outcomes of observers with those of active participants or identified factors important for the engagement of observers in simulation. We systematically searched health and education databases and reviewed journals and bibliographies for studies investigating or referring to observer roles in simulation using mannequins, simulated patients or role play simulations. A quality framework was used to rate the studies. Nine studies met the inclusion criteria. Five studies suggest learning outcomes in observer roles are as good or better than hands-on roles in simulation. Four studies document learner satisfaction in observer roles. Five studies used a tool to guide observers. Eight studies involved observers in the debrief. Learning and satisfaction in observer roles is closely associated with observer tools, learner engagement, role clarity and contribution to the debrief. Learners that valued observer roles described them as affording an overarching view, examination of details from a distance, and meaningful feedback during the debrief. Learners who did not value observer roles

Temporal Variability of Observed and Simulated Hyperspectral Earth Reflectance

NASA Technical Reports Server (NTRS)

Roberts, Yolanda; Pilewskie, Peter; Kindel, Bruce; Feldman, Daniel; Collins, William D.

2012-01-01

The Climate Absolute Radiance and Refractivity Observatory (CLARREO) is a climate observation system designed to study Earth's climate variability with unprecedented absolute radiometric accuracy and SI traceability. Observation System Simulation Experiments (OSSEs) were developed using GCM output and MODTRAN to simulate CLARREO reflectance measurements during the 21st century as a design tool for the CLARREO hyperspectral shortwave imager. With OSSE simulations of hyperspectral reflectance, Feldman et al. [2011a,b] found that shortwave reflectance is able to detect changes in climate variables during the 21st century and improve time-to-detection compared to broadband measurements. The OSSE has been a powerful tool in the design of the CLARREO imager and for understanding the effect of climate change on the spectral variability of reflectance, but it is important to evaluate how well the OSSE simulates the Earth's present-day spectral variability. For this evaluation we have used hyperspectral reflectance measurements from the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY), a shortwave spectrometer that was operational between March 2002 and April 2012. To study the spectral variability of SCIAMACHY-measured and OSSE-simulated reflectance, we used principal component analysis (PCA), a spectral decomposition technique that identifies dominant modes of variability in a multivariate data set. Using quantitative comparisons of the OSSE and SCIAMACHY PCs, we have quantified how well the OSSE captures the spectral variability of Earth?s climate system at the beginning of the 21st century relative to SCIAMACHY measurements. These results showed that the OSSE and SCIAMACHY data sets share over 99% of their total variance in 2004. Using the PCs and the temporally distributed reflectance spectra projected onto the PCs (PC scores), we can study the temporal variability of the observed and simulated reflectance spectra. Multivariate time

Interactive visualization of numerical simulation results: A tool for mission planning and data analysis

NASA Technical Reports Server (NTRS)

Berchem, J.; Raeder, J.; Walker, R. J.; Ashour-Abdalla, M.

1995-01-01

We report on the development of an interactive system for visualizing and analyzing numerical simulation results. This system is based on visualization modules which use the Application Visualization System (AVS) and the NCAR graphics packages. Examples from recent simulations are presented to illustrate how these modules can be used for displaying and manipulating simulation results to facilitate their comparison with phenomenological model results and observations.

An Observing System Simulation Experiment Approach to Meteorological Network Assessment

NASA Astrophysics Data System (ADS)

Abbasnezhadi, K.; Rasmussen, P. F.; Stadnyk, T.; Boluwade, A.

2016-12-01

A proper knowledge of the spatiotemporal distribution of rainfall is important in order to conduct a mindful investigation of water movement and storage throughout a catchment. Currently, the most accurate precipitation information available for the remote Boreal ecozones of northern Manitoba is coming from the Canadian Precipitation Analysis (CaPA) data assimilation system. Throughout the Churchill River Basin (CRB), CaPA still does not have the proper skill due to the limited number of weather stations. A new approach to experimental network design was investigated based on the concept of Observing System Simulation Experiment (OSSE). The OSSE-based network assessment procedure which simulates the CaPA system provides a scientific and hydrologically significant tool to assess the sensitivity of CaPA precipitation analysis to observation network density throughout the CRB. To simulate CaPA system, synthetic background and station data were simulated, respectively, by adding spatially uncorrelated and correlated Gaussian noises to an assumingly true daily weather field synthesized by a gridded precipitation generator which simulates CaPA data. Given the true reference field on one hand, and a set of pseudo-CaPA analyses associated with different network realizations on the other hand, a WATFLOOD hydrological model was employed to compare the modeled runoff. The simulations showed that as network density increases, the accuracy of CaPA precipitation products improves up to a certain limit beyond which adding more stations to the network does not result in further accuracy.

Simulating observations with HARMONI: the integral field spectrograph for the European Extremely Large Telescope

NASA Astrophysics Data System (ADS)

Zieleniewski, Simon; Thatte, Niranjan; Kendrew, Sarah; Houghton, Ryan; Tecza, Matthias; Clarke, Fraser; Fusco, Thierry; Swinbank, Mark

2014-07-01

With the next generation of extremely large telescopes commencing construction, there is an urgent need for detailed quantitative predictions of the scientific observations that these new telescopes will enable. Most of these new telescopes will have adaptive optics fully integrated with the telescope itself, allowing unprecedented spatial resolution combined with enormous sensitivity. However, the adaptive optics point spread function will be strongly wavelength dependent, requiring detailed simulations that accurately model these variations. We have developed a simulation pipeline for the HARMONI integral field spectrograph, a first light instrument for the European Extremely Large Telescope. The simulator takes high-resolution input data-cubes of astrophysical objects and processes them with accurate atmospheric, telescope and instrumental effects, to produce mock observed cubes for chosen observing parameters. The output cubes represent the result of a perfect data reduc- tion process, enabling a detailed analysis and comparison between input and output, showcasing HARMONI's capabilities. The simulations utilise a detailed knowledge of the telescope's wavelength dependent adaptive op- tics point spread function. We discuss the simulation pipeline and present an early example of the pipeline functionality for simulating observations of high redshift galaxies.

The Cloud Feedback Model Intercomparison Project Observational Simulator Package: Version 2

NASA Astrophysics Data System (ADS)

Swales, Dustin J.; Pincus, Robert; Bodas-Salcedo, Alejandro

2018-01-01

The Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) gathers together a collection of observation proxies or satellite simulators that translate model-simulated cloud properties to synthetic observations as would be obtained by a range of satellite observing systems. This paper introduces COSP2, an evolution focusing on more explicit and consistent separation between host model, coupling infrastructure, and individual observing proxies. Revisions also enhance flexibility by allowing for model-specific representation of sub-grid-scale cloudiness, provide greater clarity by clearly separating tasks, support greater use of shared code and data including shared inputs across simulators, and follow more uniform software standards to simplify implementation across a wide range of platforms. The complete package including a testing suite is freely available.

Simulating Cosmic Reionization and Its Observable Consequences

NASA Astrophysics Data System (ADS)

Shapiro, Paul

2017-01-01

I summarize recent progress in modelling the epoch of reionization by large- scale simulations of cosmic structure formation, radiative transfer and their interplay, which trace the ionization fronts that swept across the IGM, to predict observable signatures. Reionization by starlight from early galaxies affected their evolution, impacting reionization, itself, and imprinting the galaxies with a memory of reionization. Star formation suppression, e.g., may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. I describe CoDa (''Cosmic Dawn''), the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all the atomic-cooling galactic halos in that volume. A 90 Mpc box was simulated from a constrained realization of primordial fluctuations, chosen to reproduce present-day features of the Local Group, including the Milky Way and M31, and the local universe beyond, including the Virgo cluster. The new RAMSES-CUDATON hybrid CPU-GPU code took 11 days to perform this simulation on the Titan supercomputer at Oak Ridge National Laboratory, with 4096-cubed N-body particles for the dark matter and 4096-cubed cells for the atomic gas and ionizing radiation.

Comparison of Magnetospheric Magnetic Field Variations at Quasi-Zenith Orbit Based on Michibiki Observation and REPPU Global MHD Simulation

NASA Astrophysics Data System (ADS)

Kubota, Y.; Nagatsuma, T.; Den, M.; Nakamizo, A.; Matsumoto, H.; Tanaka, T.

2017-12-01

We are developing a numerical simulator for future space weather forecast using magnetosphere-ionosphere coupling global MHD simulation called REPPU (REProduce Plasma Universe) code. We investigate the validity of the MHD simulation result as compared with observation. In this study we simulate some events including both quiet and disturbed geomagnetic conditions using OMNIWeb solar wind data. The simulation results are compared with magnetic field observations from Michibiki satellite, which is on the quasi-zenith orbit (QZO). In quiet geomagnetic condition, magnetic field variations at QZO obtained from simulation results have good consistency as compared correspondence with those from Michibiki observation. In disturbed geomagnetic condition in which the Dst < -20 nT, however, V component of magnetic field variations from simulation results tend to deviate from observations especially at the night side. We consider that this deviation during disturbed geomagnetic condition might be due to tail and/or ring current enhancement which is already suggested by many other MHD simulation studies as compared with the magnetic field observation at geosynchronous orbit. In this presentation, we will discuss the cause of this discrepancy in more detail with studying the relationship between the magnetic field deviation and some parameters such as Dst and solar wind.

"Observing" the Circumnuclear Stars and Gas in Disk Galaxy Simulations

NASA Astrophysics Data System (ADS)

Cook, Angela; Hicks, Erin K. S.

2018-06-01

We present simulations based on theoretical models of common disk processes designed to represent potential inflow observed within the central 500 pc of local Seyfert galaxies. Mock observations of these n-body plus smoothed particle hydrodynamical simulations provide the conceptual framework in which to identify the driving inflow mechanism, for example nuclear bars, and to quantify to the inflow based on observable properties. From these mock observations the azimuthal average of the flux distribution, velocity dispersion, and velocity of both the stars and interstellar medium on scales of 50pc have been measured at a range of inclinations angles. A comparison of the simulated disk galaxies with these observed azimuthal averages in 40 Seyfert galaxies measured as part of the KONA (Keck OSIRIS Nearby AGN) survey will be presented.

Observability of ionospheric space-time structure with ISR: A simulation study

NASA Astrophysics Data System (ADS)

Swoboda, John; Semeter, Joshua; Zettergren, Matthew; Erickson, Philip J.

2017-02-01

The sources of error from electronically steerable array (ESA) incoherent scatter radar (ISR) systems are investigated both theoretically and with use of an open-source ISR simulator, developed by the authors, called Simulator for ISR (SimISR). The main sources of error incorporated in the simulator include statistical uncertainty, which arises due to nature of the measurement mechanism and the inherent space-time ambiguity from the sensor. SimISR can take a field of plasma parameters, parameterized by time and space, and create simulated ISR data at the scattered electric field (i.e., complex receiver voltage) level, subsequently processing these data to show possible reconstructions of the original parameter field. To demonstrate general utility, we show a number of simulation examples, with two cases using data from a self-consistent multifluid transport model. Results highlight the significant influence of the forward model of the ISR process and the resulting statistical uncertainty on plasma parameter measurements and the core experiment design trade-offs that must be made when planning observations. These conclusions further underscore the utility of this class of measurement simulator as a design tool for more optimal experiment design efforts using flexible ESA class ISR systems.

Comparing Simulations and Observations of Galaxy Evolution: Methods for Constraining the Nature of Stellar Feedback

NASA Astrophysics Data System (ADS)

Hummels, Cameron

Computational hydrodynamical simulations are a very useful tool for understanding how galaxies form and evolve over cosmological timescales not easily revealed through observations. However, they are only useful if they reproduce the sorts of galaxies that we see in the real universe. One of the ways in which simulations of this sort tend to fail is in the prescription of stellar feedback, the process by which nascent stars return material and energy to their immediate environments. Careful treatment of this interaction in subgrid models, so-called because they operate on scales below the resolution of the simulation, is crucial for the development of realistic galaxy models. Equally important is developing effective methods for comparing simulation data against observations to ensure galaxy models which mimic reality and inform us about natural phenomena. This thesis examines the formation and evolution of galaxies and the observable characteristics of the resulting systems. We employ extensive use of cosmological hydrodynamical simulations in order to simulate and interpret the evolution of massive spiral galaxies like our own Milky Way. First, we create a method for producing synthetic photometric images of grid-based hydrodynamical models for use in a direct comparison against observations in a variety of filter bands. We apply this method to a simulation of a cluster of galaxies to investigate the nature of the red-sequence/blue-cloud dichotomy in the galaxy color-magnitude diagram. Second, we implement several subgrid models governing the complex behavior of gas and stars on small scales in our galaxy models. Several numerical simulations are conducted with similar initial conditions, where we systematically vary the subgrid models, afterward assessing their efficacy through comparisons of their internal kinematics with observed systems. Third, we generate an additional method to compare observations with simulations, focusing on the tenuous circumgalactic

Aerosol midlatitude cyclone indirect effects in observations and high-resolution simulations

NASA Astrophysics Data System (ADS)

McCoy, Daniel T.; Field, Paul R.; Schmidt, Anja; Grosvenor, Daniel P.; Bender, Frida A.-M.; Shipway, Ben J.; Hill, Adrian A.; Wilkinson, Jonathan M.; Elsaesser, Gregory S.

2018-04-01

Aerosol-cloud interactions are a major source of uncertainty in inferring the climate sensitivity from the observational record of temperature. The adjustment of clouds to aerosol is a poorly constrained aspect of these aerosol-cloud interactions. Here, we examine the response of midlatitude cyclone cloud properties to a change in cloud droplet number concentration (CDNC). Idealized experiments in high-resolution, convection-permitting global aquaplanet simulations with constant CDNC are compared to 13 years of remote-sensing observations. Observations and idealized aquaplanet simulations agree that increased warm conveyor belt (WCB) moisture flux into cyclones is consistent with higher cyclone liquid water path (CLWP). When CDNC is increased a larger LWP is needed to give the same rain rate. The LWP adjusts to allow the rain rate to be equal to the moisture flux into the cyclone along the WCB. This results in an increased CLWP for higher CDNC at a fixed WCB moisture flux in both observations and simulations. If observed cyclones in the top and bottom tercile of CDNC are contrasted it is found that they have not only higher CLWP but also cloud cover and albedo. The difference in cyclone albedo between the cyclones in the top and bottom third of CDNC is observed by CERES to be between 0.018 and 0.032, which is consistent with a 4.6-8.3 Wm-2 in-cyclone enhancement in upwelling shortwave when scaled by annual-mean insolation. Based on a regression model to observed cyclone properties, roughly 60 % of the observed variability in CLWP can be explained by CDNC and WCB moisture flux.

Observations and Simulations of Formation of Broad Plasma Depletions Through Merging Process

NASA Technical Reports Server (NTRS)

Huang, Chao-Song; Retterer, J. M.; Beaujardiere, O. De La; Roddy, P. A.; Hunton, D.E.; Ballenthin, J. O.; Pfaff, Robert F.

2012-01-01

Broad plasma depletions in the equatorial ionosphere near dawn are region in which the plasma density is reduced by 1-3 orders of magnitude over thousands of kilometers in longitude. This phenomenon is observed repeatedly by the Communication/Navigation Outage Forecasting System (C/NOFS) satellite during deep solar minimum. The plasma flow inside the depletion region can be strongly upward. The possible causal mechanism for the formation of broad plasma depletions is that the broad depletions result from merging of multiple equatorial plasma bubbles. The purpose of this study is to demonstrate the feasibility of the merging mechanism with new observations and simulations. We present C/NOFS observations for two cases. A series of plasma bubbles is first detected by C/NOFS over a longitudinal range of 3300-3800 km around midnight. Each of the individual bubbles has a typical width of approx 100 km in longitude, and the upward ion drift velocity inside the bubbles is 200-400 m/s. The plasma bubbles rotate with the Earth to the dawn sector and become broad plasma depletions. The observations clearly show the evolution from multiple plasma bubbles to broad depletions. Large upward plasma flow occurs inside the depletion region over 3800 km in longitude and exists for approx 5 h. We also present the numerical simulations of bubble merging with the physics-based low-latitude ionospheric model. It is found that two separate plasma bubbles join together and form a single, wider bubble. The simulations show that the merging process of plasma bubbles can indeed occur in incompressible ionospheric plasma. The simulation results support the merging mechanism for the formation of broad plasma depletions.

Understanding the central kinematics of globular clusters with simulated integrated-light IFU observations

NASA Astrophysics Data System (ADS)

Bianchini, Paolo; Norris, Mark A.; van de Ven, Glenn; Schinnerer, Eva

2015-10-01

The detection of intermediate-mass black holes in the centres of globular clusters is highly controversial, as complementary observational methods often deliver significantly different results. In order to understand these discrepancies, we develop a procedure to simulate integral field unit (IFU) observations of globular clusters: Simulating IFU Star Cluster Observations (SISCO). The inputs of our software are realistic dynamical models of globular clusters that are then converted in a spectral data cube. We apply SISCO to Monte Carlo cluster simulations with a realistic number of stars and concentrations. Using independent realizations of a given simulation we are able to quantify the stochasticity intrinsic to the problem of observing a partially resolved stellar population with integrated-light spectroscopy. We show that the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in the velocity dispersion measurements up to ≃40 per cent around the expected value, preventing any sound assessment of the central kinematic and a sensible interpretation of the presence/absence of an intermediate-mass black hole. Moreover, we illustrate that, in our mock IFU observations, the average kinematic tracer has a mass of ≃0.75 M⊙, only slightly lower than the mass of the typical stars examined in studies of resolved line-of-sight velocities of giant stars. Finally, in order to recover unbiased kinematic measurements we test different masking techniques that allow us to remove the spaxels dominated by bright stars, bringing the scatter down to a level of only a few per cent. The application of SISCO will allow us to investigate state-of-the-art simulations as realistic observations.

Comparing models of star formation simulating observed interacting galaxies

NASA Astrophysics Data System (ADS)

Quiroga, L. F.; Muñoz-Cuartas, J. C.; Rodrigues, I.

2017-07-01

In this work, we make a comparison between different models of star formation to reproduce observed interacting galaxies. We use observational data to model the evolution of a pair of galaxies undergoing a minor merger. Minor mergers represent situations weakly deviated from the equilibrium configuration but significant changes in star fomation (SF) efficiency can take place, then, minor mergers provide an unique scene to study SF in galaxies in a realistic but yet simple way. Reproducing observed systems also give us the opportunity to compare the results of the simulations with observations, which at the end can be used as probes to characterize the models of SF implemented in the comparison. In this work we compare two different star formation recipes implemented in Gadget3 and GIZMO codes. Both codes share the same numerical background, and differences arise mainly in the star formation recipe they use. We use observations from Pico dos Días and GEMINI telescopes and show how we use observational data of the interacting pair in AM2229-735 to characterize the interacting pair. Later we use this information to simulate the evolution of the system to finally reproduce the observations: Mass distribution, morphology and main features of the merger-induced star formation burst. We show that both methods manage to reproduce roughly the star formation activity. We show, through a careful study, that resolution plays a major role in the reproducibility of the system. In that sense, star formation recipe implemented in GIZMO code has shown a more robust performance. Acknowledgements: This work is supported by Colciencias, Doctorado Nacional - 617 program.

Evaluating performance of risk identification methods through a large-scale simulation of observational data.

PubMed

Ryan, Patrick B; Schuemie, Martijn J

2013-10-01

There has been only limited evaluation of statistical methods for identifying safety risks of drug exposure in observational healthcare data. Simulations can support empirical evaluation, but have not been shown to adequately model the real-world phenomena that challenge observational analyses. To design and evaluate a probabilistic framework (OSIM2) for generating simulated observational healthcare data, and to use this data for evaluating the performance of methods in identifying associations between drug exposure and health outcomes of interest. Seven observational designs, including case-control, cohort, self-controlled case series, and self-controlled cohort design were applied to 399 drug-outcome scenarios in 6 simulated datasets with no effect and injected relative risks of 1.25, 1.5, 2, 4, and 10, respectively. Longitudinal data for 10 million simulated patients were generated using a model derived from an administrative claims database, with associated demographics, periods of drug exposure derived from pharmacy dispensings, and medical conditions derived from diagnoses on medical claims. Simulation validation was performed through descriptive comparison with real source data. Method performance was evaluated using Area Under ROC Curve (AUC), bias, and mean squared error. OSIM2 replicates prevalence and types of confounding observed in real claims data. When simulated data are injected with relative risks (RR) ≥ 2, all designs have good predictive accuracy (AUC > 0.90), but when RR < 2, no methods achieve 100 % predictions. Each method exhibits a different bias profile, which changes with the effect size. OSIM2 can support methodological research. Results from simulation suggest method operating characteristics are far from nominal properties.

Observational Simulation of Icing in Extreme Weather Conditions

NASA Astrophysics Data System (ADS)

Gultepe, Ismail; Heymsfield, Andrew; Agelin-Chaab, Martin; Komar, John; Elfstrom, Garry; Baumgardner, Darrel

2017-04-01

Observations and prediction of icing in extreme weather conditions are important for aviation, transportation, and shipping applications, and icing adversely affects the economy. Icing environments can be studied either in the outdoor atmosphere or in the laboratory. There have been several aircraft based in-situ studies related to weather conditions affecting aviation operations, transportation, and marine shipping that includes icing, wind, and turbulence. However, studying severe weather conditions from aircraft observations are limited due to safety and sampling issues, instrumental uncertainties, and even the possibility of aircraft producing its own physical and dynamical effects. Remote sensing based techniques (e.g. retrieval techniques) for studying severe weather conditions represent usually a volume that cannot characterize the important scales and also represents indirect observations. Therefore, laboratory simulations of atmospheric processes can help us better understand the interactions among microphysical and dynamical processes. The Climatic Wind Tunnel (CWT) in ACE at the University of Ontario Institute of Technology (UOIT) has a large semi-open jet test chamber with flow area 7-13 m2 that can precisely control temperatures down to -40°C, and up to 250 km hr-1 wind speeds, for heavy or dry snow conditions with low visibility, similar to ones observed in the Arctic and cold climate regions, or at high altitude aeronautical conditions. In this study, the ACE CWT employed a spray nozzle array suspended in its settling chamber and fed by pressurized water, creating various particle sizes from a few microns up to mm size range. This array, together with cold temperature and high wind speed, enabled simulation of severe weather conditions, including icing, visibility, strong wind and turbulence, ice fog and frost, freezing fog, heavy snow and blizzard conditions. In this study, the test results will be summarized, and their application to aircraft

Observation and simulation of the ionosphere disturbance waves triggered by rocket exhausts

NASA Astrophysics Data System (ADS)

Lin, Charles C. H.; Chen, Chia-Hung; Matsumura, Mitsuru; Lin, Jia-Ting; Kakinami, Yoshihiro

2017-08-01

Observations and theoretical modeling of the ionospheric disturbance waves generated by rocket launches are investigated. During the rocket passage, time rate change of total electron content (rTEC) enhancement with the V-shape shock wave signature is commonly observed, followed by acoustic wave disturbances and region of negative rTEC centered along the trajectory. Ten to fifteen min after the rocket passage, delayed disturbance waves appeared and propagated along direction normal to the V-shape wavefronts. These observation features appeared most prominently in the 2016 North Korea rocket launch showing a very distinct V-shape rTEC enhancement over enormous areas along the southeast flight trajectory despite that it was also appeared in the 2009 North Korea rocket launch with the eastward flight trajectory. Numerical simulations using the physical-based nonlinear and nonhydrostatic coupled model of neutral atmosphere and ionosphere reproduce promised results in qualitative agreement with the characteristics of ionospheric disturbance waves observed in the 2009 event by considering the released energy of the rocket exhaust as the disturbance source. Simulations reproduce the shock wave signature of electron density enhancement, acoustic wave disturbances, the electron density depletion due to the rocket-induced pressure bulge, and the delayed disturbance waves. The pressure bulge results in outward neutral wind flows carrying neutrals and plasma away from it and leading to electron density depletions. Simulations further show, for the first time, that the delayed disturbance waves are produced by the surface reflection of the earlier arrival acoustic wave disturbances.

Annual asymmetry in thermospheric density: Observations and simulations

NASA Astrophysics Data System (ADS)

Lei, Jiuhou; Dou, Xiankang; Burns, Alan; Wang, Wenbin; Luan, Xiaoli; Zeng, Zhen; Xu, Jiyao

2013-05-01

In this paper, the Challenging Minisatellite Payload (CHAMP) and Gravity Recovery and Climate Experiment (GRACE) observations during 2002-2010 are utilized to study the variation of the annual asymmetry in thermospheric density at 400 km under low solar activity condition (F10.7 = 80) based on the method of empirical orthogonal functions (EOFs). The derived asymmetry index (AI) in thermospheric density from the EOF analysis shows a strong latitudinal variation at night but varies a little with latitudes in daytime. Moreover, it exhibits a terdiurnal tidal signature at low to middle latitudes. The global mean value of the AI is 0.191, indicating that a 47% difference in thermosphere between the December and June solstices in the global average. In addition, the NCAR Thermosphere-Ionosphere Electrodynamics Global Circulation Model (TIEGCM) is used to explore the possible mechanisms responsible for the observed annual asymmetry in thermospheric density. It is found that the standard simulations give a lower AI and also a weaker day-to-night difference. The simulated AI shows a semidiurnal pattern in the equatorial and low-latitude regions in contrast with the terdiurnal tide signature seen in the observed AI. The daily mean AI obtained from the simulation is 0.125, corresponding to a 29% December-to-June difference in thermospheric density at 400 km. Further sensitivity simulations demonstrated that the effect of the varying Sun-Earth distance between the December and June solstices is the main process responsible for the annual asymmetry in thermospheric density, while the magnetic field configuration and tides from the lower atmosphere contribute to the temporal and spatial variations of the AI. Specifically, the simulations show that the Sun-Earth distance effect explains 93% of the difference in thermospheric density between December and June, which is mainly associated with the corresponding changes in neutral temperature. However, our calculation from the

Reconciling Simulated and Observed Views of Clouds: MODIS, ISCCP, and the Limits of Instrument Simulators in Climate Models

NASA Technical Reports Server (NTRS)

Pincus, Robert; Platnick, Steven E.; Ackerman, Steve; Hemler, Richard; Hofmann, Patrick

2011-01-01

The properties of clouds that may be observed by satellite instruments, such as optical depth and cloud top pressure, are only loosely related to the way clouds are represented in models of the atmosphere. One way to bridge this gap is through "instrument simulators," diagnostic tools that map the model representation to synthetic observations so that differences between simulator output and observations can be interpreted unambiguously as model error. But simulators may themselves be restricted by limited information available from the host model or by internal assumptions. This work examines the extent to which instrument simulators are able to capture essential differences between MODIS and ISCCP, two similar but independent estimates of cloud properties. We focus on the stark differences between MODIS and ISCCP observations of total cloudiness and the distribution of cloud optical thickness can be traced to different approaches to marginal pixels, which MODIS excludes and ISCCP treats as homogeneous. These pixels, which likely contain broken clouds, cover about 15% of the planet and contain almost all of the optically thinnest clouds observed by either instrument. Instrument simulators can not reproduce these differences because the host model does not consider unresolved spatial scales and so can not produce broken pixels. Nonetheless, MODIS and ISCCP observation are consistent for all but the optically-thinnest clouds, and models can be robustly evaluated using instrument simulators by excluding ambiguous observations.

Reliable results from stochastic simulation models

Treesearch

Donald L., Jr. Gochenour; Leonard R. Johnson

1973-01-01

Development of a computer simulation model is usually done without fully considering how long the model should run (e.g. computer time) before the results are reliable. However construction of confidence intervals (CI) about critical output parameters from the simulation model makes it possible to determine the point where model results are reliable. If the results are...

Evaluating Simulated Tropical Convective Cores using HAIC-HIWC Microphysics and Dynamics Observations

NASA Astrophysics Data System (ADS)

Stanford, M.; Varble, A.; Zipser, E. J.; Strapp, J. W.; Leroy, D.; Schwarzenboeck, A.; Korolev, A.; Potts, R.

2016-12-01

A model intercomparison study is conducted to identify biases in simulated tropical convective core microphysical properties using two popular bulk parameterization schemes (Thompson and Morrison) and the Fast Spectral Bin Microphysics (FSBM) scheme. In-situ aircraft measurements of total condensed water content (TWC) and particle size distributions are compared with output from high-resolution WRF simulations of 4 mesoscale convective system (MCS) cases during the High Altitude Ice Crystals-High Ice Water Content (HAIC-HIWC) field campaign conducted in Darwin, Australia in 2014 and Cayenne, French Guiana in 2015. Observations of TWC collected using an isokinetic evaporator probe (IKP) optimized for high IWC measurements in conjunction with particle image processing from two optical array probes aboard the Falcon-20 research aircraft were used to constrain mass-size relationships in the observational dataset. Hydrometeor mass size distributions are compared between retrievals and simulations providing insight into the well-known high bias in simulated convective radar reflectivity. For TWC > 1 g m-3 between -10 and -40°C, simulations generally produce significantly greater median mass diameters (MMDs). Observations indicate that a sharp particle size mode occurs at 300 μm for large TWC values (> 2 g m-3) regardless of temperature. All microphysics schemes fail to reproduce this feature, and relative contributions of different hydrometeor species to this size bias vary between schemes. Despite far greater sample sizes, simulations also fail to produce high TWC conditions with very little of the mass contributed by large particles for a range of temperatures, despite such conditions being observed. Considering vapor grown particles alone in comparison with observations fails to correct the bias present in all schemes. Decreasing horizontal resolution from 1 km to 333 m shifts graupel and rain size distributions to slightly smaller sizes, but increased resolution

Model-Observation "Data Cubes" for the DOE Atmospheric Radiation Measurement Program's LES ARM Symbiotic Simulation and Observation (LASSO) Workflow

NASA Astrophysics Data System (ADS)

Vogelmann, A. M.; Gustafson, W. I., Jr.; Toto, T.; Endo, S.; Cheng, X.; Li, Z.; Xiao, H.

2015-12-01

The Department of Energy's Atmospheric Radiation Measurement (ARM) Climate Research Facilities' Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) Workflow is currently being designed to provide output from routine LES to complement its extensive observations. The modeling portion of the LASSO workflow is presented by Gustafson et al., which will initially focus on shallow convection over the ARM megasite in Oklahoma, USA. This presentation describes how the LES output will be combined with observations to construct multi-dimensional and dynamically consistent "data cubes", aimed at providing the best description of the atmospheric state for use in analyses by the community. The megasite observations are used to constrain large-eddy simulations that provide a complete spatial and temporal coverage of observables and, further, the simulations also provide information on processes that cannot be observed. Statistical comparisons of model output with their observables are used to assess the quality of a given simulated realization and its associated uncertainties. A data cube is a model-observation package that provides: (1) metrics of model-observation statistical summaries to assess the simulations and the ensemble spread; (2) statistical summaries of additional model property output that cannot be or are very difficult to observe; and (3) snapshots of the 4-D simulated fields from the integration period. Searchable metrics are provided that characterize the general atmospheric state to assist users in finding cases of interest, such as categorization of daily weather conditions and their specific attributes. The data cubes will be accompanied by tools designed for easy access to cube contents from within the ARM archive and externally, the ability to compare multiple data streams within an event as well as across events, and the ability to use common grids and time sampling, where appropriate.

Evolution of Flow channels and Dipolarization Using THEMIS Observations and Global MHD Simulations

NASA Astrophysics Data System (ADS)

El-Alaoui, M.; McPherron, R. L.; Nishimura, Y.

2017-12-01

We have extensively analyzed a substorm on March 14, 2008 for which we have observations from THEMIS spacecraft located beyond 9 RE near 2100 local time. The available data include an extensive network of all sky cameras and ground magnetometers that establish the times of various auroral and magnetic events. This arrangement provided an excellent data set with which to investigate meso-scale structures in the plasma sheet. We have used a global magnetohydrodynamic simulation to investigate the structure and dynamics of the magnetotail current sheet during this substorm. Both earthward and tailward flows were found in the observations as well as the simulations. The simulation shows that the flow channels follow tortuous paths that are often reflected or deflected before arriving at the inner magnetosphere. The simulation shows a sequence of fast flows and dipolarization events similar to what is seen in the data, though not at precisely the same times or locations. We will use our simulation results combined with the observations to investigate the global convection systems and current sheet structure during this event, showing how meso-scale structures fit into the context of the overall tail dynamics during this event. Our study includes determining the location, timing and strength of several current wedges and expansion onsets during an 8-hour interval.

QUANTIFYING OBSERVATIONAL PROJECTION EFFECTS USING MOLECULAR CLOUD SIMULATIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Beaumont, Christopher N.; Offner, Stella S.R.; Shetty, Rahul

2013-11-10

The physical properties of molecular clouds are often measured using spectral-line observations, which provide the only probes of the clouds' velocity structure. It is hard, though, to assess whether and to what extent intensity features in position-position-velocity (PPV) space correspond to 'real' density structures in position-position-position (PPP) space. In this paper, we create synthetic molecular cloud spectral-line maps of simulated molecular clouds, and present a new technique for measuring the reality of individual PPV structures. Using a dendrogram algorithm, we identify hierarchical structures in both PPP and PPV space. Our procedure projects density structures identified in PPP space into correspondingmore » intensity structures in PPV space and then measures the geometric overlap of the projected structures with structures identified from the synthetic observation. The fractional overlap between a PPP and PPV structure quantifies how well the synthetic observation recovers information about the three-dimensional structure. Applying this machinery to a set of synthetic observations of CO isotopes, we measure how well spectral-line measurements recover mass, size, velocity dispersion, and virial parameter for a simulated star-forming region. By disabling various steps of our analysis, we investigate how much opacity, chemistry, and gravity affect measurements of physical properties extracted from PPV cubes. For the simulations used here, which offer a decent, but not perfect, match to the properties of a star-forming region like Perseus, our results suggest that superposition induces a ∼40% uncertainty in masses, sizes, and velocity dispersions derived from {sup 13}CO (J = 1-0). As would be expected, superposition and confusion is worst in regions where the filling factor of emitting material is large. The virial parameter is most affected by superposition, such that estimates of the virial parameter derived from PPV and PPP information typically

An observation planning algorithm applied to multi-objective astronomical observations and its simulation in COSMOS field

NASA Astrophysics Data System (ADS)

Jin, Yi; Gu, Yonggang; Zhai, Chao

2012-09-01

Multi-Object Fiber Spectroscopic sky surveys are now booming, such as LAMOST already built by China, BIGBOSS project put forward by the U.S. Lawrence Berkeley National Lab and GTC (Gran Telescopio Canarias) telescope developed by the United States, Mexico and Spain. They all use or will use this approach and each fiber can be moved within a certain area for one astrology target, so observation planning is particularly important for this Sky Surveys. One observation planning algorithm used in multi-objective astronomical observations is developed. It can avoid the collision and interference between the fiber positioning units in the focal plane during the observation in one field of view, and the interested objects can be ovserved in a limited round with the maximize efficiency. Also, the observation simulation can be made for wide field of view through multi-FOV observation. After the observation planning is built ,the simulation is made in COSMOS field using GTC telescope. Interested galaxies, stars and high-redshift LBG galaxies are selected after the removal of the mask area, which may be bright stars. Then 9 FOV simulation is completed and observation efficiency and fiber utilization ratio for every round are given. Otherwise,allocating a certain number of fibers for background sky, giving different weights for different objects and how to move the FOV to improve the overall observation efficiency are discussed.

GCM simulations of Titan's middle and lower atmosphere and comparison to observations

NASA Astrophysics Data System (ADS)

Lora, Juan M.; Lunine, Jonathan I.; Russell, Joellen L.

2015-04-01

Simulation results are presented from a new general circulation model (GCM) of Titan, the Titan Atmospheric Model (TAM), which couples the Flexible Modeling System (FMS) spectral dynamical core to a suite of external/sub-grid-scale physics. These include a new non-gray radiative transfer module that takes advantage of recent data from Cassini-Huygens, large-scale condensation and quasi-equilibrium moist convection schemes, a surface model with "bucket" hydrology, and boundary layer turbulent diffusion. The model produces a realistic temperature structure from the surface to the lower mesosphere, including a stratopause, as well as satisfactory superrotation. The latter is shown to depend on the dynamical core's ability to build up angular momentum from surface torques. Simulated latitudinal temperature contrasts are adequate, compared to observations, and polar temperature anomalies agree with observations. In the lower atmosphere, the insolation distribution is shown to strongly impact turbulent fluxes, and surface heating is maximum at mid-latitudes. Surface liquids are unstable at mid- and low-latitudes, and quickly migrate poleward. The simulated humidity profile and distribution of surface temperatures, compared to observations, corroborate the prevalence of dry conditions at low latitudes. Polar cloud activity is well represented, though the observed mid-latitude clouds remain somewhat puzzling, and some formation alternatives are suggested.

AOD trends during 2001-2010 from observations and model simulations

NASA Astrophysics Data System (ADS)

Pozzer, A.; de Meij, A.; Yoon, J.; Tost, H.; Georgoulias, A. K.; Astitha, M.

2015-05-01

The aerosol optical depth (AOD) trend between 2001 and 2010 is estimated globally and regionally from observations and results from simulations with the EMAC (ECHAM5/MESSy Atmospheric Chemistry) model. Although interannual variability is applied only to anthropogenic and biomass-burning emissions, the model is able to quantitatively reproduce the AOD trends as observed by the MODIS (Moderate Resolution Imaging Spectroradiometer) satellite sensor, while some discrepancies are found when compared to MISR (Multi-angle Imaging SpectroRadiometer) and SeaWIFS (Sea-viewing Wide Field-of-view Sensor) observations. Thanks to an additional simulation without any change in emissions, it is shown that decreasing AOD trends over the US and Europe are due to the decrease in the emissions, while over the Sahara Desert and the Middle East region, the meteorological changes play a major role. Over Southeast Asia, both meteorology and emissions changes are equally important in defining AOD trends. Additionally, decomposing the regional AOD trends into individual aerosol components reveals that the soluble components are the most dominant contributors to the total AOD, as their influence on the total AOD is enhanced by the aerosol water content.

Summarizing Simulation Results using Causally-relevant States

PubMed Central

Parikh, Nidhi; Marathe, Madhav; Swarup, Samarth

2016-01-01

As increasingly large-scale multiagent simulations are being implemented, new methods are becoming necessary to make sense of the results of these simulations. Even concisely summarizing the results of a given simulation run is a challenge. Here we pose this as the problem of simulation summarization: how to extract the causally-relevant descriptions of the trajectories of the agents in the simulation. We present a simple algorithm to compress agent trajectories through state space by identifying the state transitions which are relevant to determining the distribution of outcomes at the end of the simulation. We present a toy-example to illustrate the working of the algorithm, and then apply it to a complex simulation of a major disaster in an urban area. PMID:28042620

Coronal Mass Ejections and Dimmings: A Comparative Study using MHD Simulations and SDO Observations

NASA Astrophysics Data System (ADS)

Jin, M.; Cheung, C. M. M.; DeRosa, M. L.; Nitta, N.; Schrijver, K.

2017-12-01

Solar coronal dimmings have been observed extensively in the past two decades. Due to their close association with coronal mass ejections (CMEs), there is a critical need to improve our understanding of the physical processes that cause dimmings and determine their relationship with CMEs. In this study, we investigate coronal dimmings by combining simulation and observational efforts. By utilizing a data-driven global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we simulate coronal dimmings resulting from different CME energetics and flux rope configurations. We synthesize the emissions of different EUV spectral bands/lines and compare with SDO/AIA and EVE observations. A detailed analysis of simulation and observation data suggests that although the transient dimming / brightening patterns could relate to plasma heating processes (either by adiabatic compression or reconnection), the long-lasting "core" and "remote" (also known as "secondary") dimmings both originate from regions with open/quasi-open fields and are caused by mass loss process. The mass loss in the remote dimming region is induced by CME-driven shock. Using metrics such as dimming depth, dimming slope, and recovery time, we investigate the relationship between dimmings and CME properties (e.g., CME mass, CME speed) in the simulation. Our result suggests that coronal dimmings encode important information about CME energetics, CME-driven shock properties, and magnetic configuration of erupting flux ropes. We also discuss how our knowledge about solar coronal dimmings could be extended to the study of stellar CMEs, which may prove important for exoplanet atmospheres and habitability but which are currently not observable.

Introduction to Observing System Simulation Experiments (OSSEs)

NASA Technical Reports Server (NTRS)

Prive, Nikki C.

2014-01-01

This presentation gives a brief overview of Observing System Simulation Experiments (OSSEs), including what OSSEs are, and how and why they are performed. The intent is to educate the audience in light of the OSSE-related sections of the Forecast Improvement Act (H.R. 2413).

SARDA HITL Simulations: System Performance Results

NASA Technical Reports Server (NTRS)

Gupta, Gautam

2012-01-01

This presentation gives an overview of the 2012 SARDA human-in-the-loop simulation, and presents a summary of system performance results from the simulation, including delay, throughput and fuel consumption

Whistler Observations on DEMETER Compared with Full Electromagnetic Wave Simulations

NASA Astrophysics Data System (ADS)

Compston, A. J.; Cohen, M.; Lehtinen, N. G.; Inan, U.; Linscott, I.; Said, R.; Parrot, M.

2014-12-01

Terrestrial Very Low Frequency (VLF) electromagnetic radiation, which strongly impacts the Van Allen radiation belt electron dynamics, is injected across the ionosphere into the Earth's plasmasphere from two primary sources: man-made VLF transmitters and lightning discharges. Numerical models of trans-ionospheric propagation of such waves remain unvalidated, and early models may have overestimated the absorption, hindering a comprehensive understanding of the global impact of VLF waves in the loss of radiation belt electrons. In an attempt to remedy the problem of a lack of accurate trans-ionospheric propagation models, we have used a full electromagnetic wave method (FWM) numerical code to simulate the propagation of lightning-generated whistlers into the magnetosphere and compared the results with whistlers observed on the DEMETER satellite and paired with lightning stroke data from the National Lightning Detection Network (NLDN). We have identified over 20,000 whistlers occuring in 14 different passes of DEMETER over the central United States during the summer of 2009, and 14,000 of those occured within the 2000 km x 2000 km simulation grid we used. As shown in the attached figure, which shows a histogram of the ratio of the simulated whistler energy to the measured whistler energy for the 14,000 whistlers we compared, the simulation tends to slightly underestimate the total whistler energy injected by about 5 dB. However, the simulation underestimates the DEMETER measurements more as one gets further from the source lightning stroke, so since the signal to noise ratio of more distant whistlers will be smaller, possibly additive noise in the DEMETER measurements (which of course is not accounted for in the model) may explain some of the observed discrepancy.

Gaps, rings, and non-axisymmetric structures in protoplanetary disks. From simulations to ALMA observations

NASA Astrophysics Data System (ADS)

Flock, M.; Ruge, J. P.; Dzyurkevich, N.; Henning, Th.; Klahr, H.; Wolf, S.

2015-02-01

Aims: Recent observations by the Atacama Large Millimeter/submillimeter Array (ALMA) of disks around young stars revealed distinct asymmetries in the dust continuum emission. In this work we wish to study axisymmetric and non-axisymmetric structures that are generated by the magneto-rotational instability in the outer regions of protoplanetary disks. We combine the results of state-of-the-art numerical simulations with post-processing radiative transfer (RT) to generate synthetic maps and predictions for ALMA. Methods: We performed non-ideal global 3D magneto-hydrodynamic (MHD) stratified simulations of the dead-zone outer edge using the FARGO MHD code PLUTO. The stellar and disk parameters were taken from a parameterized disk model applied for fitting high-angular resolution multi-wavelength observations of various circumstellar disks. We considered a stellar mass of M∗ = 0.5 M⊙ and a total disk mass of about 0.085 M∗. The 2D initial temperature and density profiles were calculated consistently from a given surface density profile and Monte Carlo radiative transfer. The 2D Ohmic resistivity profile was calculated using a dust chemistry model. We considered two values for the dust-to-gas mass ratio, 10-2 and 10-4, which resulted in two different levels of magnetic coupling. The initial magnetic field was a vertical net flux field. The radiative transfer simulations were performed with the Monte Carlo-based 3D continuum RT code MC3D. The resulting dust reemission provided the basis for the simulation of observations with ALMA. Results: All models quickly turned into a turbulent state. The fiducial model with a dust-to-gas mass ratio of 10-2 developed a large gap followed by a jump in surface density located at the dead-zone outer edge. The jump in density and pressure was strong enough to stop the radial drift of particles at this location. In addition, we observed the generation of vortices by the Rossby wave instability at the jump location close to 60 AU

Asynchronous threat awareness by observer trials using crowd simulation

NASA Astrophysics Data System (ADS)

Dunau, Patrick; Huber, Samuel; Stein, Karin U.; Wellig, Peter

2016-10-01

The last few years showed that a high risk of asynchronous threats is given in every day life. Especially in large crowds a high probability of asynchronous attacks is evident. High observational abilities to detect threats are desirable. Consequently highly trained security and observation personal is needed. This paper evaluates the effectiveness of a training methodology to enhance performance of observation personnel engaging in a specific target identification task. For this purpose a crowd simulation video is utilized. The study first provides a measurement of the base performance before the training sessions. Furthermore a training procedure will be performed. Base performance will then be compared to the after training performance in order to look for a training effect. A thorough evaluation of both the training sessions as well as the overall performance will be done in this paper. A specific hypotheses based metric is used. Results will be discussed in order to provide guidelines for the design of training for observational tasks.

Quantitative Comparison of the Variability in Observed and Simulated Shortwave Reflectance

NASA Technical Reports Server (NTRS)

Roberts, Yolanda, L.; Pilewskie, P.; Kindel, B. C.; Feldman, D. R.; Collins, W. D.

2013-01-01

The Climate Absolute Radiance and Refractivity Observatory (CLARREO) is a climate observation system that has been designed to monitor the Earth's climate with unprecedented absolute radiometric accuracy and SI traceability. Climate Observation System Simulation Experiments (OSSEs) have been generated to simulate CLARREO hyperspectral shortwave imager measurements to help define the measurement characteristics needed for CLARREO to achieve its objectives. To evaluate how well the OSSE-simulated reflectance spectra reproduce the Earth s climate variability at the beginning of the 21st century, we compared the variability of the OSSE reflectance spectra to that of the reflectance spectra measured by the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY). Principal component analysis (PCA) is a multivariate decomposition technique used to represent and study the variability of hyperspectral radiation measurements. Using PCA, between 99.7%and 99.9%of the total variance the OSSE and SCIAMACHY data sets can be explained by subspaces defined by six principal components (PCs). To quantify how much information is shared between the simulated and observed data sets, we spectrally decomposed the intersection of the two data set subspaces. The results from four cases in 2004 showed that the two data sets share eight (January and October) and seven (April and July) dimensions, which correspond to about 99.9% of the total SCIAMACHY variance for each month. The spectral nature of these shared spaces, understood by examining the transformed eigenvectors calculated from the subspace intersections, exhibit similar physical characteristics to the original PCs calculated from each data set, such as water vapor absorption, vegetation reflectance, and cloud reflectance.

Atmospherical simulations of the OMEGA/MEX observations

NASA Astrophysics Data System (ADS)

Melchiorri, R.; Drossart, P.; Combes, M.; Encrenaz, T.; Fouchet, T.; Forget, F.; Bibring, J. P.; Ignatiev, N.; Moroz, V.; OMEGA Team

The modelization of the atmospheric contribution in the martian spectrum is an important step for the OMEGA data analysis.A full line by line radiative transfer calculation is made for the gas absorption; the dust opacity component, in a first approximation, is calculated as an optically thin additive component.Due to the large number of parameters needed in the calculations, the building of a huge data base to be interpolated is not envisageable, for each observed OMEGA spectrum with calculation for all the involved parameters (atmospheric pressure, water abundance, CO abundance, dust opacity and geometric angles of observation). The simulation of the observations allows us to fix all the orbital parameters and leave the unknown parameters as the only variables.Starting from the predictions of the current meteorological models of Mars we build a smaller data base corresponding on each observation. We present here a first order simulation, which consists in retrieving atmospheric contribution from the solar reflected component as a multiplicative (for gas absorption) and an additive component (for suspended dust contribution); although a fully consistent approach will require to include surface and atmosphere contributions together in synthetic calculations, this approach is sufficient for retrieving mineralogic information cleaned from atmospheric absorption at first order.First comparison to OMEGA spectra will be presented, with first order retrieval of CO2 pressure, CO and H2O abundance, and dust opacity.

An Introduction to Observing System Simulation Experiments

NASA Technical Reports Server (NTRS)

Prive, Nikki C.

2017-01-01

Observing System Simulation Experiments (OSSEs) are used to estimate the potential impact of proposed new instruments and data on numerical weather prediction. OSSEs can also be used to help design new observing platforms and to investigate the behavior of data assimilation systems. A basic overview of how to design and perform an OSSE will be given, as well as best practices and pitfalls. Some examples using the OSSE framework developed at the NASA Global Modeling and Assimilation Office will be shown.

The Aurora radiation-hydrodynamical simulations of reionization: calibration and first results

NASA Astrophysics Data System (ADS)

Pawlik, Andreas H.; Rahmati, Alireza; Schaye, Joop; Jeon, Myoungwon; Dalla Vecchia, Claudio

2017-04-01

We introduce a new suite of radiation-hydrodynamical simulations of galaxy formation and reionization called Aurora. The Aurora simulations make use of a spatially adaptive radiative transfer technique that lets us accurately capture the small-scale structure in the gas at the resolution of the hydrodynamics, in cosmological volumes. In addition to ionizing radiation, Aurora includes galactic winds driven by star formation and the enrichment of the universe with metals synthesized in the stars. Our reference simulation uses 2 × 5123 dark matter and gas particles in a box of size 25 h-1 comoving Mpc with a force softening scale of at most 0.28 h-1 kpc. It is accompanied by simulations in larger and smaller boxes and at higher and lower resolution, employing up to 2 × 10243 particles, to investigate numerical convergence. All simulations are calibrated to yield simulated star formation rate functions in close agreement with observational constraints at redshift z = 7 and to achieve reionization at z ≈ 8.3, which is consistent with the observed optical depth to reionization. We focus on the design and calibration of the simulations and present some first results. The median stellar metallicities of low-mass galaxies at z = 6 are consistent with the metallicities of dwarf galaxies in the Local Group, which are believed to have formed most of their stars at high redshifts. After reionization, the mean photoionization rate decreases systematically with increasing resolution. This coincides with a systematic increase in the abundance of neutral hydrogen absorbers in the intergalactic medium.

Comparisons of Cloud Properties over the Southern Ocean between In situ Observations and WRF Simulations

NASA Astrophysics Data System (ADS)

D'Alessandro, J.; Diao, M.; Wu, C.; Liu, X.

2017-12-01

Numerical weather models often struggle at representing clouds since small scale cloud processes must be parameterized. For example, models often utilize simple parameterizations for transitioning from liquid to ice, usually set as a function of temperature. However, supercooled liquid water (SLW) often persists at temperatures much lower than threshold values used in microphysics parameterizations. Previous observational studies of clouds over the Southern Ocean have found high frequencies of SLW (e.g., Morrison et al., 2011). Many of these studies have relied on satellite retrievals, which provide relatively low resolution observations and are often associated with large uncertainties due to assumptions of microphysical properties (e.g., particle size distributions). Recently, the NSF/NCAR O2/N2 Ratio and CO2 Airborne Southern Ocean Study (ORCAS) campaign took observations via the NSF/NCAR HIAPER research aircraft during January and February of 2016, providing in situ observations over the Southern Ocean (50°W to 92°W). We compare simulated results from the Weather Research and Forecasting (WRF) model with in situ observations from ORCAS. Differences between observations and simulations are evaluated via statistical analyses. Initial results from ORCAS reveal a high frequency of SLW at temperatures as low as -15°C, and the existence of SLW around -30°C. Recent studies have found that boundary layer clouds are underestimated by WRF in regions unaffected by cyclonic activity (Huang et al., 2014), suggesting a lack of low-level moisture due to local processes. To explore this, relative humidity distributions are examined and controlled by cloud microphysical characteristics (e.g., total water content) and relevant ambient properties (e.g., vertical velocity). A relatively low frequency of simulated SLW may in part explain the discrepancies in WRF, as cloud-top SLW results in stronger radiative cooling and turbulent motions conducive for long-lived cloud regimes

Photometric and polarimetric observations and model simulations of (216) Kleopatra

NASA Astrophysics Data System (ADS)

Takahashi, S.; Shinokawa, K.; Yoshida, F.; Mukai, T.; Ip, W. H.; Kawabata, K.

2004-10-01

We performed photometric and polarimetric observations, on November 8 and 9, 1999, of an M-type main belt asteroid, (216) Kleopatra by using the HBS spectropolarimeter installed at Dodaira observatory, National Astronomical Observatory of Japan (NAOJ). Photometric amplitude of lightcurve in the V band was 0.12 mag, and the averaged degree of polarization was -1.01±0.1%. It seems that the polarimetric data might also show a slight change in the degree of polarization ( ~0.2%) at the second minimum of the photometric lightcurve, but we could not confirm that the feature was real because of the large errors of data. With the assumption that the surface is uniform, we have carried out lightcurve simulations based on shape models by Ostro et al. (2000), Tanga et al. (2001) and Roche binary (Cellino et al., 1985). The results of simulations were compared to the configurations of lightcurves which had been obtained at different 4 geometric positions (1980, 1982, 1987 and 1999). The model by Cellino et al. (1985) reproduced almost all the data points without the 1987 observations within ~0.05 mag., which is the best result among the 3 models. The model by Tanga et al. (2001) well reproduced the lightcurves, but failed in reproducing the 1982 amplitude (difference Δdiff ~ 0 2 mag.). We also confirmed that the model by Ostro et al. (2000) could not explain the observed lightcurves.

Observed and Simulated Supercell Demise Depicted by VORTEX2 Observations

NASA Astrophysics Data System (ADS)

Letkewicz, Casey Elizabeth

technique deemed base-state substitution (BSS) was developed to incorporate a varying base-state while maintaining a degree of experimental control. BSS was designed to test the impact of a new horizontally homogeneous base-state environment on storm morphology while maintaining storm-induced perturbations. The thermodynamic and kinematic environment can be independently altered, with modifications incorporated slowly (gradual BSS) or all at once (instant BSS). Trial simulations demonstrated that model stability was maintained and domain-integrated fields were well-behaved. The simulations of the 9 June 2009 supercell were designed to isolate the effects of the changing wind profile from those of the increasingly stable boundary layer. The observed environmental modifications, when gradually incorporated into the model via BSS, were able to reasonably reproduce the observed storm morphology. Overall, the experiments demonstrated that the demise of the supercell was ultimately the result of depleting the updraft's buoyancy. The increasingly stable low-level environment predominately contributed to demise; however, changes to the wind profile also acted to notably impact storm strength and morphology. Analysis of updraft parcels in each experiment demonstrated that the mean source region became elevated over time as the low-level environment stabilized and/or shear and helicity weakened. Consequently, updraft buoyancy was eventually depleted since the drier elevated parcels (with less instability) were more negatively impacted by entrainment. The diagnostic pressure equation was utilized to examine the evolution of vertical accelerations within the simulated storm updraft and explore how storm dynamics evolved as the environment evolved. It was found that changes to nonlinear dynamic acceleration primarily influenced the evolution in total vertical acceleration, though these changes were interconnected with shifts in buoyancy and updraft strength and rotation. The findings of

Effects of experimental protocol on global vegetation model accuracy: a comparison of simulated and observed vegetation patterns for Asia

USGS Publications Warehouse

Tang, Guoping; Shafer, Sarah L.; Barlein, Patrick J.; Holman, Justin O.

2009-01-01

Prognostic vegetation models have been widely used to study the interactions between environmental change and biological systems. This study examines the sensitivity of vegetation model simulations to: (i) the selection of input climatologies representing different time periods and their associated atmospheric CO2 concentrations, (ii) the choice of observed vegetation data for evaluating the model results, and (iii) the methods used to compare simulated and observed vegetation. We use vegetation simulated for Asia by the equilibrium vegetation model BIOME4 as a typical example of vegetation model output. BIOME4 was run using 19 different climatologies and their associated atmospheric CO2 concentrations. The Kappa statistic, Fuzzy Kappa statistic and a newly developed map-comparison method, the Nomad index, were used to quantify the agreement between the biomes simulated under each scenario and the observed vegetation from three different global land- and tree-cover data sets: the global Potential Natural Vegetation data set (PNV), the Global Land Cover Characteristics data set (GLCC), and the Global Land Cover Facility data set (GLCF). The results indicate that the 30-year mean climatology (and its associated atmospheric CO2 concentration) for the time period immediately preceding the collection date of the observed vegetation data produce the most accurate vegetation simulations when compared with all three observed vegetation data sets. The study also indicates that the BIOME4-simulated vegetation for Asia more closely matches the PNV data than the other two observed vegetation data sets. Given the same observed data, the accuracy assessments of the BIOME4 simulations made using the Kappa, Fuzzy Kappa and Nomad index map-comparison methods agree well when the compared vegetation types consist of a large number of spatially continuous grid cells. The results of this analysis can assist model users in designing experimental protocols for simulating vegetation.

Simulations as a tool for higher mass resolution spectrometer: Lessons from existing observations

NASA Astrophysics Data System (ADS)

Nicolaou, Georgios; Yamauchi, Masatoshi; Nilsson, Hans; Wieser, Martin; Fedorov, Andrei

2017-04-01

Scientific requirements of each mission are crucial for the instrument's design. Ion tracing simulations of instruments can be helpful to characterize their performance, identify their limitations and improving the design for future missions. However, simulations provide the best performance in ideal case, and the actual response is determined by many other factors. Therefore, simulations should be compared with observations when possible. Characterizing the actual response of a running instrument gives valuable lessons for the future design of test instruments with the same detection principle before spending resources to build and calibrate them. In this study we use an ion tracing simulation of the Ion Composition Analyser (ICA) on board ROSETTA, in order to characterize its response and to compare it with the observations. It turned out that, due to the complicated unexpected response of the running instrument, the heavy cometary ions and molecules are sometimes difficult to be resolved. However, preliminary simulation of a slightly modified design predicts much higher mass resolution. Even after considering the complicated unexpected response, we safely expect that the modified design can resolve most abundant heavy atomic ions (e.g., O^+) and molecular ions (e.g., N_2+ and O_2^+). We show the simulation results for both designs and ICA data.

Comparisons Between TIME-GCM/MERRA Simulations and LEO Satellite Observations

NASA Astrophysics Data System (ADS)

Hagan, M. E.; Haeusler, K.; Forbes, J. M.; Zhang, X.; Doornbos, E.; Bruinsma, S.; Lu, G.

2014-12-01

We report on yearlong National Center for Atmospheric Research (NCAR) thermosphere-ionosphere-mesosphere-electrodynamics general circulation model (TIME-GCM) simulations where we utilize the recently developed lower boundary condition based on 3-hourly MERRA (Modern-Era Retrospective Analysis for Research and Application) reanalysis data to account for tropospheric waves and tides propagating upward into the model domain. The solar and geomagnetic forcing is based on prevailing geophysical conditions. The simulations show a strong day-to-day variability in the upper thermospheric neutral temperature tidal fields, which is smoothed out quickly when averaging is applied over several days, e.g. up to 50% DE3 amplitude reduction for a 10-day average. This is an important result with respect to tidal diagnostics from satellite observations where averaging over multiple days is inevitable. In order to assess TIME-GCM performance we compare the simulations with measurements from the Gravity field and steady-state Ocean Circulation Explorer (GOCE), Challenging Minisatellite Payload (CHAMP) and Gravity Recovery and Climate Experiment (GRACE) satellites.

Ice shelf basal melt rates around Antarctica from simulations and observations

NASA Astrophysics Data System (ADS)

Schodlok, M. P.; Menemenlis, D.; Rignot, E. J.

2016-02-01

We introduce an explicit representation of Antarctic ice shelf cavities in the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) ocean retrospective analysis; and compare resulting basal melt rates and patterns to independent estimates from satellite observations. Two simulations are carried out: the first is based on the original ECCO2 vertical discretization; the second has higher vertical resolution particularly at the depth range of ice shelf cavities. The original ECCO2 vertical discretization produces higher than observed melt rates and leads to a misrepresentation of Southern Ocean water mass properties and transports. In general, thicker levels at the base of the ice shelves lead to increased melting because of their larger heat capacity. This strengthens horizontal gradients and circulation within and outside the cavities and, in turn, warm water transports from the shelf break to the ice shelves. The simulation with more vertical levels produces basal melt rates (1735 ± 164 Gt/a) and patterns that are in better agreement with observations. Thinner levels in the sub-ice-shelf cavities improve the representation of a fresh/cold layer at the ice shelf base and of warm/salty water near the bottom, leading to a sharper pycnocline and reduced vertical mixing underneath the ice shelf. Improved water column properties lead to more accurate melt rates and patterns, especially for melt/freeze patterns under large cold-water ice shelves. At the 18 km grid spacing of the ECCO2 model configuration, the smaller, warm-water ice shelves cannot be properly represented, with higher than observed melt rates in both simulations.

Comparing cosmological hydrodynamic simulations with observations of high- redshift galaxy formation

NASA Astrophysics Data System (ADS)

Finlator, Kristian Markwart

We use cosmological hydrodynamic simulations to study the impact of outflows and radiative feedback on high-redshift galaxies. For outflows, we consider simulations that assume (i) no winds, (ii) a "constant-wind" model in which the mass-loading factor and outflow speed are constant, and (iii) "momentum-driven" winds in which both parameters vary smoothly with mass. In order to treat radiative feedback, we develop a moment-based radiative transfer technique that operates in both post-processing and coupled radiative hydrodynamic modes. We first ask how outflows impact the broadband spectral energy distributions (SEDs) of six observed reionization-epoch galaxies. Simulations reproduce five regardless of the outflow prescription, while the sixth suggests an unusually bursty star formation history. We conclude that (i) simulations broadly account for available constraints on reionization-epoch galaxies, (ii) individual SEDs do not constrain outflows, and (iii) SED comparisons efficiently isolate objects that challenge simulations. We next study how outflows impact the galaxy mass metallicity relation (MZR). Momentum-driven outflows uniquely reproduce observations at z = 2. In this scenario, galaxies obey two equilibria: (i) The rate at which a galaxy processes gas into stars and outflows tracks its inflow rate; and (ii) The gas enrichment rate owing to star formation balances the dilution rate owing to inflows. Combining these conditions indicates that the MZR is dominated by the (instantaneous) variation of outflows with mass, with more-massive galaxies driving less gas into outflows per unit stellar mass formed. Turning to radiative feedback, we use post-processing simulations to study the topology of reionization. Reionization begins in overdensities and then "leaks" directly into voids, with filaments reionizing last owing to their high density and low emissivity. This result conflicts with previous findings that voids ionize last. We argue that it owes to the

Trend Estimates of AERONET-Observed and Model-Simulated AOTs Between 1993 and 2013

NASA Technical Reports Server (NTRS)

Yoon, J.; Pozzer, A.; Chang, D. Y.; Lelieveld, J.; Kim, J.; Kim, M.; Lee, Y. G.; Koo, J.-H.; Lee, J.; Moon, K. J.

2015-01-01

Recently, temporal changes in Aerosol Optical Thickness (AOT) have been investigated based on model simulations, satellite and ground-based observations. Most AOT trend studies used monthly or annual arithmetic means that discard details of the generally right-skewed AOT distributions. Potentially, such results can be biased by extreme values (including outliers). This study additionally uses percentiles (i.e., the lowest 5%, 25%, 50%, 75% and 95% of the monthly cumulative distributions fitted to Aerosol Robotic Network (AERONET)-observed and ECHAM/MESSy Atmospheric Chemistry (EMAC)-model simulated AOTs) that are less affected by outliers caused by measurement error, cloud contamination and occasional extreme aerosol events. Since the limited statistical representativeness of monthly percentiles and means can lead to bias, this study adopts the number of observations as a weighting factor, which improves the statistical robustness of trend estimates. By analyzing the aerosol composition of AERONET-observed and EMAC-simulated AOTs in selected regions of interest, we distinguish the dominant aerosol types and investigate the causes of regional AOT trends. The simulated and observed trends are generally consistent with a high correlation coefficient (R = 0.89) and small bias (slope+/-2(sigma) = 0.75 +/- 0.19). A significant decrease in EMAC-decomposed AOTs by water-soluble compounds and black carbon is found over the USA and the EU due to environmental regulation. In particular, a clear reversal in the AERONET AOT trend percentiles is found over the USA, probably related to the AOT diurnal cycle and the frequency of wildfires. In most of the selected regions of interest, EMAC-simulated trends are mainly attributed to the significant changes of the dominant aerosols; e.g., significant decrease in sea salt and water soluble compounds over Central America, increase in dust over Northern Africa and Middle East, and decrease in black carbon and organic carbon over

Sedimentary basin effects in Seattle, Washington: Ground-motion observations and 3D simulations

USGS Publications Warehouse

Frankel, Arthur; Stephenson, William; Carver, David

2009-01-01

Seismograms of local earthquakes recorded in Seattle exhibit surface waves in the Seattle basin and basin-edge focusing of S waves. Spectral ratios of Swaves and later arrivals at 1 Hz for stiff-soil sites in the Seattle basin show a dependence on the direction to the earthquake, with earthquakes to the south and southwest producing higher average amplification. Earthquakes to the southwest typically produce larger basin surface waves relative to S waves than earthquakes to the north and northwest, probably because of the velocity contrast across the Seattle fault along the southern margin of the Seattle basin. S to P conversions are observed for some events and are likely converted at the bottom of the Seattle basin. We model five earthquakes, including the M 6.8 Nisqually earthquake, using 3D finite-difference simulations accurate up to 1 Hz. The simulations reproduce the observed dependence of amplification on the direction to the earthquake. The simulations generally match the timing and character of basin surface waves observed for many events. The 3D simulation for the Nisqually earth-quake produces focusing of S waves along the southern margin of the Seattle basin near the area in west Seattle that experienced increased chimney damage from the earthquake, similar to the results of the higher-frequency 2D simulation reported by Stephenson et al. (2006). Waveforms from the 3D simulations show reasonable agreement with the data at low frequencies (0.2-0.4 Hz) for the Nisqually earthquake and an M 4.8 deep earthquake west of Seattle.

Results from a limited area mesoscale numerical simulation for 10 April 1979

NASA Technical Reports Server (NTRS)

Kalb, M. W.

1985-01-01

Results are presented from a nine-hour limited area fine mesh (35-km) mesoscale model simulation initialized with SESAME-AVE I radiosonde data for Apr. 10, 1979 at 2100 GMT. Emphasis is on the diagnosis of mesoscale structure in the mass and precipitation fields. Along the Texas/Oklahoma border, independent of the short wave, convective precipitation formed several hours into the simulation and was organized into a narrow band suggestive of the observed April 10 squall line.

Retrieved Products from Simulated Hyperspectral Observations of a Hurricane

NASA Technical Reports Server (NTRS)

Susskind, Joel; Kouvaris, Louis; Iredell, Lena; Blaisdell, John

2015-01-01

Demonstrate via Observing System Simulation Experiments (OSSEs) the potential utility of flying high spatial resolution AIRS class IR sounders on future LEO and GEO missions.The study simulates and analyzes radiances for 3 sounders with AIRS spectral and radiometric properties on different orbits with different spatial resolutions: 1) Control run 13 kilometers AIRS spatial resolution at nadir on LEO in Aqua orbit; 2) 2 kilometer spatial resolution LEO sounder at nadir ARIES; 3) 5 kilometers spatial resolution sounder on a GEO orbit, radiances simulated every 72 minutes.

Thermodynamic sensitivities in observed and simulated extreme-rain-producing mesoscale convective systems

NASA Astrophysics Data System (ADS)

Schumacher, R. S.; Peters, J. M.

2015-12-01

Mesoscale convective systems (MCSs) are responsible for a large fraction of warm-season extreme rainfall events over the continental United States, as well as other midlatitude regions globally. The rainfall production in these MCSs is determined by numerous factors, including the large-scale forcing for ascent, the organization of the convection, cloud microphysical processes, and the surrounding thermodynamic and kinematic environment. Furthermore, heavy-rain-producing MCSs are most common at night, which means that well-studied mechanisms for MCS maintenance and organization such as cold pools (gravity currents) are not always at work. In this study, we use numerical model simulations and recent field observations to investigate the sensitivity of low-level MCS structures, and their influences on rainfall, to the details of the thermodynamic environment. In particular, small alterations to the initial conditions in idealized and semi-idealized simulations result in comparatively large precipitation changes, both in terms of the intensity and the spatial distribution. The uncertainties in the thermodynamic enviroments in the model simulations will be compared with high-resolution observations from the Plains Elevated Convection At Night (PECAN) field experiment in 2015. The results have implications for the paradigms of "surface-based" versus "elevated" convection, as well as for the predictability of warm-season convective rainfall.

Improve earthquake hypocenter using adaptive simulated annealing inversion in regional tectonic, volcano tectonic, and geothermal observation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ry, Rexha Verdhora, E-mail: rexha.vry@gmail.com; Nugraha, Andri Dian, E-mail: nugraha@gf.itb.ac.id

Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment.more » We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger’s method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger’s result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.« less

Observations and simulations of the western United States' hydroclimate

NASA Astrophysics Data System (ADS)

Guirguis, Kristen

While very important from an economical and societal point of view, estimating precipitation in the western United States remains an unsolved and challenging problem. This is due to difficulties in observing and modeling precipitation in complex terrain. This research examines this issue by (i) providing a systematic evaluation of precipitation observations to quantify data uncertainty; and (ii) investigating the ability of the Ocean-Land-Atmosphere Model (OLAM) to simulate the winter hydroclimate in this region. This state-of-the-art, non-hydrostatic model has the capability of simulating simultaneously all scales of motions at various resolutions. This research intercompares nine precipitation datasets commonly used in hydrometeorological research in two ways. First, using principal component analysis, a precipitation climatology is conducted for the western U.S. from which five unique precipitation climates are identified. From this analysis, data uncertainty is shown to be primarily due to differences in (i) precipitation over the Rocky Mountains, (ii) the eastward wet-to-dry precipitation gradient during the cold season, (iii) the North American Monsoon signal, and (iv) precipitation in the desert southwest during spring and summer. The second intercomparison uses these five precipitation regions to provide location-specific assessments of uncertainty, which is shown to be dependent on season, location. Long-range weather forecasts on the order of a season are important for water-scarce regions such as the western U.S. The modeling component of this research looks at the ability of the OLAM to simulate the hydroclimate in the western U.S. during the winter of 1999. Six global simulations are run, each with a different spatial resolution over the western U.S. (360 km down to 11 km). For this study, OLAM is configured as for a long-range seasonal hindcast but with observed sea surface temperatures. OLAM precipitation compares well against observations, and is

DKIST Adaptive Optics System: Simulation Results

NASA Astrophysics Data System (ADS)

Marino, Jose; Schmidt, Dirk

2016-05-01

The 4 m class Daniel K. Inouye Solar Telescope (DKIST), currently under construction, will be equipped with an ultra high order solar adaptive optics (AO) system. The requirements and capabilities of such a solar AO system are beyond those of any other solar AO system currently in operation. We must rely on solar AO simulations to estimate and quantify its performance.We present performance estimation results of the DKIST AO system obtained with a new solar AO simulation tool. This simulation tool is a flexible and fast end-to-end solar AO simulator which produces accurate solar AO simulations while taking advantage of current multi-core computer technology. It relies on full imaging simulations of the extended field Shack-Hartmann wavefront sensor (WFS), which directly includes important secondary effects such as field dependent distortions and varying contrast of the WFS sub-aperture images.

Comparison of the analytical and simulation results of the equilibrium beam profile

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Z. J.; Zhu Shaoping; Cao, L. H.

2007-10-15

The evolution of high current electron beams in dense plasmas has been investigated by using two-dimensional particle-in-cell (PIC) simulations with immobile ions. It is shown that electron beams are split into many filaments at the beginning due to the Weibel instability, and then different filamentation beams attract each other and coalesce. The profile of the filaments can be described by formulas. Hammer et al. [Phys. Fluids 13, 1831 (1970)] developed a self-consistent relativistic electron beam model that allows the propagation of relativistic electron fluxes in excess of the Alfven-Lawson critical-current limit for a fully neutralized beam. The equilibrium solution hasmore » been observed in the simulation results, but the electron distribution function assumed by Hammer et al. is different from the simulation results.« less

Modelling dust polarization observations of molecular clouds through MHD simulations

NASA Astrophysics Data System (ADS)

King, Patrick K.; Fissel, Laura M.; Chen, Che-Yu; Li, Zhi-Yun

2018-03-01

The BLASTPol observations of Vela C have provided the most detailed characterization of the polarization fraction p and dispersion in polarization angles S for a molecular cloud. We compare the observed distributions of p and S with those obtained in synthetic observations of simulations of molecular clouds, assuming homogeneous grain alignment. We find that the orientation of the mean magnetic field relative to the observer has a significant effect on the p and S distributions. These distributions for Vela C are most consistent with synthetic observations where the mean magnetic field is close to the line of sight. Our results point to apparent magnetic disorder in the Vela C molecular cloud, although it can be due to either an inclination effect (i.e. observing close to the mean field direction) or significant field tangling from strong turbulence/low magnetization. The joint correlations of p with column density and of S with column density for the synthetic observations generally agree poorly with the Vela C joint correlations, suggesting that understanding these correlations requires a more sophisticated treatment of grain alignment physics.

Observational Implications of Gamma-Ray Burst Afterglow Jet Simulations and Numerical Light Curve Calculations

NASA Astrophysics Data System (ADS)

van Eerten, Hendrik J.; MacFadyen, Andrew I.

2012-06-01

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 rad. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor of four. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late-time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late-time features such as the rise of the counterjet, detailed jet simulations remain indispensable.

Supporting ITM Missions by Observing System Simulation Experiments: Initial Design, Challenges and Perspectives

NASA Astrophysics Data System (ADS)

Yudin, V. A.; England, S.; Matsuo, T.; Wang, H.; Immel, T. J.; Eastes, R.; Akmaev, R. A.; Goncharenko, L. P.; Fuller-Rowell, T. J.; Liu, H.; Solomon, S. C.; Wu, Q.

2014-12-01

We review and discuss the capability of novel configurations of global community (WACCM-X and TIME-GCM) and planned-operational (WAM) models to support current and forthcoming space-borne missions to monitor the dynamics and composition of the Ionosphere-Thermosphere-Mesosphere (ITM) system. In the specified meteorology model configuration of WACCM-X, the lower atmosphere is constrained by operational analyses and/or short-term forecasts provided by the Goddard Earth Observing System (GEOS-5) of GMAO/NASA/GSFC. With the terrestrial weather of GEOS-5 and updated model physics, WACCM-X simulations are capable to reproduce the observed signatures of the perturbed wave dynamics and ion-neutral coupling during recent (2006-2013) stratospheric warming events, short-term, annual and year-to-year variability of prevailing flows, planetary waves, tides, and composition. With assimilation of the NWP data in the troposphere and stratosphere the planned-operational configuration of WAM can also recreate the observed features of the ITM day-to-day variability. These "terrestrial-weather" driven whole atmosphere simulations, with day-to-day variable solar and geomagnetic inputs, can provide specification of the background state (first guess) and errors for the inverse algorithms of forthcoming NASA ITM missions, such as ICON and GOLD. With two different viewing geometries (sun-synchronous, for ICON and geostationary for GOLD) these missions promise to perform complimentary global observations of temperature, winds and constituents to constrain the first-principle space weather forecast models. The paper will discuss initial designs of Observing System Simulation Experiments (OSSE) in the coupled simulations of TIME-GCM/WACCM-X/GEOS5 and WAM/GIP. As recognized, OSSE represent an excellent learning tool for designing and evaluating observing capabilities of novel sensors. The choice of assimilation schemes, forecast and observational errors will be discussed along with challenges

Multiwavelength mock observations of the WHIM in a simulated galaxy cluster

NASA Astrophysics Data System (ADS)

Planelles, Susana; Mimica, Petar; Quilis, Vicent; Cuesta-Martínez, Carlos

2018-06-01

About half of the expected total baryon budget in the local Universe is `missing'. Hydrodynamical simulations suggest that most of the missing baryons are located in a mildly overdense, warm-hot intergalactic medium (WHIM), which is difficult to be detected at most wavelengths. In this paper, we explore multiwavelength synthetic observations of a massive galaxy cluster developed in a full Eulerian-adaptive mesh refinement cosmological simulation. A novel numerical procedure is applied on the outputs of the simulation, which are post-processed with a full-radiative transfer code that can compute the change of the intensity at any frequency along the null geodesic of photons. We compare the emission from the whole intergalactic medium and from the WHIM component (defined as the gas with a temperature in the range 105-107 K) at three observational bands associated with thermal X-rays, thermal and kinematic Sunyaev-Zel'dovich effect, and radio emission. The synthetic maps produced by this procedure could be directly compared with existing observational maps and could be used as a guide for future observations with forthcoming instruments. The analysis of the different emissions associated with a high-resolution galaxy cluster is in broad agreement with previous simulated and observational estimates of both gas components.

Discrepancy between simulated and observed ethane and propane levels explained by underestimated fossil emissions

NASA Astrophysics Data System (ADS)

Dalsøren, Stig B.; Myhre, Gunnar; Hodnebrog, Øivind; Myhre, Cathrine Lund; Stohl, Andreas; Pisso, Ignacio; Schwietzke, Stefan; Höglund-Isaksson, Lena; Helmig, Detlev; Reimann, Stefan; Sauvage, Stéphane; Schmidbauer, Norbert; Read, Katie A.; Carpenter, Lucy J.; Lewis, Alastair C.; Punjabi, Shalini; Wallasch, Markus

2018-03-01

Ethane and propane are the most abundant non-methane hydrocarbons in the atmosphere. However, their emissions, atmospheric distribution, and trends in their atmospheric concentrations are insufficiently understood. Atmospheric model simulations using standard community emission inventories do not reproduce available measurements in the Northern Hemisphere. Here, we show that observations of pre-industrial and present-day ethane and propane can be reproduced in simulations with a detailed atmospheric chemistry transport model, provided that natural geologic emissions are taken into account and anthropogenic fossil fuel emissions are assumed to be two to three times higher than is indicated in current inventories. Accounting for these enhanced ethane and propane emissions results in simulated surface ozone concentrations that are 5-13% higher than previously assumed in some polluted regions in Asia. The improved correspondence with observed ethane and propane in model simulations with greater emissions suggests that the level of fossil (geologic + fossil fuel) methane emissions in current inventories may need re-evaluation.

Numerical simulation of granular flows : comparison with experimental results

NASA Astrophysics Data System (ADS)

Pirulli, M.; Mangeney-Castelnau, A.; Lajeunesse, E.; Vilotte, J.-P.; Bouchut, F.; Bristeau, M. O.; Perthame, B.

2003-04-01

Granular avalanches such as rock or debris flows regularly cause large amounts of human and material damages. Numerical simulation of granular avalanches should provide a useful tool for investigating, within realistic geological contexts, the dynamics of these flows and of their arrest phase and for improving the risk assessment of such natural hazards. Validation of debris avalanche numerical model on granular experiments over inclined plane is performed here. The comparison is performed by simulating granular flow of glass beads from a reservoir through a gate down an inclined plane. This unsteady situation evolves toward the steady state observed in the laboratory. Furthermore simulation exactly reproduces the arrest phase obtained by suddenly closing the gate of the reservoir once a thick flow has developped. The spreading of a granular mass released from rest at the top of a rough inclined plane is also investigated. The evolution of the avalanche shape, the velocity and the characteristics of the arrest phase are compared with experimental results and analysis of the involved forces are studied for various flow laws.

Understanding Southern Ocean SST Trends in Historical Simulations and Observations

NASA Astrophysics Data System (ADS)

Kostov, Yavor; Ferreira, David; Marshall, John; Armour, Kyle

2017-04-01

Historical simulations with CMIP5 global climate models do not reproduce the observed 1979-2014 Southern Ocean (SO) cooling, and most ensemble members predict gradual warming around Antarctica. In order to understand this discrepancy and the mechanisms behind the SO cooling, we analyze output from 19 CMIP5 models. For each ensemble member we estimate the characteristic responses of SO SST to step changes in greenhouse gas (GHG) forcing and in the seasonal indices of the Southern Annular Mode (SAM). Using these step-response functions and linear convolution theory, we reconstruct the original CMIP5 simulations of 1979-2014 SO SST trends. We recover the CMIP5 ensemble mean trend, capture the intermodel spread, and reproduce very well the behavior of individual models. We thus suggest that GHG forcing and the SAM are major drivers of the simulated 1979-2014 SO SST trends. In consistence with the seasonal signature of the Antarctic ozone hole, our results imply that the summer (DJF) and fall (MAM) SAM exert a particularly important effect on the SO SST. In some CMIP5 models the SO SST response to SAM partially counteracts the warming due to GHG forcing, while in other ensemble members the SAM-induced SO SST trends complement the warming effect of GHG forcing. The compensation between GHG and SAM-induced SO SST anomalies is model-dependent and is determined by multiple factors. Firstly, CMIP5 models have different characteristic SST step response functions to SAM. Kostov et al. (2016) relate these differences to biases in the models' climatological SO temperature gradients. Secondly, many CMIP5 historical simulations underestimate the observed positive trends in the DJF and MAM seasonal SAM indices. We show that this affects the models' ability to reproduce the observed SO cooling. Last but not least, CMIP5 models differ in their SO SST step response functions to GHG forcing. Understanding the diverse behavior of CMIP5 models helps shed light on the physical processes

Panchromatic spectral energy distributions of simulated galaxies: results at redshift z = 0

NASA Astrophysics Data System (ADS)

Goz, David; Monaco, Pierluigi; Granato, Gian Luigi; Murante, Giuseppe; Domínguez-Tenreiro, Rosa; Obreja, Aura; Annunziatella, Marianna; Tescari, Edoardo

2017-08-01

We present predictions of spectral energy distributions (SEDs), from the UV to the FIR, of simulated galaxies at z = 0. These were obtained by post-processing the results of an N-body+hydro simulation of a cosmological box of side 25 Mpc, which uses the Multi-Phase Particle Integrator (MUPPI) for star formation and stellar feedback, with the grasil-3d radiative transfer code that includes reprocessing of UV light by dust. Physical properties of our sample of ˜500 galaxies resemble observed ones, though with some tension at small and large stellar masses. Comparing predicted SEDs of simulated galaxies with different samples of local galaxies, we find that these resemble observed ones, when normalized at 3.6 μm. A comparison with the Herschel Reference Survey shows that the average SEDs of galaxies, divided in bins of star formation rate (SFR), are reproduced in shape and absolute normalization to within a factor of ˜2, while average SEDs of galaxies divided in bins of stellar mass show tensions that are an effect of the difference of simulated and observed galaxies in the stellar mass-SFR plane. We use our sample to investigate the correlation of IR luminosity in Spitzer and Herschel bands with several galaxy properties. SFR is the quantity that best correlates with IR light up to 160 μm, while at longer wavelengths better correlations are found with molecular mass and, at 500 μm, with dust mass. However, using the position of the FIR peak as a proxy for cold dust temperature, we assess that heating of cold dust is mostly determined by SFR, with stellar mass giving only a minor contribution. We finally show how our sample of simulated galaxies can be used as a guide to understand the physical properties and selection biases of observed samples.

Sensitivity of Stratocumulus Optical Depths to Droplet Concentrations: Satellite Observations and Large-Eddy Simulations

NASA Technical Reports Server (NTRS)

Ackerman, A. S.; Stevens, D. E.; Toon, O. B.; Coakley, J. A., Jr.; Gore, Warren J. (Technical Monitor)

2002-01-01

A number of observations and simulations have shown that increased droplet concentrations in ship tracks increase their total cross-sectional area, thereby enhancing cloud albedo and providing a negative (cooling) radiative forcing at the surface and the top of the atmosphere. In some cases cloud water has been found to be enhanced in ship tracks, which has been attributed to suppression of drizzle and implies an enhanced susceptibility of cloud albedo to droplet concentrations. However, observations from aircraft and satellite indicate that on average cloud water is instead reduced in daytime ship tracks. Such a reduction in liquid water may be attributable to cloud-burning caused by solar heating by soot within the ship exhaust, or by increased precipitation resulting from giant nuclei in the ship exhaust. We will summarize the observational evidence and present results from large-eddy simulations that evaluate these mechanisms. Along the way we will present our insights into the interpretation of satellite retrievals of cloud microphysical properties.

ON THE OBSERVATION AND SIMULATION OF SOLAR CORONAL TWIN JETS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Jiajia; Wang, Yuming; Zhang, Quanhao

We present the first observation, analysis, and modeling of solar coronal twin jets, which occurred after a preceding jet. Detailed analysis on the kinetics of the preceding jet reveals its blowout-jet nature, which resembles the one studied in Liu et al. However, the erupting process and kinetics of the twin jets appear to be different from the preceding one. Lacking detailed information on the magnetic fields in the twin jet region, we instead use a numerical simulation using a three-dimensional (3D) MHD model as described in Fang et al., and find that in the simulation a pair of twin jetsmore » form due to reconnection between the ambient open fields and a highly twisted sigmoidal magnetic flux, which is the outcome of the further evolution of the magnetic fields following the preceding blowout jet. Based on the similarity between the synthesized and observed emission, we propose this mechanism as a possible explanation for the observed twin jets. Combining our observation and simulation, we suggest that with continuous energy transport from the subsurface convection zone into the corona, solar coronal twin jets could be generated in the same fashion addressed above.« less

Teaching childbirth with high-fidelity simulation. Is it better observing the scenario during the briefing session?

PubMed

Cuerva, Marcos J; Piñel, Carlos S; Martin, Lourdes; Espinosa, Jose A; Corral, Octavio J; Mendoza, Nicolás

2018-02-12

The design of optimal courses for obstetric undergraduate teaching is a relevant question. This study evaluates two different designs of simulator-based learning activity on childbirth with regard to respect to the patient, obstetric manoeuvres, interpretation of cardiotocography tracings (CTG) and infection prevention. This randomised experimental study which differs in the content of their briefing sessions consisted of two groups of undergraduate students, who performed two simulator-based learning activities on childbirth. The first briefing session included the observations of a properly performed scenario according to Spanish clinical practice guidelines on care in normal childbirth by the teachers whereas the second group did not include the observations of a properly performed scenario, and the students observed it only after the simulation process. The group that observed a properly performed scenario after the simulation obtained worse grades during the simulation, but better grades during the debriefing and evaluation. Simulator use in childbirth may be more fruitful when the medical students observe correct performance at the completion of the scenario compared to that at the start of the scenario. Impact statement What is already known on this subject? There is a scarcity of literature about the design of optimal high-fidelity simulation training in childbirth. It is known that preparing simulator-based learning activities is a complex process. Simulator-based learning includes the following steps: briefing, simulation, debriefing and evaluation. The most important part of high-fidelity simulations is the debriefing. A good briefing and simulation are of high relevance in order to have a fruitful debriefing session. What do the results of this study add? Our study describes a full simulator-based learning activity on childbirth that can be reproduced in similar facilities. The findings of this study add that high-fidelity simulation training in

Red-light running violation prediction using observational and simulator data.

PubMed

Jahangiri, Arash; Rakha, Hesham; Dingus, Thomas A

2016-11-01

In the United States, 683 people were killed and an estimated 133,000 were injured in crashes due to running red lights in 2012. To help prevent/mitigate crashes caused by running red lights, these violations need to be identified before they occur, so both the road users (i.e., drivers, pedestrians, etc.) in potential danger and the infrastructure can be notified and actions can be taken accordingly. Two different data sets were used to assess the feasibility of developing red-light running (RLR) violation prediction models: (1) observational data and (2) driver simulator data. Both data sets included common factors, such as time to intersection (TTI), distance to intersection (DTI), and velocity at the onset of the yellow indication. However, the observational data set provided additional factors that the simulator data set did not, and vice versa. The observational data included vehicle information (e.g., speed, acceleration, etc.) for several different time frames. For each vehicle approaching an intersection in the observational data set, required data were extracted from several time frames as the vehicle drew closer to the intersection. However, since the observational data were inherently anonymous, driver factors such as age and gender were unavailable in the observational data set. Conversely, the simulator data set contained age and gender. In addition, the simulator data included a secondary (non-driving) task factor and a treatment factor (i.e., incoming/outgoing calls while driving). The simulator data only included vehicle information for certain time frames (e.g., yellow onset); the data did not provide vehicle information for several different time frames while vehicles were approaching an intersection. In this study, the random forest (RF) machine-learning technique was adopted to develop RLR violation prediction models. Factor importance was obtained for different models and different data sets to show how differently the factors influence the

Properties of Turbulence in the Reconnection Exhaust: Numerical Simulations Compared with Observations

NASA Astrophysics Data System (ADS)

Pucci, F.; Servidio, S.; Sorriso-Valvo, L.; Olshevsky, V.; Matthaeus, W. H.; Malara, F.; Goldman, M. V.; Newman, D. L.; Lapenta, G.

2017-05-01

The properties of the turbulence that develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail, investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the dissipation mechanisms in the turbulent outflows similar to the ones observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

Properties of turbulence in the reconnection exhaust: numerical simulations compared with observations

NASA Astrophysics Data System (ADS)

Pucci, Francesco; Servidio, Sergio; Sorriso-Valvo, Luca; Olshevsky, Vyacheslav; Matthaeus, William; Malara, Francesco; Goldman, Martin; Newman, David; Lapenta, Giovanni

2017-04-01

The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with observations of reconnection event in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show: the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

Properties of Turbulence in the Reconnection Exhaust: Numerical Simulations Compared with Observations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pucci, F.; Olshevsky, V.; Lapenta, G.

2017-05-20

The properties of the turbulence that develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail, investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the dissipation mechanisms in the turbulent outflows similar to the ones observed in reconnection jets inmore » the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.« less

Gas kinematics in FIRE simulated galaxies compared to spatially unresolved H I observations

NASA Astrophysics Data System (ADS)

El-Badry, Kareem; Bradford, Jeremy; Quataert, Eliot; Geha, Marla; Boylan-Kolchin, Michael; Weisz, Daniel R.; Wetzel, Andrew; Hopkins, Philip F.; Chan, T. K.; Fitts, Alex; Kereš, Dušan; Faucher-Giguère, Claude-André

2018-06-01

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotationally supported gas discs produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing V/σ and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of ˜2000 low-redshift, H I-detected galaxies with Mstar = 107-11 M⊙ and compare to the simulated galaxies. At Mstar ≳ 1010 M⊙, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotationally supported discs. Both the observed and simulated line profiles become more Gaussian like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotational to dispersion support more strongly: at Mstar = 108-10 M⊙, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion supported. Most of the lower-mass-simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotational versus dispersion support.

The Shock Dynamics of Heterogeneous YSO Jets: 3D Simulations Meet Multi-epoch Observations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hansen, E. C.; Frank, A.; Hartigan, P.

High-resolution observations of young stellar object (YSO) jets show them to be composed of many small-scale knots or clumps. In this paper, we report results of 3D numerical simulations designed to study how such clumps interact and create morphologies and kinematic patterns seen in emission line observations. Our simulations focus on clump scale dynamics by imposing velocity differences between spherical, over-dense regions, which then lead to the formation of bow shocks as faster clumps overtake slower material. We show that much of the spatial structure apparent in emission line images of jets arises from the dynamics and interactions of thesemore » bow shocks. Our simulations show a variety of time-dependent features, including bright knots associated with Mach stems where the shocks intersect, a “frothy” emission structure that arises from the presence of the Nonlinear Thin Shell Instability along the surfaces of the bow shocks, and the merging and fragmentation of clumps. Our simulations use a new non-equilibrium cooling method to produce synthetic emission maps in H α and [S ii]. These are directly compared to multi-epoch Hubble Space Telescope observations of Herbig–Haro jets. We find excellent agreement between features seen in the simulations and the observations in terms of both proper motion and morphologies. Thus we conclude that YSO jets may be dominated by heterogeneous structures and that interactions between these structures and the shocks they produce can account for many details of YSO jet evolution.« less

Simulations of VLBI observations of a geodetic satellite providing co-location in space

NASA Astrophysics Data System (ADS)

Anderson, James M.; Beyerle, Georg; Glaser, Susanne; Liu, Li; Männel, Benjamin; Nilsson, Tobias; Heinkelmann, Robert; Schuh, Harald

2018-02-01

We performed Monte Carlo simulations of very-long-baseline interferometry (VLBI) observations of Earth-orbiting satellites incorporating co-located space-geodetic instruments in order to study how well the VLBI frame and the spacecraft frame can be tied using such measurements. We simulated observations of spacecraft by VLBI observations, time-of-flight (TOF) measurements using a time-encoded signal in the spacecraft transmission, similar in concept to precise point positioning, and differential VLBI (D-VLBI) observations using angularly nearby quasar calibrators to compare their relative performance. We used the proposed European Geodetic Reference Antenna in Space (E-GRASP) mission as an initial test case for our software. We found that the standard VLBI technique is limited, in part, by the present lack of knowledge of the absolute offset of VLBI time to Coordinated Universal Time at the level of microseconds. TOF measurements are better able to overcome this problem and provide frame ties with uncertainties in translation and scale nearly a factor of three smaller than those yielded from VLBI measurements. If the absolute time offset issue can be resolved by external means, the VLBI results can be significantly improved and can come close to providing 1 mm accuracy in the frame tie parameters. D-VLBI observations with optimum performance assumptions provide roughly a factor of two higher uncertainties for the E-GRASP orbit. We additionally simulated how station and spacecraft position offsets affect the frame tie performance.

Preparation for the solar system observations with Herschel: Simulation of Jupiter observations with PACS

NASA Astrophysics Data System (ADS)

Sagawa, Hideo; Hartogh, Paul; Rengel, Miriam; de Lange, Arno; Cavalié, Thibault

2010-11-01

Observations of the water inventory as well as other chemically important species on Jupiter will be performed in the frame of the guaranteed time key project of the Herschel Space Observatory entitled "Water and related chemistry in the Solar system". Among other onboard instruments, PACS (Photodetector Array Camera and Spectrometer) will provide new data of the spectral atlas in a wide region covering the far-infrared and submillimetre domains, with an improved spectral resolution and a higher sensitivity compared to previous observations carried out by Cassini/CIRS (Composite InfraRed Spectrometer) and by ISO (Infrared Space Observatory). In order to optimise the observational plan and to prepare for the data analysis, we have simulated the expected spectra of PACS Jupiter observations. Our simulation shows that PACS will promisingly detect several H 2O emission lines. As PACS is capable of spatially resolving the Jovian disk, we will be able to discern the external oxygen sources in the giant planets by exploring the horizontal distribution of water. In addition to H 2O lines, some absorption lines due to tropospheric CH 4, HD, PH 3 and NH 3 lines will be observed with PACS. Furthermore, owing to the high sensitivity of the instrument, the current upper limit on the abundance of hydrogen halides such as HCl will be also improved.

Revisiting the stellar velocity ellipsoid-Hubble-type relation: observations versus simulations

NASA Astrophysics Data System (ADS)

Pinna, F.; Falcón-Barroso, J.; Martig, M.; Martínez-Valpuesta, I.; Méndez-Abreu, J.; van de Ven, G.; Leaman, R.; Lyubenova, M.

2018-04-01

The stellar velocity ellipsoid (SVE) in galaxies can provide important information on the processes that participate in the dynamical heating of their disc components (e.g. giant molecular clouds, mergers, spiral density waves, and bars). Earlier findings suggested a strong relation between the shape of the disc SVE and Hubble type, with later-type galaxies displaying more anisotropic ellipsoids and early types being more isotropic. In this paper, we revisit the strength of this relation using an exhaustive compilation of observational results from the literature on this issue. We find no clear correlation between the shape of the disc SVE and morphological type, and show that galaxies with the same Hubble type display a wide range of vertical-to-radial velocity dispersion ratios. The points are distributed around a mean value and scatter of σz/σR = 0.7 ± 0.2. With the aid of numerical simulations, we argue that different mechanisms might influence the shape of the SVE in the same manner and that the same process (e.g. mergers) does not have the same impact in all the galaxies. The complexity of the observational picture is confirmed by these simulations, which suggest that the vertical-to-radial axis ratio of the SVE is not a good indicator of the main source of disc heating. Our analysis of those simulations also indicates that the observed shape of the disc SVE may be affected by several processes simultaneously and that the signatures of some of them (e.g. mergers) fade over time.

Evaluation of cloud-resolving model simulations of midlatitude cirrus with ARM and A-train observations

DOE PAGES

Muhlbauer, A.; Ackerman, T. P.; Lawson, R. P.; ...

2015-07-14

Cirrus clouds are ubiquitous in the upper troposphere and still constitute one of the largest uncertainties in climate predictions. Our paper evaluates cloud-resolving model (CRM) and cloud system-resolving model (CSRM) simulations of a midlatitude cirrus case with comprehensive observations collected under the auspices of the Atmospheric Radiation Measurements (ARM) program and with spaceborne observations from the National Aeronautics and Space Administration A-train satellites. The CRM simulations are driven with periodic boundary conditions and ARM forcing data, whereas the CSRM simulations are driven by the ERA-Interim product. Vertical profiles of temperature, relative humidity, and wind speeds are reasonably well simulated bymore » the CSRM and CRM, but there are remaining biases in the temperature, wind speeds, and relative humidity, which can be mitigated through nudging the model simulations toward the observed radiosonde profiles. Simulated vertical velocities are underestimated in all simulations except in the CRM simulations with grid spacings of 500 m or finer, which suggests that turbulent vertical air motions in cirrus clouds need to be parameterized in general circulation models and in CSRM simulations with horizontal grid spacings on the order of 1 km. The simulated ice water content and ice number concentrations agree with the observations in the CSRM but are underestimated in the CRM simulations. The underestimation of ice number concentrations is consistent with the overestimation of radar reflectivity in the CRM simulations and suggests that the model produces too many large ice particles especially toward the cloud base. Simulated cloud profiles are rather insensitive to perturbations in the initial conditions or the dimensionality of the model domain, but the treatment of the forcing data has a considerable effect on the outcome of the model simulations. Despite considerable progress in observations and microphysical parameterizations

Simulating CO2 profiles using NIES TM and comparison with HIAPER Pole-to-Pole Observations

NASA Astrophysics Data System (ADS)

Song, C.; Maksyutov, S.; Belikov, D.; Takagi, H.; Shu, J.

2015-03-01

We present a study on validation of the National Institute for Environmental Studies Transport Model (NIES TM) by comparing to observed vertical profiles of atmospheric CO2. The model uses a hybrid sigma-isentropic (σ-θ) vertical coordinate that employs both terrain-following and isentropic parts switched smoothly in the stratosphere. The model transport is driven by reanalyzed meteorological fields and designed to simulate seasonal and diurnal cycles, synoptic variations, and spatial distributions of atmospheric chemical constituents in the troposphere. The model simulations were run for biosphere, fossil fuel, air-ocean exchange, biomass burning and inverse correction fluxes of carbon dioxide (CO2) by GOSAT Level 4 product. We compared the NIES TM simulated fluxes with data from the HIAPER Pole-to-Pole Observations (HIPPO) Merged 10 s Meteorology, Atmospheric Chemistry, and Aerosol Data, including HIPPO-1, HIPPO-2 and HIPPO-3 from 128.0° E to -84.0° W, and 87.0° N to -67.2° S. The simulation results were compared with CO2 observations made in January and November 2009, and March and April 2010. The analysis attests that the model is good enough to simulate vertical profiles with errors generally within 1-2 ppmv, except for the lower stratosphere in the Northern Hemisphere high latitudes.

Simulating flaring events in complex active regions driven by observed magnetograms

NASA Astrophysics Data System (ADS)

Dimitropoulou, M.; Isliker, H.; Vlahos, L.; Georgoulis, M. K.

2011-05-01

Context. We interpret solar flares as events originating in active regions that have reached the self organized critical state, by using a refined cellular automaton model with initial conditions derived from observations. Aims: We investigate whether the system, with its imposed physical elements, reaches a self organized critical state and whether well-known statistical properties of flares, such as scaling laws observed in the distribution functions of characteristic parameters, are reproduced after this state has been reached. Methods: To investigate whether the distribution functions of total energy, peak energy and event duration follow the expected scaling laws, we first applied a nonlinear force-free extrapolation that reconstructs the three-dimensional magnetic fields from two-dimensional vector magnetograms. We then locate magnetic discontinuities exceeding a threshold in the Laplacian of the magnetic field. These discontinuities are relaxed in local diffusion events, implemented in the form of cellular automaton evolution rules. Subsequent loading and relaxation steps lead the system to self organized criticality, after which the statistical properties of the simulated events are examined. Physical requirements, such as the divergence-free condition for the magnetic field vector, are approximately imposed on all elements of the model. Results: Our results show that self organized criticality is indeed reached when applying specific loading and relaxation rules. Power-law indices obtained from the distribution functions of the modeled flaring events are in good agreement with observations. Single power laws (peak and total flare energy) are obtained, as are power laws with exponential cutoff and double power laws (flare duration). The results are also compared with observational X-ray data from the GOES satellite for our active-region sample. Conclusions: We conclude that well-known statistical properties of flares are reproduced after the system has

Realistic Simulations of Coronagraphic Observations with Future Space Telescopes

NASA Astrophysics Data System (ADS)

Rizzo, M. J.; Roberge, A.; Lincowski, A. P.; Zimmerman, N. T.; Juanola-Parramon, R.; Pueyo, L.; Hu, M.; Harness, A.

2017-11-01

We present a framework to simulate realistic observations of future space-based coronagraphic instruments. This gathers state-of-the-art scientific and instrumental expertise allowing robust characterization of future instrument concepts.

The power spectrum of solar convection flows from high-resolution observations and 3D simulations

NASA Astrophysics Data System (ADS)

Yelles Chaouche, L.; Moreno-Insertis, F.; Bonet, J. A.

2014-03-01

instead of on horizontal planes. Results: A very good match between observational and simulated Fourier power spectra is obtained for the vertical velocity data for scales between 200 km and 6 Mm. Instead, a clear vertical shift is obtained when the synthetic observations are not degraded to emulate the degradation in the IMaX data. The match for the horizontal velocity data is much less impressive because of the inaccuracies of the LCT procedure. Concerning the internal comparison of the direct velocity values of the numerical boxes with those from the synthetic observations, a high correlation (0.96) is obtained for the vertical component when using the velocity values on the log τ500 = -1 surface in the box. The corresponding Fourier spectra are near each other. A lower maximum correlation (0.5) is reached (at log τ500 = 0) for the horizontal velocities as a result of the coarseness of the LCT procedure. Correspondingly, the Fourier spectra for the LCT-determined velocities is well below that from the actual velocity components. Conclusions: As measured by the Fourier spectra, realistic numerical simulations of surface magnetoconvection provide a very good match to the observational proxies for the photospheric velocity fields at least on scales from several Mm down to around 200 km. Taking into account the spatial and spectral instrumental blurring is essential for the comparison between simulations and observations. Dopplergrams are an excellent proxy for the vertical velocities on constant-τ isosurfaces, while LCT is a much less reliable method of determining the horizontal velocities.

Volcanic observation data and simulation database at NIED, Japan (Invited)

NASA Astrophysics Data System (ADS)

Fujita, E.; Ueda, H.; Kozono, T.

2009-12-01

NIED (Nat’l Res. Inst. for Earth Sci. & Disast. Prev.) has a project to develop two volcanic database systems: (1) volcanic observation database; (2) volcanic simulation database. The volcanic observation database is the data archive center obtained by the geophysical observation networks at Mt. Fuji, Miyake, Izu-Oshima, Iwo-jima and Nasu volcanoes, central Japan. The data consist of seismic (both high-sensitivity and broadband), ground deformation (tiltmeter, GPS) and those from other sensors (e.g., rain gauge, gravimeter, magnetometer, pressure gauge.) These data is originally stored in “WIN format,” the Japanese standard format, which is also at the Hi-net (High sensitivity seismic network Japan, http://www.hinet.bosai.go.jp/). NIED joins to WOVOdat and we have prepared to upload our data, via XML format. Our concept of the XML format is 1)a common format for intermediate files to upload into the WOVOdat DB, 2) for data files downloaded from the WOVOdat DB, 3) for data exchanges between observatories without the WOVOdat DB, 4) for common data files in each observatory, 5) for data communications between systems and softwares and 6)a for softwares. NIED is now preparing for (2) the volcanic simulation database. The objective of this project is to support to develop a “real-time” hazard map, i.e., the system which is effective to evaluate volcanic hazard in case of emergency, including the up-to-date conditions. Our system will include lava flow simulation (LavaSIM) and pyroclastic flow simulation (grvcrt). The database will keep many cases of assumed simulations and we can pick up the most probable case as the first evaluation in case the eruption started. The final goals of the both database will realize the volcanic eruption prediction and forecasting in real time by the combination of monitoring data and numerical simulations.

Thermal airborne multispectral aster simulator and its preliminary results

NASA Astrophysics Data System (ADS)

Mills, F.; Kannari, Y.; Watanabe, H.; Sano, M.; Chang, S. H.

1994-03-01

An Airborne ASTER Simulator (AAS) is being developed for the Japan Resources Observation System Organization (JAROS) by the Geophysical Environmental Research (GER) Corporation. The first test flights of the AAS were over Cuprite, Nevada; Long Valley, California; and Death Valley, California, in December 1991. Preliminary laboratory tests at NASA's Stennis Space Center (SSC) were completed in April 1992. The results of the these tests indicate the AAS can discriminate between silicate and non-silicate rocks. The improvements planned for the next two years may give a spectral Full-Width at Half-Maximum (FWHM) of 0.3 μm and NEΔT of 0.2 - 0.5°K. The AAS has the potential to become a good tool for airborne TIR research and can be used for simulations of future satellite-borne TIR sensors. Flight tests over Cuprite, Nevada, and Castaic Lake, California, are planned for October-December 1992.

Evaluation of Global Observations-Based Evapotranspiration Datasets and IPCC AR4 Simulations

NASA Technical Reports Server (NTRS)

Mueller, B.; Seneviratne, S. I.; Jimenez, C.; Corti, T.; Hirschi, M.; Balsamo, G.; Ciais, P.; Dirmeyer, P.; Fisher, J. B.; Guo, Z.;

Mapping and simulating systematics due to spatially-varying observing conditions in DES science verification data

DOE PAGES

Leistedt, B.; Peiris, H. V.; Elsner, F.; ...

2016-10-17

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. Wemore » illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We then analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on N(z), the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. However, they will need to be carefully characterised in upcoming phases of DES in order to avoid biasing the inferred cosmological results. The framework presented is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope, which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high

MAPPING AND SIMULATING SYSTEMATICS DUE TO SPATIALLY VARYING OBSERVING CONDITIONS IN DES SCIENCE VERIFICATION DATA

DOE Office of Scientific and Technical Information (OSTI.GOV)

Leistedt, B.; Peiris, H. V.; Elsner, F.

Spatially varying depth and the characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, particularly in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES-SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. Wemore » illustrate the complementary nature of these two approaches by comparing the SV data with BCC-UFig, a synthetic sky catalog generated by forward-modeling of the DES-SV images. We analyze the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and are well-captured by the maps of observing conditions. The combined use of the maps, the SV data, and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on N(z), the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak-lensing analyses. However, they will need to be carefully characterized in upcoming phases of DES in order to avoid biasing the inferred cosmological results. The framework presented here is relevant to all multi-epoch surveys and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope, which will require detailed null tests and realistic end-to-end image simulations to correctly interpret the deep, high

Mapping and simulating systematics due to spatially-varying observing conditions in DES science verification data

DOE Office of Scientific and Technical Information (OSTI.GOV)

Leistedt, B.; Peiris, H. V.; Elsner, F.

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. Wemore » illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We then analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on N(z), the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. However, they will need to be carefully characterised in upcoming phases of DES in order to avoid biasing the inferred cosmological results. The framework presented is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope, which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high

Mars boundary layer simulations - Comparison with Viking lander and entry observations

NASA Technical Reports Server (NTRS)

Haberle, R. M.; Houben, H. C.

1991-01-01

Diurnal variations of wind and temperature in the lower Martian atmosphere are simulated with a boundary layer model that includes radiative heating in a dusty CO2 atmosphere, turbulence generated by convection and/or shear stresses, a surface heat budget, and time varying pressure forces due to sloping terrain. Model results for early northern summer are compared with Viking lander observations to determine the model's strengths and weaknesses, and suitability as an engineering model.

GIANT MOLECULAR CLOUD FORMATION IN DISK GALAXIES: CHARACTERIZING SIMULATED VERSUS OBSERVED CLOUD CATALOGS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Benincasa, Samantha M.; Pudritz, Ralph E.; Wadsley, James

We present the results of a study of simulated giant molecular clouds (GMCs) formed in a Milky Way-type galactic disk with a flat rotation curve. This simulation, which does not include star formation or feedback, produces clouds with masses ranging between 10{sup 4} M{sub ☉} and 10{sup 7} M{sub ☉}. We compare our simulated cloud population to two observational surveys: the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey and the BIMA All-Disk Survey of M33. An analysis of the global cloud properties as well as a comparison of Larson's scaling relations is carried out. We find that simulatedmore » cloud properties agree well with the observed cloud properties, with the closest agreement occurring between the clouds at comparable resolution in M33. Our clouds are highly filamentary—a property that derives both from their formation due to gravitational instability in the sheared galactic environment, as well as to cloud-cloud gravitational encounters. We also find that the rate at which potentially star-forming gas accumulates within dense regions—wherein n{sub thresh} ≥ 10{sup 4} cm{sup –3}—is 3% per 10 Myr, in clouds of roughly 10{sup 6} M{sub ☉}. This suggests that star formation rates in observed clouds are related to the rates at which gas can be accumulated into dense subregions within GMCs via filamentary flows. The most internally well-resolved clouds are chosen for listing in a catalog of simulated GMCs—the first of its kind. The cataloged clouds are available as an extracted data set from the global simulation.« less

Partial Variance of Increments Method in Solar Wind Observations and Plasma Simulations

NASA Astrophysics Data System (ADS)

Greco, A.; Matthaeus, W. H.; Perri, S.; Osman, K. T.; Servidio, S.; Wan, M.; Dmitruk, P.

2018-02-01

The method called "PVI" (Partial Variance of Increments) has been increasingly used in analysis of spacecraft and numerical simulation data since its inception in 2008. The purpose of the method is to study the kinematics and formation of coherent structures in space plasmas, a topic that has gained considerable attention, leading the development of identification methods, observations, and associated theoretical research based on numerical simulations. This review paper will summarize key features of the method and provide a synopsis of the main results obtained by various groups using the method. This will enable new users or those considering methods of this type to find details and background collected in one place.

Observation and simulation of AGW in Space

NASA Astrophysics Data System (ADS)

Kunitsyn, Vyacheslav; Kholodov, Alexander; Andreeva, Elena; Nesterov, Ivan; Padokhin, Artem; Vorontsov, Artem

2014-05-01

Examples are presented of satellite observations and imaging of AGW and related phenomena in space travelling ionospheric disturbances (TID). The structure of AGW perturbations was reconstructed by satellite radio tomography (RT) based on the signals of Global Navigation Satellite Systems (GNSS). The experiments use different GNSS, both low-orbiting (Russian Tsikada and American Transit) and high-orbiting (GPS, GLONASS, Galileo, Beidou). The examples of RT imaging of TIDs and AGWs from anthropogenic sources such as ground explosions, rocket launching, heating the ionosphere by high-power radio waves are presented. In the latter case, the corresponding AGWs and TIDs were generated in response to the modulation in the power of the heating wave. The natural AGW-like wave disturbances are frequently observed in the atmosphere and ionosphere in the form of variations in density and electron concentration. These phenomena are caused by the influence of the near-space environment, atmosphere, and surface phenomena including long-period vibrations of the Earth's surface, earthquakes, explosions, temperature heating, seisches, tsunami waves, etc. Examples of experimental RT reconstructions of wave disturbances associated with the earthquakes and tsunami waves are presented, and RT images of TIDs caused by the variations in the corpuscular ionization are demonstrated. The results of numerical modeling of AGW generation by some surface and volume sources are discussed. The milli-Hertz AGWs generated by these sources induce perturbations with a typical scale of a few hundred of kilometers at the heights of the middle atmosphere and ionosphere. The numerical modeling is based on the solution of equations of geophysical hydrodynamics. The results of the numerical simulations agree with the observations. The authors acknowledge the support of the Russian Foundation for Basic Research (grants 14-05-00855 and 13-05-01122), grant of the President of Russian Federation MK-2670

Comparison of solar photospheric bright points between Sunrise observations and MHD simulations

NASA Astrophysics Data System (ADS)

Riethmüller, T. L.; Solanki, S. K.; Berdyugina, S. V.; Schüssler, M.; Martínez Pillet, V.; Feller, A.; Gandorfer, A.; Hirzberger, J.

2014-08-01

Bright points (BPs) in the solar photosphere are thought to be the radiative signatures (small-scale brightness enhancements) of magnetic elements described by slender flux tubes or sheets located in the darker intergranular lanes in the solar photosphere. They contribute to the ultraviolet (UV) flux variations over the solar cycle and hence may play a role in influencing the Earth's climate. Here we aim to obtain a better insight into their properties by combining high-resolution UV and spectro-polarimetric observations of BPs by the Sunrise Observatory with 3D compressible radiation magnetohydrodynamical (MHD) simulations. To this end, full spectral line syntheses are performed with the MHD data and a careful degradation is applied to take into account all relevant instrumental effects of the observations. In a first step it is demonstrated that the selected MHD simulations reproduce the measured distributions of intensity at multiple wavelengths, line-of-sight velocity, spectral line width, and polarization degree rather well. The simulated line width also displays the correct mean, but a scatter that is too small. In the second step, the properties of observed BPs are compared with synthetic ones. Again, these are found to match relatively well, except that the observations display a tail of large BPs with strong polarization signals (most likely network elements) not found in the simulations, possibly due to the small size of the simulation box. The higher spatial resolution of the simulations has a significant effect, leading to smaller and more numerous BPs. The observation that most BPs are weakly polarized is explained mainly by the spatial degradation, the stray light contamination, and the temperature sensitivity of the Fe i line at 5250.2 Å. Finally, given that the MHD simulations are highly consistent with the observations, we used the simulations to explore the properties of BPs further. The Stokes V asymmetries increase with the distance to the

Comparison of Tropical and Extratropical Gust Factors Using Observed and Simulated Data

NASA Astrophysics Data System (ADS)

Edwards, R. P.; Schroeder, J. L.

2011-12-01

Questions of whether differences exist between tropical cyclone (TC) and extratropical (ET) wind have been the subject of considerable debate. This study will focus on the behavior of the gust factor (GF), the ratio of a peak wind speed of a certain duration and a mean wind speed of a certain duration, for three types of data: TC, ET, and simulated. For this project, the Universal Spectrum, a normalized, averaged spectrum for wind, was un-normalized and used to create simulated wind speed time series at a variety of wind speeds. Additional time series were created after modifying the spectrum to simulate the additional low-frequency energy observed in the TC wind spectrum as well as the reduction of high-frequency energy caused by a mechanical anemometer. The T and ET data used for this study were collected by Texas Tech University's mobile towers as part of various field efforts since 1998. Before comparisons were made, the database was divided into four roughness regimes based on the roughness length to ensure that differences observed in the turbulence statistics are not caused by differences in upstream terrain. The mean GF for the TC data set (open roughness regime), 1.49, was slightly higher than the ET value of 1.44 (Table 1). The distributions of GFs from each data type show similarities in shape between the base-simulated and ET data sets and between the TC and modified-simulated data set (Figure 1). These similarities are expected given the spectral similarities between the TC and modified-simulated data sets, namely additional low-frequency energy relative to the ET and base-simulated data. These findings suggest that the higher amount of low-frequency energy present in the tropical wind spectrum is partially responsible for the resulting higher GF for the tropical cyclone data. However, the modest increase in GF from the base to the modified simulated data suggest that there are more factors at work.

Large-Scale Dynamics of the Magnetospheric Boundary: Comparisons between Global MHD Simulation Results and ISTP Observations

NASA Technical Reports Server (NTRS)

Berchem, J.; Raeder, J.; Ashour-Abdalla, M.; Frank, L. A.; Paterson, W. R.; Ackerson, K. L.; Kokubun, S.; Yamamoto, T.; Lepping, R. P.

1998-01-01

Understanding the large-scale dynamics of the magnetospheric boundary is an important step towards achieving the ISTP mission's broad objective of assessing the global transport of plasma and energy through the geospace environment. Our approach is based on three-dimensional global magnetohydrodynamic (MHD) simulations of the solar wind-magnetosphere- ionosphere system, and consists of using interplanetary magnetic field (IMF) and plasma parameters measured by solar wind monitors upstream of the bow shock as input to the simulations for predicting the large-scale dynamics of the magnetospheric boundary. The validity of these predictions is tested by comparing local data streams with time series measured by downstream spacecraft crossing the magnetospheric boundary. In this paper, we review results from several case studies which confirm that our MHD model reproduces very well the large-scale motion of the magnetospheric boundary. The first case illustrates the complexity of the magnetic field topology that can occur at the dayside magnetospheric boundary for periods of northward IMF with strong Bx and By components. The second comparison reviewed combines dynamic and topological aspects in an investigation of the evolution of the distant tail at 200 R(sub E) from the Earth.

Aircraft microwave observations and simulations of deep convection from 18 to 183 GHz. II - Model results

NASA Technical Reports Server (NTRS)

Yeh, Hwa-Young M.; Prasad, N.; Mack, Robert A.; Adler, Robert F.

1990-01-01

In this June 29, 1986 case study, a radiative transfer model is used to simulate the aircraft multichannel microwave brightness temperatures presented in the Adler et al. (1990) paper and to study the convective storm structure. Ground-based radar data are used to derive hydrometeor profiles of the storm, based on which the microwave upwelling brightness temperatures are calculated. Various vertical hydrometeor phase profiles and the Marshall and Palmer (M-P, 1948) and Sekhon and Srivastava (S-S, 1970) ice particle size distributions are experimented in the model. The results are compared with the aircraft radiometric data. The comparison reveals that the M-P distribution well represents the ice particle size distribution, especially in the upper tropospheric portion of the cloud; the S-S distribution appears to better simulate the ice particle size at the lower portion of the cloud, which has a greater effect on the low-frequency microwave upwelling brightness temperatures; and that, in deep convective regions, significant supercooled liquid water (about 0.5 g/cu m) may be present up to the -30 C layer, while in less convective areas, frozen hydrometeors are predominant above -10 C level.

Evaluation of NCAR CAM5 Simulated Marine Boundary Layer Cloud Properties Using a Combination of Satellite and Surface Observations

NASA Astrophysics Data System (ADS)

Zhang, Z.; Song, H.; Wang, M.; Ghan, S. J.; Dong, X.

2016-12-01

he main objective of this study is to systematically evaluate the MBL cloud properties simulated in CAM5 family models using a combination of satellite-based CloudSat/MODIS observations and ground-based observations from the ARM Azores site, with a special focus on MBL cloud microphysics and warm rain process. First, we will present a global evaluation based on satellite observations and retrievals. We will compare global cloud properties (e.g., cloud fraction, cloud vertical structure, cloud CER, COT, and LWP, as well as drizzle frequency and intensity diagnosed using the CAM5-COSP instrumental simulators) simulated in the CAM5 models with the collocated CloudSat and MODIS observations. We will also present some preliminary results from a regional evaluation based mainly on ground observations from ARM Azores site. We will compare MBL cloud properties simulated in CAM5 models over the ARM Azores site with collocated satellite (MODIS and CloudSat) and ground-based observations from the ARM site.

Ram Pressure Stripping: Observations Meet Simulations

NASA Astrophysics Data System (ADS)

Past, Matthew; Ruszkowski, Mateusz; Sharon, Keren

2017-01-01

Ram pressure stripping occurs when a galaxy falls into the potential well of a cluster, removing gas and dust as the galaxy travels through the intracluster medium. This interaction leads to filamentary gas tails stretching behind the galaxy and plays an important role in galaxy evolution. Previously, these “jellyfish” galaxies had only been observed in nearby clusters, but recently, higher redshift (z > 0.3) examples have been found from HST data imaging.Recent work has shown that cosmic rays injected by supernovae can cause galactic disks to thicken due to cosmic ray pressure. We run three-dimensional magneto-hydrodynamical simulations of ram pressure stripping including cosmic rays to compare to previous models. We study how the efficiency of the ram pressure stripping of the gas, and the morphology of the filamentary tails, depend on the magnitude of the cosmic ray pressure support. We generate mock X-ray images and radio polarization data. Simultaneously, we perform an exhaustive search of the HST archive to increase the sample of jellyfish galaxies and compare selected cases to simulations.

Retrieved Products from Simulated Hyperspectral Observations of a Hurricane

NASA Technical Reports Server (NTRS)

Susskind, Joel; Kouvaris, Louis; Iredell, Lena; Blaisdell, John

2015-01-01

Retrievals were run using the AIRS Science Team Version-6 AIRS Only retrieval algorithm, which generates a Neural-Net first guess (T(sub s))(sup 0), (T(p))(sup 0), and (q(p))(sup 0) as a function of observed AIRS radiances. AIRS Science Team Neural-Net coefficients performed very well beneath 300 mb using the simulated radiances. This means the simulated radiances are very realistic. First guess and retrieved values of T(p) above 300 mb were biased cold, but both represented the model spatial structure very well. QC'd T(p) and q(p) retrievals for all experiments had similar accuracies compared to their own truth fields, and were roughly consistent with results obtained using real data. Spatial coverage of retrievals, as well as the representativeness of the spatial structure of the storm, improved dramatically with decreasing size of the instrument's FOV. We sent QC'd values of T(p) and q(p) to Bob Atlas at AOML for use as input to OSSE Data Assimilation experiments.

Infrared Spectral Observations While Drilling into a Frozen Lunar Simulant

NASA Technical Reports Server (NTRS)

Roush, Ted L.; Colaprete, Anthony; Thompson, Sarah; Cook, Amanda; Kleinhenz, Julie

2014-01-01

Past and continuing observations indicate an enrichment of volatile materials in lunar polar regions. While these volatiles may be located near the surface, access to them will likely require subsurface sampling, during which it is desirable to monitor the volatile content. In a simulation of such activities, a multilayer lunar simulant was prepared with differing water content, and placed inside a thermal vacuum chamber at Glenn Research Center (GRC). The soil profile was cooled using liquid nitrogen. In addition to the soil, a drill and infrared (IR) spectrometer (1600-3400 nm) were also located in the GRC chamber. We report the spectral observations obtained during a sequence where the drill was repeatedly inserted and extracted, to different depths, at the same location. We observe an overall increase in the spectral signature of water ice over the duration of the test. Additionally, we observe variations in the water ice spectral signature as the drill encounters different layers.

Simulation of keratoconus observation in photorefraction

NASA Astrophysics Data System (ADS)

Chen, Ying-Ling; Tan, B.; Baker, K.; Lewis, J. W. L.; Swartz, T.; Jiang, Y.; Wang, M.

2006-11-01

In the recent years, keratoconus (KC) has increasingly gained attention due to its treatment options and to the popularity of keratorefractive surgery. This paper investigates the potential of identification of KC using photorefraction (PR), an optical technique that is similar to objective retinoscopy and is commonly used for large-scale ocular screening. Using personalized eye models of both KC and pre-LASIK patients, computer simulations were performed to achieve visualization of this ophthalmic measurement. The simulations are validated by comparing results to two sets of experimental measurements. These PR images show distinguishable differences between KC eyes and eyes that are either normal or ametropic. The simulation technique with personalized modeling can be extended to other ophthalmic instrument developments. It makes possible investigation with the least number of real human subjects. The application is also of great interest in medical training.

3-D VPIC simulation of an vortex-induced reconnection event observed by MMS

DOE Data Explorer

Nakamura, Takuma; Daughton, William

2016-01-01

The data set consists of a 3-D fully kinetic (VPIC) simulation of an in-situ observation event at the Earth's magnetopause by the NASA MMS spacecraft on September 8, 2015. The results show a turbulent development of magnetic reconnection induced by the Kelvin-Helmohltz vortex, and resulting significantly efficient plasma mixing across the magnetopause. The vortex-induced reconnection signatures are well consistent with the MMS observations. These results are published in some scientific journals such as Nature Communications. Fortran unformatted files with 1024x1536x512 cells, which have been compressed from original ones with 2048x3072x1024 cells, are archived for selected time slices of field and moment data shown in these papers.

Comparisons of Spectra from 3D Kinetic Meteor PIC Simulations with Theory and Observations

NASA Astrophysics Data System (ADS)

Oppenheim, M. M.; Tarnecki, L. K.

2017-12-01

Meteoroids smaller than a grain of sand have significant impacts on the composition, chemistry, and dynamics of the atmosphere. The processes by which they turbulently diffuse can be studied using collisional kinetic particle-in-cell (PIC) simulations. Spectral analysis is a valuable tool for comparing such simulations of turbulent, non-specular meteor trails with observations. We present three types of spectral information: full spectra along the trail in k-ω space, spectral widths at common radar frequencies, and power as a function of angle with respect to B. These properties can be compared to previously published data. Zhou et al. (2004) use radar theory to predict the power observed by a radar as a function of the angle between the meteor trail and the radar beam and the size of field-aligned irregularities (FAI) within the trail. Close et al. (2008) present observations of meteor trails from the ALTAIR radar, including power returned as a function of angle off B for a small sample of meteors. Close et al. (2008) and Zhou et al. (2004) both suggest a power drop off of 2-3 dB per degree off perpendicular to B. We compare results from our simulations with both theory and observations for a range of conditions, including trail altitude and incident neutral wind speed. For 1m waves, power fell off by 1-3 dB per degree off perpendicular to B. These comparisons help determine if small-scale simulations accurately capture the behavior of real meteors.

The role of historical forcings in simulating the observed Atlantic multidecadal oscillation

NASA Astrophysics Data System (ADS)

Murphy, Lisa N.; Bellomo, Katinka; Cane, Mark; Clement, Amy

2017-03-01

We analyze the Atlantic multidecadal oscillation (AMO) in the preindustrial (PI) and historical (HIST) simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) to assess the drivers of the observed AMO from 1865 to 2005. We draw 141 year samples from the 41 CMIP5 model's PI runs and compare the correlation and variance between the observed AMO and the simulated PI and HIST AMO. The correlation coefficients in 38 forced (HIST) models are above the 90% confidence level and explain up to 56% of the observed variance. The probability that any of the unforced (PI) models do as well is less than 3% in 31 models. Multidecadal variability is larger in 39 CMIP5 HIST simulations and in all HIST members of the Community Earth System Model Large Ensemble than their corresponding PI. We conclude that there is an essential role for external forcing in driving the observed AMO.

Low-frequency waves at comet 67P/Churyumov-Gerasimenko. Observations compared to numerical simulations

NASA Astrophysics Data System (ADS)

Koenders, C.; Perschke, C.; Goetz, C.; Richter, I.; Motschmann, U.; Glassmeier, K. H.

2016-10-01

Context. A new type of low-frequency wave was detected by the magnetometer of the Rosetta Plasma Consortium at the comet during the initial months after the arrival of the Rosetta spacecraft at comet 67P/Churyumov-Gerasimenko. This large-amplitude, nearly continuous wave activity is observed in the frequency range from 30 mHz to 80 mHz where 40 mHz to 50 mHz is the dominant frequency. This type of low frequency is not closely related to the gyrofrequency of newborn cometary ions, which differs from previous wave activity observed in the interaction region of comets with the solar wind. Aims: This work aims to reveal a global view on the wave activity region using simulations of the comet-solar wind interaction region. Parameters, such as wavelength, propagation direction, and propagation patterns, are within the focus of this study. While the Rosetta observations only provide local information, numerical simulations provide further information on the global wave properties. Methods: Standard hybrid simulations were applied to the comet-solar wind interaction scenario. In the model, the ions were described as particles, which allows us to describe kinetic processes of the ions. The electrons were described as a fluid. Results: The simulations exhibit a threefold wave structure of the interaction region. A Mach cone and a Whistler wing are observed downstream of the comet. The third kind of wave activity found are low-frequency waves at 97 mHz, which corresponds to the waves observed by Richter et al. (2015, Ann. Geophys., 33, 1031). These waves are caused by the initial pick-up of the cometary ions that are perpendicular to the solar wind flow and in the interplanetary magnetic field direction. The associated electric current becomes unstable. The simulations show that wave activity is only detectable in the + E hemisphere and that the Mach cone and whistler wings need to be distinguished from the newly found instability driven wave activity. The movie associated to

Insights into the use of time-lapse GPR data as observations for inverse multiphase flow simulations of DNAPL migration

USGS Publications Warehouse

Johnson, R.H.; Poeter, E.P.

2007-01-01

Perchloroethylene (PCE) saturations determined from GPR surveys were used as observations for inversion of multiphase flow simulations of a PCE injection experiment (Borden 9??m cell), allowing for the estimation of optimal bulk intrinsic permeability values. The resulting fit statistics and analysis of residuals (observed minus simulated PCE saturations) were used to improve the conceptual model. These improvements included adjustment of the elevation of a permeability contrast, use of the van Genuchten versus Brooks-Corey capillary pressure-saturation curve, and a weighting scheme to account for greater measurement error with larger saturation values. A limitation in determining PCE saturations through one-dimensional GPR modeling is non-uniqueness when multiple GPR parameters are unknown (i.e., permittivity, depth, and gain function). Site knowledge, fixing the gain function, and multiphase flow simulations assisted in evaluating non-unique conceptual models of PCE saturation, where depth and layering were reinterpreted to provide alternate conceptual models. Remaining bias in the residuals is attributed to the violation of assumptions in the one-dimensional GPR interpretation (which assumes flat, infinite, horizontal layering) resulting from multidimensional influences that were not included in the conceptual model. While the limitations and errors in using GPR data as observations for inverse multiphase flow simulations are frustrating and difficult to quantify, simulation results indicate that the error and bias in the PCE saturation values are small enough to still provide reasonable optimal permeability values. The effort to improve model fit and reduce residual bias decreases simulation error even for an inversion based on biased observations and provides insight into alternate GPR data interpretations. Thus, this effort is warranted and provides information on bias in the observation data when this bias is otherwise difficult to assess. ?? 2006 Elsevier B

ICE-COLA: fast simulations for weak lensing observables

NASA Astrophysics Data System (ADS)

Izard, Albert; Fosalba, Pablo; Crocce, Martin

2018-01-01

Approximate methods to full N-body simulations provide a fast and accurate solution to the development of mock catalogues for the modelling of galaxy clustering observables. In this paper we extend ICE-COLA, based on an optimized implementation of the approximate COLA method, to produce weak lensing maps and halo catalogues in the light-cone using an integrated and self-consistent approach. We show that despite the approximate dynamics, the catalogues thus produced enable an accurate modelling of weak lensing observables one decade beyond the characteristic scale where the growth becomes non-linear. In particular, we compare ICE-COLA to the MICE Grand Challenge N-body simulation for some fiducial cases representative of upcoming surveys and find that, for sources at redshift z = 1, their convergence power spectra agree to within 1 per cent up to high multipoles (i.e. of order 1000). The corresponding shear two point functions, ξ+ and ξ-, yield similar accuracy down to 2 and 20 arcmin respectively, while tangential shear around a z = 0.5 lens sample is accurate down to 4 arcmin. We show that such accuracy is stable against an increased angular resolution of the weak lensing maps. Hence, this opens the possibility of using approximate methods for the joint modelling of galaxy clustering and weak lensing observables and their covariance in ongoing and future galaxy surveys.

UAS in the NAS - Analysis Results and Recommendations for Integration of CNPC and ATC Communications Simulation Report

NASA Technical Reports Server (NTRS)

Kubat, Gregory

2016-01-01

This report addresses a deliverable to the UAS-in-the-NAS project for recommendations for integration of CNPC and ATC communications based on analysis results from modeled radio system and NAS-wide UA communication architecture simulations. For each recommendation, a brief explanation of the rationale for its consideration is provided with any supporting results obtained or observed in our simulation activity.

Evaluating Clouds in Long-Term Cloud-Resolving Model Simulations with Observational Data

NASA Technical Reports Server (NTRS)

Zeng, Xiping; Tao, Wei-Kuo; Zhang, Minghua; Peters-Lidard, Christa; Lang, Stephen; Simpson, Joanne; Kumar, Sujay; Xie, Shaocheng; Eastman, Joseph L.; Shie, Chung-Lin;

CLARREO shortwave observing system simulation experiments of the twenty-first century: Simulator design and implementation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Feldman, D.R.; Algieri, C.A.; Ong, J.R.

2011-04-01

Projected changes in the Earth system will likely be manifested in changes in reflected solar radiation. This paper introduces an operational Observational System Simulation Experiment (OSSE) to calculate the signals of future climate forcings and feedbacks in top-of-atmosphere reflectance spectra. The OSSE combines simulations from the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report for the NCAR Community Climate System Model (CCSM) with the MODTRAN radiative transfer code to calculate reflectance spectra for simulations of current and future climatic conditions over the 21st century. The OSSE produces narrowband reflectances and broadband fluxes, the latter of which have been extensivelymore » validated against archived CCSM results. The shortwave reflectance spectra contain atmospheric features including signals from water vapor, liquid and ice clouds, and aerosols. The spectra are also strongly influenced by the surface bidirectional reflectance properties of predicted snow and sea ice and the climatological seasonal cycles of vegetation. By comparing and contrasting simulated reflectance spectra based on emissions scenarios with increasing projected and fixed present-day greenhouse gas and aerosol concentrations, we find that prescribed forcings from increases in anthropogenic sulfate and carbonaceous aerosols are detectable and are spatially confined to lower latitudes. Also, changes in the intertropical convergence zone and poleward shifts in the subsidence zones and the storm tracks are all detectable along with large changes in snow cover and sea ice fraction. These findings suggest that the proposed NASA Climate Absolute Radiance and Refractivity Observatory (CLARREO) mission to measure shortwave reflectance spectra may help elucidate climate forcings, responses, and feedbacks.« less

Simulation and 'TWINS Observations of the 22 July 2009 Storm

NASA Technical Reports Server (NTRS)

Fok, Mei-Ching; Buzulukova, Natalia Y.; Chen, Sheng-Hsien; Valek, Phil; Goldstein, Jerry; McComas, David

2010-01-01

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 is the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly -80nT at 7:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the ENA emissions was seen in the midnight-to-dawn local-time sector. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in-situ ion data from THEMIS satellites. Simulation with static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always eastward than those in the observed TWINS images. On the other hand, simulation with a dynamic magnetic field gives better spatial agreements with both ENA and insitu particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Simulation and Twins Observations of the 22 July 2009 Storm

NASA Technical Reports Server (NTRS)

Fok, M.-C.; Buzulukova, N.; Chen, S.-H.; Valek, P. W.; Goldstein, J.; McComas, D. J.

2011-01-01

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 was at the time the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly .80 nT at 0700 and 1000 UT. During the main phase, and at the peak of the storm, TWINS 1 and 2 were near apogee and moving between predawn and postdawn local time. The energetic neutral atom (ENA) imagers on the two spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the high-altitude ENA emissions was seen in the midnight-to-dawn local time sector. The development of this storm has been simulated using the comprehensive ring current model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in situ ion data from a THEMIS satellite. Simulation with a static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always more eastward than those in the observed TWINS images. On the other hand, the simulation with a dynamic magnetic field gives better spatial agreement with both ENA and in situ particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Marvel-ous Dwarfs: Results from Four Heroically Large Simulated Volumes of Dwarf Galaxies

NASA Astrophysics Data System (ADS)

Munshi, Ferah; Brooks, Alyson; Weisz, Daniel; Bellovary, Jillian; Christensen, Charlotte

2018-01-01

We present results from high resolution, fully cosmological simulations of cosmic sheets that contain many dwarf galaxies. Together, they create the largest collection of simulated dwarf galaxies to date, with z=0 stellar masses comparable to the LMC or smaller. In total, we have simulated almost 100 luminous dwarf galaxies, forming a sample of simulated dwarfs which span a wide range of physical (stellar and halo mass) and evolutionary properties (merger history). We show how they can be calibrated against a wealth of observations of nearby galaxies including star formation histories, HI masses and kinematics, as well as stellar metallicities. We present preliminary results answering the following key questions: What is the slope of the stellar mass function at extremely low masses? Do halos with HI and no stars exist? What is the scatter in the stellar to halo mass relationship as a function of dwarf mass? What drives the scatter? With this large suite, we are beginning to statistically characterize dwarf galaxies and identify the types and numbers of outliers to expect.

Cardiac looping may be driven by compressive loads resulting from unequal growth of the heart and pericardial cavity. Observations on a physical simulation model

PubMed Central

Bayraktar, Meriç; Männer, Jörg

2014-01-01

The transformation of the straight embryonic heart tube into a helically wound loop is named cardiac looping. Such looping is regarded as an essential process in cardiac morphogenesis since it brings the building blocks of the developing heart into an approximation of their definitive topographical relationships. During the past two decades, a large number of genes have been identified which play important roles in cardiac looping. However, how genetic information is physically translated into the dynamic form changes of the looping heart is still poorly understood. The oldest hypothesis of cardiac looping mechanics attributes the form changes of the heart loop (ventral bending → simple helical coiling → complex helical coiling) to compressive loads resulting from growth differences between the heart and the pericardial cavity. In the present study, we have tested the physical plausibility of this hypothesis, which we call the growth-induced buckling hypothesis, for the first time. Using a physical simulation model, we show that growth-induced buckling of a straight elastic rod within the confined space of a hemispherical cavity can generate the same sequence of form changes as observed in the looping embryonic heart. Our simulation experiments have furthermore shown that, under bilaterally symmetric conditions, growth-induced buckling generates left- and right-handed helices (D-/L-loops) in a 1:1 ratio, while even subtle left- or rightward displacements of the caudal end of the elastic rod at the pre-buckling state are sufficient to direct the buckling process toward the generation of only D- or L-loops, respectively. Our data are discussed with respect to observations made in biological “models.” We conclude that compressive loads resulting from unequal growth of the heart and pericardial cavity play important roles in cardiac looping. Asymmetric positioning of the venous heart pole may direct these forces toward a biased generation of D- or L-loops. PMID

Evaluation of Cloud-Resolving Model Intercomparison Simulations Using TWP-ICE Observations: Precipitation and Cloud Structure

NASA Technical Reports Server (NTRS)

Varble, Adam; Fridlind, Ann M.; Zipser, Edward J.; Ackerman, Andrew S.; Chaboureau, Jean-Pierre; Fan, Jiwen; Hill, Adrian; McFarlane, Sally A.; Pinty, Jean-Pierre; Shipway, Ben

2011-01-01

The Tropical Warm Pool.International Cloud Experiment (TWP ]ICE) provided extensive observational data sets designed to initialize, force, and constrain atmospheric model simulations. In this first of a two ]part study, precipitation and cloud structures within nine cloud ]resolving model simulations are compared with scanning radar reflectivity and satellite infrared brightness temperature observations during an active monsoon period from 19 to 25 January 2006. Seven of nine simulations overestimate convective area by 20% or more leading to general overestimation of convective rainfall. This is balanced by underestimation of stratiform rainfall by 5% to 50% despite overestimation of stratiform area by up to 65% because of a preponderance of very low stratiform rain rates in all simulations. All simulations fail to reproduce observed radar reflectivity distributions above the melting level in convective regions and throughout the troposphere in stratiform regions. Observed precipitation ]sized ice reaches higher altitudes than simulated precipitation ]sized ice despite some simulations that predict lower than observed top ]of ]atmosphere infrared brightness temperatures. For the simulations that overestimate radar reflectivity aloft, graupel is the cause with one ]moment microphysics schemes whereas snow is the cause with two ]moment microphysics schemes. Differences in simulated radar reflectivity are more highly correlated with differences in mass mean melted diameter (Dm) than differences in ice water content. Dm is largely dependent on the mass ]dimension relationship and gamma size distribution parameters such as size intercept (N0) and shape parameter (m). Having variable density, variable N0, or m greater than zero produces radar reflectivities closest to those observed.

Arctic daily temperature and precipitation extremes: Observed and simulated physical behavior

NASA Astrophysics Data System (ADS)

Glisan, Justin Michael

Simulations using a six-member ensemble of Pan-Arctic WRF (PAW) were produced on two Arctic domains with 50-km resolution to analyze precipitation and temperature extremes for various periods. The first study used a domain developed for the Regional Arctic Climate Model (RACM). Initial simulations revealed deep atmospheric circulation biases over the northern Pacific Ocean, manifested in pressure, geopotential height, and temperature fields. Possible remedies to correct these large biases, such as modifying the physical domain or using different initial/boundary conditions, were unsuccessful. Spectral (interior) nudging was introduced as a way of constraining the model to be more consistent with observed behavior. However, such control over numerical model behavior raises concerns over how much nudging may affect unforced variability and extremes. Strong nudging may reduce or filter out extreme events, since the nudging pushes the model toward a relatively smooth, large-scale state. The question then becomes---what is the minimum spectral nudging needed to correct biases while not limiting the simulation of extreme events? To determine this, we use varying degrees of spectral nudging, using WRF's standard nudging as a reference point during January and July 2007. Results suggest that there is a marked lack of sensitivity to varying degrees of nudging. Moreover, given that nudging is an artificial forcing applied in the model, an important outcome of this work is that nudging strength apparently can be considerably smaller than WRF's standard strength and still produce reliable simulations. In the remaining studies, we used the same PAW setup to analyze daily precipitation extremes simulated over a 19-year period on the CORDEX Arctic domain for winter and summer. We defined these seasons as the three-month period leading up to and including the climatological sea ice maximum and minimum, respectively. Analysis focused on four North American regions defined using

The Photochemical Reflectance Index from Directional Cornfield Reflectances: Observations and Simulations

NASA Technical Reports Server (NTRS)

Cheng, Yen-Ben; Middleton, Elizabeth M.; Zhang, Qingyuan; Corp, Lawrence A.; Dandois, Jonathan; Kustas, William P.

2012-01-01

The two-layer Markov chain Analytical Canopy Reflectance Model (ACRM) was linked with in situ hyperspectral leaf optical properties to simulate the Photochemical Reflectance Index (PRI) for a corn crop canopy at three different growth stages. This is an extended study after a successful demonstration of PRI simulations for a cornfield previously conducted at an early vegetative growth stage. Consistent with previous in situ studies, sunlit leaves exhibited lower PRI values than shaded leaves. Since sunlit (shaded) foliage dominates the canopy in the reflectance hotspot (coldspot), the canopy PRI derived from field hyperspectral observations displayed sensitivity to both view zenith angle and relative azimuth angle at all growth stages. Consequently, sunlit and shaded canopy sectors were most differentiated when viewed along the azimuth matching the solar principal plane. These directional PRI responses associated with sunlit/shaded foliage were successfully reproduced by the ACRM. As before, the simulated PRI values from the current study were closer to in situ values when both sunlit and shaded leaves were utilized as model input data in a two-layer mode, instead of a one-layer mode with sunlit leaves only. Model performance as judged by correlation between in situ and simulated values was strongest for the mature corn crop (r = 0.87, RMSE = 0.0048), followed by the early vegetative stage (r = 0.78; RMSE = 0.0051) and the early senescent stage (r = 0.65; RMSE = 0.0104). Since the benefit of including shaded leaves in the scheme varied across different growth stages, a further analysis was conducted to investigate how variable fractions of sunlit/shaded leaves affect the canopy PRI values expected for a cornfield, with implications for 20 remote sensing monitoring options. Simulations of the sunlit to shaded canopy ratio near 50/50 +/- 10 (e.g., 60/40) matching field observations at all growth stages were examined. Our results suggest in the importance of the

Synchrony between reanalysis-driven RCM simulations and observations: variation with time scale

NASA Astrophysics Data System (ADS)

de Elía, Ramón; Laprise, René; Biner, Sébastien; Merleau, James

2017-04-01

Unlike coupled global climate models (CGCMs) that run in a stand-alone mode, nested regional climate models (RCMs) are driven by either a CGCM or a reanalysis dataset. This feature makes high correlations between the RCM simulation and its driver possible. When the driving dataset is a reanalysis, time correlations between RCM output and observations are also common and to be expected. In certain situations time correlation between driver and driven RCM is of particular interest and techniques have been developed to increase it (e.g. large-scale spectral nudging). For such cases, a question that remains open is whether aggregating in time increases the correlation between RCM output and observations. That is, although the RCM may be unable to reproduce a given daily event, whether it will still be able to satisfactorily simulate an anomaly on a monthly or annual basis. This is a preconception that the authors of this work and others in the community have held, perhaps as a natural extension of the properties of upscaling or aggregating other statistics such as the mean squared error. Here we explore analytically four particular cases that help us partially answer this question. In addition, we use observations datasets and RCM-simulated data to illustrate our findings. Results indicate that time upscaling does not necessarily increase time correlations, and that those interested in achieving high monthly or annual time correlations between RCM output and observations may have to do so by increasing correlation as much as possible at the shortest time scale. This may indicate that even when only concerned with time correlations at large temporal scale, large-scale spectral nudging acting at the time-step level may have to be used.

Virtual Observation System for Earth System Model: An Application to ACME Land Model Simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Dali; Yuan, Fengming; Hernandez, Benjamin

Investigating and evaluating physical-chemical-biological processes within an Earth system model (EMS) can be very challenging due to the complexity of both model design and software implementation. A virtual observation system (VOS) is presented to enable interactive observation of these processes during system simulation. Based on advance computing technologies, such as compiler-based software analysis, automatic code instrumentation, and high-performance data transport, the VOS provides run-time observation capability, in-situ data analytics for Earth system model simulation, model behavior adjustment opportunities through simulation steering. A VOS for a terrestrial land model simulation within the Accelerated Climate Modeling for Energy model is also presentedmore » to demonstrate the implementation details and system innovations.« less

Virtual Observation System for Earth System Model: An Application to ACME Land Model Simulations

DOE PAGES

Wang, Dali; Yuan, Fengming; Hernandez, Benjamin; ...

2017-01-01

Investigating and evaluating physical-chemical-biological processes within an Earth system model (EMS) can be very challenging due to the complexity of both model design and software implementation. A virtual observation system (VOS) is presented to enable interactive observation of these processes during system simulation. Based on advance computing technologies, such as compiler-based software analysis, automatic code instrumentation, and high-performance data transport, the VOS provides run-time observation capability, in-situ data analytics for Earth system model simulation, model behavior adjustment opportunities through simulation steering. A VOS for a terrestrial land model simulation within the Accelerated Climate Modeling for Energy model is also presentedmore » to demonstrate the implementation details and system innovations.« less

WRF-Chem Simulations of Lightning-NOx Production and Transport in Oklahoma and Colorado Thunderstorms Observed During DC3

NASA Technical Reports Server (NTRS)

Cummings, Kristin A.; Pickering, Kenneth E.; Barth, M.; Bela, M.; Li, Y.; Allen, D.; Bruning, E.; MacGorman, D.; Rutledge, S.; Basarab, B.;

Plasma and wave properties downstream of Martian bow shock: Hybrid simulations and MAVEN observations

NASA Astrophysics Data System (ADS)

Dong, Chuanfei; Winske, Dan; Cowee, Misa; Bougher, Stephen W.; Andersson, Laila; Connerney, Jack; Epley, Jared; Ergun, Robert; McFadden, James P.; Ma, Yingjuan; Toth, Gabor; Curry, Shannon; Nagy, Andrew; Jakosky, Bruce

2015-04-01

Two-dimensional hybrid simulation codes are employed to investigate the kinetic properties of plasmas and waves downstream of the Martian bow shock. The simulations are two-dimensional in space but three dimensional in field and velocity components. Simulations show that ion cyclotron waves are generated by temperature anisotropy resulting from the reflected protons around the Martian bow shock. These proton cyclotron waves could propagate downward into the Martian ionosphere and are expected to heat the O+ layer peaked from 250 to 300 km due to the wave-particle interaction. The proton cyclotron wave heating is anticipated to be a significant source of energy into the thermosphere, which impacts atmospheric escape rates. The simulation results show that the specific dayside heating altitude depends on the Martian crustal field orientations, solar cycles and seasonal variations since both the cyclotron resonance condition and the non/sub-resonant stochastic heating threshold depend on the ambient magnetic field strength. The dayside magnetic field profiles for different crustal field orientation, solar cycle and seasonal variations are adopted from the BATS-R-US Mars multi-fluid MHD model. The simulation results, however, show that the heating of O+ via proton cyclotron wave resonant interaction is not likely in the relatively weak crustal field region, based on our simplified model. This indicates that either the drift motion resulted from the transport of ionospheric O+, or the non/sub-resonant stochastic heating mechanism are important to explain the heating of Martian O+ layer. We will investigate this further by comparing the simulation results with the available MAVEN data. These simulated ion cyclotron waves are important to explain the heating of Martian O+ layer and have significant implications for future observations.

Investigating Microphysics of Intracluster Medium with Advanced Hydrodynamic Simulations and X-Ray Observations

NASA Astrophysics Data System (ADS)

Markevitch, Maxim

' -- sharp, arc-like contact discontinuities often observed in the cluster high- resolution X-ray images. Most of the observed cold fronts are very smooth, but some are visibly affected by instabilities. We will constrain the viscosity (largely independently of its exact physical nature) by including various levels of viscosity in the simulation of a strategically selected sample of cold front clusters. The forthcoming Astro-H mission opens another avenue to study the ICM viscosity and related phenomena -- by directly observing turbulence in merging and relaxed clusters. The turbulence in our merger simulations with varying viscosity could be directly compared to the levels eventually observed by Astro-H. Over the past several years, we have developed most of the machinery necessary for the above simulations by adding the appropriate code for diffusive physics into FLASH -- the high-resolution, grid-based magnetohydrodynamic code. As part of the proposed project, we will (a) add a few missing physical ingredients to the code and (b) simulate a sample of merging clusters, carefully selected from the XMM and Chandra archives to provide the strongest constraints on thermal conduction and viscosity in the ICM. The comparison with observations will result in strict limits on these microphysical quantities, adding an important element into the astrophysical foundation for the use of clusters as cosmological tools. The proposed research is directly relevant to the ROSES ATP solicitation, because it is a theoretical/numerical study of the physical processes in galaxy clusters, which will lead to predictions that can be tested with observations by the NASA space astrophysics missions XMM, Chandra and Astro-H.

The role of historical forcings in simulating the observed Atlantic Multidecadal Oscillation

NASA Astrophysics Data System (ADS)

Goes, L. M.; Cane, M. A.; Bellomo, K.; Clement, A. C.

2016-12-01

The variation in basin-wide North Atlantic sea surface temperatures (SST), known as the Atlantic multidecadal oscillation (AMO), affects climate throughout the Northern Hemisphere and tropics, yet the forcing mechanisms are not fully understood. Here, we analyze the AMO in the Coupled Model Intercomparison Project phase 5 (CMIP5) Pre-industrial (PI) and Historical (HIST) simulations to determine the role of historical climate forcings in producing the observed 20th century shifts in the AMO (OBS, 1865-2005). We evaluate whether the agreement between models and observations is better with historical forcings or without forcing - i.e. due to processes internal to the climate system, such as the Atlantic Meridional Overturning Circulation (AMOC). To do this we draw 141-year samples from 38 CMIP5 PI runs and compare the correlation between the PI and HIST AMO to the observed AMO. We find that in the majority of models (24 out of 38), it is very unlikely (less than 10% chance) that the unforced simulations produce agreement with observations that are as high as the forced simulations. We also compare the amplitude of the simulated AMO and find that 87% of models produce multi-decadal variance in the AMO with historical forcings that is very likely higher than without forcing, but most models underestimate the variance of the observed AMO. This indicates that over the 20th century external rather than internal forcing was crucial in setting the pace, phase and amplitude of the AMO.

The Arctic clouds from model simulations and long-term observations at Barrow, Alaska

NASA Astrophysics Data System (ADS)

Zhao, Ming

The Arctic is a region that is very sensitive to global climate change while also experiencing significant changes in its surface air temperature, sea-ice cover, atmospheric circulation, precipitation, snowfall, biogeochemical cycling, and land surface. Although previous studies have shown that the arctic clouds play an important role in the arctic climate changes, the arctic clouds are poorly understood and simulated in climate model due to limited observations. Furthermore, most of the studies were based on short-term experiments and typically only cover the warm seasons, which do not provide a full understanding of the seasonal cycle of arctic clouds. To address the above concerns and to improve our understanding of arctic clouds, six years of observational and retrieval data from 1999 to 2004 at the Atmospheric Radiation Management (ARM) Climate Research Facility (ACRF) North Slope of Alaska (NSA) Barrow site are used to understand the arctic clouds and related radiative processes. In particular, we focus on the liquid-ice mass partition in the mixed-phase cloud layer. Statistical results show that aerosol type and concentration are important factors that impact the mixed-phase stratus (MPS) cloud microphysical properties: liquid water path (LWP) and liquid water fraction (LWF) decrease with the increase of cloud condensation nuclei (CCN) number concentration; the high dust loading and dust occurrence in the spring are possible reasons for the much lower LWF than the other seasons. The importance of liquid-ice mass partition on surface radiation budgets was analyzed by comparing cloud longwave radiative forcings under the same LWP but different ice water path (IWP) ranges. Results show the ice phase enhance the surface cloud longwave (LW) forcing by 8˜9 W m-2 in the moderately thin MPS. This result provides an observational evidence on the aerosol glaciation effect in the moderately thin MPS, which is largely unknown so far. The above new insights are

A Comparison of Observed and Simulated 1990 – 2010 U.S. Ozone Trends

EPA Science Inventory

In this study, we analyze ozone concentrations from long-term (1990 – 2010) WRF-CMAQ simulations driven by year specific meteorology and emissions. These simulations allow us to compare observed and simulated ozone trends in order to evaluate the model’s ability to pr...

Evaluation of Cloud Microphysics in JMA-NHM Simulations Using Bin or Bulk Microphysical Schemes through Comparison with Cloud Radar Observations

NASA Technical Reports Server (NTRS)

Iguchi, Takamichi; Nakajima, Teruyuki; Khain, Alexander P.; Saito, Kazuo; Takemura, Toshihiko; Okamoto, Hajime; Nishizawa, Tomoaki; Tao, Wei-Kuo

2012-01-01

Numerical weather prediction (NWP) simulations using the Japan Meteorological Agency NonhydrostaticModel (JMA-NHM) are conducted for three precipitation events observed by shipborne or spaceborneW-band cloud radars. Spectral bin and single-moment bulk cloud microphysics schemes are employed separatelyfor an intercomparative study. A radar product simulator that is compatible with both microphysicsschemes is developed to enable a direct comparison between simulation and observation with respect to theequivalent radar reflectivity factor Ze, Doppler velocity (DV), and path-integrated attenuation (PIA). Ingeneral, the bin model simulation shows better agreement with the observed data than the bulk modelsimulation. The correction of the terminal fall velocities of snowflakes using those of hail further improves theresult of the bin model simulation. The results indicate that there are substantial uncertainties in the masssizeand sizeterminal fall velocity relations of snowflakes or in the calculation of terminal fall velocity of snowaloft. For the bulk microphysics, the overestimation of Ze is observed as a result of a significant predominanceof snow over cloud ice due to substantial deposition growth directly to snow. The DV comparison shows thata correction for the fall velocity of hydrometeors considering a change of particle size should be introducedeven in single-moment bulk cloud microphysics.

Improve observation-based ground-level ozone spatial distribution by compositing satellite and surface observations: A simulation experiment

NASA Astrophysics Data System (ADS)

Zhang, Yuzhong; Wang, Yuhang; Crawford, James; Cheng, Ye; Li, Jianfeng

2018-05-01

Obtaining the full spatial coverage of daily surface ozone fields is challenging because of the sparsity of the surface monitoring network and the difficulty in direct satellite retrievals of surface ozone. We propose an indirect satellite retrieval framework to utilize the information from satellite-measured column densities of tropospheric NO2 and CH2O, which are sensitive to the lower troposphere, to derive surface ozone fields. The method is applicable to upcoming geostationary satellites with high-quality NO2 and CH2O measurements. To prove the concept, we conduct a simulation experiment using a 3-D chemical transport model for July 2011 over the eastern US. The results show that a second order regression using both NO2 and CH2O column densities can be an effective predictor for daily maximum 8-h average ozone. Furthermore, this indirect retrieval approach is shown to be complementary to spatial interpolation of surface observations, especially in regions where the surface sites are sparse. Combining column observations of NO2 and CH2O with surface site measurements leads to an improved representation of surface ozone over simple kriging, increasing the R2 value from 0.53 to 0.64 at a surface site distance of 252 km. The improvements are even more significant with larger surface site distances. The simulation experiment suggests that the indirect satellite retrieval technique can potentially be a useful tool to derive the full spatial coverage of daily surface ozone fields if satellite observation uncertainty is moderate.

Simulated observations of young gravitationally unstable protoplanetary discs

NASA Astrophysics Data System (ADS)

Douglas, T. A.; Caselli, P.; Ilee, J. D.; Boley, A. C.; Hartquist, T. W.; Durisen, R. H.; Rawlings, J. M. C.

2013-08-01

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. Atacama Large Millimetre/sub-millimetre Array (ALMA) will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here, we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M⊙) self-gravitating disc embedded in a 10 M⊙ dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3→2, HCO+ 3→2, OCS 26→25 and H2CO 404→303 lines can be used to probe the disc structure and kinematics at all scales.

Rip current evidence by hydrodynamic simulations, bathymetric surveys and UAV observation

NASA Astrophysics Data System (ADS)

Benassai, Guido; Aucelli, Pietro; Budillon, Giorgio; De Stefano, Massimo; Di Luccio, Diana; Di Paola, Gianluigi; Montella, Raffaele; Mucerino, Luigi; Sica, Mario; Pennetta, Micla

2017-09-01

The prediction of the formation, spacing and location of rip currents is a scientific challenge that can be achieved by means of different complementary methods. In this paper the analysis of numerical and experimental data, including RPAS (remotely piloted aircraft systems) observations, allowed us to detect the presence of rip currents and rip channels at the mouth of Sele River, in the Gulf of Salerno, southern Italy. The dataset used to analyze these phenomena consisted of two different bathymetric surveys, a detailed sediment analysis and a set of high-resolution wave numerical simulations, completed with Google EarthTM images and RPAS observations. The grain size trend analysis and the numerical simulations allowed us to identify the rip current occurrence, forced by topographically constrained channels incised on the seabed, which were compared with observations.

Trend estimates of AERONET-observed and model-simulated AOT percentiles between 1993 and 2013

NASA Astrophysics Data System (ADS)

Yoon, Jongmin; Pozzer, Andrea; Chang, Dong Yeong; Lelieveld, Jos

2016-04-01

Recent Aerosol Optical thickness (AOT) trend studies used monthly or annual arithmetic means that discard details of the generally right-skewed AOT distributions. Potentially, such results can be biased by extreme values (including outliers). This study additionally uses percentiles (i.e., the lowest 5%, 25%, 50%, 75% and 95% of the monthly cumulative distributions fitted to Aerosol Robotic Network (AERONET)-observed and ECHAM/MESSy Atmospheric Chemistry (EMAC)-model simulated AOTs) that are less affected by outliers caused by measurement error, cloud contamination and occasional extreme aerosol events. Since the limited statistical representativeness of monthly percentiles and means can lead to bias, this study adopts the number of observations as a weighting factor, which improves the statistical robustness of trend estimates. By analyzing the aerosol composition of AERONET-observed and EMAC-simulated AOTs in selected regions of interest, we distinguish the dominant aerosol types and investigate the causes of regional AOT trends. The simulated and observed trends are generally consistent with a high correlation coefficient (R = 0.89) and small bias (slope±2σ = 0.75 ± 0.19). A significant decrease in EMAC-decomposed AOTs by water-soluble compounds and black carbon is found over the USA and the EU due to environmental regulation. In particular, a clear reversal in the AERONET AOT trend percentiles is found over the USA, probably related to the AOT diurnal cycle and the frequency of wildfires.

AOD trends during 2001-2010 from observations and model simulations

NASA Astrophysics Data System (ADS)

Pozzer, Andrea; de Meij, Alexander; Yoon, Jongmin; Astitha, Marina

2016-04-01

The trend of aerosol optical depth (AOD) between 2001 and 2010 is estimated globally and regionally from remote sensed observations by the MODIS (Moderate Resolution Imaging Spectroradiometer), MISR (Multi-angle Imaging SpectroRadiometer) and SeaWIFS (Sea-viewing Wide Field-of-view Sensor) satellite sensor. The resulting trends have been compared to model results from the EMAC (ECHAM5/MESSy Atmospheric Chemistry {[1]}), model. Although interannual variability is applied only to anthropogenic and biomass-burning emissions, the model is able to quantitatively reproduce the AOD trends as observed by MODIS, while some discrepancies are found when compared to MISR and SeaWIFS. An additional numerical simulation with the same model was performed, neglecting any temporal change in the emissions, i.e. with no interannual variability for any emission source. It is shown that decreasing AOD trends over the US and Europe are due to the decrease in the (anthropogenic) emissions. On contrary over the Sahara Desert and the Middle East region, the meteorological/dynamical changes in the last decade play a major role in driving the AOD trends. Further, over Southeast Asia, both meteorology and emissions changes are equally important in defining AOD trends {[2]}. Finally, decomposing the regional AOD trends into individual aerosol components reveals that the soluble components are the most dominant contributors to the total AOD, as their influence on the total AOD is enhanced by the aerosol water content. {[1]}: Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H., Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717-752, doi:10.5194/gmd-3-717-2010, 2010. {[2]}: Pozzer, A., de Meij, A., Yoon, J., Tost, H., Georgoulias, A. K., and Astitha, M.: AOD trends during 2001-2010 from observations and model simulations, Atmos. Chem. Phys., 15, 5521-5535, doi:10.5194/acp-15-5521-2015, 2015.

Congested traffic states in empirical observations and microscopic simulations

NASA Astrophysics Data System (ADS)

Treiber, Martin; Hennecke, Ansgar; Helbing, Dirk

2000-08-01

We present data from several German freeways showing different kinds of congested traffic forming near road inhomogeneities, specifically lane closings, intersections, or uphill gradients. The states are localized or extended, homogeneous or oscillating. Combined states are observed as well, like the coexistence of moving localized clusters and clusters pinned at road inhomogeneities, or regions of oscillating congested traffic upstream of nearly homogeneous congested traffic. The experimental findings are consistent with a recently proposed theoretical phase diagram for traffic near on-ramps [D. Helbing, A. Hennecke, and M. Treiber, Phys. Rev. Lett. 82, 4360 (1999)]. We simulate these situations with a continuous microscopic single-lane model, the ``intelligent driver model,'' using empirical boundary conditions. All observations, including the coexistence of states, are qualitatively reproduced by describing inhomogeneities with local variations of one model parameter. We show that the results of the microscopic model can be understood by formulating the theoretical phase diagram for bottlenecks in a more general way. In particular, a local drop of the road capacity induced by parameter variations has essentially the same effect as an on-ramp.

Arctic sea-ice diffusion from observed and simulated Lagrangian trajectories

NASA Astrophysics Data System (ADS)

Rampal, Pierre; Bouillon, Sylvain; Bergh, Jon; Ólason, Einar

2016-07-01

We characterize sea-ice drift by applying a Lagrangian diffusion analysis to buoy trajectories from the International Arctic Buoy Programme (IABP) dataset and from two different models: the standalone Lagrangian sea-ice model neXtSIM and the Eulerian coupled ice-ocean model used for the TOPAZ reanalysis. By applying the diffusion analysis to the IABP buoy trajectories over the period 1979-2011, we confirm that sea-ice diffusion follows two distinct regimes (ballistic and Brownian) and we provide accurate values for the diffusivity and integral timescale that could be used in Eulerian or Lagrangian passive tracers models to simulate the transport and diffusion of particles moving with the ice. We discuss how these values are linked to the evolution of the fluctuating displacements variance and how this information could be used to define the size of the search area around the position predicted by the mean drift. By comparing observed and simulated sea-ice trajectories for three consecutive winter seasons (2007-2011), we show how the characteristics of the simulated motion may differ from or agree well with observations. This comparison illustrates the usefulness of first applying a diffusion analysis to evaluate the output of modeling systems that include a sea-ice model before using these in, e.g., oil spill trajectory models or, more generally, to simulate the transport of passive tracers in sea ice.

Observed and Simulated Eddy Diffusivity Upstream of the Drake Passage

NASA Astrophysics Data System (ADS)

Tulloch, R.; Ferrari, R. M.; Marshall, J.

2012-12-01

Estimates of eddy diffusivity in the Southern Ocean are poorly constrained due to lack of observations. We compare the first direct estimate of isopycnal eddy diffusivity upstream of the Drake Passage (from Ledwell et al. 2011) with a numerical simulation. The estimate is computed from a point tracer release as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). We find that the observational diffusivity estimate of about 500m^2/s at 1500m depth is close to that computed in a data-constrained, 1/20th of a degree simulation of the Drake Passage region. This tracer estimate also agrees with Lagrangian float calculations in the model. The role of mean flow suppression of eddy diffusivity at shallower depths will also be discussed.

Constraining ammonia dairy emissions during NASA DISCOVER-AQ California: surface and airborne observation comparisons with CMAQ simulations

NASA Astrophysics Data System (ADS)

Miller, D. J.; Liu, Z.; Sun, K.; Tao, L.; Nowak, J. B.; Bambha, R.; Michelsen, H. A.; Zondlo, M. A.

2014-12-01

Agricultural ammonia (NH3) emissions are highly uncertain in current bottom-up inventories. Ammonium nitrate is a dominant component of fine aerosols in agricultural regions such as the Central Valley of California, especially during winter. Recent high resolution regional modeling efforts in this region have found significant ammonium nitrate and gas-phase NH3 biases during summer. We compare spatially-resolved surface and boundary layer gas-phase NH3 observations during NASA DISCOVER-AQ California with Community Multi-Scale Air Quality (CMAQ) regional model simulations driven by the EPA NEI 2008 inventory to constrain wintertime NH3 model biases. We evaluate model performance with respect to aerosol partitioning, mixing and deposition to constrain contributions to modeled NH3 concentration biases in the Central Valley Tulare dairy region. Ammonia measurements performed with an open-path mobile platform on a vehicle are gridded to 4 km resolution hourly background concentrations. A peak detection algorithm is applied to remove local feedlot emission peaks. Aircraft NH3, NH4+ and NO3- observations are also compared with simulations extracted along the flight tracks. We find NH3 background concentrations in the dairy region are underestimated by three to five times during winter and NH3 simulations are moderately correlated with observations (r = 0.36). Although model simulations capture NH3 enhancements in the dairy region, these simulations are biased low by 30-60 ppbv NH3. Aerosol NH4+ and NO3- are also biased low in CMAQ by three and four times respectively. Unlike gas-phase NH3, CMAQ simulations do not capture typical NH4+ or NO3- enhancements observed in the dairy region. In contrast, boundary layer height simulations agree well with observations within 13%. We also address observational constraints on simulated NH3 deposition fluxes. These comparisons suggest that NEI 2008 wintertime dairy emissions are underestimated by a factor of three to five. We test

Understanding radio polarimetry. V. Making matrix self-calibration work: processing of a simulated observation

NASA Astrophysics Data System (ADS)

Hamaker, J. P.

2006-09-01

Context: .This is Paper V in a series on polarimetric aperture synthesis based on the algebra of 2×2 matrices. Aims: .It validates the matrix self-calibration theory of the preceding Paper IV and outlines the algorithmic methods that had to be developed for its application. Methods: .New avenues of polarimetric self-calibration opened up in Paper IV are explored by processing a simulated observation. To focus on the polarimetric issues, it is set up so as to sidestep some of the common complications of aperture synthesis, yet properly represent physical conditions. In addition to a representative collection of observing errors, the simulated instrument includes strongly varying Faraday rotation and antennas with unequal feeds. The selfcal procedure is described in detail, including aspects in which it differs from the scalar case, and its effects are demonstrated with a number of intermediate image results. Results: .The simulation's outcome is in full agreement with the theory. The nonlinear matrix equations for instrumental parameters are readily solved by iteration; a convergence problem is easily remedied with a new ancillary algorithm. Instrumental effects are cleanly separated from source properties without reference to changes in parallactic rotation during the observation. Polarimetric images of high purity and dynamic range result. As theory predicts, polarimetric errors that are common to all sources inevitably remain; prior knowledge of the statistics of linear and circular polarization in a typical observed field can be applied to eliminate most of them. Conclusions: .The paper conclusively demonstrates that matrix selfcal per se is a viable method that may foster substantial advancement in the art of radio polarimetry. For its application in real observations, a number of issues must be resolved that matrix selfcal has in common with its scalar sibling, such as the treatment of extended sources and the familiar sampling and aliasing problems. The

Comparison of Malaria Simulations Driven by Meteorological Observations and Reanalysis Products in Senegal

PubMed Central

Diouf, Ibrahima; Rodriguez-Fonseca, Belen; Deme, Abdoulaye; Caminade, Cyril; Morse, Andrew P.; Cisse, Moustapha; Sy, Ibrahima; Dia, Ibrahima; Ermert, Volker; Ndione, Jacques-André; Gaye, Amadou Thierno

2017-01-01

The analysis of the spatial and temporal variability of climate parameters is crucial to study the impact of climate-sensitive vector-borne diseases such as malaria. The use of malaria models is an alternative way of producing potential malaria historical data for Senegal due to the lack of reliable observations for malaria outbreaks over a long time period. Consequently, here we use the Liverpool Malaria Model (LMM), driven by different climatic datasets, in order to study and validate simulated malaria parameters over Senegal. The findings confirm that the risk of malaria transmission is mainly linked to climate variables such as rainfall and temperature as well as specific landscape characteristics. For the whole of Senegal, a lag of two months is generally observed between the peak of rainfall in August and the maximum number of reported malaria cases in October. The malaria transmission season usually takes place from September to November, corresponding to the second peak of temperature occurring in October. Observed malaria data from the Programme National de Lutte contre le Paludisme (PNLP, National Malaria control Programme in Senegal) and outputs from the meteorological data used in this study were compared. The malaria model outputs present some consistencies with observed malaria dynamics over Senegal, and further allow the exploration of simulations performed with reanalysis data sets over a longer time period. The simulated malaria risk significantly decreased during the 1970s and 1980s over Senegal. This result is consistent with the observed decrease of malaria vectors and malaria cases reported by field entomologists and clinicians in the literature. The main differences between model outputs and observations regard amplitude, but can be related not only to reanalysis deficiencies but also to other environmental and socio-economic factors that are not included in this mechanistic malaria model framework. The present study can be considered as a

Comparison of Malaria Simulations Driven by Meteorological Observations and Reanalysis Products in Senegal.

PubMed

Diouf, Ibrahima; Rodriguez-Fonseca, Belen; Deme, Abdoulaye; Caminade, Cyril; Morse, Andrew P; Cisse, Moustapha; Sy, Ibrahima; Dia, Ibrahima; Ermert, Volker; Ndione, Jacques-André; Gaye, Amadou Thierno

2017-09-25

The analysis of the spatial and temporal variability of climate parameters is crucial to study the impact of climate-sensitive vector-borne diseases such as malaria. The use of malaria models is an alternative way of producing potential malaria historical data for Senegal due to the lack of reliable observations for malaria outbreaks over a long time period. Consequently, here we use the Liverpool Malaria Model (LMM), driven by different climatic datasets, in order to study and validate simulated malaria parameters over Senegal. The findings confirm that the risk of malaria transmission is mainly linked to climate variables such as rainfall and temperature as well as specific landscape characteristics. For the whole of Senegal, a lag of two months is generally observed between the peak of rainfall in August and the maximum number of reported malaria cases in October. The malaria transmission season usually takes place from September to November, corresponding to the second peak of temperature occurring in October. Observed malaria data from the Programme National de Lutte contre le Paludisme (PNLP, National Malaria control Programme in Senegal) and outputs from the meteorological data used in this study were compared. The malaria model outputs present some consistencies with observed malaria dynamics over Senegal, and further allow the exploration of simulations performed with reanalysis data sets over a longer time period. The simulated malaria risk significantly decreased during the 1970s and 1980s over Senegal. This result is consistent with the observed decrease of malaria vectors and malaria cases reported by field entomologists and clinicians in the literature. The main differences between model outputs and observations regard amplitude, but can be related not only to reanalysis deficiencies but also to other environmental and socio-economic factors that are not included in this mechanistic malaria model framework. The present study can be considered as a

Exploring the observational constraints on the simulation of brown carbon

NASA Astrophysics Data System (ADS)

Wang, X.; Heald, C. L.; Liu, J.; Weber, R. J.; Campuzano-Jost, P.; Jimenez, J. L.; Schwarz, J. P.; Perring, A. E.

2017-12-01

Brown carbon (BrC) is the component of organic aerosols (OA) which strongly absorbs solar radiation in the near-UV range of the spectrum. However the sources, evolution, and optical properties of BrC remain highly uncertain, and therefore constitute a large source of uncertainty in estimating the global direct radiative effect (DRE) of aerosols. Previous modeling studies of BrC optical properties and DRE have been unable to fully evaluate the skill of their simulations, given the lack of direct measurements of organic aerosol absorption. In this study, we develop a global model simulation (GEOS-Chem) of BrC and test it against BrC absorption measurements from two aircraft campaigns in the U.S. (SEAC4RS and DC3). To our knowledge, this is the first study to compare simulated BrC absorption with direct, continuous ambient measurements. We show that the laboratory-based BrC absorption properties from biomass burning overestimate the aircraft measurements of ambient BrC. In addition, applying a photochemical whitening scheme to simulated BrC is better able to represent the observed BrC absorption. These observations are consistent with a mass absorption coefficient (MAC) of freshly emitted biomass burning OA of 0.57m2g-1. Using the RRTMG model integrated with GEOS-Chem, we estimate that the all-sky top-of-atmosphere direct radiative effect (DRE) of OA is -0.350 Wm-2, 10% higher than that without consideration of BrC absorption. Therefore, our best estimate of the absorption DRE of BrC is +0.042 Wm-2. We suggest that the DRE of BrC has been overestimated previously due to the lack of observational constraints from direct measurements as well as neglect of the effects of photochemical whitening.

Ionospheric Simulation System for Satellite Observations and Global Assimilative Model Experiments - ISOGAME

NASA Technical Reports Server (NTRS)

Pi, Xiaoqing; Mannucci, Anthony J.; Verkhoglyadova, Olga; Stephens, Philip; Iijima, Bryron A.

2013-01-01

Modeling and imaging the Earth's ionosphere as well as understanding its structures, inhomogeneities, and disturbances is a key part of NASA's Heliophysics Directorate science roadmap. This invention provides a design tool for scientific missions focused on the ionosphere. It is a scientifically important and technologically challenging task to assess the impact of a new observation system quantitatively on our capability of imaging and modeling the ionosphere. This question is often raised whenever a new satellite system is proposed, a new type of data is emerging, or a new modeling technique is developed. The proposed constellation would be part of a new observation system with more low-Earth orbiters tracking more radio occultation signals broadcast by Global Navigation Satellite System (GNSS) than those offered by the current GPS and COSMIC observation system. A simulation system was developed to fulfill this task. The system is composed of a suite of software that combines the Global Assimilative Ionospheric Model (GAIM) including first-principles and empirical ionospheric models, a multiple- dipole geomagnetic field model, data assimilation modules, observation simulator, visualization software, and orbit design, simulation, and optimization software.

Features of the accretion in the EX Hydrae system: Results of numerical simulation

NASA Astrophysics Data System (ADS)

Isakova, P. B.; Zhilkin, A. G.; Bisikalo, D. V.; Semena, A. N.; Revnivtsev, M. G.

2017-07-01

A two-dimensional numerical model in the axisymmetric approximation that describes the flow structure in the magnetosphere of the white dwarf in the EX Hya system has been developed. Results of simulations show that the accretion in EX Hya proceeds via accretion columns, which are not closed and have curtain-like shapes. The thickness of the accretion curtains depends only weakly on the thickness of the accretion disk. This thickness developed in the simulations does not agree with observations. It is concluded that the main reason for the formation of thick accretion curtains in the model is the assumption that the magnetic field penetrates fully into the plasma of the disk. An analysis based on simple estimates shows that a diamagnetic disk that fully or partially shields the magnetic field of the star may be a more attractive explanation for the observed features of the accretion in EX Hya.

Simulations of the observation of clouds and aerosols with the Experimental Lidar in Space Equipment system.

PubMed

Liu, Z; Voelger, P; Sugimoto, N

2000-06-20

We carried out a simulation study for the observation of clouds and aerosols with the Japanese Experimental Lidar in Space Equipment (ELISE), which is a two-wavelength backscatter lidar with three detection channels. The National Space Development Agency of Japan plans to launch the ELISE on the Mission Demonstrate Satellite 2 (MDS-2). In the simulations, the lidar return signals for the ELISE are calculated for an artificial, two-dimensional atmospheric model including different types of clouds and aerosols. The signal detection processes are simulated realistically by inclusion of various sources of noise. The lidar signals that are generated are then used as input for simulations of data analysis with inversion algorithms to investigate retrieval of the optical properties of clouds and aerosols. The results demonstrate that the ELISE can provide global data on the structures and optical properties of clouds and aerosols. We also conducted an analysis of the effects of cloud inhomogeneity on retrievals from averaged lidar profiles. We show that the effects are significant for space lidar observations of optically thick broken clouds.

Presenting simulation results in a nested loop plot.

PubMed

Rücker, Gerta; Schwarzer, Guido

2014-12-12

Statisticians investigate new methods in simulations to evaluate their properties for future real data applications. Results are often presented in a number of figures, e.g., Trellis plots. We had conducted a simulation study on six statistical methods for estimating the treatment effect in binary outcome meta-analyses, where selection bias (e.g., publication bias) was suspected because of apparent funnel plot asymmetry. We varied five simulation parameters: true treatment effect, extent of selection, event proportion in control group, heterogeneity parameter, and number of studies in meta-analysis. In combination, this yielded a total number of 768 scenarios. To present all results using Trellis plots, 12 figures were needed. Choosing bias as criterion of interest, we present a 'nested loop plot', a diagram type that aims to have all simulation results in one plot. The idea was to bring all scenarios into a lexicographical order and arrange them consecutively on the horizontal axis of a plot, whereas the treatment effect estimate is presented on the vertical axis. The plot illustrates how parameters simultaneously influenced the estimate. It can be combined with a Trellis plot in a so-called hybrid plot. Nested loop plots may also be applied to other criteria such as the variance of estimation. The nested loop plot, similar to a time series graph, summarizes all information about the results of a simulation study with respect to a chosen criterion in one picture and provides a suitable alternative or an addition to Trellis plots.

Multi-spacecraft observations and transport simulations of solar energetic particles for the May 17th 2012 event

NASA Astrophysics Data System (ADS)

Battarbee, M.; Guo, J.; Dalla, S.; Wimmer-Schweingruber, R.; Swalwell, B.; Lawrence, D. J.

2018-05-01

Context. The injection, propagation and arrival of solar energetic particles (SEPs) during eruptive solar events is an important and current research topic of heliospheric physics. During the largest solar events, particles may have energies up to a few GeVs and sometimes even trigger ground-level enhancements (GLEs) at Earth. These large SEP events are best investigated through multi-spacecraft observations. Aims: We aim to study the first GLE-event of solar cycle 24, from 17th May 2012, using data from multiple spacecraft (SOHO, GOES, MSL, STEREO-A, STEREO-B and MESSENGER). These spacecraft are located throughout the inner heliosphere, at heliocentric distances between 0.34 and 1.5 astronomical units (au), covering nearly the whole range of heliospheric longitudes. Methods: We present and investigate sub-GeV proton time profiles for the event at several energy channels, obtained via different instruments aboard the above spacecraft. We investigated issues caused by magnetic connectivity, and present results of three-dimensional SEP propagation simulations. We gathered virtual time profiles and perform qualitative and quantitative comparisons with observations, assessed longitudinal injection and transport effects as well as peak intensities. Results: We distinguish different time profile shapes for well-connected and weakly connected observers, and find our onset time analysis to agree with this distinction. At select observers, we identify an additional low-energy component of Energetic Storm Particles (ESPs). Using well-connected observers for normalisation, our simulations are able to accurately recreate both time profile shapes and peak intensities at multiple observer locations. Conclusions: This synergetic approach combining numerical modelling with multi-spacecraft observations is crucial for understanding the propagation of SEPs within the interplanetary magnetic field. Our novel analysis provides valuable proof of the ability to simulate SEP propagation

Results of the Simulation and Assimilation of Doppler Wind Lidar Observations in Preparation for European Space Agency's Aeolus Mission

NASA Technical Reports Server (NTRS)

McCarty, Will

2011-01-01

With the launch of the European Space Agency's Aeolus Mission in 2013, direct spaceborne measurements of vertical wind profiles are imminent via Doppler wind lidar technology. Part of the preparedness for such missions is the development of the proper data assimilation methodology for handling such observations. Since no heritage measurements exist in space, the Joint Observing System Simulation Experiment (Joint OSSE) framework has been utilized to generate a realistic proxy dataset as a precursor to flight. These data are being used for the development of the Gridpoint Statistical Interpolation (GSI) data assimilation system utilized at a number of centers through the United States including the Global Modeling and Assimilation Office (GMAO) at NASA/Goddard Space Flight Center and at the National Centers for Environmental Prediction (NOAA/NWS/NCEP) as an activity through the Joint Center for Satellite Data Assimilation. An update of this ongoing effort will be presented, including the methodology of proxy data generation, the limitations of the proxy data, the handling of line-of-sight wind measurements within the GSI, and the impact on both analyses and forecasts with the addition of the new data type.

Skill of ship-following large-eddy simulations in reproducing MAGIC observations across the northeast Pacific stratocumulus to cumulus transition region

NASA Astrophysics Data System (ADS)

McGibbon, J.; Bretherton, C. S.

2017-06-01

During the Marine ARM GPCI Investigation of Clouds (MAGIC) in October 2011 to September 2012, a container ship making periodic cruises between Los Angeles, CA, and Honolulu, HI, was instrumented with surface meteorological, aerosol and radiation instruments, a cloud radar and ceilometer, and radiosondes. Here large-eddy simulation (LES) is performed in a ship-following frame of reference for 13 four day transects from the MAGIC field campaign. The goal is to assess if LES can skillfully simulate the broad range of observed cloud characteristics and boundary layer structure across the subtropical stratocumulus to cumulus transition region sampled during different seasons and meteorological conditions. Results from Leg 15A, which sampled a particularly well-defined stratocumulus to cumulus transition, demonstrate the approach. The LES reproduces the observed timing of decoupling and transition from stratocumulus to cumulus and matches the observed evolution of boundary layer structure, cloud fraction, liquid water path, and precipitation statistics remarkably well. Considering the simulations of all 13 cruises, the LES skillfully simulates the mean diurnal variation of key measured quantities, including liquid water path (LWP), cloud fraction, measures of decoupling, and cloud radar-derived precipitation. The daily mean quantities are well represented, and daily mean LWP and cloud fraction show the expected correlation with estimated inversion strength. There is a -0.6 K low bias in LES near-surface air temperature that results in a high bias of 5.6 W m-2 in sensible heat flux (SHF). Overall, these results build confidence in the ability of LES to represent the northeast Pacific stratocumulus to trade cumulus transition region.